🧲 Refactor croco 3

documentation/scripts/plotting_scripts/2nd_gen_rebco_critical_surface.py

@@ -22,6 +22,8 @@
         (
             j_scaling[i, j],
             _,
+            _,
+            _,
         ) = superconductors.jcrit_rebco(temp_grid[i, j], b_grid[i, j])
         # Convert from A/m² to kA/mm² (1 A/m² = 1e-6 A/mm²)
         j_scaling[i, j] *= 1e-9

documentation/source/eng-models/superconductors.md

@@ -383,7 +383,7 @@ The geometry of a single CroCo cable is calculated as follows:
 1. The diameter of the circular internal tape region is given by the outer copper diameter minus its thickness:
 
     $$
-    \overbrace{D_{\text{cable,internal}}}^{\texttt{dia_croco_strand_tape_region}} = \overbrace{D_{\text{cable}}}^{\texttt{dia_croco_strand}} - \overbrace{dx_{\text{cable,copper}}}^{\texttt{dx_croco_strand_copper}}
+    \overbrace{D_{\text{cable,internal}}}^{\texttt{dia_croco_strand_tape_region}} = \overbrace{D_{\text{cable}}}^{\texttt{dia_croco_strand}} - 2 \times \overbrace{dx_{\text{cable,copper}}}^{\texttt{dx_croco_strand_copper}}
     $$
 
 2. The total thickness of the HTS tape is found:

documentation/source/eng-models/tf-coil-superconducting.md

@@ -652,6 +652,29 @@ WIP
 
 WIP
 
-### Superconductor properties | `supercon_croco()`
+### Cable space dimensions | `tf_turn_croco_cable_space_properties()`
+
+This function calculates the dimensions of the circular cable space in the CroCo turn and that of the cables required to properly pack into the cable space.
+
+The required diameter of a single cable element in the turn is given by:
+
+$$
+\overbrace{D_{\text{CroCo cable}}}^{\texttt{dia_tf_turn_croco_cable}} = \frac{dx_{\text{turn,jacket}}}{3} -\left(\overbrace{dx_{\text{steel}}}^{\texttt{dx_tf_turn_steel}} \times \frac{2}{3}\right)
+$$
+
+The full area of the circular cable space is:
+
+$$
+\overbrace{A_{\text{cable space, full}}}^{\texttt{a_tf_turn_cable_space_no_void}} = \frac{9\pi}{4} - D_{\text{CroCo cable}}^2
+$$
+
+The effective cable space is just the total cable space minus the central full copper cable in the middle. 
+
+$$
+\overbrace{A_{\text{cable space, eff}}}^{\texttt{a_tf_turn_cable_space_effective}} = A_{\text{cable space, full}} - \left(\frac{\pi}{4} \times D_{\text{CroCo cable}}^2\right)
+$$
+
+
+### Superconductor properties | `tf_croco_superconductor_properties()`
 
 ### Quench voltage | `croco_voltage()`

process/core/init.py

@@ -357,7 +357,7 @@ def check_process(inputs, data):  # noqa: ARG001
             data.impurity_radiation.f_nd_impurity_electrons[imp]
         )
 
-    # Stop the run if oacdcp is used as an optimisation variable
+    # Stop the run if j_tf_coil_full_area is used as an optimisation variable
     # As the current density is now calculated from b_plasma_toroidal_on_axis without constraint 10
 
     if (data_structure.numerics.ixc[: data_structure.numerics.nvar] == 12).any():

process/core/input.py

@@ -251,7 +251,7 @@ def __post_init__(self):
     "copper_rrr": InputVariable(
         data_structure.rebco_variables, float, range=(1.0, 10000.0)
     ),
-    "dx_hts_tape_copper": InputVariable(
+    "dx_tf_hts_tape_copper": InputVariable(
         data_structure.rebco_variables, float, range=(0.0, 0.001)
     ),
     "copperaoh_m2": InputVariable(
@@ -282,7 +282,7 @@ def __post_init__(self):
     "crane_arm_h": InputVariable("buildings", float, range=(1.0, 100.0)),
     "crane_clrnc_h": InputVariable("buildings", float, range=(0.0, 10.0)),
     "crane_clrnc_v": InputVariable("buildings", float, range=(0.0, 10.0)),
-    "dx_croco_strand_copper": InputVariable(
+    "dx_tf_croco_strand_copper": InputVariable(
         data_structure.rebco_variables, float, range=(0.001, 0.1)
     ),
     "cryomag_h": InputVariable("buildings", float, range=(1.0, 100.0)),
@@ -540,7 +540,7 @@ def __post_init__(self):
     "gas_buildings_w": InputVariable("buildings", float, range=(10.0, 1000.0)),
     "ground_clrnc": InputVariable("buildings", float, range=(0.0, 10.0)),
     "n_ecrh_harmonic": InputVariable("current_drive", float, range=(1.0, 10.0)),
-    "dx_hts_tape_hastelloy": InputVariable(
+    "dx_tf_hts_tape_hastelloy": InputVariable(
         data_structure.rebco_variables, float, range=(1e-08, 0.001)
     ),
     "hccl": InputVariable("buildings", float, range=(0.0, 10.0)),
@@ -639,7 +639,7 @@ def __post_init__(self):
     "nd_plasma_pedestal_electron": InputVariable("physics", float, range=(0.0, 1e21)),
     "nd_plasma_separatrix_electron": InputVariable("physics", float, range=(0.0, 1e21)),
     "nflutfmax": InputVariable("constraints", float, range=(0.0, 1e24)),
-    "oacdcp": InputVariable(
+    "j_tf_coil_full_area": InputVariable(
         data_structure.tfcoil_variables, float, range=(10000.0, 1000000000.0)
     ),
     "f_a_cs_turn_steel": InputVariable("pf_coil", float, range=(0.001, 0.999)),
@@ -733,7 +733,7 @@ def __post_init__(self):
     "reactor_hall_w": InputVariable("buildings", float, range=(10.0, 1000.0)),
     "reactor_roof_thk": InputVariable("buildings", float, range=(0.25, 25.0)),
     "reactor_wall_thk": InputVariable("buildings", float, range=(0.25, 25.0)),
-    "dx_hts_tape_rebco": InputVariable(
+    "dx_tf_hts_tape_rebco": InputVariable(
         data_structure.rebco_variables,
         float,
         range=(1e-08, 0.0001),
@@ -858,10 +858,10 @@ def __post_init__(self):
     "t_turn_tf_max": InputVariable(
         data_structure.tfcoil_variables, float, range=(0.0, 1.0)
     ),
-    "dx_hts_tape_total": InputVariable(
+    "dx_tf_hts_tape_total": InputVariable(
         data_structure.rebco_variables, float, range=(0.0, 0.1)
     ),
-    "dr_hts_tape": InputVariable(
+    "dr_tf_hts_tape": InputVariable(
         data_structure.rebco_variables, float, range=(0.0, 0.1)
     ),
     "tauee_in": InputVariable("physics", float, range=(0.0, 100.0)),
@@ -1072,8 +1072,8 @@ def __post_init__(self):
     "ccls_ma": InputVariable("pf_coil", float, array=True),
     "cfind": InputVariable("costs", float, array=True),
     "i_blkt_coolant_type": InputVariable("fwbs", int, choices=[1, 2]),
-    "coppera_m2_max": InputVariable(
-        data_structure.rebco_variables, float, range=(1.0e6, 1.0e10)
+    "tf_coppera_m2_max": InputVariable(
+        data_structure.superconducting_tf_coil_variables, float, range=(1.0e6, 1.0e10)
     ),
     "cost_model": InputVariable("costs", int, choices=[0, 1, 2]),
     "i_vac_pump_dwell": InputVariable("vacuum", int, choices=[0, 1, 2]),

process/core/io/obsolete_vars.py

@@ -420,12 +420,12 @@
     "t_structural_vertical": "dz_cs_turn_conduit",
     "t_cable_tf": "dx_tf_turn_cable_space_general",
     "t_turn_tf": "dx_tf_turn_general",
-    "copper_thick": "dx_hts_tape_copper",
-    "croco_thick": "dx_croco_strand_copper",
-    "hastelloy_thickness": "dx_hts_tape_hastelloy",
-    "rebco_thickness": "dx_hts_tape_rebco",
-    "tape_thickness": "dx_hts_tape_total",
-    "tape_width": "dr_hts_tape",
+    "copper_thick": "dx_tf_hts_tape_copper",
+    "croco_thick": "dx_tf_croco_strand_copper",
+    "hastelloy_thickness": "dx_tf_hts_tape_hastelloy",
+    "rebco_thickness": "dx_tf_hts_tape_rebco",
+    "tape_thickness": "dx_tf_hts_tape_total",
+    "tape_width": "dr_tf_hts_tape",
     "beta": "beta_total_vol_avg",
     "beta_max": "beta_vol_avg_max",
     "beta_min": "beta_vol_avg_min",
@@ -448,6 +448,14 @@
     "vcool": "vel_cp_coolant_midplane",
     "rcool": "radius_cp_coolant_channel",
     "fl_h_threshold": None,
+    "dx_hts_tape_copper": "dx_tf_hts_tape_copper",
+    "dx_croco_strand_copper": "dx_tf_croco_strand_copper",
+    "dx_hts_tape_hastelloy": "dx_tf_hts_tape_hastelloy",
+    "dx_hts_tape_rebco": "dx_tf_hts_tape_rebco",
+    "oacdcp": "j_tf_coil_full_area",
+    "dx_hts_tape_total": "dx_tf_hts_tape_total",
+    "dr_hts_tape": "dr_tf_hts_tape",
+    "coppera_m2_max": "tf_coppera_m2_max",
 }
 
 OBS_VARS_HELP = {

process/core/io/plot/scans.py

@@ -82,7 +82,7 @@ def plot_scan(
         2: "pflux_div_heat_load_max_mw",
         3: "p_plant_electric_net_mw",
         4: "hfact",
-        5: "oacdcp",
+        5: "j_tf_coil_full_area",
         6: "pflux_fw_neutron_max_mw",
         7: "beamfus0",
         8: "Obsolete",  # OBSOLETE