Add blanket arc radius calculations and update plotting to visualize arc points

Author: chris-ashe

process/blanket_library.py

@@ -107,6 +107,9 @@ def component_volumes(self):
                 build_variables.a_blkt_inboard_surface_full_coverage,
                 build_variables.a_blkt_outboard_surface_full_coverage,
                 build_variables.a_blkt_total_surface_full_coverage,
+                fwbs_variables.r_blkt_arc_1,
+                fwbs_variables.r_blkt_arc_2,
+                fwbs_variables.r_blkt_arc_3,
             ) = self.calculate_elliptical_blkt_areas(
                 rmajor=physics_variables.rmajor,
                 rminor=physics_variables.rminor,
@@ -394,7 +397,7 @@ def calculate_elliptical_blkt_areas(
         dr_shld_outboard: float,
         dr_blkt_outboard: float,
         dz_blkt_half: float,
-    ) -> tuple[float, float, float]:
+    ) -> tuple[float, float, float, float, float, float]:
         """Calculate elliptical blanket surface areas.
 
         :param rmajor: Major radius of the plasma (m)
@@ -418,7 +421,8 @@ def calculate_elliptical_blkt_areas(
         :param dz_blkt_half: Half-height of the blanket (m)
         :type dz_blkt_half: float
 
-        :return: Tuple containing inboard blanket surface area (m²), outboard blanket surface area (m²), and total blanket surface area (m²)
+        :return: Tuple containing inboard blanket surface area (m²), outboard blanket surface area (m²), and total blanket surface area (m²),
+        and three blanket arc radii (m)
 
         """
         # Major radius to centre of inboard and outboard ellipses (m)
@@ -445,7 +449,14 @@ def calculate_elliptical_blkt_areas(
             a_blkt_total_surface,
         ) = eshellarea(r1, r2, r3, dz_blkt_half)
 
-        return a_blkt_inboard_surface, a_blkt_outboard_surface, a_blkt_total_surface
+        return (
+            a_blkt_inboard_surface,
+            a_blkt_outboard_surface,
+            a_blkt_total_surface,
+            r1,
+            r2,
+            r3,
+        )
 
     @staticmethod
     def calculate_elliptical_blkt_volumes(
@@ -707,6 +718,27 @@ def output_blkt_volumes_and_areas(self):
             build_variables.a_blkt_total_surface_full_coverage,
         )
 
+        po.oblnkl(self.outfile)
+
+        po.ovarst(
+            self.outfile,
+            "First radial blanket arc point (m)",
+            "(r_blkt_arc_1)",
+            fwbs_variables.r_blkt_arc_1,
+        )
+        po.ovarst(
+            self.outfile,
+            "Second radial blanket arc point (m)",
+            "(r_blkt_arc_2)",
+            fwbs_variables.r_blkt_arc_2,
+        )
+        po.ovarst(
+            self.outfile,
+            "Third radial blanket arc point (m)",
+            "(r_blkt_arc_3)",
+            fwbs_variables.r_blkt_arc_3,
+        )
+
     def primary_coolant_properties(self, output: bool):
         """Calculates the fluid properties of the Primary Coolant in the FW and BZ.
         Uses middle value of input and output temperatures of coolant.

process/io/plot_proc.py

@@ -12817,6 +12817,21 @@ def plot_blkt_structure(
         label="Midplane",
     )
 
+    # Plot the points as black dots, number them, and connect them with lines
+    xs = [m_file.get(f"r_blkt_arc_{i}", scan=scan) for i in range(1, 4)]
+    for i, (x,) in enumerate(zip(xs, strict=False), 1):
+        ax.plot(x, 0.0, "ko", markersize=8)
+        ax.text(
+            x,
+            0.0,
+            str(i),
+            color="red",
+            fontsize=5,
+            ha="center",
+            va="center",
+            fontweight="bold",
+        )
+
     textstr_blkt_areas = (
         f"$\\mathbf{{Blanket \\ Areas:}}$\n\n"
         f"Inboard blanket, with holes and gaps: {m_file.get('a_blkt_inboard_surface', scan=scan):,.3f} $\\text{{m}}^2$ \n"