PLUTO
4.4-patch2
|
Shock finding algorithm. More...
#include "pluto.h"
Functions | |
void | FlagShock (const Data *d, Grid *grid) |
Search and flag computational zones lying in a shock wave. The flagging strategy is based on two switches designed to detect the presence of compressive motion or shock waves in the fluid:
where sets the shock strength. At the discrete level we replace the two conditions by
where is a vector of integer numbers giving the position of a computational zone, while
is a unit vector in the direction given by
d
. Once a zone has been tagged as lying in a shock, different flags may be switched on or off to control the update strategy in these critical regions.
This function can be called called when:
SHOCK_FLATTENING
has been set to MULTID:
in this case shocked zones are tagged with FLAG_MINMOD
and FLAG_HLL
that will later be used to force the reconstruction with the minmod limiter and the Riemann solver to HLL.ENTROPY_SWITCH
has been turned on: this flag will be checked later in the ConsToPrim() functions in order to recover pressure from the entropy density rather than from the total energy density. The update process is: {V,s}
where s
is the entropy;{V}
;{s}
from {V}
; {U,S}
to primitive: Reference