wxm::apps::mhd_two_temperature::imhd Namespace Reference

Classes
class	CylSource
	Implements the source terms in the Two Temperature Ideal MHD equations arising from a cylindrical geometry. More...
class	ElectricField
	Class that calculates the electric field from Ohm's Law using ideal terms E = -u x B. More...
class	Ideal_MHD
	Class used implement Ideal MHD for the Two-Temperature MHD model. More...
class	InjectionSource
	Plasma source injection source for Two-Temperature Ideal MHD. More...

Functions
real	calculate_max_wave_speed (const real *q, const real gas_gamma, const real *norm, const real rho_min, const real press_min)
	Calculate the maximum wave speed for Ideal MHD at an element face.

Function Documentation

calculate_max_wave_speed()

real wxm::apps::mhd_two_temperature::imhd::calculate_max_wave_speed	(	const real *	q,
		const real	gas_gamma,
		const real *	norm,
		const real	rho_min,
		const real	press_min
	)

Calculate the maximum wave speed for Ideal MHD at an element face.

\[ c = c_{f} + \left|u\right|, \]

where \( c_{f}\) is the fast magnetosonic wave speed and \( u \) is the x-directed fluid speed in the direction of the face normal.

\[ c_{f} = \sqrt\left[\frac{1}{2} \left(v_{A}^{2} + c_{s}^{2}\right) + \frac{1}{2}\sqrt\left\{\\left(v_{A}^{2} + c_{s}^{2}\right)^{2} -4 v_{Ax}^{2}c_{s}^{2}right\} \right] \]

where \( v_{A}^{2} = \frac{B^{2}}{\rho} \), \( c_{s}^{2} = \frac{\gamma p}{\rho} \), and \( v_{Ax}^{2} = \frac{\left(\boldsymbol{B}\cdot \boldsymbol{n}\right)^{2}}{\rho} \)

Parameters

q	The raw MHD variables.
gas_gamma	Ratio of specific heats.
norm	normal vector on element face.
rho_min	Normalized minimum mass density.
press_min	Normalized minimum pressure.

Returns

c maximum wavespeed