I think that beta_n is off by a factor of 3/2. It is being computed from the volume average of W_tot instead of P (where Pressure per m^3 is n_eT_e+n_iT_i, but W_tot per m^3 is 3/2 (n_eT_e+n_iT_i), but the docstring in the W_tot def. says they're equal.) Fixing this makes my popcon runs match ITER and MANTA beta_n values much better. Let me know if I'm missing something, and if not, I'm happy to implement a fix @matt-pharr.
Beta_n def:
|
def BetaN(s, T_i_keV, n_e_20) -> float: |
W_tot def:
|
def _W_tot_prof(s, rho, T_i_keV: float, n_e_20: float): |
I think that beta_n is off by a factor of 3/2. It is being computed from the volume average of W_tot instead of P (where Pressure per m^3 is n_eT_e+n_iT_i, but W_tot per m^3 is 3/2 (n_eT_e+n_iT_i), but the docstring in the W_tot def. says they're equal.) Fixing this makes my popcon runs match ITER and MANTA beta_n values much better. Let me know if I'm missing something, and if not, I'm happy to implement a fix @matt-pharr.
Beta_n def:
OpenPOPCON/openpopcon/core.py
Line 436 in 8d6ab81
W_tot def:
OpenPOPCON/openpopcon/core.py
Line 512 in 8d6ab81