Athena: integrators/integrate_3d_ctu.c File Reference

Integrate MHD equations using 3D version of the directionally unsplit CTU integrator of Colella (1990). More...

Go to the source code of this file.

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.
void	integrate_3d_ctu (DomainS *pD)
	3D CTU integrator for MHD using 6-solve method
void	integrate_init_3d (MeshS *pM)
	Allocate temporary integration arrays.
void	integrate_destruct_3d (void)
	Free temporary integration arrays.
Variables
static Cons1DS ***	Ul_x1Face = NULL
static Cons1DS ***	Ur_x1Face = NULL
static Cons1DS ***	Ul_x2Face = NULL
static Cons1DS ***	Ur_x2Face = NULL
static Cons1DS ***	Ul_x3Face = NULL
static Cons1DS ***	Ur_x3Face = NULL
Cons1DS ***	x1Flux = NULL
Cons1DS ***	x2Flux = NULL
Cons1DS ***	x3Flux = NULL
static Real ***	B1_x1Face = NULL
static Real ***	B2_x2Face = NULL
static Real ***	B3_x3Face = NULL
Real ***	emf1 = NULL
Real ***	emf2 = NULL
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 *	W = NULL
static Prim1DS *	Wl = NULL
static Prim1DS *	Wr = NULL
static Cons1DS *	U1d = NULL
static Real ***	dhalf = NULL
static Real ***	phalf = NULL
Real	etah
static Real ***	eta1 = NULL
static Real ***	eta2 = NULL
static Real ***	eta3 = NULL
static Real **	remapEyiib = NULL
static Real **	remapEyoib = NULL
static Real ***	geom_src = NULL

Detailed Description

Integrate MHD equations using 3D version of the directionally unsplit CTU integrator of Colella (1990).

PURPOSE: Integrate MHD equations using 3D version of the directionally unsplit CTU integrator of Colella (1990). 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, optically thin cooling, shearing box source terms, and the H-correction of Sanders et al.
For adb hydro, requires (9*Cons1DS + 3*Real) = 48 3D arrays
For adb mhd, requires (9*Cons1DS + 10*Real) = 73 3D arrays The H-correction of Sanders et al. adds another 3 arrays.

REFERENCES:

P. Colella, "Multidimensional upwind methods for hyperbolic conservation laws", JCP, 87, 171 (1990)

T. Gardiner & J.M. Stone, "An unsplit Godunov method for ideal MHD via constrained transport in three dimensions", JCP, 227, 4123 (2008)

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

J.M. Stone et al., "Athena: A new code for astrophysical MHD", ApJS, 178, 137 (2008)

CONTAINS PUBLIC FUNCTIONS:

Definition in file integrate_3d_ctu.c.

Function Documentation

void integrate_3d_ctu ( DomainS * pD )

3D CTU integrator for MHD using 6-solve method

Definition at line 110 of file integrate_3d_ctu.c.

References B1_x1Face, ConsS::B1c, GridS::B1i, B2_x2Face, ConsS::B2c, GridS::B2i, B3_x3Face, ConsS::B3c, GridS::B3i, Bxc, Bxi, Prim1DS::By, Cons1DS::By, Prim1DS::Bz, Bz(), Cons1DS::Bz, cc_pos(), cfast(), GridS::CGrid, Cons1D_to_Prim1D(), CoolingFunc, GridS::Coup, d(), Prim1DS::d, ConsS::d, Cons1DS::d, d1, dhalf, dim, GridS::dt, GridS::dx1, GridS::dx2, GridS::dx3, E, ConsS::E, Cons1DS::E, GPCouple::Eloss, emf1, emf1_cc, emf2, emf2_cc, emf3, emf3_cc, eta1, eta2, eta3, etah, exchange_feedback(), GPCouple::fb1, GPCouple::fb2, GPCouple::fb3, feedback_predictor(), fluxes(), four_pi_G, Gamma, Gamma_1, geom_src, get_myGridIndex(), grav_mean_rho, DomainS::Grid, GPCouple::grid_cs, GPCouple::grid_d, GPCouple::grid_v1, GPCouple::grid_v2, GPCouple::grid_v3, GridS::ie, GridOvrlpS::ijke, GridOvrlpS::ijks, il, integrate_emf1_corner(), integrate_emf2_corner(), integrate_emf3_corner(), Integrate_Particles(), GridS::is, Iso_csound2, iu, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, lr_states(), ConsS::M1, ConsS::M2, ConsS::M3, Cons1DS::Mx, Cons1DS::My, my_iproc, my_jproc, my_kproc, GridOvrlpS::myEMF1, GridOvrlpS::myEMF2, GridOvrlpS::myEMF3, GridOvrlpS::myFlx, myID_Comm_world, Cons1DS::Mz, n, GridS::NCGrid, DomainS::NGrid, GridS::NPGrid, Omega_0, OrbitalProfile, Prim1DS::P, Cons1DS::Pflux, GridS::PGrid, phalf, GridS::Phi, Prim1D_to_Cons1D(), qshear, GridS::r, r, RemapEy_ix1(), RemapEy_ox1(), remapEyiib, remapEyoib, GridS::ri, ri, Cons1DS::s, ConsS::s, s, ShearingBoxPot, StaticGravPot, GridS::U, Ur_x1Face, Ur_x2Face, Ur_x3Face, Prim1DS::Vx, Prim1DS::Vy, Prim1DS::Vz, Wl, Wr, GridS::x1MassFlux, x1vc(), x2Flux, GridS::x2MassFlux, x3Flux, and GridS::x3MassFlux.

Here is the call graph for this function:

void integrate_destruct_3d ( void )

Free temporary integration arrays.

Definition at line 3498 of file integrate_3d_ctu.c.

Referenced by integrate_destruct().

Here is the caller graph for this function:

static void integrate_emf1_corner ( const GridS * pG ) [static]

Integrates face centered B-fluxes to compute corner EMFs.

Note:

x1Flux.By = VxBy - BxVy = v1*b2-b1*v2 = -EMFZ
x1Flux.Bz = VxBz - BxVz = v1*b3-b1*v3 = EMFY
x2Flux.By = VxBy - BxVy = v2*b3-b2*v3 = -EMFX
x2Flux.Bz = VxBz - BxVz = v2*b1-b2*v1 = EMFZ
x3Flux.By = VxBy - BxVy = v3*b1-b3*v1 = -EMFY
x3Flux.Bz = VxBz - BxVz = v3*b2-b3*v2 = EMFX

Definition at line 3567 of file integrate_3d_ctu.c.

References Cons1DS::By, Cons1DS::Bz, d(), emf1, emf1_cc, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, x2Flux, and x3Flux.

Referenced by integrate_3d_ctu().

Here is the call graph for this function:

Here is the caller graph for this function:

static void integrate_emf2_corner ( const GridS * pG ) [static]

Integrates face centered B-fluxes to compute corner EMFs.

Note:

x1Flux.By = VxBy - BxVy = v1*b2-b1*v2 = -EMFZ
x1Flux.Bz = VxBz - BxVz = v1*b3-b1*v3 = EMFY
x2Flux.By = VxBy - BxVy = v2*b3-b2*v3 = -EMFX
x2Flux.Bz = VxBz - BxVz = v2*b1-b2*v1 = EMFZ
x3Flux.By = VxBy - BxVy = v3*b1-b3*v1 = -EMFY
x3Flux.Bz = VxBz - BxVz = v3*b2-b3*v2 = EMFX

Definition at line 3636 of file integrate_3d_ctu.c.

References Cons1DS::By, Cons1DS::Bz, d(), emf2, emf2_cc, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, and x3Flux.

Referenced by integrate_3d_ctu().

Here is the call graph for this function:

Here is the caller graph for this function:

static void integrate_emf3_corner ( const GridS * pG ) [static]

Integrates face centered B-fluxes to compute corner EMFs.

Note:

x1Flux.By = VxBy - BxVy = v1*b2-b1*v2 = -EMFZ
x1Flux.Bz = VxBz - BxVz = v1*b3-b1*v3 = EMFY
x2Flux.By = VxBy - BxVy = v2*b3-b2*v3 = -EMFX
x2Flux.Bz = VxBz - BxVz = v2*b1-b2*v1 = EMFZ
x3Flux.By = VxBy - BxVy = v3*b1-b3*v1 = -EMFY
x3Flux.Bz = VxBz - BxVz = v3*b2-b3*v2 = EMFX

Definition at line 3705 of file integrate_3d_ctu.c.

References Cons1DS::By, Bz(), Cons1DS::Bz, d(), emf3, emf3_cc, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, GridS::r, GridS::ri, and x2Flux.

Referenced by integrate_3d_ctu().

Here is the call graph for this function:

Here is the caller graph for this function:

void integrate_init_3d ( MeshS * pM )

Allocate temporary integration arrays.

Definition at line 3374 of file integrate_3d_ctu.c.

Referenced by integrate_init().

Here is the caller graph for this function:

Variable Documentation

Real*** B1_x1Face = NULL [static]

Definition at line 63 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real *** B2_x2Face = NULL [static]

Definition at line 63 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real *** B3_x3Face = NULL [static]

Definition at line 63 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real* Bxc = NULL [static]

Definition at line 69 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real * Bxi = NULL [static]

Definition at line 69 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real*** dhalf = NULL [static]

Definition at line 74 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real*** emf1 = NULL

Definition at line 64 of file integrate_3d_ctu.c.

Referenced by FixCell(), integrate_3d_ctu(), integrate_3d_vl(), integrate_destruct_3d(), integrate_emf1_corner(), and integrate_init_3d().

Real*** emf1_cc = NULL [static]

Definition at line 65 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), integrate_emf1_corner(), and integrate_init_3d().

Real *** emf2 = NULL

Definition at line 64 of file integrate_3d_ctu.c.

Referenced by FixCell(), integrate_3d_ctu(), integrate_3d_vl(), integrate_destruct_3d(), integrate_emf2_corner(), and integrate_init_3d().

Real *** emf2_cc = NULL [static]

Definition at line 65 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), integrate_emf2_corner(), and integrate_init_3d().

Real *** emf3 = NULL

Definition at line 64 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), integrate_emf3_corner(), and integrate_init_3d().

Real *** emf3_cc = NULL [static]

Definition at line 65 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), integrate_emf3_corner(), and integrate_init_3d().

Real*** eta1 = NULL [static]

Definition at line 79 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real *** eta2 = NULL [static]

Definition at line 79 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real *** eta3 = NULL [static]

Definition at line 79 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Definition at line 33 of file roe.c.

Real*** geom_src = NULL [static]

Definition at line 89 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real *** phalf = NULL [static]

Definition at line 74 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real** remapEyiib = NULL [static]

Definition at line 84 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Real ** remapEyoib = NULL [static]

Definition at line 84 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Cons1DS* U1d = NULL [static]

Definition at line 71 of file integrate_3d_ctu.c.

Cons1DS*** Ul_x1Face = NULL [static]

Definition at line 56 of file integrate_3d_ctu.c.

Cons1DS*** Ul_x2Face = NULL [static]

Definition at line 57 of file integrate_3d_ctu.c.

Cons1DS*** Ul_x3Face = NULL [static]

Definition at line 58 of file integrate_3d_ctu.c.

Referenced by integrate_3d_vl().

Cons1DS *** Ur_x1Face = NULL [static]

Definition at line 56 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Cons1DS *** Ur_x2Face = NULL [static]

Definition at line 57 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Cons1DS *** Ur_x3Face = NULL [static]

Definition at line 58 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_3d_vl(), integrate_destruct_3d(), and integrate_init_3d().

Prim1DS* W = NULL [static]

Definition at line 70 of file integrate_3d_ctu.c.

Prim1DS * Wl = NULL [static]

Definition at line 70 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Prim1DS * Wr = NULL [static]

Definition at line 70 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), and integrate_init_3d().

Cons1DS*** x1Flux = NULL

Definition at line 59 of file integrate_3d_ctu.c.

Cons1DS *** x2Flux = NULL

Definition at line 59 of file integrate_3d_ctu.c.

Referenced by integrate_3d_ctu(), integrate_destruct_3d(), integrate_emf1_corner(), integrate_emf3_corner(), and integrate_init_3d().

Cons1DS *** x3Flux = NULL

Definition at line 59 of file integrate_3d_ctu.c.

Referenced by FixCell(), integrate_3d_ctu(), integrate_3d_vl(), integrate_destruct_3d(), integrate_emf1_corner(), integrate_emf2_corner(), integrate_init_3d(), viscosity(), viscosity_destruct(), viscosity_init(), ViscStress_aniso(), and ViscStress_iso().