Problem generator for non-linear streaming instability in stratified disks. More...

Functions
double	ran2 (long int *idum)
double	Normal (long int *idum)
	Normal distribution random number generator.
Real	Erf (Real z)
	Error function.
void	MultiNSH (int n, Real tstop, Real mratio, Real etavk, Real uxNSH, Real uyNSH, Real wxNSH, Real wyNSH)
	Multi-species NSH equilibrium.
static Real	hst_rho_Vx_dVy (const Grid *pG, const int i, const int j, const int k)
static void	close_ix2 (Grid *pGrid)
static void	close_ox2 (Grid *pGrid)
static Real	ShearingBoxPot (const Real x1, const Real x2, const Real x3)
static int	property_limit (const Grain gr, const GrainAux grsub)
	User defined particle selection function (1: true; 0: false).
static int	property_trace (const Grain gr, const GrainAux grsub)
	User defined particle selection function (1: true; 0: false).
static int	property_type (const Grain gr, const GrainAux grsub)
	User defined particle selection function (1: true; 0: false).
static int	property_dense (const Grain gr, const GrainAux grsub)
	User defined particle selection function (1: true; 0: false).
Real	expr_dpar (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle density.
Real	expr_V1par (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle 1-velocity.
Real	expr_V2par (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle 2-velocity.
Real	expr_V3par (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle 3-velocity.
Real	expr_V2 (const Grid *pG, const int i, const int j, const int k)
void	problem (Grid pGrid, Domain pDomain)
void	problem_write_restart (Grid pG, Domain pD, FILE *fp)
void	problem_read_restart (Grid pG, Domain pD, FILE *fp)
Gasfun_t	get_usr_expr (const char *expr)
	Sets pointer to expression for special output data.
VGFunout_t	get_usr_out_fun (const char *name)
	Returns a user defined output function pointer.
PropFun_t	get_usr_par_prop (const char *name)
	Returns a user defined particle selection function.
void	gasvshift (const Real x1, const Real x2, const Real x3, Real u1, Real u2, Real *u3)
	Gas velocity shift.
void	Userforce_particle (Vector *ft, const Real x1, const Real x2, const Real x3, const Real v1, const Real v2, const Real v3)
void	Userwork_in_loop (Grid pGrid, Domain pDomain)
void	Userwork_after_loop (Grid pGrid, Domain pDomain)
Variables
Real	vsc1
Real	vsc2
int	ipert
Real	x1min
Real	x1max
Real	x2min
Real	x2max
Real	Lx
Real	Lz
Real	Lg
long	Npar
long	ntrack
long	nlis
int	mytype
Real	dpar_thresh

Detailed Description

Problem generator for non-linear streaming instability in stratified disks.

PURPOSE: Problem generator for non-linear streaming instability in stratified disks. This code works in 2D ONLY. Isothermal eos is assumed, and the value etavk/iso_sound is fixed. MPI domain decomposition in x is allowed, but not in z.

Perturbation modes:

ipert = 0: multi-nsh equilibtium
ipert = 1: white noise within the entire grid
ipert = 2: non-nsh velocity

Must be configured using --enable-shearing-box and --with-eos=isothermal. FARGO is recommended.

Reference: Johansen & Youdin, 2007, ApJ, 662, 627 Bai & Stone, 2009, in preparation

Definition in file par_strat2d.c.

Function Documentation

static void close_ix2 ( Grid * pGrid ) [static]

Definition at line 418 of file par_strat2d.c.

References il, and iu.

Referenced by problem().

Here is the caller graph for this function:

static void close_ox2 ( Grid * pGrid ) [static]

Definition at line 437 of file par_strat2d.c.

References il, and iu.

Referenced by problem().

Here is the caller graph for this function:

Real Erf ( Real z )

Error function.

Definition at line 736 of file par_strat2d.c.

References il, and iu.

Real expr_dpar	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)

Wrapper for particle density.

Definition at line 291 of file output_particle.c.

Real expr_V1par	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)

Wrapper for particle 1-velocity.

Definition at line 312 of file output_particle.c.

Real expr_V2	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)

Real expr_V2par	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)

Wrapper for particle 2-velocity.

Definition at line 319 of file output_particle.c.

Real expr_V3par	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)

Wrapper for particle 3-velocity.

Definition at line 326 of file output_particle.c.

void gasvshift	(	const Real	x1,
		const Real	x2,
		const Real	x3,
		Real *	u1,
		Real *	u2,
		Real *	u3
	)

Gas velocity shift.

Definition at line 475 of file par_strat2d.c.

Gasfun_t get_usr_expr ( const char * expr )

Sets pointer to expression for special output data.

Definition at line 456 of file par_strat2d.c.

VGFunout_t get_usr_out_fun ( const char * name )

Returns a user defined output function pointer.

Definition at line 461 of file par_strat2d.c.

PropFun_t get_usr_par_prop ( const char * name )

Returns a user defined particle selection function.

Definition at line 466 of file par_strat2d.c.

static Real hst_rho_Vx_dVy	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)			`[static]`

Definition at line 405 of file par_strat2d.c.

References cc_pos(), Omega_0, and qshear.

Referenced by problem(), and problem_read_restart().

Here is the call graph for this function:

Here is the caller graph for this function:

void MultiNSH	(	int	n,
		Real *	tstop,
		Real *	mratio,
		Real	etavk,
		Real *	uxNSH,
		Real *	uyNSH,
		Real *	wxNSH,
		Real *	wyNSH
	)

Multi-species NSH equilibrium.

Input: # of particle types (n), dust stopping time and mass ratio array, and drift speed etavk. Output: gas NSH equlibrium velocity (u), and dust NSH equilibrium velocity array (w).

Definition at line 599 of file par_strat2d.c.

References calloc_1d_array(), calloc_2d_array(), free_2d_array(), InverseMatrix(), and MatrixMult().

Referenced by problem().

Here is the call graph for this function: