🐛 Update minimum beta constraint to work on the total beta and not th…

…ermal
ukaea · Dec 19, 2024 · 9510a25 · 9510a25
1 parent 43eba4d
commit 9510a25
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/documentation/proc-pages/physics-models/plasma_beta/plasma_beta.md b/documentation/proc-pages/physics-models/plasma_beta/plasma_beta.md
@@ -191,6 +191,8 @@ The value of `beta_poloidal_max` can be set to the desired maximum poloidal beta
 
 This constraint can be activated by stating `icc = 84` in the input file.
 
+The value of `beta_max` can be set to the desired minimum total beta. The scaling value `fbeta_min` can be varied also.
+
 [^1]: N.A. Uckan and ITER Physics Group, 'ITER Physics Design Guidelines: 1989',
 
 [^2]: D.J. Ward, 'PROCESS Fast Alpha Pressure', Work File Note F/PL/PJK/PROCESS/CODE/050

diff --git a/source/fortran/constraint_equations.f90 b/source/fortran/constraint_equations.f90
@@ -3137,9 +3137,8 @@ subroutine constraint_eqn_084(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units)
       !! fbeta_min : input real : f-value for constraint beta-beta_fast_alpha > beta_min
       !! beta_min : input real :  Lower limit for beta
       !! beta : input real :  plasma beta
-      !! beta_fast_alpha : input real : Alpha particle beta
 
-      use physics_variables, only: beta_min, beta, beta_fast_alpha
+      use physics_variables, only: beta_min, beta
       use constraint_variables, only: fbeta_min
       implicit none
             real(dp), intent(out) :: tmp_cc
@@ -3149,9 +3148,9 @@ subroutine constraint_eqn_084(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units)
       character(len=10), intent(out) :: tmp_units
 
 
-      tmp_cc = 1.0D0 - fbeta_min * (beta-beta_fast_alpha)/beta_min
+      tmp_cc = 1.0D0 - fbeta_min * (beta)/beta_min
       tmp_con = beta_min * (1.0D0 - tmp_cc)
-      tmp_err = (beta-beta_fast_alpha) * tmp_cc
+      tmp_err = (beta) * tmp_cc
       tmp_symbol = '>'
       tmp_units = ''