PLUTO
4.4-patch2
|
HLLEM Riemann solver for RMHD. More...
#include "pluto.h"
Functions | |
void | HLLEM_Solver (const Sweep *sweep, int beg, int end, double *cmax, Grid *grid) |
Reference:
restart; B[x] := 0; B[y] := 0; B[z] := 0;
lambda := delta + v[x]: # We use Lambda = lambda-v[x] to have a simpler
delta := gamma*(lambda - v[x]): vB := v[x]*B[x] + v[y]*B[y] + v[z]*B[z]: b[x] := B[x]/gamma + gamma*v[x]*vB: b[0] := gamma*vB: B2 := B[x]^2 + B[y]^2 + B[z]^2: BB := b[x] - lambda*b[0]: b2 := B2/gamma^2 + vB^2: fMB := w*(1 - c[s]^2)*a^4 - (1 - lambda^2)*((b2 + w*c[s]^2)*a^2 - c[s]^2*BB^2): fMB := simplify(fMB); c[4] := coeff(fMB,delta,4); c[3] := coeff(fMB,delta,3); c[2] := coeff(fMB,delta,2); c[1] := coeff(fMB,delta,1); c[0] := coeff(fMB,delta,0);