PLUTO
4.4-patch2
|
Lax-Friedrechs (Rusanov) Riemann solver for HD. More...
#include "pluto.h"
Functions | |
void | LF_Solver (const Sweep *sweep, int beg, int end, double *cmax, Grid *grid) |
Solve the Riemann problem for the adiabatic and isothermal HD equations using the Lax-Friedrichs Rusanov Riemann solver with local maximum characteristic speed:
where is a function of the arithmetic average between the left and the right states.
On input, this function takes left and right primitive state vectors stateL->v
and stateR->v
at zone edge i+1/2; On output, return flux and pressure vectors at the same interface i+1/2
(note that the i
refers to i+1/2
).
Also during this step, compute maximum wave propagation speed (cmax) for explicit time step computation.
Reference: