Create specific turn tf classes

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

@@ -629,3 +629,29 @@ It is also worth mentioning that the the ripple must be evaluated layer-by-layer
 used to quantify the SC cross-section area per layer. Finally, resistive
 coils do not suffer from on-coil ripple as there is no radial case present.
 
+-----------------------------
+
+## Cable in Conduit TF coil class | `CICCSuperconductingTFCoil(SuperconductingTFCoil)`
+
+WIP
+
+### Averaged turn geometry | `tf_cable_in_conduit_averaged_turn_geometry()`
+
+### Integer turn geometry | `tf_cable_in_conduit_integer_turn_geometry()`
+
+### Superconductor length | `calculate_cable_in_conduit_superconductor_length()`
+
+### Superconductor strand count | `calculate_cable_in_conduit_strand_count()`
+
+### Superconductor properties | `tf_cable_in_conduit_superconductor_properties()`
+
+
+-----------------------------
+
+## Cross Conductor TF coil class | `CROCOSuperconductingTFCoil(SuperconductingTFCoil)`
+
+WIP
+
+### Superconductor properties | `supercon_croco()`
+
+### Quench voltage | `croco_voltage()`

process/core/caller.py

@@ -16,6 +16,7 @@
 from process.core.solver.iteration_variables import set_scaled_iteration_variable
 from process.core.solver.objectives import objective_function
 from process.models.tfcoil.base import TFConductorModel
+from process.models.tfcoil.superconducting import SuperconductingTFTurnType
 
 if TYPE_CHECKING:
     from process.core.model import DataStructure
@@ -296,7 +297,20 @@ def _call_models_once(self, xc: np.ndarray):
 
         # Toroidal field coil superconductor model
         if data_structure.tfcoil_variables.i_tf_sup == TFConductorModel.SUPERCONDUCTING:
-            self.models.sctfcoil.run()
+            if (
+                SuperconductingTFTurnType(
+                    data_structure.superconducting_tf_coil_variables.i_tf_turn_type
+                )
+                == SuperconductingTFTurnType.CABLE_IN_CONDUIT
+            ):
+                self.models.cicc_sctfcoil.run()
+            elif (
+                SuperconductingTFTurnType(
+                    data_structure.superconducting_tf_coil_variables.i_tf_turn_type
+                )
+                == SuperconductingTFTurnType.CROSS_CONDUCTOR
+            ):
+                self.models.croco_sctfcoil.run()
 
         if (
             data_structure.tfcoil_variables.i_tf_sup

process/core/input.py

@@ -1963,6 +1963,9 @@ def __post_init__(self):
     "n_tf_wp_layers": InputVariable(
         data_structure.tfcoil_variables, int, range=(1, 100)
     ),
+    "i_tf_turn_type": InputVariable(
+        data_structure.superconducting_tf_coil_variables, int, choices=[1, 2]
+    ),
     "n_liq_recirc": InputVariable(data_structure.fwbs_variables, int, range=(1, 50)),
     "n_tf_wp_pancakes": InputVariable(
         data_structure.tfcoil_variables, int, range=(1, 100)

process/core/output.py

@@ -1,6 +1,9 @@
 from process import data_structure
 from process.core.log import logging_model_handler
 from process.models.tfcoil.base import TFConductorModel
+from process.models.tfcoil.superconducting import (
+    SuperconductingTFTurnType,
+)
 
 
 def write(models, _outfile):
@@ -70,7 +73,18 @@ def write(models, _outfile):
 
     # Toroidal field coil superconductor model
     if data_structure.tfcoil_variables.i_tf_sup == TFConductorModel.SUPERCONDUCTING:
-        models.sctfcoil.output()
+        tf_turn_type = SuperconductingTFTurnType(
+            data_structure.superconducting_tf_coil_variables.i_tf_turn_type
+        )
+        if tf_turn_type == SuperconductingTFTurnType.CABLE_IN_CONDUIT:
+            models.cicc_sctfcoil.output()
+        elif tf_turn_type == SuperconductingTFTurnType.CROSS_CONDUCTOR:
+            models.croco_sctfcoil.output()
+        else:
+            raise ValueError(
+                "Unsupported superconducting TF turn type: "
+                f"{data_structure.superconducting_tf_coil_variables.i_tf_turn_type}"
+            )
 
     # Toroidal field coil aluminium model
     if (

process/data_structure/superconducting_tf_coil_variables.py

@@ -216,6 +216,9 @@
 j_tf_superconductor: float = None
 """Current density in the superconducting cable [A/m^2]"""
 
+i_tf_turn_type: int = None
+"""Switch for TF turn geometry type"""
+
 # Vacuum Vessel stress on TF coil quench
 
 vv_stress_quench: float = None
@@ -321,6 +324,7 @@ def init_superconducting_tf_coil_variables():
         f_a_tf_turn_cable_space_cooling, \
         c_tf_turn_cables_critical, \
         j_tf_superconductor, \
+        i_tf_turn_type, \
         vv_stress_quench
 
     is_leg_cp_temp_same = 0
@@ -382,4 +386,5 @@ def init_superconducting_tf_coil_variables():
     f_a_tf_turn_cable_space_cooling = 0.0
     c_tf_turn_cables_critical = 0.0
     j_tf_superconductor = 0.0
+    i_tf_turn_type = 1
     vv_stress_quench = 0.0

process/main.py

@@ -107,7 +107,11 @@
     CopperTFCoil,
     ResistiveTFCoil,
 )
-from process.models.tfcoil.superconducting import SuperconductingTFCoil
+from process.models.tfcoil.superconducting import (
+    CICCSuperconductingTFCoil,
+    CROCOSuperconductingTFCoil,
+    SuperconductingTFCoil,
+)
 from process.models.vacuum import Vacuum, VacuumVessel
 from process.models.water_use import WaterUse
 
@@ -602,6 +606,8 @@ def __init__(self, data: DataStructure):
         self.cryostat = Cryostat()
         self.build = Build()
         self.sctfcoil = SuperconductingTFCoil()
+        self.cicc_sctfcoil = CICCSuperconductingTFCoil()
+        self.croco_sctfcoil = CROCOSuperconductingTFCoil()
         self.tfcoil = TFCoil(build=self.build)
         self.resistive_tf_coil = ResistiveTFCoil()
         self.copper_tf_coil = CopperTFCoil()

process/models/tfcoil/base.py

@@ -565,6 +565,12 @@ def outtf(self):
             "(i_tf_sup)",
             tfcoil_variables.i_tf_sup,
         )
+        po.ovarin(
+            self.outfile,
+            "Superconducting TF coil turn type",
+            "(i_tf_turn_type)",
+            superconducting_tf_coil_variables.i_tf_turn_type,
+        )
 
         if tfcoil_variables.i_tf_sup == TFConductorModel.WATER_COOLED_COPPER:
             po.ocmmnt(