add tally to c5g7, set up smr-mg

Author: ilhamv

05-C5G7-original/part1/input.py

@@ -1,5 +1,5 @@
 import numpy as np
-import h5py, mcdc
+import h5py, mcdc, math
 
 # =============================================================================
 # Materials
@@ -284,8 +284,24 @@ def set_mat(mat):
 # Set tally, setting, and run mcdc
 # =============================================================================
 
+# Tally
+mcdc.tally.mesh_tally(
+    scores=["flux"],
+    x=np.linspace(0.0, pitch * 17 * 3, 17 * 3 + 1),
+    y=np.linspace(-pitch * 17 * 3, 0.0, 17 * 3 + 1),
+    z=np.linspace(-(core_height / 2 + reflector_thickness), (core_height / 2 + reflector_thickness), int(math.ceil((core_height + 2.0 * reflector_thickness) / pitch)) + 1),
+    g=np.array([-0.5, 4.5, 6.5])
+)
+mcdc.tally.mesh_tally(
+    scores=["fission"],
+    x=np.linspace(0.0, pitch * 17 * 2, 17 * 2 + 1),
+    y=np.linspace(-pitch * 17 * 2, 0.0, 17 * 2 + 1),
+    z=np.linspace(-(core_height / 2), (core_height / 2), int(math.ceil(core_height / pitch)) + 1),
+)
+
 # Setting
-mcdc.setting(N_particle=1e5)
+#mcdc.setting(N_particle=1e5)
+mcdc.setting(N_particle=1e3)
 mcdc.eigenmode(N_inactive=50, N_active=150, gyration_radius="all", save_particle=True)
 mcdc.population_control()
 

08-SMR_MG-4_phase/kinetic/input.py

@@ -9,32 +9,51 @@
 lib = h5py.File("../../data/MGXS-C5G7.h5", "r")
 
 # Setter
-def set_mat(mat):
+def set_mat(mat, frac):
+    capture = mat[0]['capture'][:] * frac[0]
+    scatter = mat[0]['scatter'][:] * frac[0]
+    fission = mat[0]['fission'][:] * frac[0]
+    nu_p = mat[0]['nu_p'][:] * frac[0]
+    nu_d = mat[0]['nu_d'][:] * frac[0]
+    chi_p = mat[0]['chi_p'][:] * frac[0]
+    chi_d = mat[0]['chi_d'][:] * frac[0]
+    speed = mat[0]['speed'][:] * frac[0]
+    decay = mat[0]['decay'][:] * frac[0]
+    for i in range(1, len(mat)):
+        capture += mat[i]['capture'][:] * frac[i]
+        scatter += mat[i]['scatter'][:] * frac[i]
+        fission += mat[i]['fission'][:] * frac[i]
+        nu_p += mat[i]['nu_p'][:] * frac[i]
+        nu_d += mat[i]['nu_d'][:] * frac[i]
+        chi_p += mat[i]['chi_p'][:] * frac[i]
+        chi_d += mat[i]['chi_d'][:] * frac[i]
+        speed += mat[i]['speed'][:] * frac[i]
+        decay += mat[i]['decay'][:] * frac[i]
     return mcdc.material(
-        capture=mat["capture"][:],
-        scatter=mat["scatter"][:],
-        fission=mat["fission"][:],
-        nu_p=mat["nu_p"][:],
-        nu_d=mat["nu_d"][:],
-        chi_p=mat["chi_p"][:],
-        chi_d=mat["chi_d"][:],
-        speed=mat["speed"],
-        decay=mat["decay"],
+        capture=capture,
+        scatter=scatter,
+        fission=fission,
+        nu_p=nu_p,
+        nu_d=nu_d,
+        chi_p=chi_p,
+        chi_d=chi_d,
+        speed=speed,
+        decay=decay,
     )
 
 # Materials
-m_helium = set_mat(lib["gt"])
-m_inconel = set_mat(lib["gt"])
-m_ss302 = set_mat(lib["gt"])
-m_ss304 = set_mat(lib["gt"])
-m_csteel = set_mat(lib["gt"])
-m_zr4 = set_mat(lib["gt"])
-m_m5 = set_mat(lib["gt"])
-m_water = set_mat(lib["mod"])
-m_fuel31 = set_mat(lib["mox87"])
-m_fuel24 = set_mat(lib["mox7"])
-m_fuel16 = set_mat(lib["mox43"])
-m_cr = set_mat(lib["cr"])
+m_helium = set_mat([lib["gt"]], [1.0])
+m_inconel = set_mat([lib["gt"]], [1.0])
+m_ss302 = set_mat([lib["gt"]], [1.0])
+m_ss304 = set_mat([lib["gt"]], [1.0])
+m_csteel = set_mat([lib["gt"]], [1.0])
+m_zr4 = set_mat([lib["gt"]], [1.0])
+m_m5 = set_mat([lib["gt"]], [1.0])
+m_water = set_mat([lib["mod"], lib['cr']], [0.97, 0.03])
+m_fuel31 = set_mat([lib["uo2"]], [1.0])
+m_fuel24 = set_mat([lib["mox87"]], [1.0])
+m_fuel16 = set_mat([lib["mox43"]], [1.0])
+m_cr = set_mat([lib["cr"]], [2.0])
 
 # ======================================================================================
 # Geometry
@@ -152,11 +171,11 @@ def set_mat(mat):
 
 # Problem spec.
 config['time'] = 0.0
-config['frac_r'] = 0.0
-config['frac_rx'] = 0.0
-config['frac_s'] = 0.0
-config['frac_sx2'] = 0.0
-config['frac_sx3'] = 0.0
+config['frac_r'] = 1.0
+config['frac_rx'] = 1.0
+config['frac_s'] = 1.0
+config['frac_sx2'] = 0.65
+config['frac_sx3'] = 1.0
 cases = ['r', 'rx', 's', 'sx2', 'sx3']
 
 # All cases
@@ -893,7 +912,7 @@ def set_mat(mat):
 # =============================================================================
 
 # Setting
-mcdc.setting(N_particle=1e3)
+mcdc.setting(N_particle=1e5)
 mcdc.eigenmode(N_inactive=50, N_active=150, gyration_radius="all")
 mcdc.population_control()