📝 Enhance documentation for resistive plasma heating and clarify refe…

…rences to ITER Physics Design Guidelines

documentation/proc-pages/physics-models/plasma_current/plasma_resistive_heating.md

@@ -1 +1,25 @@
-# Resistive Plasma Heating
+# Resistive Plasma Heating
+
+The ohmic component of the plasma heating is given by that from the ITER 1989 Physics Design Guidelines[^1]
+
+Using the resistive loop voltage for a reference profile of parabolic shape with:
+
+$$
+\alpha_n \approx 0.5, \alpha_T \approx 1.0,  \alpha_J \approx 1.5
+$$
+
+
+$$
+U \  [\text{V}] \approx 2.15 \times 10^{-3} Z_{\text{eff}}\langle \gamma_{\text{NC}} \rangle \frac{I_{\text{ind}}R_0}{\kappa a^2} \frac{1}{T_{10}^{1.5}}
+$$
+
+The neoclassical (avergae) resisitivity enhancement factor $\left(\langle \gamma_{\text{NC}} \rangle \right)$ is given by an empirical fit:
+
+$$
+\langle \gamma_{\text{NC}} \rangle = 4.3 -0.6A
+$$
+
+where $A$ is valid in the range of 2.5 - 4.0.
+
+
+[^1]: N.A. Uckan and ITER Physics Group, 'ITER Physics Design Guidelines: 1989',

process/physics.py

@@ -2937,7 +2937,7 @@ def plasma_ohmic_heating(
         # Density weighted electron temperature in 10 keV units
         t10 = ten / 10.0
 
-        # Plasma resistance, from loop voltage calculation in IPDG89
+        # Plasma resistance, from loop voltage calculation in ITER Physics Design Guidelines: 1989
         res_plasma = (
             physics_variables.plasma_res_factor
             * 2.15e-9
@@ -2947,8 +2947,8 @@ def plasma_ohmic_heating(
         )
 
         # Neo-classical resistivity enhancement factor
-        # Taken from  N. A. Uckan et al, Fusion Technology 13 (1988) p.411.
-        # The expression is valid for aspect ratios in the range 2.5--4.
+        # Taken from ITER Physics Design Guidelines: 1989
+        # The expression is valid for aspect ratios in the range 2.5 to 4.0
         rpfac = 4.3 - 0.6 * rmajor / rminor
         res_plasma = res_plasma * rpfac