Athena: prob/shk_cloud.c File Reference

Problem generator for shock-cloud problem; planar shock impacting a single spherical cloud. More...

Go to the source code of this file.

Functions
void	shk_cloud_iib (GridS *pGrid)
	Sets boundary condition on left X boundary (iib).
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.
static Real	color (const GridS *pG, const int i, const int j, const int k)
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	dl
static Real	pl
static Real	ul
static Real	bxl
static Real	byl
static Real	bzl

Detailed Description

Problem generator for shock-cloud problem; planar shock impacting a single spherical cloud.

PURPOSE: Problem generator for shock-cloud problem; planar shock impacting a single spherical cloud. Input parameters are:

problem/Mach = Mach number of incident shock
problem/drat = density ratio of cloud to ambient
problem/beta = ratio of Pgas/Pmag
problem/iprob = integer flag to determine problem

The cloud radius is fixed at 1.0. The center of the coordinate system defines the center of the cloud, and should be in the middle of the cloud. The shock is initially at x1=-2.0. A typical grid domain should span x1 in [-3.0,7.0] , y and z in [-2.5,2.5] (see input file in /tst) Various test cases are possible:

(iprob=1): B parallel to shock normal
(iprob=2): B perpendicular to shock normal -- NOT YET IMPLEMENTED

If the code is configured with nscalars>0, the cloud material is labeled with U[k][j][i].s[0]=1.

PRIVATE FUNCTION PROTOTYPES: shk_cloud_iib() - fixes BCs on L-x1 (left edge) of grid to postshock flow.

Definition in file shk_cloud.c.

Function Documentation

static Real color	(	const GridS *	pG,
		const int	i,
		const int	j,
		const int	k
	)			`[static]`

Definition at line 250 of file shk_cloud.c.

References ConsS::d, ConsS::s, and GridS::U.

ConsFun_t get_usr_expr ( const char * expr )

Sets pointer to expression for special output data.

Definition at line 256 of file shk_cloud.c.

VOutFun_t get_usr_out_fun ( const char * name )

Returns a user defined output function pointer.

Definition at line 265 of file shk_cloud.c.

void problem ( DomainS * pDomain )

Problem initial conditions, edited for each problem.

Definition at line 54 of file shk_cloud.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, beta, bvals_mhd_fun(), bxl, byl, bzl, cc_pos(), ConsS::d, DomainS::Disp, dl, dr, ConsS::E, Gamma, Gamma_1, DomainS::Grid, GridS::ie, iprob, GridS::is, Iso_csound, GridS::je, GridS::js, GridS::ke, GridS::ks, left_x1, ConsS::M1, ConsS::M2, ConsS::M3, Mach, par_getd(), par_geti(), pl, pr, rad, ConsS::s, shk_cloud_iib(), GridS::U, ul, and ur.

Here is the call graph for this function:

void problem_read_restart	(	MeshS *	pM,
		FILE *	fp
	)

Problem specific subroutine for use with reading restart files, Edited for each problem.

Definition at line 244 of file shk_cloud.c.

void problem_write_restart	(	MeshS *	pM,
		FILE *	fp
	)

Problem specific subroutine for use with writing restart files, Edited for each problem.

Definition at line 239 of file shk_cloud.c.

void shk_cloud_iib ( GridS * pGrid )

Sets boundary condition on left X boundary (iib).

Note quantities at this boundary are held fixed at the downstream state

Definition at line 288 of file shk_cloud.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, bxl, byl, bzl, ConsS::d, dl, ConsS::E, Gamma_1, GridS::je, GridS::js, GridS::ke, GridS::ks, ConsS::M1, ConsS::M2, ConsS::M3, pl, ConsS::s, GridS::U, and ul.

Referenced by problem().

Here is the caller graph for this function:

void Userwork_after_loop ( MeshS * pM )

Problem work executed after the main Athena loop, see main.c Edited for each problem.

Definition at line 274 of file shk_cloud.c.

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 269 of file shk_cloud.c.

Variable Documentation

Real bxl [static]

Definition at line 40 of file shk_cloud.c.

Referenced by problem(), and shk_cloud_iib().

Real byl [static]

Definition at line 40 of file shk_cloud.c.

Referenced by problem(), and shk_cloud_iib().

Real bzl [static]

Definition at line 40 of file shk_cloud.c.

Referenced by problem(), and shk_cloud_iib().

Real dl [static]

Definition at line 38 of file shk_cloud.c.

Referenced by problem(), and shk_cloud_iib().

Real pl [static]

Definition at line 38 of file shk_cloud.c.

Referenced by problem(), and shk_cloud_iib().

Real ul [static]

Definition at line 38 of file shk_cloud.c.

Referenced by problem(), and shk_cloud_iib().