Integrate MHD equations using 3D version of the directionally unsplit MUSCL-Hancock (VL) integrator. More...

Functions
static void	integrate_emf1_corner (const GridS *pG)
	Integrates face centered B-fluxes to compute corner EMFs.
static void	integrate_emf2_corner (const GridS *pG)
	Integrates face centered B-fluxes to compute corner EMFs.
static void	integrate_emf3_corner (const GridS *pG)
	Integrates face centered B-fluxes to compute corner EMFs.
static void	FixCell (GridS *pG, Int3Vect)
	Uses first order fluxes to fix negative d,P or superluminal v.
void	integrate_3d_vl (DomainS *pD)
	3D van Leer unsplit integrator for MHD.
void	integrate_init_3d (MeshS *pM)
	Allocate temporary integration arrays.
void	integrate_destruct_3d (void)
	Free temporary integration arrays.
Variables
static Prim1DS ***	Wl_x1Face = NULL
static Prim1DS ***	Wr_x1Face = NULL
static Prim1DS ***	Wl_x2Face = NULL
static Prim1DS ***	Wr_x2Face = NULL
static Prim1DS ***	Wl_x3Face = NULL
static Prim1DS ***	Wr_x3Face = NULL
static Cons1DS ***	x1Flux = NULL
static Cons1DS ***	x2Flux = NULL
static Cons1DS ***	x3Flux = NULL
static Cons1DS ***	x1FluxP = NULL
static Cons1DS ***	x2FluxP = NULL
static Cons1DS ***	x3FluxP = NULL
static Real ***	emf1P = NULL
static Real ***	emf2P = NULL
static Real ***	emf3P = NULL
static Real ***	S = NULL
static Real ***	Shalf = NULL
static Real ***	x1FluxS = NULL
static Real ***	x2FluxS = NULL
static Real ***	x3FluxS = NULL
static Real ***	x1FluxSP = NULL
static Real ***	x2FluxSP = NULL
static Real ***	x3FluxSP = NULL
static Real ***	B1_x1Face = NULL
static Real ***	B2_x2Face = NULL
static Real ***	B3_x3Face = NULL
static Real ***	emf1 = NULL
static Real ***	emf2 = NULL
static Real ***	emf3 = NULL
static Real ***	emf1_cc = NULL
static Real ***	emf2_cc = NULL
static Real ***	emf3_cc = NULL
static Real *	Bxc = NULL
static Real *	Bxi = NULL
static Prim1DS *	W1d = NULL
static Prim1DS *	Wl = NULL
static Prim1DS *	Wr = NULL
static Cons1DS *	U1d = NULL
static Cons1DS *	Ul = NULL
static Cons1DS *	Ur = NULL
static PrimS ***	W = NULL
static ConsS ***	Uhalf = NULL
static PrimS ***	Whalf = NULL
Real	etah
static Real ***	eta1 = NULL
static Real ***	eta2 = NULL
static Real ***	eta3 = NULL

Detailed Description

Integrate MHD equations using 3D version of the directionally unsplit MUSCL-Hancock (VL) integrator.

PURPOSE: Integrate MHD equations using 3D version of the directionally unsplit MUSCL-Hancock (VL) integrator. The variables updated are:

U.[d,M1,M2,M3,E,B1c,B2c,B3c,s] -- where U is of type ConsS
B1i, B2i, B3i -- interface magnetic field Also adds gravitational source terms, self-gravity, and the H-correction of Sanders et al.

For adb hydro, requires (9*Cons1DS + 3*Real + 1*ConsS) = 53 3D arrays

For adb mhd, requires (9*Cons1DS + 9*Real + 1*ConsS) = 80 3D arrays

REFERENCE:

J.M Stone & T.A. Gardiner, "A simple, unsplit Godunov method for multidimensional MHD", NewA 14, 139 (2009)

R. Sanders, E. Morano, & M.-C. Druguet, "Multidimensinal dissipation for upwind schemes: stability and applications to gas dynamics", JCP, 145, 511 (1998)

CONTAINS PUBLIC FUNCTIONS:

Definition in file integrate_3d_vl_sr.c.