PLUTO
4.4-patch2
|
Compute the resistive MHD flux for cell-centered field. More...
#include "pluto.h"
Functions | |
void | ResistiveFlux (Data_Arr V, Data_Arr curlB, double **res_flux, double **dcoeff, int beg, int end, Grid *grid) |
Compute the resistive fluxes for the induction and energy equations to update cell-centered fields. Fluxes at zone faces by averaging the current available at cell edges. Cell-centered fields are updated also when the CT formalism is employed.
In the induction equation, fluxes are computed by explicitly writing the curl operator in components. In Cartesian components, for instance, one has
where is the resistive diagonal tensor and
is the current density. The corresponding contribution to the energy equation is
The sign of the flux is given by writing each equation as
(flux is positive on the right hand side).
References
void ResistiveFlux | ( | Data_Arr | V, |
Data_Arr | curlB, | ||
double ** | res_flux, | ||
double ** | dcoeff, | ||
int | beg, | ||
int | end, | ||
Grid * | grid | ||
) |
[in] | V | 3D data array of primitive variables |
[out] | curlB | the current |
[out] | dcoeff | the diffusion coefficients evaluated at cell interfaces |
[in] | beg | initial index of computation |
[in] | end | final index of computation |
[in] | grid | pointer to an array of Grid structures. |