PLUTO
4.4-patch2
|
Thermal conduction (TC) module header file. More...
Go to the source code of this file.
Functions | |
void | TC_Flux (double ***, const Sweep *, double *, int, int, Grid *) |
void | TC_kappa (double *, double, double, double, double *, double *, double *) |
void | TC_RHS (const Data *, Data_Arr, double *, double **, double, int, int, Grid *) |
void | GetGradient (double ***, double **, int, int, Grid *) |
Contains prototypes for the thermal conduction module.
void GetGradient | ( | double *** | T, |
double ** | gradT, | ||
int | beg, | ||
int | end, | ||
Grid * | grid | ||
) |
Compute the gradient of a 3D scalar quantity T in the direction given by g_dir. Return a 1D array (dT/dx, dT/dy, dT/dz) along that direction computed at cell interfaces, e.g.
if g_dir == IDIR –> compute
[ dT/dl1, dT/dl2, dT/dl3 ] at interface (i+1/2,j,k)
if g_dir == JDIR –> compute
[ dT/dl1, dT/dl2, dT/dl3 ] at interface (i,j+1/2,k)
if g_dir == KDIR –> compute
[ dT/dl1, dT/dl2, dT/dl3 ] at interface (i,j,k+1/2)
Here dl1, dl2 and dl3 are the line element in the thre directions:
Cartesian: {dl1, dl2, dl3} = {dx, dy, dz} Cylindrical: {dl1, dl2, dl3} = {dr, dz, - } Polar: {dl1, dl2, dl3} = {dr, r.dphi, dz} Spherical: {dl1, dl2, dl3} = {dr, r.dtheta, r.sin(theta).dphi}
LAST MODIFIED
14 Apr 2011 by T. Matsakos, A. Mignone
Compute the thermal conduction flux, sweep->par_flx.
[in] | T | 3D array containing the dimensionless temperature |
[in,out] | sweep | pointer to a Sweep structure |
[out] | dcoeff | the diffusion coefficient needed for computing the time step. |
[in] | beg | initial index of computation |
[in] | end | final index of computation |
[in] | grid | pointer to an array of Grid structures |
void TC_kappa | ( | double * | v, |
double | x1, | ||
double | x2, | ||
double | x3, | ||
double * | kpar, | ||
double * | knor, | ||
double * | phi | ||
) |
Compute thermal conduction coefficients.
[in] | v | array of primitive variables |
[in] | x1 | coordinate in the X1 direction |
[in] | x2 | coordinate in the X2 direction |
[in] | x3 | coordinate in the X3 direction |
[out] | kpar | pointer to the conduction coefficient ![]() |
[out] | knor | pointer to the conduction coefficient ![]() |
[out] | phi | pointer to the parameter ![]() |
void TC_RHS | ( | const Data * | d, |
Data_Arr | dU, | ||
double * | dcoeff, | ||
double ** | aflux, | ||
double | dt, | ||
int | beg, | ||
int | end, | ||
Grid * | grid | ||
) |
[in] | d | pointer to PLUTO Data structure |
[out] | dU | a 4D array containing conservative variables increment |
[out] | dcoeff | 1D array of diffusion coefficients |
[out] | aflux | pointer to 2D array for AMR re-fluxing operations |
[in] | dt | the current time-step |
[in] | beg,end | initial and final interface indices |
[in] | grid | pointer to Grid structure. |