fiqus 2025.11.0__py3-none-any.whl → 2026.1.0__py3-none-any.whl

fiqus/MainFiQuS.py

@@ -25,6 +25,9 @@ from fiqus.mains.MainCCT import MainCCT
 from fiqus.mains.MainMultipole import MainMultipole
 from fiqus.mains.MainPancake3D import MainPancake3D
 from fiqus.mains.MainConductorAC_Strand import MainConductorAC_Strand
+from fiqus.mains.MainHomogenizedConductor import MainHomogenizedConductor
+from fiqus.mains.MainConductorAC_Rutherford import MainConductorAC_Rutherford
+from fiqus.mains.MainConductorAC_CC import MainConductorAC_CC
 
 class MainFiQuS:
     """
@@ -90,6 +93,12 @@ class MainFiQuS:
             self.main_magnet = MainPancake3D(fdm=self.fdm, verbose=verbose)
         elif self.fdm.magnet.type == "CACStrand":
             self.main_magnet = MainConductorAC_Strand(fdm=self.fdm, inputs_folder_path=pathlib.Path(input_file_path).parent, outputs_folder_path=model_folder, verbose=verbose)
+        elif self.fdm.magnet.type == "CACCC":
+            self.main_magnet = MainConductorAC_CC(fdm=self.fdm, inputs_folder_path=pathlib.Path(input_file_path).parent, verbose=verbose)
+        elif self.fdm.magnet.type == "HomogenizedConductor":
+            self.main_magnet = MainHomogenizedConductor(fdm=self.fdm, inputs_folder_path=pathlib.Path(input_file_path).parent, outputs_folder_path=model_folder, verbose=verbose)
+        elif self.fdm.magnet.type == "CACRutherford":
+            self.main_magnet = MainConductorAC_Rutherford(fdm=self.fdm, inputs_folder_path=pathlib.Path(input_file_path).parent, verbose=verbose)
         elif self.fdm.magnet.type == "multipole":
             self.file_name = os.path.basename(input_file_path)[:-5]
             if not self.fdm.magnet.geometry.geom_file_path:

fiqus/data/DataConductor.py

@@ -116,6 +116,27 @@ class ProDefined(BaseModel):
     v: Optional[float] = None  # [-]
     B0: Optional[float] = None  # [-]
 
+class Succi_fixed(BaseModel):
+    """
+    Level 3: Class for cable Succi's YBCO fit
+    """
+
+    type: Literal["Succi_fixed"]
+    Jc_factor: Optional[float] = Field(gt=0.0,
+        default = 1.0,
+        description = "This factor multiplies the Jc returned by the function."
+    )
+    # all other parameters of the Succi fit are hardcoded
+class Fujikura(BaseModel):
+    """
+    Level 3: Class for cable Fujikura's fit
+    """
+
+    type: Literal["Fujikura"]
+    Jc_factor: Optional[float] = Field(gt=0.0,
+        default = 1.0,
+        description = "This factor multiplies the Jc returned by the function."
+    )
 
 class BSCCO_2212_LBNL(BaseModel):
     """
@@ -223,6 +244,31 @@ class Ribbon(BaseModel):
     material_outer_voids: Optional[str] = None
     material_core: Optional[str] = None
 
+class TSTC(BaseModel):
+    """
+    Twisted Stacked-Tape Cable (TSTC) type:
+    """
+    type: Literal["TSTC"]
+    stack_layout: Optional[List[Literal[-1,0,1]]] = Field(default=None,description="2D: Tape stack layout ordered TOP->BOTTOM. The numbers represent: 1 = a CC tape, -1 = a flipped CC tape, 0 = a shunt.")
+    nb_tapes: Optional[int] = Field(default=None, description="3D: Number of tapes in the stack")
+    tape_width: Optional[float] = Field(default=None, description="3D and 2D: Width of each tape")
+    tape_thickness: Optional[float] = Field(default=None, description="3D and 2D: Thickness of each tape")
+    twist_pitch: Optional[float] = Field(default=None, description="3D: Length over which tapes are twisted by full rotation")
+    pitch_fraction: Optional[float] = Field(default=1.0, description="3D: Fraction of the twist pitch to be modelled (1.0 = full pitch, 0.5 = half pitch, etc.)")
+    bare_cable_width: Optional[float] = Field(default=None, description="Cable width, typically the same as CC width")
+    bare_cable_height_low: Optional[float] = Field(default=None, description="Narrow end (if applicable) cable height (thickness), typically the same as (CC thickness + Cu stabilizer thickness) * number of tapes.")
+    bare_cable_height_high: Optional[float] = Field(default=None, description="Wide end (if applicable) cable height (thickness), typically the same as (CC thickness + Cu stabilizer thickness) * number of tapes.")
+    bare_cable_height_mean: Optional[float] = Field(default=None, description="Average (if applicable) cable height (thickness), typically the same as (CC thickness + Cu stabilizer thickness) * number of tapes.")
+    th_insulation_along_width: Optional[float] = Field(default=None, description="Insulation thickness along the width ")
+    th_insulation_along_height: Optional[float] = Field(default=None, description="Insulation thickness along the height ")
+    f_superconductor: Optional[float] = Field(default=None, description="Fraction of superconductor related to the total area of the cable (winding cell)")
+    f_stabilizer: Optional[float] = Field(default=None, description="Fraction of stabilizer related to the total area of the cable (winding cell)")
+    f_silver: Optional[float] = Field(default=None, description="Fraction of silver related to the total area of the cable (winding cell)")
+    f_substrate: Optional[float] = Field(default=None, description="Fraction of substrate (including buffer layers and silver overlay) related to the total area of the cable (winding cell)")
+    f_shunt: Optional[float] = Field(default=None, description="Fraction of substrate (including buffer layers and silver overlay) related to the total area of the cable (winding cell)")
+    f_insulation: Optional[float] = Field(default=None, description="Fraction of cable insulation related to the total area of the cable (winding cell)")
+    f_inner_voids: Optional[float] = Field(default=None, description="Fraction of additional material (typically insulation) related to the total area of the cable (winding cell)")
+    f_outer_voids: Optional[float] = Field(default=None, description="Fraction of additional material (typically helium impregnating the windings) related to the total area of the cable (winding cell)")
 
 # ------------------- Conductors ---------------------------#
 
@@ -328,6 +374,69 @@ class Rectangular(BaseModel):
     # superconductor: MaterialSuperconductor = MaterialSuperconductor()
     # stabilizer: MaterialStabilizer = MaterialStabilizer()
     
+class Copper_thickness(BaseModel):
+    left: Optional[float] = Field(default=None, description="On the left side.")
+    right: Optional[float] = Field(default=None, description="On the right side.")
+    top: Optional[float] = Field(default=None, description="On the top side.")
+    bottom: Optional[float] = Field(default=None, description="On the bottom side.")
+
+class Silver_thickness(BaseModel):
+    top: Optional[float] = Field(default=None, description="On the top side.")
+    bottom: Optional[float] = Field(default=None, description="On the bottom side.")
+
+class CC(BaseModel):
+    """
+    Level 2: Class for coated conductor parameters
+    """
+    type: Literal["CC"]
+
+    # Core layer sizes
+    HTS_thickness: Optional[float] = Field(default=None, description="HTS thickness in meters.")
+    HTS_width: Optional[float] = Field(default=None, description="HTS width in meters.")
+    number_of_filaments: Optional[int] = Field(default=1, ge=1, description="Number of HTS filaments. If 1, no striation case")
+    gap_between_filaments: Optional[float] = Field(default=None, description="Gap between HTS filaments in meters. Only applies when number_of_filaments > 1.")
+    substrate_thickness: Optional[float] = Field(default=None, description="Substrate layer thickness in meters.")
+
+    # Plating/stabilizer
+    copper_thickness: Copper_thickness = Field(default=Copper_thickness(), description="Copper thickness in meters")
+    silver_thickness: Silver_thickness = Field(default=Silver_thickness(), description="Silver thickness in meters")
+
+    # -- Superconductor parameters -- #
+    material_superconductor: Optional[str] = Field(default=None, description="Material of the superconductor. E.g. NbTi, Nb3Sn, etc.")
+    n_value_superconductor: Optional[float] = Field(default=None, description="n value of the superconductor (for power law fit).")
+    ec_superconductor: Optional[float] = Field(default=None, description="Critical electric field of the superconductor.")
+    minimum_jc_fraction: Optional[float] = Field(gt=0, le=1, default=None, description="Fraction of Jc(minimum_jc_field, T) to use as minimum Jc for the power law"
+                                                                                       " fit to avoid division by zero when Jc(B_local, T) decreases to zero."
+                                                                           "Typical value would be 0.001 (so the Jc_minimum is 0.1% of Jc(minimum_jc_field, T))"
+                                                                            "This fraction is only allowed to be greater than 0.0 and less than or equal to 1.0")
+    minimum_jc_field: Optional[float] = Field(default=None, description="Magnetic flux density in tesla used for calculation of Jc(minimum_jc_field, T)."
+                                                                        "This gets multiplied by minimum_jc_fraction and used as minimum Jc for the power law")
+    k_material_superconductor: Optional[Union[str, float]] = Field(default=None, description="Thermal conductivity of the superconductor.")
+    Cv_material_superconductor: Optional[Union[str, float]] = Field(default=None, description="Material function for specific heat of the superconductor.")
+    k_material_stabilizer: Optional[Union[str, float]] = Field(default=None, description="Thermal conductivity of the stabilizer, typically copper.")
+    Cv_material_stabilizer: Optional[Union[str, float]] = Field(default=None, description="Material function for specific heat of the stabilizer, typically copper.")
+    rho_material_stabilizer: Optional[Union[str, float]] = Field(default=None, description="Material function for resistivity of the stabilizer. Constant resistivity can be given as float.")
+    RRR: Optional[Union[float, List[float]]] = Field(default=None, description="Residual resistivity ratio of the stabilizer. If a list of RRR is provided it needs to match in length the number of matrix regions in the geometry (typically 3)")
+    T_ref_RRR_high: Optional[float] = Field(default=None, description="Upper reference temperature for RRR measurements.")
+    T_ref_RRR_low: Optional[float] = Field(default=None, description="Lower reference temperature for RRR measurements.")
+    k_material_silver: Optional[Union[str, float]] = Field(default=None, description="Thermal conductivity of the silver")
+    Cv_material_silver: Optional[Union[str, float]] = Field(default=None, description="Material function for specific heat of the silver")
+    rho_material_silver: Optional[Union[str, float]] = Field(default=None, description="Material function for resistivity of the silver. Constant resistivity can be given as float.")
+    RRR_silver: Optional[Union[float, List[float]]] = Field(default=None, description="Residual resistivity ratio of the silver. If a list of RRR is provided it needs to match in length the number of matrix regions in the geometry (typically 3)")
+    T_ref_RRR_high_silver: Optional[float] = Field(default=None, description="Upper reference temperature for RRR measurements for silver.")
+    T_ref_RRR_low_silver: Optional[float] = Field(default=None, description="Lower reference temperature for RRR measurements for silver.")
+    rho_material_substrate: Optional[Union[str, float]] = Field(default=None, description="Material function for resistivity of the substrate. Constant resistivity can be given as float.")
+    k_material_substrate: Optional[Union[str, float]] = Field(default=None, description="Thermal conductivity of the substrate.")
+    Cv_material_substrate: Optional[Union[str, float]] = Field(default=None, description="Material function for specific heat of the substrate.")
+
+class Homogenized(BaseModel):
+    """
+    Level 2: Class for homogenized strand parameters, to be used in the Rutherford cable model
+    """
+    type: Literal["Homogenized"]
+
+    # Strand diameter (used in the geometry step)
+    diameter: Optional[float] = Field(default=None, description="Undeformed round strand diameter. Used in the geometry step if keep_strand_area==true, the strand is deformed while preserving its surface area. Not used otherwise.")
 
 
 # ------------------- Conductors ---------------------------#
@@ -341,10 +450,10 @@ class Conductor(BaseModel):
     version: Optional[str] = None
     case: Optional[str] = None
     state: Optional[str] = None
-    cable: Union[Rutherford, Mono, Ribbon] = {
+    cable: Union[Rutherford, Mono, Ribbon, TSTC] = {
         "type": "Rutherford"
     }  # TODO: Busbar, Rope, Roebel, CORC, TSTC, CICC
-    strand: Union[Round, Rectangular] = {"type": "Round"}  # TODO: Tape, WIC
-    Jc_fit: Union[ConstantJc, Bottura, CUDI1, CUDI3, Summers, Bordini, Nb3Sn_HFM, BSCCO_2212_LBNL, Ic_A_NbTi, ProDefined] = {
+    strand: Union[Round, Rectangular, CC, Homogenized] = {"type": "Round"}  # TODO: Tape, WIC
+    Jc_fit: Union[ConstantJc, Bottura, CUDI1, CUDI3, Summers, Bordini, Nb3Sn_HFM, BSCCO_2212_LBNL, Ic_A_NbTi, ProDefined, Succi_fixed, Fujikura] = {
         "type": "CUDI1"
     }  # TODO: CUDI other numbers? , Roxie?

fiqus/data/DataFiQuS.py

@@ -7,8 +7,9 @@ from fiqus.data.DataFiQuSCCT import CCT
 from fiqus.data.DataFiQuSMultipole import Multipole
 from fiqus.data.DataFiQuSPancake3D import Pancake3D
 from fiqus.data.DataFiQuSConductorAC_Strand import CACStrand
-
-
+from fiqus.data.DataFiQuSHomogenizedConductor import HomogenizedConductor
+from fiqus.data.DataFiQuSConductorAC_Rutherford import CACRutherford
+from fiqus.data.DataFiQuSConductorAC_CC import CACCC
 
 class RunFiQuS(BaseModel):
     """
@@ -144,7 +145,7 @@ class FDM(BaseModel):
 
     general: GeneralFiQuS = GeneralFiQuS()
     run: RunFiQuS = RunFiQuS()
-    magnet: Union[Multipole, CCT, Pancake3D, CACStrand] = Field(
+    magnet: Union[Multipole, CCT, Pancake3D, CACStrand, HomogenizedConductor, CACRutherford, CACCC] = Field(
         default=Multipole(), discriminator="type"
     )
     circuit: Circuit_Class = Circuit_Class()

fiqus/data/DataFiQuSConductorAC_CC.py

@@ -0,0 +1,345 @@
+
+"""DataFiQuSConductorAC_CC.py"""
+
+from pydantic import BaseModel, Field, field_validator
+from typing import Literal, Optional, Union, List
+
+
+
+
+class CACCCGeometry(BaseModel):
+    """
+    Level 2: Geometry for CACCC.
+    """
+    air_radius: Optional[float] = Field(
+        default=None, 
+        description = "Radius of air region."
+    )
+
+
+class CACCCGeneralparameters(BaseModel):
+    """
+    Level 3: Class for general parameters
+    """
+    temperature: float = Field(
+        default=1.9, 
+        description = "Temperature (K)."
+    )
+    noOfMPITasks: Optional[Union[bool, int]] = Field(
+        default=False, title = "No. of tasks for MPI parallel run of GetDP",
+        description = "If integer, GetDP will be run in parallel using MPI. This is only valid"
+                      " if MPI is installed on the system and an MPI-enabled GetDP is used." 
+                      " If False, GetDP will be run in serial without invoking mpiexec."
+    )
+
+
+
+class CACCCSolveInitialconditions(BaseModel):
+    """
+    Level 3: Class for initial conditions
+    """
+    init_type: Optional[Literal['virgin', 'pos_file', 'uniform_field']] = Field(
+        default='virgin', 
+        description = "Type of initialization for the simulation. (i) 'virgin' is the default type, the initial magnetic field is zero,"
+                      "(ii) 'pos_file' is to initialize from the solution of another solution, given by the solution_to_init_from entry,"
+        " and (iii) 'uniform_field' is to initialize at a uniform field, which will be the applied field at the initial time of the simulation."
+        " Note that the uniform_field option does not allow any non-zero transport current (initialization from pos_file is needed for this)."
+    )
+    solution_to_init_from: Optional[Union[int, str]] = Field(
+        default=None,
+        description = "Name xxx of the solution from which the simulation should be initialized. " 
+                      "The file last_magnetic_field.pos of folder Solution_xxx will be used for the initial solution."
+                      "It must be in the Geometry_.../Mesh_.../ folder in which the Solution_xxx will be saved.",
+    )
+
+
+
+class CACCCSolveSourceparametersSine(BaseModel):
+    """
+    Level 4: Class for Sine source parameters
+    """
+    frequency: Optional[float] = Field(
+        default=None, 
+        description = "Frequency of the sine source (Hz)."
+    )
+    field_amplitude: Optional[float] = Field(
+        default=None, 
+        description = "Amplitude of the sine field (T)."
+    )
+    current_amplitude: Optional[float] = Field(
+        default=None, 
+        description = "Amplitude of the sine current (A)."
+    )
+
+
+
+class CACCCSolveSourceparametersPiecewise(BaseModel):
+    """
+    Level 4: Class for piecewise (linear) source parameters
+    """
+    source_csv_file: Optional[str] = Field(
+        default=None, 
+        description = "File name for the from_file source type defining the time evolution of current and field (in-phase)."
+                      "Multipliers are used for each of them." 
+                      "The file should contain two columns: 'time' (s) and 'value' (field/current (T/A)), with these headers."
+                      "If this field is set, times, applied_fields_relative and transport_currents_relative are ignored."
+    )
+    times: Optional[List[float]] = Field(
+        default=None, 
+        description = "Time instants (s) defining the piecewise linear sources." 
+                      "Used only if source_csv_file is not set." 
+                      "Can be scaled by time_multiplier."
+    )
+    applied_fields_relative: Optional[List[float]] = Field(
+        default=None, 
+        description = "Applied fields relative to multiplier applied_field_multiplier at the time instants 'times'."
+                      "Used only if source_csv_file is not set."
+    )
+    transport_currents_relative: Optional[List[float]] = Field(
+        default=None, 
+        description = "Transport currents relative to multiplier transport_current_multiplier at the time instants 'times'." 
+                      "Used only if source_csv_file is not set."
+    )
+    time_multiplier: Optional[float] = Field(
+        default=None, 
+        description = "Multiplier for the time values in times (scales the time values)." 
+                      "Also used for the time values in the source_csv_file."
+    )
+    applied_field_multiplier: Optional[float] = Field(
+        default=None, 
+        description = "Multiplier for the applied fields in applied_fields_relative."
+                      "Also used for the values in the source_csv_file."
+    )
+    transport_current_multiplier: Optional[float] = Field(
+        default=None, 
+        description = "Multiplier for the transport currents in transport_currents_relative."
+                      "Also used for the values in the source_csv_file."
+    )
+
+
+
+class CACCCSolveSourceparameters(BaseModel):
+    """
+    Level 3: Class for material properties
+    """
+    source_type: Literal['sine', 'piecewise'] = Field(
+        default='sine',
+        description = "Time evolution of applied current and magnetic field. " 
+                      "Supported options are: sine, piecewise."
+    )
+    sine: CACCCSolveSourceparametersSine = CACCCSolveSourceparametersSine()
+    piecewise: CACCCSolveSourceparametersPiecewise = CACCCSolveSourceparametersPiecewise()
+    field_angle_with_respect_to_normal_direction: Optional[float] = Field(
+        default=None, 
+        description = "Angle of the source magnetic field with respect to the y-axis (normal to the tape) (degrees)."
+    )
+
+
+
+class CACCCSolveNumericalparametersSine(BaseModel):
+    """ 
+    Level 4: Numerical parameters corresponding to the sine source
+    """
+    timesteps_per_period: Optional[float] = Field(
+        default=None, 
+        description = "Initial value for number of time steps (-) per period for the sine source." 
+                      "Determines the initial time step size."
+    )
+    number_of_periods_to_simulate: Optional[float] = Field(
+        default=None, 
+        description = "Number of periods (-) to simulate for the sine source."
+    )
+
+
+
+class CACCCSolveNumericalparametersPiecewise(BaseModel):
+    """
+    Level 4: Numerical parameters corresponding to the piecewise source
+    """
+    time_to_simulate: Optional[float] = Field(
+        default=None, 
+        description = "Total time to simulate (s). Used for the piecewise source."
+    )
+    timesteps_per_time_to_simulate: Optional[float] = Field(
+        default=None, 
+        description = "If variable_max_timestep is False. Number of time steps (-) per period for the piecewise source."
+    )
+    force_stepping_at_times_piecewise_linear: bool = Field(
+        default=False, 
+        description = "If True, time-stepping will contain exactly the time instants that are in" 
+                      "the times_source_piecewise_linear list (to avoid truncation maximum applied field/current values)."
+    )
+    variable_max_timestep: bool = Field(
+        default=False, 
+        description = "If False, the maximum time step is kept constant through the simulation. " 
+                      "If True, it varies according to the piecewise definition."
+    )
+    times_max_timestep_piecewise_linear: Optional[List[float]] = Field(
+        default=None, 
+        description = "Time instants (s) defining the piecewise linear maximum time step."
+    )
+    max_timestep_piecewise_linear: Optional[List[float]] = Field(
+        default=None, 
+        description = "Maximum time steps (s) at the times_max_timestep_piecewise_linear. " 
+                      "Above the limits, linear extrapolation of the last two values."
+    )
+
+
+
+class CACCCSolveNumericalparameters(BaseModel):
+    """
+    Level 3: Class for numerical parameters
+    """
+    relative_tolerance: Optional[float] = Field(default=1e-6, description="Tolerance on the relative change of the power indicator for the convergence criterion (1e-6 is usually a safe choice).")
+    voltage_per_meter_stopping_criterion: Optional[float] = Field(default=None, description="If a non-zero value is given, the simulation will stop if the transport voltage per meter reaches this value (in absolute value).")
+    relaxation_factors: Optional[bool] = Field(default=True, description="Use of relaxation factors to help convergence (automatic selection based on the lowest residual).")
+    sine: CACCCSolveNumericalparametersSine = CACCCSolveNumericalparametersSine()
+    piecewise: CACCCSolveNumericalparametersPiecewise = CACCCSolveNumericalparametersPiecewise()
+
+
+
+class CACCCSolve(BaseModel):
+    """
+    Level 2: Solve block for CACCC
+    """
+    pro_template: Optional[str] = Field(
+        default='CAC_CC_template.pro', 
+        description = "Name of the .pro template file."
+    )
+    conductor_name: Optional[str] = Field(
+        default=None, 
+        description = "Name of the conductor. Must match a conductor name in "
+                      "the conductors section of the input YAML-file."
+    )
+    general_parameters: CACCCGeneralparameters = (
+        CACCCGeneralparameters()
+    )
+    initial_conditions: CACCCSolveInitialconditions = (
+        CACCCSolveInitialconditions()
+    )
+    source_parameters: CACCCSolveSourceparameters = (
+        CACCCSolveSourceparameters()
+    )
+    numerical_parameters: CACCCSolveNumericalparameters = (
+        CACCCSolveNumericalparameters()
+    )
+
+
+
+class CACCCMesh(BaseModel):
+    """
+    Level 2: Mesh parameters for CACCC.
+    """
+    HTS_n_elem_width: Optional[int] = Field(
+        default=None,
+        description = "Number of elements along HTS width (x-direction)."
+    )
+    HTS_n_elem_thickness: Optional[int] = Field(
+        default=None,
+        description = "Number of elements through HTS thickness (y-direction)."
+    )
+    substrate_elem_scale: Optional[float] = Field(
+        default=None,
+        description = "Element-count scale factor for substrate layer."
+    )
+    substrate_side_progression: Optional[float] = Field(
+        default=None,
+        description = "Progression factor for substrate vertical sides near the HTS side."
+    )
+    silver_elem_scale: Optional[float] = Field(
+        default=None,
+        description = "Element-count scale factor for silver layers."
+    )
+    copper_elem_scale: Optional[float] = Field(
+        default=None,
+        description = "Element-count scale factor for copper layers."
+    )
+    air_boundary_mesh_size_ratio: Optional[float] = Field(
+        default=None,
+        description = "Ratio of air outer-boundary mesh size to the HTS base size."
+    )
+    scaling_global: Optional[float] = Field(
+        default=None,
+        description = "Global refinement factor."
+    )
+    bump_coef: Optional[float] = Field(
+        default=None,
+        description = "Unified bump coefficient for transfinite horizontal edges. "
+                      "Used for both HTS and SilverTop when applying 'Bump' distributions. "
+                      "Values < 1 cluster nodes toward the edges; values > 1 cluster toward the center."
+    )
+
+
+
+
+class CACCCPostprocPosFiles(BaseModel):
+    """
+    Level 3: Class for post-pro .pos file requests
+    """
+    quantities: Optional[List[str]] = Field(
+        default=None, description = "List of GetDP postprocessing quantities to write to .pos file. "
+                                    "Examples of valid entry is: phi, h, b, j, jz, power"
+    )
+    regions: Optional[List[str]] = Field(
+        default=None, 
+        description = "List of GetDP regions to write to .pos file postprocessing for. " 
+                      "Examples of a valid entry is: " 
+                      "Matrix, Filaments, Omega (full domain), OmegaC (conducting domain), OmegaCC (non conducting domain)"
+    )
+
+    @field_validator("quantities", 'regions', mode="before")
+    @classmethod
+    def convert_none_to_empty_list(cls, v):
+        pass
+        if v is None:
+            return []   # if None default to empty list. Jinja relies on list for looping in the pro template. The default is not a list to be consistent with SDK entries.
+        return v
+
+    model_config = {
+        "validate_default": True
+    }
+
+
+
+class CWSStrandPostprocCleanup(BaseModel):
+    """
+    Level 3: Class for cleanup settings
+    """
+    remove_pre_file: bool = Field(
+        default=False,
+        description = "Set True to remove the .pre-file after post-processing, to save disk space."
+    )
+    remove_res_file: bool = Field(
+        default=False,
+        description = "Set True to remove the .res-file after post-processing, to save disk space."
+    )
+    remove_msh_file: bool = Field(
+        default=False,
+        description = "Set True to remove the .msh-file after post-processing, to save disk space."
+    )
+
+
+
+class CACCCPostproc(BaseModel):
+    """
+    Post-processing options for CACCC
+    """
+    pos_files: CACCCPostprocPosFiles = Field(
+        default=CACCCPostprocPosFiles(),
+        description = "Entries controlling output of .pos files." 
+                      "If None or empty lists are given, no .pos files are written." 
+                      "Note that not all combinations of quantities and regions make sense."
+    )
+    cleanup: CWSStrandPostprocCleanup = CWSStrandPostprocCleanup()
+
+
+class CACCC(BaseModel):
+    """
+    Level 1: Class for FiQuS CACCC
+    """
+    type: Literal["CACCC"]
+    geometry: CACCCGeometry = CACCCGeometry()
+    mesh: CACCCMesh = CACCCMesh()
+    solve: CACCCSolve = CACCCSolve()
+    postproc: CACCCPostproc = CACCCPostproc()
+