Skip to content

tfcoil_variables

Module containing global variables relating to the toroidal field coil systems

References

  • ITER Magnets design description document DDD11-2 v2 2 (2009)

N_RADIAL_ARRAY = 50 module-attribute

Size of the radial distribution arrays per layers used for stress, strain and displacement distibution

a_tf_coil_inboard_case = None module-attribute

external case area per coil (inboard leg) (m2)

a_tf_coil_outboard_case = None module-attribute

external case area per coil (outboard leg) (m2)

a_tf_turn_steel = None module-attribute

area of the cable conduit (m2)

a_tf_wp_conductor = None module-attribute

Winding pack conductor area [m2] Does not include the area of voids and central helium channel

a_res_tf_coil_conductor = None module-attribute

Area of resistive conductor in resistive TF coil [m2]

a_tf_turn_cable_space_no_void = None module-attribute

Cable space area (per turn) [m2] Includes the area of voids and central helium channel

a_tf_turn_insulation = None module-attribute

single turn insulation area (m2)

a_tf_coil_wp_turn_insulation = None module-attribute

winding pack turn insulation area per coil (m2)

sig_tf_case_max = None module-attribute

Allowable maximum shear stress (Tresca criterion) in TF coil case (Pa)

sig_tf_wp_max = None module-attribute

Allowable maximum shear stress (Tresca criterion) in TF coil conduit (Pa)

a_tf_leg_outboard = None module-attribute

outboard TF leg area (m2)

a_tf_wp_steel = None module-attribute

Total area of all winding pack steel (sum of each conduit steel from each turn) (m2)

a_tf_wp_extra_void = None module-attribute

winding pack void (He coolant) area (m2)

a_tf_wp_coolant_channels = None module-attribute

winding pack He coil area (m2)

bcritsc = None module-attribute

upper critical field (T) for Nb3Sn superconductor at zero temperature and strain (i_tf_sc_mat=4, =bc20m)

b_tf_inboard_peak_symmetric = None module-attribute

mean peak field at TF coil (T)

b_tf_inboard_peak_with_ripple = None module-attribute

peak field at TF conductor with ripple (T)

casestr = None module-attribute

case strain

dr_tf_plasma_case = None module-attribute

inboard TF coil case plasma side thickness (m) (calculated for stellarators)

f_dr_tf_plasma_case = None module-attribute

inboard TF coil case plasma side thickness as a fraction of dr_tf_inboard

i_f_dr_tf_plasma_case = None module-attribute

logical switch to make dr_tf_plasma_case a fraction of TF coil thickness (f_dr_tf_plasma_case)

dx_tf_side_case_min = None module-attribute

inboard TF coil minimum sidewall case thickness (m) (calculated for stellarators)

dx_tf_side_case_peak = None module-attribute

inboard TF coil peak sidewall case thickness (m) (calculated for stellarators)

casths_fraction = None module-attribute

inboard TF coil sidewall case thickness as a fraction of dx_tf_inboard_out_toroidal

tfc_sidewall_is_fraction = None module-attribute

logical switch to make dx_tf_side_case_min a fraction of TF coil thickness (casths_fraction)

t_conductor = None module-attribute

Conductor (cable + steel conduit) area averaged dimension [m]

dx_tf_turn_general = None module-attribute

TF coil turn edge length including turn insulation [m] If the turn is not a square (i_tf_turns_integer = 1) a squared turn of equivelent size is use to calculated this quantity If the dx_tf_turn_general is non zero, c_tf_turn is calculated

i_dx_tf_turn_general_input = None module-attribute

Boolean switch to activated when the user set the TF coil turn dimensions Not an input

t_turn_tf_max = None module-attribute

TF turn edge length including turn insulation upper limit [m] If the turn is not a square (i_tf_turns_integer = 1) a squared turn of equivelent size is use for this constraint constraint equation icc = 86

dx_tf_turn_cable_space_general = None module-attribute

TF coil superconducting cable squared/rounded dimensions [m] If the turn is not a square (i_tf_turns_integer = 1) a squared cable of equivelent size is use to calculated this quantity If the dx_tf_turn_cable_space_general is non zero, c_tf_turn is calculated

i_dx_tf_turn_cable_space_general_input = None module-attribute

Boolean switch to activated when the user set the TF coil cable dimensions Not an input

acs = None module-attribute

Area of space inside conductor (m2)

cdtfleg = None module-attribute

TF outboard leg current density (A/m2) (resistive coils only)

cforce = None module-attribute

centering force on inboard leg (per coil) (N/m)

cplen = None module-attribute

length of TF coil inboard leg ('centrepost') (i_tf_sup = 1)

c_tf_turn = None module-attribute

TF coil current per turn (A). (calculated for stellarators) (calculated for integer-turn TF coils i_tf_turns_integer=1) (iteration variable 60)

c_tf_turn_max = None module-attribute

Max TF coil current per turn [A]. (for stellarators and i_tf_turns_integer=1) (constraint equation 77)

den_tf_coil_case = None module-attribute

density of coil case (kg/m3)

dcond = None module-attribute

density of superconductor type given by i_tf_sc_mat/i_cs_superconductor/i_pf_superconductor (kg/m3)

den_tf_wp_turn_insulation = None module-attribute

density of conduit + ground-wall insulation (kg/m3)

dia_tf_turn_coolant_channel = None module-attribute

diameter of central helium channel in TF winding (m)

e_tf_magnetic_stored_total_gj = None module-attribute

total magnetic stored energy in the toroidal field coils (GJ)

e_tf_coil_magnetic_stored = None module-attribute

Stored magnetic energy in a single TF coil (J)

b_crit_upper_nbti = None module-attribute

upper critical field of GL_nbti

t_crit_nbti = None module-attribute

critical temperature of GL_nbti

max_force_density = None module-attribute

Maximal (WP averaged) force density in TF coils at 1 point. (MN/m3)

f_a_tf_turn_cable_copper = None module-attribute

copper fraction of cable conductor (TF coils) (iteration variable 59)

fhts = None module-attribute

technology adjustment factor for critical current density fit for isumat..=2 Bi-2212 superconductor, to describe the level of technology assumed (i.e. to account for stress, fatigue, radiation, AC losses, joints or manufacturing variations; 1.0 would be very optimistic)

insstrain = None module-attribute

Radial strain in insulator

i_tf_stress_model = None module-attribute

Switch for the TF coil stress model 0 : Generalized plane strain formulation, Issues #977 and #991, O(n^3) 1 : Old plane stress model (only for SC) 2 : Axisymmetric extended plane strain, Issues #1414 and #998, O(n)

i_tf_tresca = None module-attribute

Switch for TF coil conduit Tresca stress criterion: 0 : Tresca (no adjustment); 1 : Tresca with CEA adjustment factors (radial+2%, vertical+60%)

i_tf_wp_geom = None module-attribute

Switch for TF WP geometry selection 0 : Rectangular geometry 1 : Double rectangular geometry 2 : Trapezoidal geometry (constant lateral casing thickness) Default setting for backward compatibility if i_tf_turns_integer = 0 : Double rectangular if i_tf_turns_integer = 1 : Rectangular

i_tf_case_geom = None module-attribute

Switch for TF case geometry selection 0 : Circular front case (ITER design) 1 : Straight front case

i_tf_turns_integer = None module-attribute

Switch for TF coil integer/non-integer turns: 0 : non-integer turns 1 : integer turns

i_tf_sc_mat = None module-attribute

Switch for superconductor material in TF coils: - =1 ITER Nb3Sn critical surface model with standard ITER parameters - =2 Bi-2212 high temperature superconductor (range of validity T < 20K, adjusted field b < 104 T, B > 6 T) - =3 NbTi - =4 ITER Nb3Sn model with user-specified parameters - =5 WST Nb3Sn parameterisation - =6 REBCO HTS tape in CroCo strand - =7 Durham Ginzburg-Landau critical surface model for Nb-Ti - =8 Durham Ginzburg-Landau critical surface model for REBCO - =9 Hazelton experimental data + Zhai conceptual model for REBCO

i_tf_sup = None module-attribute

Switch for TF coil conductor model: - =0 copper - =1 superconductor - =2 Cryogenic aluminium

i_tf_shape = None module-attribute

Switch for TF coil toroidal shape: - =0 Default value : Picture frame coil for TART / PROCESS D-shape for non itart - =1 PROCESS D-shape : parametrise with 2 arcs - =2 Picture frame coils

i_tf_cond_eyoung_axial = None module-attribute

Switch for the behavior of the TF coil conductor elastic axial properties - =0 Young's modulus is set to zero, and the conductor is not considered in the stress calculation. This corresponds to the case that the conductor is much less stiff than the conduit, or the case that the conductor is prevented (isolated) from taking axial loads. - =1 Elastic properties are set by user input, using the variable eyoung_cond_axial - =2 Elastic properties are set to reasonable defaults taking into account the superconducting material i_tf_sc_mat

i_tf_cond_eyoung_trans = None module-attribute

Switch for the behavior of the elastic properties of the TF coil conductorin the transverse direction. Only active if i_tf_cond_eyoung_axial == 2 - =0 Cable not potted in solder. Transverse Young's modulus set to zero. - =1 Cable potted in solder. If i_tf_cond_eyoung_axial == 2, the transverse Young's modulus of the conductor is equal to the axial, which is set to a sensible material-dependent default.

n_tf_wp_pancakes = None module-attribute

Number of pancakes in TF coil. Only used if i_tf_turns_integer=1

n_tf_wp_layers = None module-attribute

Number of layers in TF coil. Only used if i_tf_turns_integer=1

n_rad_per_layer = None module-attribute

Size of the arrays per layers storing the radial dependent stress quantities (stresses, strain displacement etc..)

i_tf_bucking = None module-attribute

Switch for TF inboard suport structure design: Default setting for backward compatibility - if copper resistive TF (i_tf_sup = 0) : Free standing TF without bucking structure - if Superconducting TF (i_tf_sup = 1) : Free standing TF with a steel casing - if aluminium TF (i_tf_sup = 2) : Free standing TF with a bucking structure Rem : the case is a bucking structure - =0 : Free standing TF without case/bucking cyliner (only a conductor layer) - =1 : Free standing TF with a case/bucking cylinder made of - if copper resistive TF (i_tf_sup = 0) : used defined bucking cylinder - if Superconducting TF (i_tf_sup = 1) : Steel casing - if aluminium resisitive TF (i_tf_sup = 2) : used defined bucking cylinder - =2 : The TF is in contact with the CS : "bucked and wedged design" Fast version : thin TF-CS interface neglected in the stress calculations (3 layers) The CS is frictionally decoupled from the TF, does not carry axial tension - =3 : The TF is in contact with the CS : "bucked and wedged design" Full version : thin TF-CS Kapton interface introduced in the stress calculations (4 layers) The CS and kaptop are frictionally decoupled from the TF, do not carry axial tension

n_tf_graded_layers = None module-attribute

Number of layers of different stress properties in the WP. If n_tf_graded_layers > 1, a graded coil is condidered

n_tf_stress_layers = None module-attribute

Number of layers considered for the inboard TF stress calculations set in initial.f90 from i_tf_bucking and n_tf_graded_layers

n_tf_wp_stress_layers = None module-attribute

Maximum number of layers that can be considered in the TF coil composited/smeared stress analysis. This is the layers of one turn, not the entire WP. Default: 5. void, conductor, copper, conduit, insulation.

j_tf_bus = None module-attribute

bussing current density (A/m2)

j_crit_str_tf = None module-attribute

j_crit_str : superconductor strand critical current density under operating conditions (A/m2). Necessary for the cost calculation in $/kAm

j_crit_str_0 = None module-attribute

j_crit_str_pf_0 : superconductor strand critical current density at 6 T and 4.2 K (A/m2) Necessary for the cost calculation in $/kAm

j_tf_wp_critical = None module-attribute

critical current density for winding pack (A/m2)

j_tf_wp_quench_heat_max = None module-attribute

allowable TF coil winding pack current density, for dump temperature rise protection (A/m2)

j_tf_wp = None module-attribute

winding pack engineering current density (A/m2)

oacdcp = None module-attribute

Overall current density in TF coil inboard legs midplane (A/m2) Rem SK : Not used in tfcoil to set the current any more. Should not be used as iteration variable 12 any more. It is now calculated.

eyoung_ins = None module-attribute

Insulator Young's modulus [Pa]. Default value (1.0D8) setup the following values - SC TF, eyoung_ins = 20 Gpa (default value from DDD11-2 v2 2 (2009)) - Al TF, eyoung_ins = 2.5 GPa (Kapton polymer)

eyoung_steel = None module-attribute

Steel case Young's modulus (Pa) (default value from DDD11-2 v2 2 (2009))

eyoung_cond_axial = None module-attribute

SC TF coil conductor Young's modulus in the parallel (along the wire/tape) direction [Pa] Set by user input only if i_tf_cond_eyoung_axial == 1; otherwise set by the behavior of that switch.

eyoung_cond_trans = None module-attribute

SC TF coil conductor Young's modulus in the transverse direction [Pa] Set by user input only if i_tf_cond_eyoung_axial == 1; otherwise set by the behavior of that switch.

eyoung_res_tf_buck = None module-attribute

Resistive TF magnets bucking cylinder young modulus (Pa)

eyoung_copper = None module-attribute

Copper young modulus. Default value taken from wikipedia

eyoung_al = None module-attribute

Aluminium young modulus. Default value taken from wikipedia

poisson_steel = None module-attribute

poisson_copper = None module-attribute

poisson_al = None module-attribute

poisson_ins = None module-attribute

Insulation Poisson's ratio. Default: Kapton. Source : DuPont™ Kapton® HN datasheet.

poisson_cond_axial = None module-attribute

SC TF coil conductor Poisson's ratio in the parallel-transverse direction

poisson_cond_trans = None module-attribute

SC TF coil conductor Poisson's ratio in the transverse-transverse direction

r_b_tf_inboard_peak = None module-attribute

Radius of maximum TF B-field (m)

res_tf_leg = None module-attribute

TF coil leg resistance (ohm)

toroidalgap = None module-attribute

Minimal distance between two toroidal coils. (m)

ripple_b_tf_plasma_edge_max = None module-attribute

maximum allowable toroidal field ripple amplitude at plasma edge (%)

ripple_b_tf_plasma_edge = None module-attribute

peak/average toroidal field ripple at plasma edge (%)

c_tf_total = None module-attribute

total (summed) current in TF coils (A)

radial_array = None module-attribute

Array refining the radii of the stress calculations arrays

sig_tf_r = None module-attribute

TF Inboard leg radial stress in steel r distribution at mid-plane [Pa]

sig_tf_t = None module-attribute

TF Inboard leg tangential stress in steel r distribution at mid-plane [Pa]

deflect = None module-attribute

TF coil radial deflection (displacement) radial distribution [m]

sig_tf_z = None module-attribute

TF Inboard leg vertical tensile stress in steel at mid-plane [Pa]

sig_tf_vmises = None module-attribute

TF Inboard leg Von-Mises stress in steel r distribution at mid-plane [Pa]

s_shear_tf = None module-attribute

TF Inboard leg maximum shear stress (Tresca criterion) in steel r distribution at mid-plane [Pa]

sig_tf_cs_bucked = None module-attribute

sig_tf_case = None module-attribute

Maximum shear stress (Tresca criterion) in TF casing steel structures (Pa)

sig_tf_wp = None module-attribute

str_cs_con_res = None module-attribute

Residual manufacturing strain in CS superconductor material

str_pf_con_res = None module-attribute

Residual manufacturing strain in PF superconductor material

str_tf_con_res = None module-attribute

Residual manufacturing strain in TF superconductor material If i_str_wp == 0, used to compute the critical surface. Otherwise, the self-consistent winding pack str_wp is used.

str_wp = None module-attribute

Axial (vertical) strain in the TF coil winding pack found by self-consistent stress/strain calculation. if i_str_wp == 1, used to compute the critical surface. Otherwise, the input value str_tf_con_res is used. Constrain the absolute value using constraint equation 88 You can't have constraint 88 and i_str_wp = 0 at the same time

str_wp_max = None module-attribute

Maximum allowed absolute value of the strain in the TF coil (Constraint equation 88)

i_str_wp = None module-attribute

Switch for the behavior of the TF strain used to compute the strain-dependent critical surface: - =0 str_tf_con_res is used - =1 str_wp is used

quench_model = None module-attribute

switch for TF coil quench model (Only applies to REBCO magnet at present, issue #522): - ='exponential' exponential quench with constant discharge resistor - ='linear' quench with constant voltage

time1 = None module-attribute

Time at which TF quench is detected (s)

tcritsc = None module-attribute

critical temperature (K) for superconductor at zero field and strain (i_tf_sc_mat=4, =tc0m)

t_tf_superconductor_quench = None module-attribute

fast discharge time for TF coil in event of quench (s) (iteration variable 56) For REBCO model, meaning depends on quench_model: - exponential quench : e-folding time (s)` - linear quench : discharge time (s)

a_tf_inboard_total = None module-attribute

Area of inboard midplane TF legs (m2)

len_tf_bus = None module-attribute

TF coil bus length (m)

m_tf_bus = None module-attribute

TF coil bus mass (kg)

tfckw = None module-attribute

available DC power for charging the TF coils (kW)

tfcmw = None module-attribute

Peak power per TF power supply (MW)

p_cp_resistive_mw = None module-attribute

Peak resistive TF coil inboard leg power (MW)

p_tf_joints_resistive_mw = None module-attribute

TF joints resistive power losses (MW)

tfcryoarea = None module-attribute

surface area of toroidal shells covering TF coils (m2)

tficrn = None module-attribute

TF coil half-width - inner dr_bore (m)

ind_tf_coil = None module-attribute

TF coil inductance (H)

dx_tf_wp_insertion_gap = None module-attribute

TF coil WP insertion gap (m)

p_tf_leg_resistive_mw = None module-attribute

TF coil outboard leg resistive power (MW)

rho_cp = None module-attribute

TF coil inboard leg resistivity [Ohm-m]. If itart=0, this variable is the average resistivity over the whole magnet

rho_tf_leg = None module-attribute

Resistivity of a TF coil leg (Ohm-m)

rho_tf_bus = None module-attribute

Resistivity of a TF coil bus (Ohm-m). Default values is for that of GLIDCOP AL-15 (C15715) at 293K

frhocp = None module-attribute

Centrepost resistivity enhancement factor. For itart=0, this factor is used for the whole magnet

frholeg = None module-attribute

Ouboard legs resistivity enhancement factor. Only used for itart=1.

i_cp_joints = None module-attribute

Switch for CP demoutable joints type -= 0 : Clampled joints -= 1 : Sliding joints Default value (-1) choses : Sliding joints for resistive magnets (i_tf_sup = 0, 2) Clampled joints for superconducting magents (i_tf_sup = 1)

rho_tf_joints = None module-attribute

TF joints surfacic resistivity [ohm.m]. Feldmetal joints assumed.

n_tf_joints_contact = None module-attribute

Number of contact per turn

n_tf_joints = None module-attribute

Number of joints Ex: n_tf_joints = 2 for top and bottom CP joints

th_joint_contact = None module-attribute

TF sliding joints contact pad width [m]

p_tf_joints_resistive = None module-attribute

Calculated TF joints resistive power losses [W]

len_tf_coil = None module-attribute

TF coil circumference (m)

eff_tf_cryo = None module-attribute

TF cryoplant efficiency (compared to pefect Carnot cycle). Using -1 set the default value depending on magnet technology: - i_tf_sup = 1 : SC magnet, eff_tf_cryo = 0.13 (ITER design) - i_tf_sup = 2 : Cryo-aluminium, eff_tf_cryo = 0.4

n_tf_coils = None module-attribute

Number of TF coils (default = 50 for stellarators). Number of TF coils outer legs for ST

tfocrn = None module-attribute

TF coil half-width - outer dr_bore (m)

tfsai = None module-attribute

area of the inboard TF coil legs (m2)

tfsao = None module-attribute

area of the outboard TF coil legs (m2)

tftmp = None module-attribute

peak helium coolant temperature in TF coils and PF coils (K)

dx_tf_inboard_out_toroidal = None module-attribute

TF coil toroidal thickness (m)

dx_tf_turn_insulation = None module-attribute

conduit insulation thickness (m)

layer_ins = None module-attribute

Additional insulation thickness between layers (m)

dr_tf_nose_case = None module-attribute

inboard TF coil case outer (non-plasma side) thickness (m) (iteration variable 57) (calculated for stellarators)

dr_tf_wp_with_insulation = None module-attribute

radial thickness of winding pack (m) (iteration variable 140) (issue #514)

dx_tf_turn_steel = None module-attribute

TF coil turn steel conduit case thickness (m) (iteration variable 58)

dx_tf_wp_insulation = None module-attribute

Thickness of the ground insulation layer surrounding (m) - Superconductor TF (i_tf_sup == 1) : The TF coil Winding packs - Resistive magnets (i_tf_sup /= 1) : The TF coil wedges Rem : Thickness calculated for stellarators.

temp_tf_superconductor_margin_min = None module-attribute

minimum allowable temperature margin : TF coils (K)

temp_cs_superconductor_margin_min = None module-attribute

minimum allowable temperature margin : CS (K)

tmargmin = None module-attribute

minimum allowable temperature margin : TFC AND CS (K)

temp_margin = None module-attribute

temperature margin (K)

temp_tf_superconductor_margin = None module-attribute

TF coil superconductor temperature margin (K)

temp_tf_conductor_quench_max = None module-attribute

maximum temp during a quench for protection (K)

temp_croco_quench_max = None module-attribute

CroCo strand: maximum permitted temp during a quench (K)

temp_croco_quench = None module-attribute

CroCo strand: Actual temp reached during a quench (K)

temp_tf_cryo = None module-attribute

coil temperature for cryogenic plant power calculation (K)

n_tf_coil_turns = None module-attribute

number of turns per TF coil

v_tf_coil_dump_quench_max_kv = None module-attribute

max voltage across TF coil during quench (kV) (iteration variable 52)

vforce = None module-attribute

vertical tension on inboard leg/coil (N)

f_vforce_inboard = None module-attribute

Fraction of the total vertical force taken by the TF inboard leg tension Not used for resistive itart=1 (sliding joints)

vforce_outboard = None module-attribute

Vertical tension on outboard leg/coil (N)

f_a_tf_turn_cable_space_extra_void = None module-attribute

coolant fraction of TFC 'cable' (i_tf_sup=1), or of TFC leg (i_tf_ssup=0)

voltfleg = None module-attribute

volume of each TF coil outboard leg (m3)

vtfkv = None module-attribute

TF coil voltage for resistive coil including bus (kV)

v_tf_coil_dump_quench_kv = None module-attribute

voltage across a TF coil during quench (kV)

m_tf_coil_case = None module-attribute

mass per coil of external case (kg)

m_tf_coil_conductor = None module-attribute

TF coil conductor mass per coil (kg/coil). For itart=1, coil is return limb plus centrepost/n_tf_coils

m_tf_coil_copper = None module-attribute

copper mass in TF coil conductor (kg/coil). For itart=1, coil is return limb plus centrepost/n_tf_coils

whtconal = None module-attribute

Aluminium mass in TF coil conductor (kg/coil). For itart=1, coil is return limb plus centrepost/n_tf_coils

m_tf_coil_wp_turn_insulation = None module-attribute

conduit insulation mass in TF coil conductor (kg/coil)

m_tf_coil_superconductor = None module-attribute

superconductor mass in TF coil cable (kg/coil)

m_tf_wp_steel_conduit = None module-attribute

steel conduit mass in TF coil conductor (kg/coil)

m_tf_coil_wp_insulation = None module-attribute

mass of ground-wall insulation layer per coil (kg/coil)

m_tf_coils_total = None module-attribute

total mass of the TF coils (kg)

dx_tf_wp_primary_toroidal = None module-attribute

width of first step of winding pack (m)

dx_tf_wp_secondary_toroidal = None module-attribute

width of second step of winding pack (m)

dthet = None module-attribute

angle of arc i (rad)

radctf = None module-attribute

radius of arc i (m)

r_tf_arc = None module-attribute

x location of arc point i on surface (m)

xctfc = None module-attribute

x location of arc centre i (m)

z_tf_arc = None module-attribute

y location of arc point i on surface (m)

yctfc = None module-attribute

y location of arc centre i (m)

tfa = None module-attribute

Horizontal radius of inside edge of TF coil (m)

tfb = None module-attribute

Vertical radius of inside edge of TF coil (m)

drtop = None module-attribute

centrepost taper maximum radius adjustment (m)

dztop = None module-attribute

centrepost taper height adjustment (m)

etapump = None module-attribute

centrepost coolant pump efficiency

fcoolcp = None module-attribute

coolant fraction of TF coil inboard legs (iteration variable 23)

f_a_tf_cool_outboard = None module-attribute

coolant fraction of TF coil outboard legs

a_cp_cool = None module-attribute

Centrepost cooling area toroidal cross-section (constant over the whole CP)

n_cp_coolant_channels_total = None module-attribute

number of centrepost coolant tubes

p_cp_coolant_pump_elec = None module-attribute

centrepost coolant pump power (W)

p_cp_resistive = None module-attribute

resistive power in the centrepost (itart=1) [W]. If itart=0, this variable is the ressitive power on the whole magnet

p_tf_leg_resistive = None module-attribute

Summed resistive power in the TF coil legs [W]. Remain 0 if itart=0.

temp_cp_max = None module-attribute

maximum peak centrepost temperature (K) (constraint equation 44)

radius_cp_coolant_channel = None module-attribute

average radius of coolant channel (m) (iteration variable 69)

temp_cp_coolant_inlet = None module-attribute

centrepost coolant inlet temperature (K)

dtemp_cp_coolant = None module-attribute

inlet / outlet TF coil coolant temperature rise (K)

temp_cp_average = None module-attribute

Average temperature of centrepost called CP (K). Only used for resistive coils to compute the resisitive heating. Must be an iteration variable for ST (itart=1) (iteration variable 20)

tcpav2 = None module-attribute

Computed centrepost average temperature (K) (for consistency)

temp_tf_legs_outboard = None module-attribute

Average temperature of the TF outboard legs [K]. If temp_tf_legs_outboard=-1.0, the ouboard legs and CP temperatures are the same. Fixed for now, should use a contraints eq like temp_cp_average

temp_cp_peak = None module-attribute

peak centrepost temperature (K)

vel_cp_coolant_midplane = None module-attribute

inlet centrepost coolant flow speed at midplane (m/s) (iteration variable 70)

vol_cond_cp = None module-attribute

Exact conductor volume in the centrepost (m3)

whtcp = None module-attribute

mass of TF coil inboard legs (kg)

whttflgs = None module-attribute

mass of the TF coil legs (kg)

cryo_cool_req = None module-attribute

Cryo cooling requirement at helium temp 4.5K (kW)

theta1_coil = None module-attribute

The angle of the outboard arc forming the TF coil current center line [deg]

theta1_vv = None module-attribute

The angle of the outboard arc forming the Vacuum Vessel current center line [deg]

max_vv_stress = None module-attribute

The allowable peak maximum shear stress in the vacuum vessel due to quench and fast discharge of the TF coils [Pa]

t_tf_quench_detection = None module-attribute

TF coil quench detection time (s). Only used for TF coil quench protection.

rrr_tf_cu = None module-attribute

TF coil copper residual-resistance-ratio (RRR). Only used for quench protection.

init_tfcoil_variables()

Source code in process/data_structure/tfcoil_variables.py
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369
1370
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
1411
1412
1413
1414
1415
1416
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
1503
1504
1505
1506
1507
1508
1509
1510
1511
1512
1513
1514
1515
1516
1517
1518
1519
1520
1521
1522
1523
1524
1525
1526
def init_tfcoil_variables():
    global \
        a_tf_coil_inboard_case, \
        a_tf_coil_outboard_case, \
        a_tf_turn_steel, \
        a_tf_wp_conductor, \
        a_res_tf_coil_conductor, \
        a_tf_turn_cable_space_no_void, \
        a_tf_turn_insulation, \
        a_tf_coil_wp_turn_insulation, \
        sig_tf_case_max, \
        sig_tf_wp_max, \
        a_tf_leg_outboard, \
        a_tf_wp_steel, \
        a_tf_wp_extra_void, \
        a_tf_wp_coolant_channels, \
        bcritsc, \
        b_tf_inboard_peak_symmetric, \
        b_tf_inboard_peak_with_ripple, \
        casestr, \
        dr_tf_plasma_case, \
        f_dr_tf_plasma_case, \
        i_f_dr_tf_plasma_case, \
        dx_tf_side_case_min, \
        dx_tf_side_case_peak, \
        casths_fraction, \
        tfc_sidewall_is_fraction, \
        t_conductor, \
        dx_tf_turn_general, \
        i_dx_tf_turn_general_input, \
        t_turn_tf_max, \
        dx_tf_turn_cable_space_general, \
        i_dx_tf_turn_cable_space_general_input, \
        acs, \
        cdtfleg, \
        cforce, \
        cplen, \
        c_tf_turn, \
        c_tf_turn_max, \
        den_tf_coil_case, \
        dcond, \
        den_tf_wp_turn_insulation, \
        dia_tf_turn_coolant_channel, \
        e_tf_magnetic_stored_total_gj, \
        e_tf_coil_magnetic_stored, \
        b_crit_upper_nbti, \
        t_crit_nbti, \
        max_force_density, \
        f_a_tf_turn_cable_copper, \
        fhts, \
        insstrain, \
        i_tf_stress_model, \
        i_tf_tresca, \
        i_tf_wp_geom, \
        i_tf_case_geom, \
        i_tf_turns_integer, \
        i_tf_sc_mat, \
        i_tf_sup, \
        i_tf_shape, \
        i_tf_cond_eyoung_axial, \
        i_tf_cond_eyoung_trans, \
        n_tf_wp_pancakes, \
        n_tf_wp_layers, \
        n_rad_per_layer, \
        i_tf_bucking, \
        n_tf_graded_layers, \
        n_tf_stress_layers, \
        n_tf_wp_stress_layers, \
        j_tf_bus, \
        j_crit_str_tf, \
        j_crit_str_0, \
        j_tf_wp_critical, \
        j_tf_wp_quench_heat_max, \
        j_tf_wp, \
        oacdcp, \
        eyoung_ins, \
        eyoung_steel, \
        eyoung_cond_axial, \
        eyoung_cond_trans, \
        eyoung_res_tf_buck, \
        eyoung_copper, \
        eyoung_al, \
        poisson_steel, \
        poisson_copper, \
        poisson_al, \
        poisson_ins, \
        poisson_cond_axial, \
        poisson_cond_trans, \
        r_b_tf_inboard_peak, \
        res_tf_leg, \
        toroidalgap, \
        ripple_b_tf_plasma_edge_max, \
        ripple_b_tf_plasma_edge, \
        c_tf_total, \
        radial_array, \
        sig_tf_r, \
        sig_tf_t, \
        deflect, \
        sig_tf_z, \
        sig_tf_vmises, \
        s_shear_tf, \
        sig_tf_cs_bucked, \
        sig_tf_case, \
        sig_tf_wp, \
        str_cs_con_res, \
        str_pf_con_res, \
        str_tf_con_res, \
        str_wp, \
        str_wp_max, \
        i_str_wp, \
        quench_model, \
        time1, \
        tcritsc, \
        t_tf_superconductor_quench, \
        a_tf_inboard_total, \
        len_tf_bus, \
        m_tf_bus, \
        tfckw, \
        tfcmw, \
        p_cp_resistive_mw, \
        p_tf_joints_resistive_mw, \
        tfcryoarea, \
        tficrn, \
        ind_tf_coil, \
        dx_tf_wp_insertion_gap, \
        p_tf_leg_resistive_mw, \
        rho_cp, \
        rho_tf_leg, \
        rho_tf_bus, \
        frhocp, \
        frholeg, \
        i_cp_joints, \
        rho_tf_joints, \
        n_tf_joints_contact, \
        n_tf_joints, \
        th_joint_contact, \
        p_tf_joints_resistive, \
        len_tf_coil, \
        eff_tf_cryo, \
        n_tf_coils, \
        tfocrn, \
        tfsai, \
        tfsao, \
        tftmp, \
        dx_tf_inboard_out_toroidal, \
        dx_tf_turn_insulation, \
        layer_ins, \
        dr_tf_nose_case, \
        dr_tf_wp_with_insulation, \
        dx_tf_turn_steel, \
        dx_tf_wp_insulation, \
        temp_tf_superconductor_margin_min, \
        temp_cs_superconductor_margin_min, \
        tmargmin, \
        temp_margin, \
        temp_tf_superconductor_margin, \
        temp_tf_conductor_quench_max, \
        temp_croco_quench_max, \
        temp_croco_quench, \
        temp_tf_cryo, \
        n_tf_coil_turns, \
        v_tf_coil_dump_quench_max_kv, \
        vforce, \
        f_vforce_inboard, \
        vforce_outboard, \
        f_a_tf_turn_cable_space_extra_void, \
        voltfleg, \
        vtfkv, \
        v_tf_coil_dump_quench_kv, \
        m_tf_coil_case, \
        m_tf_coil_conductor, \
        m_tf_coil_copper, \
        whtconal, \
        m_tf_coil_wp_turn_insulation, \
        m_tf_coil_superconductor, \
        m_tf_wp_steel_conduit, \
        m_tf_coil_wp_insulation, \
        m_tf_coils_total, \
        dx_tf_wp_primary_toroidal, \
        dx_tf_wp_secondary_toroidal, \
        dthet, \
        radctf, \
        r_tf_arc, \
        xctfc, \
        z_tf_arc, \
        yctfc, \
        tfa, \
        tfb, \
        drtop, \
        dztop, \
        etapump, \
        fcoolcp, \
        f_a_tf_cool_outboard, \
        a_cp_cool, \
        n_cp_coolant_channels_total, \
        p_cp_coolant_pump_elec, \
        p_cp_resistive, \
        p_tf_leg_resistive, \
        temp_cp_max, \
        radius_cp_coolant_channel, \
        temp_cp_coolant_inlet, \
        dtemp_cp_coolant, \
        temp_cp_average, \
        tcpav2, \
        temp_tf_legs_outboard, \
        temp_cp_peak, \
        vel_cp_coolant_midplane, \
        vol_cond_cp, \
        whtcp, \
        whttflgs, \
        cryo_cool_req, \
        theta1_coil, \
        theta1_vv, \
        max_vv_stress, \
        t_tf_quench_detection, \
        rrr_tf_cu

    a_tf_coil_inboard_case = 0.0
    a_tf_coil_outboard_case = 0.0
    a_tf_turn_steel = 0.0
    a_tf_wp_conductor = 0.0
    a_res_tf_coil_conductor = 0.0
    a_tf_turn_cable_space_no_void = 0.0
    a_tf_turn_insulation = 0.0
    a_tf_coil_wp_turn_insulation = 0.0
    sig_tf_case_max = 6.0e8
    sig_tf_wp_max = 6.0e8
    a_tf_leg_outboard = 0.0
    a_tf_wp_steel = 0.0
    a_tf_wp_extra_void = 0.0
    a_tf_wp_coolant_channels = 0.0
    bcritsc = 24.0
    b_tf_inboard_peak_symmetric = 0.0
    b_tf_inboard_peak_with_ripple = 0.0
    casestr = 0.0
    dr_tf_plasma_case = 0.0
    f_dr_tf_plasma_case = 0.05
    i_f_dr_tf_plasma_case = False
    dx_tf_side_case_min = 0.0
    dx_tf_side_case_peak = 0.0
    casths_fraction = 0.06
    t_conductor = 0.0
    dx_tf_turn_cable_space_general = 0.0
    i_dx_tf_turn_cable_space_general_input = False
    dx_tf_turn_general = 0.0
    i_dx_tf_turn_general_input = False
    t_turn_tf_max = 0.05
    acs = 0.0
    cdtfleg = 0.0
    cforce = 0.0
    cplen = 0.0
    c_tf_turn = 7.0e4
    c_tf_turn_max = 9.0e4
    den_tf_coil_case = 8000.0
    dcond = np.array([
        6080.0,
        6080.0,
        6070.0,
        6080.0,
        6080.0,
        8500.0,
        6070.0,
        8500.0,
        8500.0,
    ])
    den_tf_wp_turn_insulation = 1800.0
    dia_tf_turn_coolant_channel = 0.005
    e_tf_magnetic_stored_total_gj = 0.0
    e_tf_coil_magnetic_stored = 0.0
    b_crit_upper_nbti = 14.86
    t_crit_nbti = 9.04
    max_force_density = 0.0
    f_a_tf_turn_cable_copper = 0.69
    fhts = 0.5
    insstrain = 0.0
    i_tf_stress_model = 1
    i_tf_tresca = 0
    i_tf_wp_geom = -1
    i_tf_case_geom = 0
    i_tf_turns_integer = 0
    i_tf_sc_mat = 1
    i_tf_sup = 1
    i_tf_shape = 0
    i_tf_cond_eyoung_axial = 0
    i_tf_cond_eyoung_trans = 1
    n_tf_wp_pancakes = 10
    n_tf_wp_layers = 20
    n_rad_per_layer = 100
    i_tf_bucking = -1
    n_tf_graded_layers = 1
    n_tf_stress_layers = 0
    n_tf_wp_stress_layers = 5
    j_tf_bus = 1.25e6
    j_crit_str_tf = 0.0
    j_crit_str_0 = np.array([
        596905475.80390120,
        1925501534.8512938,
        724544682.96063495,
        549858624.45072436,
        669284509.85818779,
        0.0,
        898964415.36996782,
        1158752995.2559297,
        865652122.9071957,
    ])
    j_tf_wp_critical = 0.0
    j_tf_wp_quench_heat_max = 0.0
    j_tf_wp = 0.0
    oacdcp = 0.0
    eyoung_ins = 1.0e8
    eyoung_steel = 2.05e11
    eyoung_cond_axial = 6.6e8
    eyoung_cond_trans = 0.0
    eyoung_res_tf_buck = 150.0e9
    eyoung_copper = 117.0e9
    eyoung_al = 69.0e9
    poisson_steel = 0.3
    poisson_copper = 0.35
    poisson_al = 0.35
    poisson_ins = 0.34
    poisson_cond_axial = 0.3
    poisson_cond_trans = 0.3
    r_b_tf_inboard_peak = 0.0
    res_tf_leg = 0.0
    toroidalgap = 1.0  # [m]
    ripple_b_tf_plasma_edge_max = 1.0
    ripple_b_tf_plasma_edge = 0.0
    c_tf_total = 0.0
    radial_array = np.zeros(N_RADIAL_ARRAY)
    sig_tf_r = np.zeros(N_RADIAL_ARRAY)
    sig_tf_t = np.zeros(N_RADIAL_ARRAY)
    deflect = np.zeros(N_RADIAL_ARRAY)
    sig_tf_z = 0.0
    sig_tf_vmises = np.zeros(N_RADIAL_ARRAY)
    s_shear_tf = np.zeros(N_RADIAL_ARRAY)
    sig_tf_cs_bucked = 0.0
    sig_tf_case = 0.0
    sig_tf_wp = 0.0
    str_cs_con_res = -0.005
    str_pf_con_res = -0.005
    str_tf_con_res = -0.005
    str_wp = 0.0
    str_wp_max = 0.7e-2
    i_str_wp = 1
    quench_model = "exponential"
    time1 = 0
    tcritsc = 16.0
    t_tf_superconductor_quench = 10.0
    a_tf_inboard_total = 0.0
    len_tf_bus = 300.0
    m_tf_bus = 0.0
    tfckw = 0.0
    tfcmw = 0.0
    p_cp_resistive_mw = 0.0
    p_tf_joints_resistive_mw = 0.0
    tfcryoarea = 0.0
    tficrn = 0.0
    ind_tf_coil = 0.0
    dx_tf_wp_insertion_gap = 0.01
    p_tf_leg_resistive_mw = 0.0
    rho_cp = 0.0
    rho_tf_leg = 0.0
    rho_tf_bus = 1.86e-8
    frhocp = 1.0
    frholeg = 1.0
    rho_tf_joints = 2.5e-10
    n_tf_joints_contact = 6
    n_tf_joints = 4
    th_joint_contact = 0.03
    p_tf_joints_resistive = 0.0
    len_tf_coil = 0.0
    eff_tf_cryo = -1.0
    n_tf_coils = 16.0
    tfocrn = 0.0
    tfsai = 0.0
    tfsao = 0.0
    tftmp = 4.5
    dx_tf_inboard_out_toroidal = 1.0
    dx_tf_turn_insulation = 8e-4
    layer_ins = 0.0
    dr_tf_nose_case = 0.3
    dr_tf_wp_with_insulation = 0.0
    dx_tf_turn_steel = 8e-3
    dx_tf_wp_insulation = 0.018
    temp_tf_superconductor_margin_min = 0.0
    temp_cs_superconductor_margin_min = 0.0
    tmargmin = 0.0
    temp_margin = 0.0
    temp_tf_superconductor_margin = 0.0
    temp_tf_conductor_quench_max = 150.0
    temp_croco_quench_max = 200.0
    temp_croco_quench = 0.0
    temp_tf_cryo = 4.5
    n_tf_coil_turns = 0.0
    v_tf_coil_dump_quench_max_kv = 20.0
    vforce = 0.0
    f_vforce_inboard = 0.5
    vforce_outboard = 0.0
    f_a_tf_turn_cable_space_extra_void = 0.4
    voltfleg = 0.0
    vtfkv = 0.0
    v_tf_coil_dump_quench_kv = 0.0
    m_tf_coil_case = 0.0
    m_tf_coil_conductor = 0.0
    m_tf_coil_copper = 0.0
    whtconal = 0.0
    m_tf_coil_wp_turn_insulation = 0.0
    m_tf_coil_superconductor = 0.0
    m_tf_wp_steel_conduit = 0.0
    m_tf_coil_wp_insulation = 0.0
    m_tf_coils_total = 0.0
    dx_tf_wp_primary_toroidal = 0.0
    dx_tf_wp_secondary_toroidal = 0.0
    dthet = np.zeros(4)
    radctf = np.zeros(4)
    r_tf_arc = np.zeros(5)
    xctfc = np.zeros(4)
    z_tf_arc = np.zeros(5)
    yctfc = np.zeros(4)
    tfa = np.zeros(4)
    tfb = np.zeros(4)
    drtop = 0.0
    dztop = 0.0
    etapump = 0.8
    fcoolcp = 0.3
    f_a_tf_cool_outboard = 0.2
    a_cp_cool = 0.0
    n_cp_coolant_channels_total = 0.0
    p_cp_coolant_pump_elec = 0.0
    p_cp_resistive = 0.0
    p_tf_leg_resistive = 0.0
    temp_cp_max = 473.15  # 200 C
    radius_cp_coolant_channel = 0.005
    temp_cp_coolant_inlet = 313.15  # 40 C
    dtemp_cp_coolant = 0.0
    temp_cp_average = 373.15  # 100 C
    tcpav2 = 0.0
    temp_tf_legs_outboard = -1.0
    temp_cp_peak = 0.0
    vel_cp_coolant_midplane = 20.0
    vol_cond_cp = 0.0
    whtcp = 0.0
    whttflgs = 0.0
    tfc_sidewall_is_fraction = False
    i_cp_joints = -1
    cryo_cool_req = 0.0
    theta1_coil = 45.0
    theta1_vv = 1.0  # 1 Deg
    max_vv_stress = 143.0e6
    rrr_tf_cu = 100.0
    t_tf_quench_detection = 3.0