Compute the right hand side of the conservative HD/MHD equations. More...
#include "pluto.h"Functions | |
| static void | TotalFlux (const State_1D *, double *, int, int, Grid *) |
| void | RightHandSide (const State_1D *state, Time_Step *Dts, int beg, int end, double dt, Grid *grid) |
This function constructs the one-dimensional right hand side of the conservative MHD or HD equations in the direction given by g_dir in different geometries. The right hand side is computed as a two-point flux difference term plus a source term:
where
: interface areas
: cell volume
: interface fluxes
: time step
: source term including geometrical terms and body forces.See also RHS_page The right hand side is assembled through the following steps:
Compute the total flux in order to enforce conservation of angular momentum and energy in presence of FARGO source terms, rotation or gravitational potential.
1.8.2