Merge pull request #845 from ukaea/adding_gap_between_divertor_and_tf

added gap size

Author: shimwell

docs/source/paramak.parametric_reactors.rst

@@ -19,7 +19,7 @@ are red, PF coil cases are yellow, the center column shielding is dark green,
 the blanket is light green, the divertor is orange, the firstwall is grey
 and the rear wall of the blanket is teal.
 
-.. image:: https://user-images.githubusercontent.com/8583900/116444066-41367180-a84c-11eb-8e1b-dfa4815f560b.png
+.. image:: https://user-images.githubusercontent.com/8583900/119011892-6458ca80-b98d-11eb-92b1-bbc370cb9c84.png
    :width: 450
    :align: center
 

paramak/parametric_reactors/ball_reactor.py

@@ -35,6 +35,8 @@ class BallReactor(paramak.Reactor):
         plasma_gap_vertical_thickness: the vertical thickness of the gap
             between the plasma and firstwall (cm). If left as None then the
             outer_plasma_gap_radial_thickness is used.
+        divertor_to_tf_gap_vertical_thickness: the vertical thickness of the
+            gap between the divertor and the TF coils.
         number_of_tf_coils: the number of tf coils
         pf_coil_radial_thicknesses: the radial
             thickness of each poloidal field coil.
@@ -76,6 +78,7 @@ def __init__(
             elongation: float,
             triangularity: float,
             plasma_gap_vertical_thickness: Optional[float] = None,
+            divertor_to_tf_gap_vertical_thickness: Optional[float] = 0,
             number_of_tf_coils: Optional[int] = 12,
             rear_blanket_to_tf_gap: Optional[float] = None,
             pf_coil_radial_thicknesses: Optional[Union[float, List[float]]] = None,
@@ -123,6 +126,7 @@ def __init__(
         self.rotation_angle = rotation_angle
 
         self.plasma_gap_vertical_thickness = plasma_gap_vertical_thickness
+        self.divertor_to_tf_gap_vertical_thickness = divertor_to_tf_gap_vertical_thickness
         if self.plasma_gap_vertical_thickness is None:
             self.plasma_gap_vertical_thickness = \
                 self.outer_plasma_gap_radial_thickness
@@ -315,23 +319,22 @@ def _make_vertical_build(self):
         # in a similar manner to the radial build
 
         self._firstwall_start_height = (
-            self._plasma.high_point[1] + self.plasma_gap_vertical_thickness
-        )
+            self._plasma.high_point[1] +
+            self.plasma_gap_vertical_thickness)
         self._firstwall_end_height = self._firstwall_start_height + \
             self.firstwall_radial_thickness
 
         self._blanket_start_height = self._firstwall_end_height
-        self._blanket_end_height = \
-            self._blanket_start_height + self.blanket_radial_thickness
+        self._blanket_end_height = self._blanket_start_height + self.blanket_radial_thickness
 
         self._blanket_rear_wall_start_height = self._blanket_end_height
-        self._blanket_rear_wall_end_height = (
-            self._blanket_rear_wall_start_height +
-            self.blanket_rear_wall_radial_thickness)
+        self._blanket_rear_wall_end_height = self._blanket_rear_wall_start_height + \
+            self.blanket_rear_wall_radial_thickness
+
+        self._tf_coil_start_height = self._blanket_rear_wall_end_height + \
+            self.divertor_to_tf_gap_vertical_thickness
 
-        self._tf_coil_height = self._blanket_rear_wall_end_height
-        self._center_column_shield_height = \
-            self._blanket_rear_wall_end_height * 2
+        self._center_column_shield_height = self._blanket_rear_wall_end_height * 2
 
         if self.rear_blanket_to_tf_gap is not None:
             self._tf_coil_start_radius = self._blanket_rear_wall_end_radius + \
@@ -343,7 +346,7 @@ def _make_vertical_build(self):
     def _make_inboard_tf_coils(self):
 
         self._inboard_tf_coils = paramak.CenterColumnShieldCylinder(
-            height=self._tf_coil_height * 2,
+            height=self._tf_coil_start_height * 2,
             inner_radius=self._inboard_tf_coils_start_radius,
             outer_radius=self._inboard_tf_coils_end_radius,
             rotation_angle=self.rotation_angle,
@@ -534,7 +537,7 @@ def _make_tf_coils(self):
                 with_inner_leg=False,
                 horizontal_start_point=(
                     self._inboard_tf_coils_start_radius,
-                    self._tf_coil_height),
+                    self._tf_coil_start_height),
                 vertical_mid_point=(
                     self._tf_coil_start_radius, 0),
                 thickness=self.outboard_tf_coil_radial_thickness,