PLUTO
4.4-patch2
|
Compute time-centered interface states using characteristic tracing. More...
#include "pluto.h"
Macros | |
#define | CHTR_REF_STATE 2 |
Functions | |
void | CharTracingStep (const Sweep *sweep, int beg, int end, Grid *grid) |
Advance 1-D left and right interface states, previously computed with any of the States functions, to the half time level (n+1/2) by extrapolating characteristic variables. This is done using an upwind selection rule that discards waves not reaching the interface in dt/2.
References:
#define CHTR_REF_STATE 2 |
Flag to control the choice of the reference state in the upwind selection rule. Set CHTR_REF_STATE to 1,2,3 to use
Compute interface states using characteristic tracing step.
[in] | sweep | pointer to a Sweep structure |
[in] | beg | initial index of computation |
[in] | end | final index of computation |
[in] | grid | pointer to Grid structure |
Tasks are numbered below.
2a) Define characteristic cfl coefficients.
2b) Obtain characteristic variable increments dwp and dwm.
2c) Initialize vp and vm to the reference state. Since this is somewhat arbitrary we use the value of CHTR_REF_STATE to select one of the following cases:
Passive scalars use always CHTR_REF_STATE == 2.
2d) Compute left and right states in primitive variables. This step also depends on the value of CHTR_REF_STATE and include:
2e) Add source term to L/R states
2f) Repeat construction for passive scalars