PLUTO
4.4-patch2
|
Compute the outermost wave speeds for HLL-based solvers. More...
#include "pluto.h"
Functions | |
void | HLL_Speed (const State *stateL, const State *stateR, double *SL, double *SR, int beg, int end) |
HLL_Speed() computes an estimate to the leftmost and rightmost wave signal speeds bounding the Riemann fan based on the input states ::stateL->v, ::stateR->v. Depending on the estimate, several variants are possible.
void HLL_Speed | ( | const State * | stateL, |
const State * | stateR, | ||
double * | SL, | ||
double * | SR, | ||
int | beg, | ||
int | end | ||
) |
Compute leftmost (SL) and rightmost (SR) speeds for the Riemann fan.
[in] | stateL | pointer to a state structure for the left state |
[in] | stateR | pointer to a state structure for the right state |
[out] | SL | the (estimated) leftmost speed of the Riemann fan |
[out] | SR | the (estimated) rightmost speed of the Riemann fan |
[in] | beg | starting index of computation |
[in] | end | final index of computation |
Switches:
ROE_ESTIMATE (YES/NO), DAVIS_ESTIMATE (YES/NO). TVD_ESTIMATE (YES/NO) JAN_HLL (YES/NO)
These switches set how the wave speed estimates are going to be computed. Only one can be set to 'YES', and the rest of them must be set to 'NO'
ROE_ESTIMATE: b_m = (0, (u_R - c_R, u_L - c_L, u_{roe} - c_{roe})) b_m = (0, (u_R + c_R, u_L + c_L, u_{roe} + c_{roe}))
where u_{roe} and c_{roe} are computed using Roe averages.
DAVIS_ESTIMATE: b_m = (0, (u_R - c_R, u_L - c_L)) b_m = (0, (u_R + c_R, u_L + c_L))