microphysics/conduction.c File Reference

Adds explicit thermal conduction term to the energy equation, dE/dt = Div(Q). More...

Go to the source code of this file.

Functions
void	HeatFlux_iso (DomainS *pD)
	Calculate heat fluxes with isotropic conduction.
void	HeatFlux_aniso (DomainS *pD)
	Calculate heat fluxes with anisotropic conduction.
void	conduction (DomainS *pD)
	Explicit thermal conduction.
void	conduction_init (MeshS *pM)
	Allocate temporary arrays.
void	conduction_destruct (void)
	Free temporary arrays.
Variables
static Real ***	Temp = NULL
static Real3Vect ***	Q = NULL

---

Detailed Description

Adds explicit thermal conduction term to the energy equation, dE/dt = Div(Q).

where

• Q = kappa_iso Grad(T) + kappa_aniso([b Dot Grad(T)]b) = heat flux
• T = (P/d)*(mbar/k_B) = temperature
• b = magnetic field unit vector

Here

• kappa_iso is the coeffcient for isotropic conduction
• kappa_aniso is the coeffcient for anisotropic conduction

The heat flux Q is calculated by calls to HeatFlux_* functions.

CONTAINS PUBLIC FUNCTIONS:

• conduction() - updates energy equation with thermal conduction
• conduction_init() - allocates memory needed
• conduction_destruct() - frees memory used

Definition in file conduction.c.

---

Function Documentation

void conduction

(

DomainS *

pD

)

Explicit thermal conduction.

Definition at line 56 of file conduction.c.

References ConsS::B1c, ConsS::B2c, ConsS::B3c, ConsS::d, GridS::dt, GridS::dx1, GridS::dx2, GridS::dx3, ConsS::E, Gamma_1, DomainS::Grid, HeatFlux_aniso(), HeatFlux_iso(), GridS::ie, GridS::is, GridS::je, jl, GridS::js, ju, kappa_aniso, kappa_iso, GridS::ke, kl, GridS::ks, ku, ConsS::M1, ConsS::M2, ConsS::M3, GridS::Nx, Temp, GridS::U, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by integrate_diff().

Here is the call graph for this function:

Here is the caller graph for this function:

void conduction_destruct

(

void

)

Free temporary arrays.

Definition at line 429 of file conduction.c.

References free_3d_array(), and Temp.

Referenced by conduction_init(), and integrate_diff_destruct().

Here is the call graph for this function:

Here is the caller graph for this function:

void conduction_init

(

MeshS *

pM

)

Allocate temporary arrays.

Definition at line 379 of file conduction.c.

References ath_error(), calloc_3d_array(), conduction_destruct(), MeshS::Domain, MeshS::DomainsPerLevel, DomainS::Grid, MeshS::NLevels, MeshS::Nx, GridS::Nx, and Temp.

Referenced by integrate_diff_init().

Here is the call graph for this function:

Here is the caller graph for this function:

void HeatFlux_aniso

(

DomainS *

pD

)

Calculate heat fluxes with anisotropic conduction.

Definition at line 195 of file conduction.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, GridS::dx1, GridS::dx2, GridS::dx3, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, kappa_aniso, GridS::ke, GridS::ks, DomainS::Nx, Temp, GridS::U, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by conduction().

Here is the caller graph for this function:

void HeatFlux_iso

(

DomainS *

pD

)

Calculate heat fluxes with isotropic conduction.

Definition at line 149 of file conduction.c.

References GridS::dx1, GridS::dx2, GridS::dx3, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, kappa_iso, GridS::ke, GridS::ks, GridS::Nx, Temp, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by conduction().

Here is the caller graph for this function:

---

Variable Documentation

Real3Vect*** Q = NULL [static]

Definition at line 40 of file conduction.c.

Referenced by getVChar_pluto().

Real*** Temp = NULL [static]

Definition at line 39 of file conduction.c.

Referenced by conduction(), conduction_destruct(), conduction_init(), HeatFlux_aniso(), and HeatFlux_iso().