Added cladding circumference, 52-element (and alt 51-element) Apollonian bunddles

1 files changed, 5 insertions, 0 deletions

diff --git a/candu-37.py b/candu-37.py
index e9642e7..b58dd68 100644
--- a/candu-37.py
+++ b/candu-37.py
@@ -108,9 +108,12 @@ def make_ring(n: int, r: float, alpha: float) -> None:
         theta = (i*(360/n) + alpha)*np.pi/180
         x = r*np.cos(theta)
         y = r*np.sin(theta)
+        #print((float(x),float(y)))
 
         make_pin(rf=r_fuel, rc=r_clad, x0=x, y0=y)
 
+    return
+
 
 #> Start with the innermost pin
 make_pin(rf=r_fuel, rc=r_clad, x0=0.0, y0=0.0)
@@ -184,6 +187,7 @@ settings.export_to_xml()
 
 V_fuel = len(fuel_region_list)*np.pi*r_fuel**2
 V_clad = len(clad_region_list)*np.pi*(r_clad**2 - r_fuel**2)
+C_clad = len(clad_region_list)*2*np.pi*r_clad
 V_mod = np.pi*pt_inner.r**2 - (V_fuel + V_clad) + (28.575**2 - np.pi*ct_outer.r**2)
 
 d_fuel = mat_fuel.get_nuclide_atom_densities()
@@ -194,6 +198,7 @@ DTU_ratio = (V_mod * N_mod) / (V_fuel * N_fuel)
 
 #> Print these for the final table
 print(f"Flow area, cm² = {V_mod}")
+print(f"Cladding ciricumference, cm = {C_clad}")
 print(f"Fuel mass per length, g/cm = {V_fuel*mat_fuel.density}")
 print(f"Cladding mass per length, g/cm = {V_clad*mat_zircaloy_4.density}")
 print(f"DTU ratio = {DTU_ratio}")