constraint_variables Module

This module contains global variables relating to the constraint equations (f-values, limits, etc.).

References

-



Contents


Variables

TypeVisibility AttributesNameInitial
real(kind=dp), public :: auxmin

minimum auxiliary power (MW) (constraint equation 40)

real(kind=dp), public :: betpmx

maximum poloidal beta (constraint equation 48)

real(kind=dp), public :: bigqmin

minimum fusion gain Q (constraint equation 28)

real(kind=dp), public :: bmxlim

maximum peak toroidal field (T) (constraint equation 25)

real(kind=dp), public :: fauxmn

f-value for minimum auxiliary power (constraint equation 40, iteration variable 64)

real(kind=dp), public :: fbeta

f-value for epsilon beta-poloidal (constraint equation 6, iteration variable 8)

real(kind=dp), public :: fbetap

f-value for poloidal beta (constraint equation 48, iteration variable 79)

real(kind=dp), public :: fbetatry

f-value for beta limit (constraint equation 24, iteration variable 36)

real(kind=dp), public :: fbetatry_lower

f-value for (lower) beta limit (constraint equation 84, iteration variable 173)

real(kind=dp), public :: fcpttf

f-value for TF coil current per turn upper limit (constraint equation 77, iteration variable 146)

real(kind=dp), public :: fcwr

f-value for conducting wall radius / rminor limit (constraint equation 23, iteration variable 104)

real(kind=dp), public :: fdene

f-value for density limit (constraint equation 5, iteration variable 9) (invalid if ipedestal=3)

real(kind=dp), public :: fdivcol

f-value for divertor collisionality (constraint equation 22, iteration variable 34)

real(kind=dp), public :: fdtmp

f-value for first wall coolant temperature rise (constraint equation 38, iteration variable 62)

real(kind=dp), public :: fecrh_ignition

f-value for ecrh ignition constraint (constraint equation 91, iteration variable 168)

real(kind=dp), public :: fflutf

f-value for neutron fluence on TF coil (constraint equation 53, iteration variable 92)

real(kind=dp), public :: ffuspow

f-value for maximum fusion power (constraint equation 9, iteration variable 26)

real(kind=dp), public :: fgamcd

f-value for current drive gamma (constraint equation 37, iteration variable 40)

real(kind=dp), public :: fhldiv

f-value for divertor heat load (constraint equation 18, iteration variable 27)

real(kind=dp), public :: fiooic

f-value for TF coil operating current / critical current ratio (constraint equation 33, iteration variable 50)

real(kind=dp), public :: fipir

f-value for Ip/Irod upper limit constraint equation icc = 46 iteration variable ixc = 72

real(kind=dp), public :: fjohc

f-value for central solenoid current at end-of-flattop (constraint equation 26, iteration variable 38)

real(kind=dp), public :: fjohc0

f-value for central solenoid current at beginning of pulse (constraint equation 27, iteration variable 39)

real(kind=dp), public :: fjprot

f-value for TF coil winding pack current density (constraint equation 35, iteration variable 53)

real(kind=dp), public :: flhthresh

f-value for L-H power threshold (constraint equation 15, iteration variable 103)

real(kind=dp), public :: fmva

f-value for maximum MVA (constraint equation 19, iteration variable 30)

real(kind=dp), public :: fnbshinef

f-value for maximum neutral beam shine-through fraction (constraint equation 59, iteration variable 105)

real(kind=dp), public :: fncycle

f-value for minimum CS coil stress load cycles (constraint equation 90, iteration variable 167)

real(kind=dp), public :: fnesep

f-value for Eich critical separatrix density (constraint equation 76, iteration variable 144)

real(kind=dp), public :: foh_stress

f-value for Tresca yield criterion in Central Solenoid (constraint equation 72, iteration variable 123)

real(kind=dp), public :: fpeakb

f-value for maximum toroidal field (constraint equation 25, iteration variable 35)

real(kind=dp), public :: fpinj

f-value for injection power (constraint equation 30, iteration variable 46)

real(kind=dp), public :: fpnetel

f-value for net electric power (constraint equation 16, iteration variable 25)

real(kind=dp), public :: fportsz

f-value for neutral beam tangency radius limit (constraint equation 20, iteration variable 33)

real(kind=dp), public :: fpsepbqar

f-value for maximum Psep*Bt/qAR limit (constraint equation 68, iteration variable 117)

real(kind=dp), public :: fpsepr

f-value for maximum Psep/R limit (constraint equation 56, iteration variable 97)

real(kind=dp), public :: fptemp

f-value for peak centrepost temperature (constraint equation 44, iteration variable 68)

real(kind=dp), public :: fptfnuc

f-value for maximum TF coil nuclear heating (constraint equation 54, iteration variable 95)

real(kind=dp), public :: fq

f-value for edge safety factor (constraint equation 45, iteration variable 71)

real(kind=dp), public :: fqval

f-value for Q (constraint equation 28, iteration variable 45)

real(kind=dp), public :: fradpwr

f-value for core radiation power limit (constraint equation 17, iteration variable 28)

real(kind=dp), public :: fradwall

f-value for upper limit on radiation wall load (constr. equ. 67, iteration variable 116)

real(kind=dp), public :: freinke

f-value for Reinke detachment criterion (constr. equ. 78, iteration variable 147)

real(kind=dp), public :: frminor

f-value for minor radius limit (constraint equation 21, iteration variable 32)

real(kind=dp), public :: fstrcase

f-value for maximum TF coil case Tresca yield criterion (constraint equation 31, iteration variable 48)

real(kind=dp), public :: fstrcond

f-value for maxiumum TF coil conduit Tresca yield criterion (constraint equation 32, iteration variable 49)

real(kind=dp), public :: fstr_wp

f-value for maxiumum TF coil strain absolute value (constraint equation 88, iteration variable 165)

real(kind=dp), public :: fmaxvvstress

f-value for maximum permitted stress of the VV (constraint equation 65, iteration variable 113)

real(kind=dp), public :: ftbr

f-value for minimum tritium breeding ratio (constraint equation 52, iteration variable 89)

real(kind=dp), public :: ft_burn

f-value for minimum burn time (constraint equation 13, iteration variable 21)

real(kind=dp), public :: ftcycl

f-value for cycle time (constraint equation 42, iteration variable 67)

real(kind=dp), public :: ftmargoh

f-value for central solenoid temperature margin (constraint equation 60, iteration variable 106)

real(kind=dp), public :: ftmargtf

f-value for TF coil temperature margin (constraint equation 36, iteration variable 54)

real(kind=dp), public :: ft_current_ramp_up

f-value for plasma current ramp-up time (constraint equation 41, iteration variable 66)

real(kind=dp), public :: ftpeak

f-value for first wall peak temperature (constraint equation 39, iteration variable 63)

real(kind=dp), public :: fvdump

f-value for dump voltage (constraint equation 34, iteration variable 51)

real(kind=dp), public :: fvs

f-value for flux-swing (V-s) requirement (STEADY STATE) (constraint equation 12, iteration variable 15)

real(kind=dp), public :: fvvhe

f-value for vacuum vessel He concentration limit (iblanket = 2) (constraint equation 55, iteration variable 96)

real(kind=dp), public :: fwalld

f-value for maximum wall load (constraint equation 8, iteration variable 14)

real(kind=dp), public :: fzeffmax

f-value for maximum zeff (constraint equation 64, iteration variable 112)

real(kind=dp), public :: gammax

maximum current drive gamma (constraint equation 37)

real(kind=dp), public :: maxradwallload

Maximum permitted radiation wall load (MW/m^2) (constraint equation 67)

real(kind=dp), public :: mvalim

maximum MVA limit (constraint equation 19)

real(kind=dp), public :: nbshinefmax

maximum neutral beam shine-through fraction (constraint equation 59)

real(kind=dp), public :: nflutfmax

max fast neutron fluence on TF coil (n/m2) (blktmodel>0) (constraint equation 53) Also used for demontable magnets (itart = 1) and superconducting coils (i_tf_sup = 1) To set the CP lifetime (constraint equation 85)

real(kind=dp), public :: pdivtlim

Minimum pdivt [MW] (constraint equation 80)

real(kind=dp), public :: peakfactrad

peaking factor for radiation wall load (constraint equation 67)

real(kind=dp), public :: peakradwallload

Peak radiation wall load (MW/m^2) (constraint equation 67)

real(kind=dp), public :: pnetelin

required net electric power (MW) (constraint equation 16)

real(kind=dp), public :: powfmax

maximum fusion power (MW) (constraint equation 9)

real(kind=dp), public :: psepbqarmax

maximum ratio of Psep*Bt/qAR (MWT/m) (constraint equation 68)

real(kind=dp), public :: pseprmax

maximum ratio of power crossing the separatrix to plasma major radius (Psep/R) (MW/m) (constraint equation 56)

real(kind=dp), public :: ptfnucmax

maximum nuclear heating in TF coil (MW/m3) (constraint equation 54)

real(kind=dp), public :: tbrmin

minimum tritium breeding ratio (constraint equation 52)

real(kind=dp), public :: t_burn_min

minimum burn time (s) (KE - no longer itv., see issue #706)

real(kind=dp), public :: tcycmn

minimum cycle time (s) (constraint equation 42)

real(kind=dp), public :: t_current_ramp_up_min

minimum plasma current ramp-up time (s) (constraint equation 41)

real(kind=dp), public :: vvhealw

allowed maximum helium concentration in vacuum vessel at end of plant life (appm) (iblanket =2) (constraint equation 55)

real(kind=dp), public :: walalw

allowable neutron wall-load (MW/m2) (constraint equation 8)

real(kind=dp), public :: taulimit

Lower limit on taup/taueff the ratio of alpha particle to energy confinement times (constraint equation 62)

real(kind=dp), public :: ftaulimit

f-value for lower limit on taup/taueff the ratio of alpha particle to energy confinement times (constraint equation 62, iteration variable 110)

real(kind=dp), public :: fniterpump

f-value for constraint that number of pumps < tfno (constraint equation 63, iteration variable 111)

real(kind=dp), public :: zeffmax

maximum value for Zeff (constraint equation 64)

real(kind=dp), public :: fpoloidalpower

f-value for constraint on rate of change of energy in poloidal field (constraint equation 66, iteration variable 115)

real(kind=dp), public :: fpsep

f-value to ensure separatrix power is less than value from Kallenbach divertor (Not required as constraint 69 is an equality)

real(kind=dp), public :: fcqt

TF coil quench temparature remains below tmax_croco (constraint equation 74, iteration variable 141)


Subroutines

public subroutine init_constraint_variables()

Initialise module variables

Arguments

None