Add general plasma functions 1

documentation/physics-models/detailed_physics.md

@@ -0,0 +1,85 @@
+# Detailed Plasma Physics
+
+It can sometimes be useful to calculate rough values for key plasma paramters that are normally used in higher fidelity codes. The `DetailedPhysics()` class stores functions that are called and the end of the run to show rough values for key plasma behavior parameters. The calculation is done at the end as no other methods currently depend on these values.
+
+## Detailed Plasma Physics | `DetailedPhysics()`
+
+
+------------------
+
+### Debye length | `calculate_debye_length()`
+
+Calculates the Debye lenght given by:
+
+$$
+\lambda_{D} = \sqrt{\frac{\epsilon_0 k_B T_e}{n e^2}}
+$$
+
+-------------------
+
+### Relativistic particle speed | `calculate_relativistic_particle_speed()`
+
+$$
+v = c \times \sqrt{\left(1- \frac{1}{\left(1+\frac{E}{mc^2}\right)^2}\right)}
+$$
+
+------------------
+
+### Coulomb Logarithm | `calculate_coulomb_log_from_impact()`
+
+Calculates the Coulomb logarithm assuming a straight line Landau-Spitzer method
+
+$$
+\ln \Lambda = \ln{\left(\frac{b_{\text{max}}}{b_{\text{min}}}\right)}
+$$
+
+The maximum impact parameter is given by the Debye length calculated by [`calculate_debye_length()`](#debye-length--calculate_debye_length)
+$$
+b_{\text{max}} = \lambda_{\text{Debye}}
+$$
+
+The minimum impact paramter is the largest of either the classical distance of closest approach or the Debye length. 
+
+$$
+\begin{split}b_{\text{min}} ≡
+\left\{
+    \begin{array}{ll}
+               λ_{\text{de Broglie}} & \mbox{if } λ_{\text{de Broglie}} ≥ ρ_⟂ \\
+               ρ_⟂         & \mbox{if } ρ_⟂ ≥ λ_{\text{de Broglie}}
+    \end{array}
+\right.\end{split}
+$$
+
+$ρ_⟂$ is the classical distance of closest approach calculated by [`calculate_classical_distance_of_closest_approach()`](#classical-distance-of-closest-approach----calculate_classical_distance_of_closest_approach)
+
+------------------
+
+### Classical distance of closest approach |   `calculate_classical_distance_of_closest_approach()`
+
+$$
+\frac{Z_1Z_2e^2}{4\pi \epsilon_0 E_{\text{kinetic}}}
+$$
+
+---------------------
+
+### DeBroglie Wavelength | `calculate_debroglie_wavelength()`
+
+$$
+\lambda_{\text{DeBroglie}} = \frac{h}{2\pi m v}
+$$
+
+----------------------
+
+### Plasma Frequency | `calculate_plasma_frequency()`
+
+$$
+\omega_p = \sqrt{\frac{n_ie^2}{\epsilon_0 m_i}}
+$$
+
+---------------------
+
+### Larmor Frequency | `calculate_larmor_frequency()`
+
+$$
+f_{\text{Larmor}} = \frac{Z_ieB}{2\pi m_i}
+$$

mkdocs.yml

@@ -70,6 +70,7 @@ nav:
               - Confinement time: physics-models/plasma_confinement.md
               - L-H transition: physics-models/plasma_h_mode.md
               - Plasma Core Power Balance: physics-models/plasma_power_balance.md
+              - Detailed Plasma Physics: physics-models/detailed_physics.md
           - Pulsed Plant Operation: physics-models/pulsed-plant.md
       - Engineering Models:
           - Machine Build: eng-models/machine-build.md

process/constants.py

@@ -189,12 +189,15 @@
 https://physics.nist.gov/cgi-bin/Compositions/stand_alone.pl?ele=W
 """
 
-SPEED_LIGHT = 299792458
+SPEED_LIGHT = 299792458.0
 """Speed of light in vacuum (c) [m/s]
 Reference: National Institute of Standards and Technology (NIST)
 https://physics.nist.gov/cgi-bin/cuu/Value?c|search_for=light
 """
 
+PLANCK_CONSTANT = 6.62607015e-34
+"""Planck's constant [J.s]"""
+
 D_T_ENERGY = (
     (DEUTERON_MASS + TRITON_MASS) - (ALPHA_MASS + NEUTRON_MASS)
 ) * SPEED_LIGHT**2

process/data_structure/physics_variables.py

@@ -1320,6 +1320,24 @@
 n_charge_plasma_effective_mass_weighted_vol_avg: float = None
 """Plasma mass-weighted volume averaged plasma effective charge"""
 
+len_plasma_debye_electron_profile: list[float] = None
+"""Profile of electron Debye length in plasma (m)"""
+
+len_plasma_debye_electron_vol_avg: float = None
+"""Volume averaged electron Debye length in plasma (m)"""
+
+vel_plasma_electron_profile: list[float] = None
+"""Profile of electron thermal velocity in plasma (m/s)"""
+
+plasma_coulomb_log_electron_electron_profile: list[float] = None
+"""Profile of electron-electron Coulomb logarithm in plasma"""
+
+freq_plasma_electron_profile: list[float] = None
+"""Electron plasma frequency profile (Hz)"""
+
+freq_plasma_larmor_toroidal_electron_profile: list[float] = None
+"""Profile of electron Larmor frequency in plasma due to toroidal magnetic field (Hz)"""
+
 
 def init_physics_module():
     """Initialise the physics module"""
@@ -1639,7 +1657,13 @@ def init_physics_variables():
         a_plasma_poloidal, \
         n_charge_plasma_effective_vol_avg, \
         n_charge_plasma_effective_profile, \
-        n_charge_plasma_effective_mass_weighted_vol_avg
+        n_charge_plasma_effective_mass_weighted_vol_avg, \
+        len_plasma_debye_electron_profile, \
+        len_plasma_debye_electron_vol_avg, \
+        vel_plasma_electron_profile, \
+        plasma_coulomb_log_electron_electron_profile, \
+        freq_plasma_electron_profile, \
+        freq_plasma_larmor_toroidal_electron_profile
 
     m_beam_amu = 0.0
     m_fuel_amu = 0.0
@@ -1900,3 +1924,9 @@ def init_physics_variables():
     n_charge_plasma_effective_vol_avg = 0.0
     n_charge_plasma_effective_profile = []
     n_charge_plasma_effective_mass_weighted_vol_avg = 0.0
+    len_plasma_debye_electron_profile = []
+    len_plasma_debye_electron_vol_avg = 0.0
+    vel_plasma_electron_profile = []
+    plasma_coulomb_log_electron_electron_profile = []
+    freq_plasma_electron_profile = []
+    freq_plasma_larmor_toroidal_electron_profile = []