diff options
| author | Connor Moore <connor.moore@psi.ch> | 2026-06-29 13:44:59 +0200 |
|---|---|---|
| committer | Connor Moore <connor.moore@psi.ch> | 2026-06-29 13:44:59 +0200 |
| commit | e5b2ecd4df2db10f7835388554975987ffccbf69 (patch) | |
| tree | e410808e5a37084cfe07bae96a0816ee73fae6ea /hpmr.py | |
| parent | e07688afc0e63eb1aaa4379e2d8de3b8c9a55977 (diff) | |
Added TRISO stacking printout for packing
Diffstat (limited to 'hpmr.py')
| -rwxr-xr-x | hpmr.py | 15 |
1 files changed, 10 insertions, 5 deletions
@@ -12,8 +12,8 @@ cli_parser = argparse.ArgumentParser( description="Program to investigate various homogenization techniques for TRISO compacts in the HP-MR benchmark." ) -cli_parser.add_argument("-t","--technique",default="vwh",choices=["none","vwh","rpt","rrpt"],help="What homogenization technique to use (if any) [none/vwh/rpt/rrpt], default VWH") -cli_parser.add_argument("-r","--radius",default="0.8",type=float,help="Radius for RPT or RRPT homogenization [cm], default 0.5 cm. In the case of RRPT, this is the inner radius.") +cli_parser.add_argument("-t","--technique",default="none",choices=["none","vwh","rpt","rrpt"],help="What homogenization technique to use (if any) [none/vwh/rpt/rrpt], default none") +cli_parser.add_argument("-r","--radius",default="0.8",type=float,help="Radius for RPT or RRPT homogenization [cm], default 0.8 cm. In the case of RRPT, this is the inner radius.") cli_args = cli_parser.parse_args() ############ INTPUT-OUTPUT MAPPING ############ @@ -473,6 +473,11 @@ geo_vwh_fracs = [] geo_rpt_fracs = [] geo_triso_rpt_region = -openmc.ZCylinder(r=cli_args.radius,x0=0.0,y0=0.0)&+geo_constant_z_min&-geo_constant_z_max geo_n_triso = len(geo_triso_compact_packing) +print("\n TRISO Stacking Information from Packing Calculation") +print(f" -> Total TRISO pebbles packed: {geo_n_triso}") +print(f" -> Total TRISO volume in compact: {geo_n_triso*4/3*math.pi*geo_triso_spheres[4].r**3:.4f}") +print(f" -> Minimum RPT cylinder radius: {(geo_n_triso*4/3*geo_triso_spheres[4].r**3/(2*geo_constant_height))**0.5:.4f}") +print(" TRISO Stacking calculation complete!\n") # Get volume fractions for mixing geo_triso_compact_vwh_v_total = math.pi*geo_constant_triso_radius**2*2*geo_constant_height @@ -626,9 +631,9 @@ geometry = openmc.Geometry(geo_assembly_universe) ############ SETTINGS ############ settings = openmc.Settings() -settings.particles = 10000 -settings.batches = 400 -settings.inactive = 80 +settings.particles = 100000 +settings.batches = 200 +settings.inactive = 50 # Source sampling set_source_pts = [] |