PLUTO
4.4-patch2
|
GLM module implementation. More...
#include "pluto.h"
Functions | |
void | GLM_Solve (const Sweep *sweep, int beg, int end, Grid *grid) |
void | GLM_Source (const Data *data, double dt, Grid *grid) |
void | GLM_ExtendedSource (const Sweep *sweep, double dt, int beg, int end, Grid *grid) |
void | GLM_Init (const Data *d, const timeStep *Dts, Grid *grid) |
void | GLM_ComputeDivE (const Sweep *sweep, Grid *grid) |
Variables | |
double | glm_ch = -1.0 |
Contains functions for the GLM module.
Reference
Compute the divergence of E using Godunov fluxes previously obtained at cell interfaces. This function may be used in Resistive RMHD.
Add source terms to the right hand side of the conservative equations, momentum and energy equations only. This yields the extended GLM equations given by Eq. (24a)–(24c) in
"Hyperbolic Divergence cleaning for the MHD Equations" Dedner et al. (2002), JcP, 175, 645
Initialize the maximum propagation speed glm_ch at the beginning of integration cycle.
Solve the 2x2 linear hyperbolic GLM-MHD system given by the divergence cleaning approach. Modify inteface states (Bx and psi components) for input to full Riemann problem. We use Eq. (42) of Dedner et al (2002)
[in,out] | sweep | pointer to a Sweep structure |
[in] | beg | starting index of computation |
[in] | end | final index of computation |
[in] | grid | pointer to Grid structure |
The purpose of this function is two-fold:
The following MAPLE script has been used
Include the damping source term of the Lagrangian multiplier equation in a split fashion for the mixed GLM formulation. Ref. Mignone & Tzeferacos, JCP (2010) 229, 2117, Equation (27).
double glm_ch = -1.0 |
The propagation speed of divergence error.