Athena: prob/carbuncle.c File Reference

Problem generator for carbuncle instability. More...

Go to the source code of this file.

Functions
static void	initialize_states (void)
	Uses Rankine Hugoniot relations for adiabatic gas to shock jump conditions.
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	Mach
static Real	dl
static Real	pl
static Real	ul
static Real	dr
static Real	pr
static Real	ur

Detailed Description

Problem generator for carbuncle instability.

PURPOSE: Problem generator for carbuncle instability. Sets up a planar shock propagating in the x1-direction with arbitrary Mach number input from problem file. Adds perturbations to the transverse velocity of arbitrary amplitude A in a strip of one zone ahead of shock. If perturbation amplitude is zero code keeps shock exactly planar. If A/Cs = 10^{-4}, shock completely disintegrates without H-correction. H-correction completely fixes problem. Runs two problems: shk_flag = 0 - standing shock in middle of grid (obc_x1=2 in input file) shk_flag = 1 - flow at Ux=Mach into wall (obc_x1=1 in input file)

PRIVATE FUNCTION PROTOTYPES:

initialize_states() - sets shock jumps given Mach number

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

Definition in file carbuncle.c.

Function Documentation

ConsFun_t get_usr_expr ( const char * expr )

Sets pointer to expression for special output data.

Definition at line 158 of file carbuncle.c.

VOutFun_t get_usr_out_fun ( const char * name )

Returns a user defined output function pointer.

Definition at line 163 of file carbuncle.c.

static void initialize_states ( void ) [static]

Uses Rankine Hugoniot relations for adiabatic gas to shock jump conditions.

Definition at line 183 of file carbuncle.c.

References dl, dr, Gamma, Gamma_1, Mach, pl, pr, ul, and ur.

Referenced by problem().

Here is the caller graph for this function:

void problem ( DomainS * pDomain )

Problem initial conditions, edited for each problem.

Definition at line 46 of file carbuncle.c.

References amp, cc_pos(), ConsS::d, dl, dr, ConsS::E, Gamma, Gamma_1, DomainS::Grid, GridS::ie, initialize_states(), GridS::is, Iso_csound, GridS::je, GridS::js, GridS::ke, GridS::ks, ConsS::M1, ConsS::M2, ConsS::M3, Mach, par_getd(), pl, pr, DomainS::RootMaxX, DomainS::RootMinX, 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 153 of file carbuncle.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 148 of file carbuncle.c.

void Userwork_after_loop ( MeshS * pM )

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

Definition at line 171 of file carbuncle.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 167 of file carbuncle.c.

Variable Documentation

Real dl [static]

Definition at line 33 of file carbuncle.c.

Referenced by fluxes(), initialize_states(), and problem().

Real dr [static]

Definition at line 33 of file carbuncle.c.

Referenced by fluxes(), initialize_states(), problem(), and Userwork_in_loop().

Real Mach [static]

Definition at line 33 of file carbuncle.c.

Referenced by initialize_states(), and problem().

Real pl [static]

Definition at line 33 of file carbuncle.c.

Referenced by initialize_states(), and problem().

Real pr [static]

Definition at line 33 of file carbuncle.c.

Referenced by initialize_states(), and problem().

Real ul [static]

Definition at line 33 of file carbuncle.c.

Referenced by initialize_states(), and problem().

Real ur [static]

Definition at line 33 of file carbuncle.c.

Referenced by initialize_states(), and problem().