The cylindrical analogue of the Bondi accretion (Parker wind) problem. More...
Go to the source code of this file.
Functions | |
| Real | grav_pot (const Real x1, const Real x2, const Real x3) |
| Gravitational potential. | |
| Real | grav_acc (const Real x1, const Real x2, const Real x3) |
| Gravitational acceleration. | |
| Real | myfunc (const Real x, const Real v) |
| This funciton is used to calculate velocity v as a function of position x using lambda_c, the critical value of the dimensionless mass wind/accretion rate. | |
| 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 | b0 |
| static Real | lambda_s |
| static int | iprob |
| static ConsS *** | RootSoln = NULL |
Detailed Description
The cylindrical analogue of the Bondi accretion (Parker wind) problem.
Hydrodynamic, rotationless, axisymmetric.
- REFERENCE: F. Shu, "The Physics of Astrophysics, Vol. II: Gas Dynamics", 1992. ISBN 0935702652
- REFERENCE: L. Spitzer, "Physical Processes in the Interstellar Medium", 1998. ISBN 0471293350
Definition in file cylwind.c.
Function Documentation
| ConsFun_t get_usr_expr | ( | const char * | expr | ) |
| VOutFun_t get_usr_out_fun | ( | const char * | name | ) |
This funciton is used to calculate velocity v as a function of position x using lambda_c, the critical value of the dimensionless mass wind/accretion rate.
This function is used to calculate y (ie.
This function is used to calculate v as a function of x, gamma, and lambda using the bisection method.
Standard bisection is used to find the root(s).
rho) as a function of x (ie. R), gamma, eta, theta, omega, and E using the bisection method.
Definition at line 200 of file cylwind.c.
References Gamma_1, and lambda_s.
Referenced by problem().
| void problem | ( | DomainS * | pDomain | ) |
Problem initial conditions, edited for each problem.
Definition at line 41 of file cylwind.c.
References a, ath_error(), b0, ConsS::B1c, GridS::B1i, beta, bisection(), bvals_mhd_fun(), calloc_3d_array(), cc_pos(), ConsS::d, do_nothing_bc(), GridS::dx1, ConsS::E, Gamma, Gamma_1, DomainS::Grid, GridS::ie, il, iprob, GridS::is, iu, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, lambda_s, left_x1, ConsS::M1, myfunc(), nx1, nx2, nx3, par_getd(), par_geti(), pgas, pgas0, right_x1, Soln, StaticGravPot, GridS::U, and x1GravAcc.
| 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 173 of file cylwind.c.
References compute_l1_error().
| void Userwork_in_loop | ( | MeshS * | pM | ) |
Variable Documentation
Definition at line 31 of file cylwind.c.
Referenced by getVChar_echo(), getVChar_pluto(), int_par_fulimp(), and problem().
