diff options
Diffstat (limited to 'candu-37.py')
| -rw-r--r-- | candu-37.py | 5 |
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}") |