diff options
| -rwxr-xr-x | pin.py | 51 |
1 files changed, 32 insertions, 19 deletions
@@ -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 } |