PLUTO
4.4-patch2
|
Computes viscous fluxes and source terms for the HD/MHD equations. More...
#include "pluto.h"
Functions | |
void | ViscousFlux (const Data *d, double **ViF, double **ViS, double *dcoeff, int beg, int end, Grid *grid) |
Compute the stress tensor components (at the i+1/2 face of each cell) and adds explicit viscous terms to the energy and momentum equation. It is called in the during the sweep integrators. The stress tensor is given by
where is the dynamic viscosity and
is the unit tensor.
In order to compute the viscous stress tensor, we report here the expression for the gradient of a vector in cylindrical and spherical coordinates (see the book "I do like CFD, vol. I" by K. Masatsuka, Sec. 1.5.5, Eq. 1.5.35):
References
void ViscousFlux | ( | const Data * | d, |
double ** | ViF, | ||
double ** | ViS, | ||
double * | dcoeff, | ||
int | beg, | ||
int | end, | ||
Grid * | grid | ||
) |
[in] | V | data array containing cell-centered quantities |
[in,out] | ViF | pointer to viscous fluxes |
[in,out] | ViS | pointer to viscous source terms |
[in,out] | dcoeff | pointer to diffusion coefficient for dt calculation |
[in] | beg | integer, index for loop beg |
[in] | end | integer, index for loop end |
[in] | grid | pointer to array of Grid structures |