PlasmaProperties.jl

PlasmaProperties, a simple plasma calculator module in Julia

The module contains basic functions to compute common plasma properties, such as characteristic frequencies and lengths.

It builds on Unitful, which provides a simple unit system for physical quantities.

P = Plasma()

Mutable type holding plasma properties like density, temperature, etc. The plasma is quasineutral, cold-ion. The show() method prints a table of the properties and derived quantities.

A name can be specified to quickly select the ion mass and charge. If units are not given, the default is mi in kg, n in 1/m^3, Te in eV, B in Gauss.

CoulombLog(ne,Te)

Approximate the Coulomb logarithm.

c_s(mi,Zi,Te)
c_s(P::QuasineutralPlasma)

Compute the sound speed of ions in m/s. This ignores numerical factors and gammae.

cth(m,Z,T)

Compute the thermal speed of a charged species in m/s. This ignores numerical factors and gammae.

c_the(Te)
c_the(P::QuasineutralPlasma)

Compute the thermal speed of electrons in km/s. This ignores numerical factors and gammae.

checkorassume(value,units)

If value is not a Unitful quantity, assume it is in the given units.

f_c(m,Z,B)

Compute the (unsigned) cyclotron frequency of a charged species in MHz. If m does not have units, it is assumed to be in kg. If B does not have Unitful units, it is assumed to be in Tesla.

f_ce(B)
f_ce(P::QuasineutralPlasma)

Compute the (unsigned) cyclotron frequency of electrons in GHz. If B does not have Unitful units, it is assumed to be in Tesla.

f_ci(mi,Zi,B)
f_ci(P::QuasineutralPlasma)

Compute the (unsigned) cyclotron frequency of ions of mass mi in kHz. If B does not have Unitful units, it is assumed to be in Tesla.

f_pe(m,Z,n)

Compute the plasma frequency of a charged species in MHz. If m does not have units, it is assumed to be in kg. If n does not have Unitful units, it is assumed to be in m^-3.

f_pe(ne)
f_pe(P::QuasineutralPlasma)

Compute the plasma frequency of electrons in GHz. If ne does not have Unitful units, it is assumed to be in m^-3.

f_pi(mi,Z,ni)
f_pi(P::QuasineutralPlasma)

Compute the plasma frequency of ions of mass mi in MHz. If ni does not have Unitful units, it is assumed to be in m^-3.

force_energy_units!(value)

If value is in temperature units, convert it to energy units.

rL(m,Z,B,vperp)

Compute the gyro radius of a charged species in m. If m does not have units, it is assumed to be in kg. If B does not have Unitful units, it is assumed to be in Tesla. If vperp does not have Unitful units, it is assumed to be in m/s.

r_Le(B,vperpe)
r_Le(P::QuasineutralPlasma)

Compute the electron gyro radius in mm. If B does not have Unitful units, it is assumed to be in Tesla. If vperpe does not have Unitful units, it is assumed to be in m/s.

r_Lsi(mi,Z,B,Te)
r_Lsi(P::QuasineutralPlasma)

Compute the sonic ion gyro radius based on the sound speed, in cm. If mi does not have units, it is assumed to be in kg. If B does not have Unitful units, it is assumed to be in Tesla. If Te does not have Unitful units, it is assumed to be in eV.

λ_De(ne,Te)
λ_De(P::QuasineutralPlasma)

Compute the Debye length of the plasm in mm. If ne does not have Unitful units, it is assumed to be in m^-3. If Te does not have Unitful units, it is assumed to be in eV.

ν_ee(ne,Te)
ν_ee(P::QuasineutralPlasma)

Electron-electron momentum collision frequency in MHz.

ν_ei(Zi,np,Te)
ν_ei(P::QuasineutralPlasma)

Electron-ion momentum collision frequency in MHz.

ν_ie(mi,Zi,np,Te)
ν_ie(P::QuasineutralPlasma)

Ion-electron momentum collision frequency in MHz.

χ(Zi,np,Te,B)
χ(P::QuasineutralPlasma)

Hall parameter based on the electron-ion collision frequency only.