Added TSL for materials

1 files changed, 9 insertions, 0 deletions

diff --git a/pin.py b/pin.py
index 597a4e7..0382712 100755
--- a/pin.py
+++ b/pin.py
@@ -43,6 +43,7 @@ mat_graphite_matrix.add_nuclide("C12",0.9999997)
 mat_graphite_matrix.add_nuclide("B10",3e-7)
 mat_graphite_matrix.set_density("g/cc",1.806)
 mat_graphite_matrix.temperature = mat_temp_global
+mat_graphite_matrix.add_s_alpha_beta("c_Graphite")
 mat_list.append(mat_graphite_matrix)
 
 ## 200s - TRISO Pins ##
@@ -50,6 +51,8 @@ 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.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")
@@ -59,6 +62,8 @@ mat_triso_uco_10wo.add_nuclide("C12",0.13793)
 mat_triso_uco_10wo.add_nuclide("O16",0.51724)
 mat_triso_uco_10wo.set_density("g/cc",10.744)
 mat_triso_uco_10wo.temperature = mat_temp_fuel
+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")
@@ -68,6 +73,8 @@ 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")
@@ -78,6 +85,8 @@ mat_triso_un.add_nuclide("N14",2.8971e-4)
 mat_triso_un.add_nuclide("N15",2.8681e-2)
 mat_triso_un.set_density("g/cc",14.32)
 mat_triso_un.temperature = mat_temp_fuel
+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")