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
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.