PLUTO
4.4-patch2
|
Compute source terms for Powell formulation. More...
#include "pluto.h"
Functions | |
void | Roe_DivBSource (const Sweep *sweep, int beg, int end, Grid *grid) |
void | HLL_DivBSource (const Sweep *sweep, double **Uhll, int beg, int end, Grid *grid) |
This file contains two implementations of Powell's source term in the 8-wave formulation for MHD:
where B
is the magnetic field in the normal direction (g_dir), Q
is a fluid quantity, A
is the area and the denominator is the cell volume. The first function, Roe_DivBSource() is called by Roe_Solver() and TVDLF_Solver() and computes the normal component using arithmetic average of the left and right states: . The second implementation contained in HLL_DivBSource() computes the term using upwinding:
Reference:
Include div.B source term to momentum, induction and energy equation. Used in conjunction with an HLL-type Riemann solver.