Added UCO 15 w/o fuel, disabled incorrect TSL on RPT homogenized material

1 files changed, 21 insertions, 9 deletions

diff --git a/pin.py b/pin.py
index 8ebc961..766b906 100755
--- a/pin.py
+++ b/pin.py
@@ -63,8 +63,8 @@ mat_triso_uco_5wo.add_nuclide("C12",0.13793)
 mat_triso_uco_5wo.add_nuclide("O16",0.51724)
 mat_triso_uco_5wo.set_density("g/cc",10.744)
 mat_triso_uco_5wo.temperature = mat_temp_fuel
-mat_triso_uco_5wo.add_s_alpha_beta("c_U_in_UC_10p")
-mat_triso_uco_5wo.add_s_alpha_beta("c_C_in_UC_10p")
+mat_triso_uco_5wo.add_s_alpha_beta("c_U_in_UC_5p")
+mat_triso_uco_5wo.add_s_alpha_beta("c_C_in_UC_5p")
 mat_list.append(mat_triso_uco_5wo)
 
 mat_triso_uco_10wo = openmc.Material(202,name="TRISO Pin 10% UCO Kernel")
@@ -78,7 +78,18 @@ mat_triso_uco_10wo.add_s_alpha_beta("c_U_in_UC_10p")
 mat_triso_uco_10wo.add_s_alpha_beta("c_C_in_UC_10p")
 mat_list.append(mat_triso_uco_10wo)
 
-mat_triso_uco_20wo = openmc.Material(203,name="TRISO Pin 20% UCO Kernel")
+mat_triso_uco_15wo = openmc.Material(203,name="TRISO Pin 15.5% UCO Kernel")
+mat_triso_uco_15wo.add_nuclide("U235",0.051732)
+mat_triso_uco_15wo.add_nuclide("U238",0.293200)
+mat_triso_uco_15wo.add_nuclide("C12",0.13793)
+mat_triso_uco_15wo.add_nuclide("O16",0.51724)
+mat_triso_uco_15wo.set_density("g/cc",10.5)
+mat_triso_uco_15wo.temperature = mat_temp_fuel
+mat_triso_uco_15wo.add_s_alpha_beta("c_U_in_UC_10p")
+mat_triso_uco_15wo.add_s_alpha_beta("c_C_in_UC_10p")
+mat_list.append(mat_triso_uco_15wo)
+
+mat_triso_uco_20wo = openmc.Material(204,name="TRISO Pin 20% UCO Kernel")
 mat_triso_uco_20wo.add_nuclide("U235",0.068794)
 mat_triso_uco_20wo.add_nuclide("U238",0.27604)
 mat_triso_uco_20wo.add_nuclide("C12",0.13793)
@@ -89,7 +100,7 @@ mat_triso_uco_20wo.add_s_alpha_beta("c_U_in_UC_HALEU")
 mat_triso_uco_20wo.add_s_alpha_beta("c_C_in_UC_HALEU")
 mat_list.append(mat_triso_uco_20wo)
 
-mat_triso_un = openmc.Material(204,name="TRISO Pin UN Kernel")
+mat_triso_un = openmc.Material(205,name="TRISO Pin UN Kernel")
 mat_triso_un.add_nuclide("U235",2.1125e-4)
 mat_triso_un.add_nuclide("U238",2.8758e-2)
 mat_triso_un.add_nuclide("U234",1.6795e-6)
@@ -101,21 +112,21 @@ mat_triso_un.add_s_alpha_beta("c_U_in_UN")
 mat_triso_un.add_s_alpha_beta("c_N_in_UN")
 mat_list.append(mat_triso_un)
 
-mat_triso_buffer = openmc.Material(205,name="TRISO Pin Buffer")
+mat_triso_buffer = openmc.Material(206,name="TRISO Pin Buffer")
 mat_triso_buffer.add_nuclide("C12",1.0)
 mat_triso_buffer.set_density("g/cc",1.04)
 mat_triso_buffer.temperature = mat_temp_fuel
 mat_triso_buffer.add_s_alpha_beta("c_Graphite")
 mat_list.append(mat_triso_buffer)
 
-mat_triso_pyc1 = openmc.Material(206,name="TRISO Pin PyC1")
+mat_triso_pyc1 = openmc.Material(207,name="TRISO Pin PyC1")
 mat_triso_pyc1.add_nuclide("C12",1.0)
 mat_triso_pyc1.set_density("g/cc",1.882)
 mat_triso_pyc1.temperature = mat_temp_fuel
 mat_triso_pyc1.add_s_alpha_beta("c_Graphite")
 mat_list.append(mat_triso_pyc1)
 
-mat_triso_sic = openmc.Material(207,name="TRISO Pin SiC")
+mat_triso_sic = openmc.Material(208,name="TRISO Pin SiC")
 mat_triso_sic.add_nuclide("Si28",0.4611)
 mat_triso_sic.add_nuclide("Si29",0.0234)
 mat_triso_sic.add_nuclide("Si30",0.0154)
@@ -126,7 +137,7 @@ mat_triso_sic.add_s_alpha_beta("c_C_in_SiC")
 mat_triso_sic.add_s_alpha_beta("c_Si_in_SiC")
 mat_list.append(mat_triso_sic)
 
-mat_triso_pyc2 = openmc.Material(208,name="TRISO Pin PyC2")
+mat_triso_pyc2 = openmc.Material(209,name="TRISO Pin PyC2")
 mat_triso_pyc2.add_nuclide("C12",1.0)
 mat_triso_pyc2.add_s_alpha_beta("c_Graphite")
 mat_triso_pyc2.temperature = mat_temp_fuel
@@ -138,6 +149,7 @@ mat_fuel_map = {
         "uo2":mat_triso_uo2,
         "uco5":mat_triso_uco_5wo,
         "uco10":mat_triso_uco_10wo,
+        "uco15":mat_triso_uco_15wo,
         "uco20":mat_triso_uco_20wo,
         "un":mat_triso_un
         }
@@ -252,7 +264,7 @@ if cli_args.technique=="rpt": # There is a minimum viable radius for RPT, so onl
             percent_type="vo",
             name="RPT Homogenized TRISO Compact"
             )
-    mat_triso_rpt.add_s_alpha_beta("c_Graphite")
+    #mat_triso_rpt.add_s_alpha_beta("c_Graphite")
 else:
     mat_triso_rpt = openmc.Material(name="PLACEHOLDER//DO NOT USE")
     mat_triso_rpt.add_nuclide("U235",1.0)