A dynamic test of force balance using a B_R-only, time-dependent, non-axisymmetric magnetic field. More...
Go to the source code of this file.
Functions | |
Real | d (const Real x1, const Real x2, const Real x3) |
Real | M2 (const Real x1, const Real x2, const Real x3) |
Real | B1 (const Real x1, const Real x2, const Real x3) |
1-component of B field | |
Real | B1i (const Real x2) |
1-component of B field at interface | |
Real | Pgas (const Real x1, const Real x2, const Real x3) |
gas pressure | |
static Real | grav_pot (const Real x1, const Real x2, const Real x3) |
Gravitational potential. | |
static Real | grav_acc (const Real x1, const Real x2, const Real x3) |
Gravitational acceleration. | |
void | cylbr_ix1 (GridS *pG) |
Inner-R boundary conditions. | |
void | cylbr_ox1 (GridS *pG) |
B_R = B_0/R boundary conditions, Outer x1 boundary. | |
void | problem (DomainS *pDomain) |
Problem initial conditions, edited for each problem. | |
void | problem_write_restart (MeshS *pM, FILE *fp) |
Problem specific subroutine for use with writing restart files, Edited for each problem. | |
void | problem_read_restart (MeshS *pM, FILE *fp) |
Problem specific subroutine for use with reading restart files, Edited for each problem. | |
ConsFun_t | get_usr_expr (const char *expr) |
Sets pointer to expression for special output data. | |
VOutFun_t | get_usr_out_fun (const char *name) |
Returns a user defined output function pointer. | |
void | Userwork_in_loop (MeshS *pM) |
Problem work executed during each step in the main Athena loop, see main.c. | |
void | Userwork_after_loop (MeshS *pM) |
Problem work executed after the main Athena loop, see main.c Edited for each problem. | |
Variables | |
static Real | br0 |
static Real | omega0 |
static Real | vz0 |
static Real | rho0 |
static Real | pgas0 |
static Real | a |
static Real | x2min |
static Real | x2max |
static Real | phi0 |
static Real | x1save |
static Real | x3save |
static ConsS *** | RootSoln = NULL |
A dynamic test of force balance using a B_R-only, time-dependent, non-axisymmetric magnetic field.
Definition in file cylbr.c.
1-component of B field
Definition at line 190 of file cylbr.c.
Referenced by B1i(), cylbr_ix1(), cylbr_ox1(), and problem().
1-component of B field at interface
Definition at line 196 of file cylbr.c.
References B1(), x1save, and x3save.
Referenced by cylbr_ix1(), cylbr_ox1(), and problem().
void cylbr_ix1 | ( | GridS * | pG | ) |
Inner-R boundary conditions.
d, M2, B1, B1i, and P are all functions of R, phi, and t.
Definition at line 212 of file cylbr.c.
References a, avg2d(), B1(), ConsS::B1c, B1i(), GridS::B1i, ConsS::B2c, ConsS::B3c, cc_pos(), d(), ConsS::d, GridS::dx1, GridS::dx2, ConsS::E, Gamma_1, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, ConsS::M1, M2(), ConsS::M2, ConsS::M3, omega0, Pgas(), phi0, qsimp(), GridS::time, GridS::U, vz0, x1save, x2min, and x3save.
Referenced by problem().
void cylbr_ox1 | ( | GridS * | pG | ) |
B_R = B_0/R boundary conditions, Outer x1 boundary.
Definition at line 253 of file cylbr.c.
References a, avg2d(), B1(), ConsS::B1c, B1i(), GridS::B1i, ConsS::B2c, ConsS::B3c, cc_pos(), d(), ConsS::d, GridS::dx1, GridS::dx2, ConsS::E, Gamma_1, GridS::ie, GridS::je, GridS::js, GridS::ke, GridS::ks, ConsS::M1, M2(), ConsS::M2, ConsS::M3, omega0, Pgas(), phi0, qsimp(), GridS::time, GridS::U, vz0, x1save, x2min, and x3save.
Referenced by problem().
ConsFun_t get_usr_expr | ( | const char * | expr | ) |
VOutFun_t get_usr_out_fun | ( | const char * | name | ) |
gas pressure
Definition at line 202 of file cylbr.c.
References a, br0, pgas0, phi0, and rho0.
Referenced by cylbr_ix1(), cylbr_ox1(), problem(), and ScaleToBeta().
void problem | ( | DomainS * | pDomain | ) |
Problem initial conditions, edited for each problem.
Definition at line 46 of file cylbr.c.
References a, ath_error(), avg2d(), B1(), ConsS::B1c, B1i(), GridS::B1i, ConsS::B2c, ConsS::B3c, br0, bvals_mhd_fun(), calloc_3d_array(), cylbr_ix1(), cylbr_ox1(), d(), ConsS::d, GridS::dx1, GridS::dx2, ConsS::E, Gamma_1, DomainS::Grid, GridS::ie, il, GridS::is, iu, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, left_x1, ConsS::M1, M2(), ConsS::M2, ConsS::M3, nx1, nx2, nx3, omega0, par_getd(), Pgas(), pgas0, phi0, qsimp(), rho0, right_x1, StaticGravPot, GridS::time, GridS::U, vz0, x1GravAcc, x1save, x2max, x2min, and x3save.
void problem_read_restart | ( | MeshS * | pM, | |
FILE * | fp | |||
) |
void problem_write_restart | ( | MeshS * | pM, | |
FILE * | fp | |||
) |
void Userwork_after_loop | ( | MeshS * | pM | ) |
Problem work executed after the main Athena loop, see main.c Edited for each problem.
Definition at line 157 of file cylbr.c.
References compute_l1_error().
void Userwork_in_loop | ( | MeshS * | pM | ) |
Problem work executed during each step in the main Athena loop, see main.c.
Edited for each problem.
Definition at line 152 of file cylbr.c.
References compute_div_b(), MeshS::Domain, and DomainS::Grid.
Definition at line 31 of file cylbr.c.
Referenced by B1(), bisection(), cylbr_ix1(), cylbr_ox1(), d(), esys_prim_adb_hyd(), esys_prim_adb_mhd(), esys_roe_adb_hyd(), fluxes(), get_ts_general(), getVlim_RS(), getwei_linear(), getwei_QP(), getwei_TSC(), M2(), Mollifier(), PFunc(), PFuncDeriv(), Pgas(), problem(), and property_mybin().
Definition at line 31 of file cylbr.c.
Referenced by B1(), grav_acc(), grav_pot(), Pgas(), and problem().
Definition at line 31 of file cylbr.c.
Referenced by cylbr_ix1(), cylbr_ox1(), grav_acc(), grav_pot(), M2(), and problem().
Definition at line 31 of file cylbr.c.
Referenced by B1(), cylbr_ix1(), cylbr_ox1(), d(), M2(), Pgas(), and problem().
Definition at line 31 of file cylbr.c.
Referenced by cylbr_ix1(), cylbr_ox1(), and problem().
Definition at line 31 of file cylbr.c.
Referenced by B1i(), cylbr_ix1(), cylbr_ox1(), and problem().
Definition at line 31 of file cylbr.c.
Referenced by B1i(), cylbr_ix1(), cylbr_ox1(), and problem().