Integrate MHD equations in 2D using the directionally unsplit CTU method of Colella (1990). More...
Go to the source code of this file.
Functions | |
| static void | integrate_emf3_corner (GridS *pG) |
| Corner EMF terms in the upwind CT method. | |
| void | integrate_2d_ctu (DomainS *pD) |
| CTU integrator in 2D. | |
| void | integrate_init_2d (MeshS *pM) |
| Allocate temporary integration arrays. | |
| void | integrate_destruct_2d (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 ** | x1Flux = NULL |
| static Cons1DS ** | x2Flux = NULL |
| static Real ** | B1_x1Face = NULL |
| static Real ** | B2_x2Face = NULL |
| static Real ** | emf3 = 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 Cons1DS * | Ul = NULL |
| static Cons1DS * | Ur = NULL |
| static Real ** | dhalf = NULL |
| static Real ** | phalf = NULL |
| Real | etah |
| static Real ** | eta1 = NULL |
| static Real ** | eta2 = NULL |
| static Real ** | geom_src = NULL |
Detailed Description
Integrate MHD equations in 2D using the directionally unsplit CTU method of Colella (1990).
PURPOSE: Integrate MHD equations in 2D using the directionally unsplit CTU method 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 -- interface magnetic field Also adds gravitational source terms, self-gravity, optically-thin cooling, and H-correction of Sanders et al.
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", JCP, 205, 509 (2005)
- 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_2d_ctu.c.
Function Documentation
| void integrate_2d_ctu | ( | DomainS * | pD | ) |
CTU integrator in 2D.
The numbering of steps follows the numbering in the 3D version. NOT ALL STEPS ARE NEEDED IN 2D.
Definition at line 97 of file integrate_2d_ctu.c.
References B1_x1Face, ConsS::B1c, GridS::B1i, B2_x2Face, ConsS::B2c, GridS::B2i, 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, dby, dhalf, dim, GridS::dt, GridS::dx1, GridS::dx2, E, ConsS::E, Cons1DS::E, GPCouple::Eloss, emf3, emf3_cc, eta1, eta2, etah, exchange_feedback(), GPCouple::fb1, GPCouple::fb2, GPCouple::fb3, feedback_predictor(), fluxes(), four_pi_G, Gamma, Gamma_1, geom_src, 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_emf3_corner(), Integrate_Particles(), GridS::is, Iso_csound2, iu, GridS::je, jl, GridS::js, ju, GridS::ks, lr_states(), ConsS::M1, ConsS::M2, ConsS::M3, Cons1DS::Mx, Cons1DS::My, GridOvrlpS::myEMF3, GridOvrlpS::myFlx, Cons1DS::Mz, n, GridS::NCGrid, GridS::NPGrid, Omega_0, Prim1DS::P, Cons1DS::Pflux, GridS::PGrid, phalf, GridS::Phi, Prim1D_to_Cons1D(), qshear, GridS::r, r, GridS::ri, ri, Cons1DS::s, ConsS::s, s, ShBoxCoord, ShearingBoxPot, StaticGravPot, GridS::U, Ur_x1Face, Ur_x2Face, Prim1DS::Vx, Prim1DS::Vy, Prim1DS::Vz, Wl, Wr, GridS::x1MassFlux, x1vc(), x2Flux, and GridS::x2MassFlux.
| void integrate_destruct_2d | ( | void | ) |
Free temporary integration arrays.
Definition at line 1981 of file integrate_2d_ctu.c.
Referenced by integrate_destruct().
| static void integrate_emf3_corner | ( | GridS * | pG | ) | [static] |
Corner EMF terms in the upwind CT method.
Definition at line 2028 of file integrate_2d_ctu.c.
References Cons1DS::By, Bz(), Cons1DS::Bz, d(), emf3, emf3_cc, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::r, GridS::ri, and x2Flux.
Referenced by integrate_2d_ctu().
| void integrate_init_2d | ( | MeshS * | pM | ) |
Allocate temporary integration arrays.
Definition at line 1885 of file integrate_2d_ctu.c.
Referenced by integrate_init().
Variable Documentation
Definition at line 56 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 56 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 61 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 61 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 66 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 57 of file integrate_2d_ctu.c.
Referenced by FixCell(), integrate_2d_ctu(), integrate_2d_vl(), integrate_3d_vl(), integrate_destruct_2d(), integrate_destruct_3d(), integrate_emf3_corner(), integrate_init_2d(), and integrate_init_3d().
Definition at line 57 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), integrate_emf3_corner(), and integrate_init_2d().
Definition at line 71 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 71 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 33 of file roe.c.
Referenced by fluxes(), integrate_2d_ctu(), integrate_2d_vl(), integrate_3d_ctu(), and integrate_3d_vl().
Definition at line 76 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 66 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 63 of file integrate_2d_ctu.c.
Definition at line 63 of file integrate_2d_ctu.c.
Referenced by integrate_destruct_2d(), and integrate_init_2d().
Definition at line 50 of file integrate_2d_ctu.c.
Definition at line 51 of file integrate_2d_ctu.c.
Referenced by integrate_2d_vl(), and integrate_3d_vl().
Definition at line 63 of file integrate_2d_ctu.c.
Referenced by integrate_destruct_2d(), and integrate_init_2d().
Definition at line 50 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 51 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_2d_vl(), integrate_3d_vl(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 62 of file integrate_2d_ctu.c.
Definition at line 62 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 62 of file integrate_2d_ctu.c.
Referenced by integrate_2d_ctu(), integrate_destruct_2d(), and integrate_init_2d().
Definition at line 52 of file integrate_2d_ctu.c.
Definition at line 52 of file integrate_2d_ctu.c.
Referenced by FixCell(), integrate_2d_ctu(), integrate_2d_vl(), integrate_3d_vl(), integrate_destruct_2d(), integrate_destruct_3d(), integrate_emf1_corner(), integrate_emf3_corner(), integrate_init_2d(), integrate_init_3d(), viscosity(), viscosity_destruct(), viscosity_init(), ViscStress_aniso(), and ViscStress_iso().
