Merging MSRE Xe Poisoning Model

doc/content/bib/vtb.bib

@@ -1551,6 +1551,17 @@ @techreport{steffy1970
             year = 1970
 }
 
+@techreport{osti_2562160,
+  author       = {Shahbazi, Shayan and Fei, Tingzhou and Fang, Jun and Jeong, Yeongshin and Tano, Mauricio and de Oliveira, Rodrigo and Walker, Samuel and Hartanto, Donny and Davidson, Eva},
+  title        = {Towards the Development of Molten Salt Reactor Multiphysics Modeling Capabilities: Depletion, Thermochemistry, Reactor Dynamics, and Species Tracking},
+  institution  = {Argonne National Laboratory},
+  number       = {ANL/NSE--24/62},
+  doi          = {10.2172/2562160},
+  url          = {https://www.osti.gov/biblio/2562160},
+  place        = {United States},
+  year         = {2024},
+  month        = {09}}
+
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % references for the MSFR Thermochimica multiphysics core model
 

doc/content/msr/msre/index.md

@@ -6,6 +6,8 @@
 
 [Multiphysics RZ Model](multiphysics_rz_model/index.md)
 
+[Modeling Xenon poisoning](msre/msre_xe.md)
+
 [Cardinal model of pipe thermal striping](cardinal_thermal_striping.md)
 
 [nekRS CFD model of lower plenum](msre/lp_nekrs_model.md)

msr/msre/xe_poisoning/.gitattributes

@@ -0,0 +1 @@
+neu_xe_init.e filter=lfs diff=lfs merge=lfs -text

msr/msre/xe_poisoning/gold/neu_xe_out.csv.gold

@@ -0,0 +1,4 @@
+time,Fuel_avg_Temp,Fuel_core_avg_Temp,Fuel_core_max_Temp,Fuel_max_Temp,Leakage,Mod_avg_Temp,Mod_max_Temp,absorption_RR,eigenvalue,fission_RR,ngamma_RR,nufission_RR,power_scaling,total_power
+0,923,923,923,923,0,923,923,0,1,0,0,0,3.456287690711e+15,5000000
+1,923,923,923,923,37.712728282201,923,923,70.041277574854,1.0123676840583,44.597005994381,24.868374260478,108.36777298497,3.6167661946162e+15,5000000
+2,923,923,923,923,37.712728282201,923,923,70.041277574854,1.0123676840583,44.597005994381,24.868374260478,108.36777298497,3.6167661946162e+15,5000000

msr/msre/xe_poisoning/gold/th_xe_out.csv.gold

@@ -0,0 +1,3 @@
+time,Xe_avg,area_pp_downcomer_inlet,inlet_T,inlet_p,outlet_T,outlet_p,return_line_fraction,vfr_downcomer,vfr_pump,vfr_return_line
+0,0,0.19677333595724,0,0,0,0,0,0,0,0
+0.1,2.6314560213783e-27,0.19677333595724,923.00465965952,-413648.28395474,909.53681689842,-480719.92049981,-0.18272592545552,-0.028604114493781,0.029091320728937,0.0053157385029184

msr/msre/xe_poisoning/mesh/.gitattributes

@@ -0,0 +1 @@
+strip_mesh_in.e filter=lfs diff=lfs merge=lfs -text

msr/msre/xe_poisoning/mesh/strip_mesh.i

@@ -0,0 +1,160 @@
+fuel_pipe_D = '${fparse 10*0.0254}'
+fuel_pipe_R = '${fparse fuel_pipe_D/2.0}'
+core_internal_R = '${fparse 0.496855718 + 2.5*0.0254 - fuel_pipe_R}'
+core_outer_gap = 0.098795213
+core_barrel_thickness = 0.005
+downcomer_thickness = 0.045589414
+height_retrun_elbow = 0.05
+
+external_pipe_length = '${fparse 6.71 * 2.0}'
+external_piping_volume = '${fparse pi*((2.5*0.0254)^2)*external_pipe_length}'
+height_pump = '${fparse fuel_pipe_R / 1.5}'
+outer_R = '${fparse fuel_pipe_R +core_internal_R + core_outer_gap + core_barrel_thickness + downcomer_thickness}'
+top_disc_volume = '${fparse pi*outer_R^2*height_pump}'
+remaining_volume = '${fparse external_piping_volume-top_disc_volume}'
+#core_riser_flow_area = '${fparse pi*fuel_pipe_R^2}'
+core_downcomer_flow_area = '${fparse pi*(outer_R^2 - (outer_R-downcomer_thickness)^2)}'
+piping_height = '${fparse remaining_volume/(core_downcomer_flow_area + core_downcomer_flow_area)}'
+
+[Mesh]
+  coord_type = 'RZ'
+  type = MeshGeneratorMesh
+  block_id = '1 2 3 4 6 7 8 9 10 11 14 15'
+  block_name = 'core lower_plenum  upper_plenum  down_comer  core_barrel riser pump elbow void bypass_line strip return_line'
+  uniform_refine = 1
+
+  [cartesian_mesh]
+    type = CartesianMeshGenerator
+    dim = 2
+    dx = '${height_retrun_elbow}
+          ${fparse fuel_pipe_R - height_retrun_elbow}
+          ${height_retrun_elbow}
+          ${fparse core_internal_R - height_retrun_elbow}
+          ${core_outer_gap}
+          ${core_barrel_thickness}
+          ${downcomer_thickness}'
+    ix = '2
+          4
+          2
+          8
+          2
+          1
+          1'
+    dy = '0.1715   0.100   0.100   0.246   0.246   0.246
+          0.246   0.246   0.100   0.100   0.1715 ${piping_height} ${height_pump}
+          ${height_retrun_elbow} ${height_retrun_elbow}'
+    iy = '6  4  4  10  10  10  10  10  4  4  6  4 4 4 4'
+    subdomain_id = ' 2  2  2  2  2  2  2
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     1  1  1  1  1  6  4
+                     3  3  3  3  3  6  4
+                     7  7  10 10 10 10 9
+                     8  8  9  9  9  9  9
+                     15 13 11 12 12 12 12
+                     15 14 11 12 12 12 12'
+  []
+  [core_downcomer_boundary]
+    input = cartesian_mesh
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '6'
+    paired_block = '1 4 2 3'
+    new_boundary = core_downcomer_boundary
+  []
+  [loop_boundary]
+    input = core_downcomer_boundary
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '9 7 3'
+    paired_block = '10'
+    new_boundary = 'loop_boundary'
+  []
+  [bypass_line_boundary]
+    input = loop_boundary
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '11'
+    paired_block = '12'
+    new_boundary = 'bypass_line_boundary'
+  []
+  [adding_top_to_return_line]
+    input = bypass_line_boundary
+    type = RenameBoundaryGenerator
+    old_boundary = 'top'
+    new_boundary = 'extraction'
+  []
+  [strip_line_boundary]
+    input = adding_top_to_return_line
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '14'
+    paired_block = '15 11'
+    new_boundary = 'strip_line_boundary'
+  []
+  [return_line_boundary]
+    input = strip_line_boundary
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '11 14 15 8'
+    paired_block = '13'
+    new_boundary = 'return_line_boundary'
+  []
+  [internal_pump_sideset]
+    input = return_line_boundary
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '7'
+    paired_block = '8'
+    new_boundary = 'pump_outlet'
+  []
+  [internal_downcomer_sideset]
+    input = internal_pump_sideset
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '4'
+    paired_block = '9'
+    new_boundary = 'downcomer_inlet'
+  []
+  [top_core_barrel]
+    input = internal_downcomer_sideset
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '6'
+    paired_block = '10'
+    new_boundary = 'top_core_barrel'
+  []
+  [top_boundary_return_line]
+    input = top_core_barrel
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '9'
+    paired_block = '12'
+    new_boundary = 'top'
+  []
+  [return_line_internal_sideset]
+    input = top_boundary_return_line
+    type = SideSetsBetweenSubdomainsGenerator
+    primary_block = '15'
+    paired_block = '8'
+    new_boundary = 'return_line_sideset'
+  []
+  # [missing_boundary]
+  #   input = return_line_internal_sideset
+  #   type = SideSetsBetweenSubdomainsGenerator
+  #   primary_block = '13'
+  #   paired_block = '8'
+  #   new_boundary = 'top_of_pump'
+  # []
+  [block_delete]
+    input = return_line_internal_sideset
+    type = BlockDeletionGenerator
+    block = '10 12 13'
+  []
+  #[block_delete]
+  #  type = BlockDeletionGenerator
+  #  input = top_core_barrel
+  #  block = '10'
+  #[]
+[]
+
+[Executioner]
+  type = Steady
+[]

msr/msre/xe_poisoning/mgxs/.gitattributes

@@ -0,0 +1 @@
+xs_msre_micro_xe.xml filter=lfs diff=lfs merge=lfs -text