Added maximum RRPT radius calculation, driving script for parametric studies

1 files changed, 5 insertions, 1 deletions

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()