Third order (piecewise parabolic) spatial reconstruction using characteristic interpolation in the primitive variables. More...

Functions
void	lr_states (const GridS *pG, const Prim1DS W[], const Real Bxc[], const Real dt, const Real dx, const int il, const int iu, Prim1DS Wl[], Prim1DS Wr[], const int dir)
	Computes L/R states Input Arguments: W = PRIMITIVE variables at cell centers along 1-D slice Bxc = B in direction of slice at cell center dtodx = dt/dx il,iu = lower and upper indices of zone centers in slice W and Bxc must be initialized over [il-2:iu+2].
void	lr_states_init (MeshS *pM)
	NoOp for first order, but included for compatibility with integrator (needed for 2nd and 3rd order).
void	lr_states_destruct (void)
	NoOp for first order, but included for compatibility with integrator (needed for 2nd and 3rd order).
Variables
static Real **	pW = NULL
static Real **	dWm = NULL
static Real **	Wim1h = NULL

Detailed Description

Third order (piecewise parabolic) spatial reconstruction using characteristic interpolation in the primitive variables.

PURPOSE: Third order (piecewise parabolic) spatial reconstruction using characteristic interpolation in the primitive variables. With the CTU integrator, a time-evolution (characteristic tracing) step is used to interpolate interface values to the half time level {n+1/2}

NOTATION:

W_{L,i-1/2} is reconstructed value on the left-side of interface at i-1/2
W_{R,i-1/2} is reconstructed value on the right-side of interface at i-1/2

The L- and R-states at the left-interface in each cell are indexed i.

W_{L,i-1/2} is denoted by Wl[i ]; W_{R,i-1/2} is denoted by Wr[i ]
W_{L,i+1/2} is denoted by Wl[i+1]; W_{R,i+1/2} is denoted by Wr[i+1]

Internally, in this routine, Wlv and Wrv are the reconstructed values on the left-and right-side of cell center. Thus (see Step 18),

W_{L,i-1/2} = Wrv(i-1); W_{R,i-1/2} = Wlv(i)

REFERENCE: P. Colella & P. Woodward, "The piecewise parabolic method (PPM) for gas-dynamical simulations", JCP, 54, 174 (1984).

CONTAINS PUBLIC FUNCTIONS:

lr_states() - computes L/R states
lr_states_init() - initializes memory for static global arrays
lr_states_destruct() - frees memory for static global arrays

Definition in file lr_states_ppm.c.

Function Documentation

void lr_states	(	const GridS *	pG,
		const Prim1DS	W[],
		const Real	Bxc[],
		const Real	dt,
		const Real	dx,
		const int	il,
		const int	iu,
		Prim1DS	Wl[],
		Prim1DS	Wr[],
		const int	dir
	)

Computes L/R states Input Arguments: W = PRIMITIVE variables at cell centers along 1-D slice Bxc = B in direction of slice at cell center dtodx = dt/dx il,iu = lower and upper indices of zone centers in slice W and Bxc must be initialized over [il-2:iu+2].

Computes L/R states Input Arguments:

W = PRIMITIVE variables at cell centers along 1-D slice
Bxc = B in direction of slice at cell center
dtodx = dt/dx
il,iu = lower and upper indices of zone centers in slice W and Bxc must be initialized over [il-2:iu+2].

Computes L/R states.

Output Arguments: Wl,Wr = L/R-states of PRIMITIVE variables at interfaces over [il:iu+1]

Input Arguments:

W = PRIMITIVE variables at cell centers along 1-D slice
Bxc = B in direction of slice at cell center
dtodx = dt/dx
il,iu = lower and upper indices of zone centers in slice W and Bxc must be initialized over [il-3:iu+3]

Output Arguments: