From e4cd463d4274a642e8cfe739b39e6ab5721a9066 Mon Sep 17 00:00:00 2001 From: Connor Moore Date: Mon, 29 Jun 2026 14:51:42 +0200 Subject: Added maximum RRPT radius calculation, driving script for parametric studies --- hpmr.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) (limited to 'hpmr.py') diff --git a/hpmr.py b/hpmr.py index 4075538..20ea672 100755 --- a/hpmr.py +++ b/hpmr.py @@ -6,6 +6,9 @@ import hashlib import sys import time +# For LCLRS +openmc.config["cross_sections"] = "/libs/endfb81_official/cross_sections.xml" + ############ COMMAND LINE ARGUMENTS ############ cli_parser = argparse.ArgumentParser( prog="HP-RM TRISO Homogenization Study", @@ -477,6 +480,7 @@ 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(f" -> Maximum RRPT ring radius: {(1.15**2 - 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 @@ -660,4 +664,4 @@ geometry.export_to_xml() settings.export_to_xml() # Run! -openmc.run() +#openmc.run() -- cgit v1.2.3