PLUTO  4.4-patch2
Functions
hllem.c File Reference

HLLEM Riemann solver for RMHD. More...

#include "pluto.h"

Functions

void HLLEM_Solver (const Sweep *sweep, int beg, int end, double *cmax, Grid *grid)
 

Detailed Description

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

form of the equations.

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);

Authors
G. Mattia A. Mignone
Date
Dec 30, 2020

Function Documentation

◆ HLLEM_Solver()

void HLLEM_Solver ( const Sweep sweep,
int  beg,
int  end,
double *  cmax,
Grid grid 
)

Solve Riemann problem for the RMHD equation using HLLEM solver.

Parameters
[in,out]sweeppointer to Sweep structure
[in]beginitial grid index
[out]endfinal grid index
[out]cmax1D array of maximum characteristic speeds
[in]gridpointer to array of Grid structures.