Athena: microphysics/resistivity.c File Reference

Adds explicit resistivity terms to the induction and energy eqns. More...

Go to the source code of this file.

Functions
void	EField_Ohm (DomainS *pD)
	Resistive EMF from Ohmic dissipation.
void	EField_Hall (DomainS *pD)
	Resistive EMF from Hall effect.
void	EField_AD (DomainS *pD)
	Resistive EMF from ambipolar diffusion.
void	EField_Hall_sub (DomainS pD, Real3Vect Bs, Real3Vect Js, Real3Vect **emfs, int noff)
	Update the Hall emfs (subroutine for the Hall calculation).
void	hyper_diffusion6 (DomainS *pD)
	Calculate the higher-order derivatives of J -- 6th order diffusion.
void	resistivity (DomainS *pD)
	Explicit resistivity.
void	resistivity_init (MeshS *pM)
	Allocate temporary arrays.
void	resistivity_destruct ()
	Free temporary arrays.
Variables
static Real3Vect ***	J = NULL
static Real3Vect ***	emf = NULL
static Real3Vect ***	EnerFlux = NULL
static Real3Vect ***	Jcor = NULL
static Real3Vect ***	Bcor = NULL
static Real3Vect ***	emfp = NULL

Detailed Description

Adds explicit resistivity terms to the induction and energy eqns.

PURPOSE: Adds explicit resistivity terms to the induction and energy eqns,

dB/dt = -Curl(E)
dE/dt = Div(B X E) where
E = eta_Ohm J + eta_Hall(J X B)/B + eta_AD J_perp = (emf)
J = Curl(B) = current density
eta_Ohm = Ohmic resistivity
eta_Hall = Hall diffusion coefficient
eta_AD = ambipolar diffusion coefficient The induction equation is updated using CT to keep div(B)=0. The total electric field (resistive EMF) is computed from calls to the EField_* functions.

CONTAINS PUBLIC FUNCTIONS:

resistivity() - updates induction and energy eqns with resistive term.
resistivity_init() - allocates memory needed
resistivity_destruct() - frees memory used

Definition in file resistivity.c.

Function Documentation

void EField_AD ( DomainS * pD )

Resistive EMF from ambipolar diffusion.

E = Q_AD (J X B) X B

Definition at line 687 of file resistivity.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, emf, GridS::eta_AD, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, GridS::Nx, GridS::U, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by resistivity().

Here is the caller graph for this function:

void EField_Hall ( DomainS * pD )

Resistive EMF from Hall effect.

E = Q_H (J X B)

Definition at line 431 of file resistivity.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, Bcor, GridS::dt, GridS::dx1, GridS::dx2, GridS::dx3, EField_Hall_sub(), emf, emfp, GridS::eta_Hall, DomainS::Grid, hyper_diffusion6(), GridS::ie, il, GridS::is, iu, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, GridS::Nx, GridS::U, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by resistivity().

Here is the call graph for this function:

Here is the caller graph for this function:

void EField_Hall_sub	(	DomainS *	pD,
		Real3Vect ***	Bs,
		Real3Vect ***	Js,
		Real3Vect ***	emfs,
		int	noff
	)

Update the Hall emfs (subroutine for the Hall calculation).

Definition at line 565 of file resistivity.c.

References GridS::eta_Hall, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, GridS::Nx, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by EField_Hall().

Here is the caller graph for this function:

void EField_Ohm ( DomainS * pD )

Resistive EMF from Ohmic dissipation.

E = J

Definition at line 353 of file resistivity.c.

References emf, GridS::eta_Ohm, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, GridS::Nx, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by resistivity().

Here is the caller graph for this function:

void hyper_diffusion6 ( DomainS * pD )

Calculate the higher-order derivatives of J -- 6th order diffusion.

Definition at line 867 of file resistivity.c.

References GridS::dt, GridS::dx1, GridS::dx2, GridS::dx3, emf, GridS::eta_Hall, fac, DomainS::Grid, GridS::ie, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, GridS::Nx, Real3Vect::x, Real3Vect::y, and Real3Vect::z.

Referenced by EField_Hall().

Here is the caller graph for this function:

void resistivity ( DomainS * pD )

Explicit resistivity.

Definition at line 67 of file resistivity.c.

References ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, GridS::dt, GridS::dx1, GridS::dx2, GridS::dx3, ConsS::E, EField_AD(), EField_Hall(), EField_Ohm(), emf, emfp, EnerFlux, eta_Ohm, DomainS::Grid, GridS::ie, GridS::is, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, GridS::Nx, Q_AD, Q_Hall, 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 resistivity_destruct ( )

Free temporary arrays.

Definition at line 1073 of file resistivity.c.

References Bcor, emf, emfp, EnerFlux, free_3d_array(), and get_myeta.

Referenced by integrate_diff_destruct(), and resistivity_init().

Here is the call graph for this function:

Here is the caller graph for this function:

void resistivity_init ( MeshS * pM )

Allocate temporary arrays.

Definition at line 991 of file resistivity.c.

References ath_error(), Bcor, calloc_3d_array(), MeshS::Domain, MeshS::DomainsPerLevel, emf, emfp, EnerFlux, get_myeta, DomainS::Grid, MeshS::NLevels, MeshS::Nx, GridS::Nx, par_geti_def(), Q_Hall, and resistivity_destruct().

Referenced by integrate_diff_init().

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

Real3Vect *** Bcor = NULL [static]

Definition at line 43 of file resistivity.c.

Referenced by EField_Hall(), resistivity_destruct(), and resistivity_init().

Real3Vect *** emf = NULL [static]

Definition at line 40 of file resistivity.c.

Referenced by EField_AD(), EField_Hall(), EField_Ohm(), hyper_diffusion6(), resistivity(), resistivity_destruct(), and resistivity_init().

Real3Vect *** emfp = NULL [static]

Definition at line 43 of file resistivity.c.

Referenced by EField_Hall(), resistivity(), resistivity_destruct(), and resistivity_init().

Real3Vect *** EnerFlux = NULL [static]

Definition at line 40 of file resistivity.c.

Referenced by resistivity(), resistivity_destruct(), and resistivity_init().

Real3Vect*** J = NULL [static]

Definition at line 40 of file resistivity.c.

Referenced by getShockVars(), getVb_Shock(), getVRel_2S(), and getVRel_RS().

Real3Vect*** Jcor = NULL [static]

Definition at line 43 of file resistivity.c.