Added 5w/o UCO fuel and changed 1975 UCO to 20 UCO

1 files changed, 32 insertions, 19 deletions

diff --git a/pin.py b/pin.py
index 0382712..f0b666d 100755
--- a/pin.py
+++ b/pin.py
@@ -15,7 +15,7 @@ cli_parser = argparse.ArgumentParser(
 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.5",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("-p","--packing-factor",default=0.29,type=float,help="Packing factor for TRISO packing (0.0, 1.0), default 0.29")
-cli_parser.add_argument("-f","--fuel",default="uo2",choices=["uo2","uco10","uco20","un"],help="What type of fuel to use [uo2/uco10/uco20/un], default uo2.")
+cli_parser.add_argument("-f","--fuel",default="uo2",choices=["uo2","uco5","uco10","uco20","un"],help="What type of fuel to use [uo2/uco10/uco20/un], default uo2.")
 cli_parser.add_argument("-e","--execute",default=False,action="store_true",help="Whether or not to automatically run OpenMC, default false.")
 cli_args = cli_parser.parse_args()
 
@@ -50,12 +50,24 @@ mat_list.append(mat_graphite_matrix)
 mat_triso_uo2 = openmc.Material(200,name="TRISO Pin UO2 Kernel")
 mat_triso_uo2.add_element("U",1.0,enrichment=5.0)
 mat_triso_uo2.add_element("O",2.0)
+mat_triso_uo2.set_density("g/cc",10.97)
 mat_triso_uo2.temperature = mat_temp_fuel
 mat_triso_uo2.add_s_alpha_beta("c_U_in_UO2")
 mat_triso_uo2.add_s_alpha_beta("c_O_in_UO2")
 mat_list.append(mat_triso_uo2)
 
-mat_triso_uco_10wo = openmc.Material(201,name="TRISO Pin UCO Kernel")
+mat_triso_uco_5wo = openmc.Material(201,name="TRISO Pin 5% UCO Kernel")
+mat_triso_uco_5wo.add_nuclide("U235",0.0172417)
+mat_triso_uco_5wo.add_nuclide("U238",0.3275923)
+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_list.append(mat_triso_uco_5wo)
+
+mat_triso_uco_10wo = openmc.Material(202,name="TRISO Pin 10% UCO Kernel")
 mat_triso_uco_10wo.add_nuclide("U235",0.034483)
 mat_triso_uco_10wo.add_nuclide("U238",0.310400)
 mat_triso_uco_10wo.add_nuclide("C12",0.13793)
@@ -66,18 +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_1975wo = openmc.Material(202,name="TRISO Pin UCO Kernel")
-mat_triso_uco_1975wo.add_nuclide("U235",0.068794)
-mat_triso_uco_1975wo.add_nuclide("U238",0.27604)
-mat_triso_uco_1975wo.add_nuclide("C12",0.13793)
-mat_triso_uco_1975wo.add_nuclide("O16",0.51724)
-mat_triso_uco_1975wo.set_density("g/cc",10.744)
-mat_triso_uco_1975wo.temperature = mat_temp_fuel
-mat_triso_uco_1975wo.add_s_alpha_beta("c_U_in_UC_HALEU")
-mat_triso_uco_1975wo.add_s_alpha_beta("c_C_in_UC_HALEU")
-mat_list.append(mat_triso_uco_1975wo)
-
-mat_triso_un = openmc.Material(203,name="TRISO Pin UN Kernel")
+mat_triso_uco_20wo = openmc.Material(203,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)
+mat_triso_uco_20wo.add_nuclide("O16",0.51724)
+mat_triso_uco_20wo.set_density("g/cc",10.744)
+mat_triso_uco_20wo.temperature = mat_temp_fuel
+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.add_nuclide("U235",2.1125e-4)
 mat_triso_un.add_nuclide("U238",2.8758e-2)
 mat_triso_un.add_nuclide("U234",1.6795e-6)
@@ -89,21 +101,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(204,name="TRISO Pin Buffer")
+mat_triso_buffer = openmc.Material(205,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(205,name="TRISO Pin PyC1")
+mat_triso_pyc1 = openmc.Material(206,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(206,name="TRISO Pin SiC")
+mat_triso_sic = openmc.Material(207,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)
@@ -114,7 +126,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(207,name="TRISO Pin PyC2")
+mat_triso_pyc2 = openmc.Material(208,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
@@ -124,8 +136,9 @@ mat_list.append(mat_triso_pyc2)
 ## Dynamic Fuel Material ##
 mat_fuel_map = {
         "uo2":mat_triso_uo2,
+        "uco5":mat_triso_uco_5wo,
         "uco10":mat_triso_uco_10wo,
-        "uco20":mat_triso_uco_1975wo,
+        "uco20":mat_triso_uco_20wo,
         "un":mat_triso_un
         }