microphysics/viscosity.c File Reference

Adds explicit viscosity terms to the momentum and energy equations. More...

Go to the source code of this file.

Classes
struct	ViscFluxS
	Structure to contain 4-components of the viscous fluxes. More...
Functions
void	ViscStress_iso (DomainS *pD)
	Calculate viscous stresses with isotropic (NS) viscosity.
void	ViscStress_aniso (DomainS *pD)
	Calculate viscous stresses with anisotropic (Braginskii) viscosity.
void	viscosity (DomainS *pD)
	Adds explicit viscosity terms to the momentum and energy equations.
void	viscosity_init (MeshS *pM)
	Allocate temporary arrays.
void	viscosity_destruct (void)
	Free temporary arrays.
Variables
static ViscFluxS ***	x1Flux = NULL
static ViscFluxS ***	x2Flux = NULL
static ViscFluxS ***	x3Flux = NULL
static Real3Vect ***	Vel = NULL
static Real ***	divv = NULL

---

Detailed Description

Adds explicit viscosity terms to the momentum and energy equations.

PURPOSE: Adds explicit viscosity terms to the momentum and energy equations,

• dM/dt = Div(T)
• dE/dt = Div(v.T) where
• T = nu_iso Grad(V) + T_Brag = TOTAl viscous stress tensor

Note T contains contributions from both isotropic (Navier-Stokes) and anisotropic (Braginskii) viscosity. These contributions are computed in calls to ViscStress_* functions.

CONTAINS PUBLIC FUNCTIONS:

• viscosity() - updates momentum and energy equations with viscous terms
• viscosity_init() - allocates memory needed
• viscosity_destruct() - frees memory used

Definition in file viscosity.c.

---

Function Documentation

void viscosity

(

DomainS *

pD

)

Adds explicit viscosity terms to the momentum and energy equations.

Definition at line 62 of file viscosity.c.

References cc_pos(), ConsS::d, divv, GridS::dt, GridS::dx1, GridS::dx2, GridS::dx3, ConsS::E, Cons1DS::E, ViscFluxS::E, DomainS::Grid, GridS::ie, GridS::is, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, ConsS::M1, ConsS::M2, ConsS::M3, Cons1DS::Mx, ViscFluxS::Mx, Cons1DS::My, ViscFluxS::My, Cons1DS::Mz, ViscFluxS::Mz, nu_aniso, nu_iso, GridS::Nx, Omega_0, qshear, GridS::U, ViscStress_aniso(), ViscStress_iso(), Real3Vect::x, x2Flux, x3Flux, 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 viscosity_destruct

(

void

)

Free temporary arrays.

Definition at line 812 of file viscosity.c.

References divv, free_3d_array(), x2Flux, and x3Flux.

Referenced by integrate_diff_destruct(), and viscosity_init().

Here is the call graph for this function:

Here is the caller graph for this function:

void viscosity_init

(

MeshS *

pM

)

Allocate temporary arrays.

Definition at line 754 of file viscosity.c.

References ath_error(), calloc_3d_array(), divv, MeshS::Domain, MeshS::DomainsPerLevel, DomainS::Grid, MeshS::NLevels, MeshS::Nx, GridS::Nx, viscosity_destruct(), x2Flux, and x3Flux.

Referenced by integrate_diff_init().

Here is the call graph for this function:

Here is the caller graph for this function:

void ViscStress_aniso

(

DomainS *

pD

)

Calculate viscous stresses with anisotropic (Braginskii) viscosity.

Definition at line 326 of file viscosity.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, ConsS::d, GridS::dx1, GridS::dx2, GridS::dx3, Cons1DS::E, ViscFluxS::E, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, Cons1DS::Mx, ViscFluxS::Mx, Cons1DS::My, ViscFluxS::My, Cons1DS::Mz, ViscFluxS::Mz, nu_aniso, DomainS::Nx, GridS::U, Real3Vect::x, x2Flux, x3Flux, Real3Vect::y, and Real3Vect::z.

Referenced by viscosity().

Here is the caller graph for this function:

void ViscStress_iso

(

DomainS *

pD

)

Calculate viscous stresses with isotropic (NS) viscosity.

Definition at line 209 of file viscosity.c.

References ConsS::d, divv, GridS::dx1, GridS::dx2, GridS::dx3, Cons1DS::E, ViscFluxS::E, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, Cons1DS::Mx, ViscFluxS::Mx, Cons1DS::My, ViscFluxS::My, Cons1DS::Mz, ViscFluxS::Mz, nu_iso, GridS::Nx, GridS::U, Real3Vect::x, x2Flux, x3Flux, Real3Vect::y, and Real3Vect::z.

Referenced by viscosity().

Here is the caller graph for this function:

---

Variable Documentation

Real*** divv = NULL [static]

Definition at line 45 of file viscosity.c.

Referenced by viscosity(), viscosity_destruct(), viscosity_init(), and ViscStress_iso().

Real3Vect*** Vel = NULL [static]

Definition at line 44 of file viscosity.c.

ViscFluxS*** x1Flux = NULL [static]

Definition at line 43 of file viscosity.c.

ViscFluxS *** x2Flux = NULL [static]

Definition at line 43 of file viscosity.c.

ViscFluxS *** x3Flux = NULL [static]

Definition at line 43 of file viscosity.c.