SubChannel Requirements Traceability Matrix

[blockedtest]
  type = Exodiff
  input = 'blocked_test.i'
  exodiff = 'blocked_test_out.e'
  recover = false
  design = 'SCMBlockedMassFlowRateAux.md'
  issues = '#29237'
  requirement = 'The system shall apply inlet boundary massflow conditions for a blocked/partially blocked inlet'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  design = 'SCMMassFlowRateAux.md'
  issues = '#29511'
  requirement = 'The system shall apply inlet boundary massflow conditions'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  requirement = 'The system shall calculate linear heat rate on fuel pin surface'
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  recover = False
  requirement = 'The system shall calculate a uniform massflow rate inlet boundary'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = False
  requirement = 'The system will calculate the geometric parameters for FCTF'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = False
  requirement = 'The system shall calculate the flow area and wetted perimeter for a MARVEL-type microreactor core modeled in SCM'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = False
  requirement = 'The system will apply initial massflow rate'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = False
  requirement = 'The system will calculate the subchannel flow area for quadrilateral assemblies'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = False
  requirement = 'The system will calculate the linear heat rate on subchannels/fuel-pins in a quadrilateral lattice'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = False
  requirement = 'The system will calculate the wetted perimeter of subchannels in a quadrilateral lattice'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = false
  requirement = 'The system will calculate the subchannel flow area for triangular assemblies'
[]

[test_pin_power]
  type = Exodiff
  input = 'test_pin_power.i'
  exodiff = 'test_pin_power_out.e'
  recover = false
  requirement = 'The system will calculate the linear heat rate on subchannels/fuel-pins in a triangular lattice'
[]

[test_axial_shape]
  type = Exodiff
  input = 'test_axial_shape.i'
  exodiff = 'test_axial_shape_out.e'
  recover = false
  requirement = 'The system will calculate the linear heat rate on subchannels/fuel-pins in a triangular lattice, but with an axial power profile'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  recover = false
  requirement = 'The system will calculate the wetted perimeter of subchannels in a triangular lattice'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a detailed mesh of the inter-wrapper for quadrilateral assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a detailed mesh of the fuel pins for quadrilateral assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_3x3.e
  cli_args = "--mesh-only coords_3x3.e"
  recover = false
  requirement = 'The system will create a detailed mesh of the subchannels for quadrilateral assemblies (3x3)'
[]

[1x2]
  type = Exodiff
  input = coords.i
  exodiff = coords_1x2.e
  cli_args = "--mesh-only coords_1x2.e Mesh/subchannel/nx=1 Mesh/subchannel/ny=2"
  recover = false
  requirement = 'The system will create a detailed mesh of the subchannels for quadrilateral assemblies (1x2)'
[]

[1x3]
  type = Exodiff
  input = coords.i
  exodiff = coords_1x3.e
  cli_args = "--mesh-only coords_1x3.e Mesh/subchannel/nx=1 Mesh/subchannel/ny=3"
  recover = false
  requirement = 'The system will create a detailed mesh of the subchannels for quadrilateral assemblies (1x3)'
[]

[3x1]
  type = Exodiff
  input = coords.i
  exodiff = coords_3x1.e
  cli_args = "--mesh-only coords_3x1.e Mesh/subchannel/nx=3 Mesh/subchannel/ny=1"
  recover = false
  requirement = 'The system will create a detailed mesh of the subchannels for quadrilateral assemblies (3x1)'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a detailed mesh of the fuel pins for triangular assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a detailed mesh of the subchannels in triangular assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a mesh of the inter-wrapper for quadrilateral assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  issues = '#29241'
  requirement = 'The system will create a mesh of the subchannels for quadrilateral assemblies'
[]

[subchannel_number]
  type = 'RunException'
  input = 'subchannel_number.i'
  expect_err = "The number of subchannels cannot be less than 2 in both directions \(x and y\)\. Smallest assembly allowed is either 2X1 or 1X2\."
  issues = '#29224'
  requirement = 'The system shall check the minimum number of subchannels'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a mesh of the fuel pins for quadrilateral assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a mesh of the inter-wrapper for triangualar assemblies'
[]

[coords]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  mesh_mode = 'replicated'
  requirement = 'The system will create a hexagonal duct mesh and a pin mesh for triangular assemblies'
[]

[test]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  mesh_mode = 'replicated'
  requirement = 'The system will create a hexagonal duct mesh for triangular assemblies'
[]

[tricoords3]
  type = Exodiff
  input = tricoords3.i
  exodiff = tricoords3_in.e
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a subchannel mesh for triangular assemblies'
[]

[check]
  type = CSVDiff
  input = 'main_app.i'
  csvdiff = 'main_app_out_sc0_sub_0001.csv'
  requirement = "The system shall be able to solve a subchannel problem as a child application."
  recover = false
  max_parallel = 1
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'Temp_out_matrix.csv'
  recover = False
  max_parallel = 1
  requirement = 'The system will print a variable on a plane in a quadrilateral assembly'
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  recover = False
  requirement = 'The system will print a variable on a specific subchannel and height in quadrilateral assemblies'
[]

[subchannel_height]
  type = 'RunException'
  input = 'test_error.i'
  expect_err = ": No element located at given height."
  requirement = 'The system will test the user provided height'
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  recover = False
  requirement = 'The system will print a variable on a specific subchannel and height in triangular assemblies'
[]

[test_LBE]
  type = CSVDiff
  input = test_LBE-19pin.i
  csvdiff = test_LBE-19pin_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  max_parallel = 1
  requirement = 'The system will examine the tri subchannel solver for lead bismuth eutetic coolant'
[]

[test_LEAD]
  type = CSVDiff
  input = test_LEAD-19pin.i
  csvdiff = test_LEAD-19pin_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  max_parallel = 1
  requirement = 'The system will examine the tri subchannel solver for lead coolant'
[]

[SHRT-17_SS_test]
  type = CSVDiff
  input = XX09_SS_SHRT17.i
  csvdiff = XX09_SS_SHRT17_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the subchannel solver for the EBR-II case'
[]

[test_explicit]
  type = CSVDiff
  input = test19_explicit.i
  csvdiff = test19_explicit_out.csv
  valgrind = NONE
  recover = false
  allow_warnings = true
  max_parallel = 1
  requirement = 'The system will examine the subchannel explicit solver for a 19pin liquid sodium cooled assembly'
[]

[test_implicit]
  type = CSVDiff
  input = test19_implicit.i
  csvdiff = test19_implicit_out.csv
  valgrind = NONE
  recover = false
  allow_warnings = true
  max_parallel = 1
  requirement = 'The system will examine the subchannel implicit solver for a 19pin liquid sodium cooled assembly'
[]

[test_monolithic]
  type = CSVDiff
  input = test19_monolithic.i
  csvdiff = test19_monolithic_out.csv
  abs_zero = 1e-7
  valgrind = NONE
  recover = false
  max_parallel = 1
  requirement = 'The system will examine the subchannel monolithic solver for a 19pin liquid sodium cooled assembly'
[]

[test_full_monolithic]
  type = CSVDiff
  input = test19_full_monolithic.i
  csvdiff = test19_full_monolithic_out.csv
  rel_err = 5e-4
  abs_zero = 1e-7
  valgrind = NONE
  recover = false
  max_parallel = 1
  requirement = 'The system will examine the subchannel full monolithic solver for a 19pin liquid sodium cooled assembly'
[]

[coupling_test]
  type = CSVDiff
  input = main.i
  csvdiff = "main_out_viz0_sample_center_pin_0001.csv
              main_out_sub0_sample_center_line_0001.csv"
  recover = false
  mesh_mode = 'replicated'
  allow_warnings = true
  max_parallel = 1
  requirement = 'The system will examine the coupling using SCM as main app'
[]

[quad_interwrapper]
  type = Exodiff
  input = quad_interwrapper.i
  exodiff = quad_interwrapper_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the explicit interwrapper solver for quadrilateral assemblies'
[]

[quad_interwrapper_monolithic]
  type = Exodiff
  input = quad_interwrapper_monolithic.i
  exodiff = quad_interwrapper_monolithic_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the monolithic interwrapper solver for quadrilateral assemblies'
[]

[quad_interwrapper_monolithic_staggered]
  type = Exodiff
  input = quad_interwrapper_monolithic_staggered.i
  exodiff = quad_interwrapper_monolithic_staggered_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the monolitich staggered grid interwrapper solver for quadrilateral assemblies'
[]

[tri_interwrapper]
  type = Exodiff
  input = tri_interwrapper.i
  exodiff = tri_interwrapper_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the explicit interwrapper solver for triangular assemblies'
[]

[tri_interwrapper_tight_side]
  type = Exodiff
  input = tri_interwrapper.i
  exodiff = tri_interwrapper_tight_side.e
  cli_args = "TriInterWrapperMesh/sub_channel/tight_side_bypass=true Outputs/file_base=tri_interwrapper_tight_side"
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the explicit interwrapper solver for triangular assemblies and side bypass'
[]

[tri_interwrapper_implicit]
  type = Exodiff
  input = tri_interwrapper_implicit.i
  exodiff = tri_interwrapper_implicit_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the implicit interwrapper solver for triangular assemblies'
[]

[tri_interwrapper_monolithic]
  type = Exodiff
  input = tri_interwrapper_monolithic.i
  exodiff = tri_interwrapper_monolithic_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the monolithic interwrapper solver for triangular assemblies'
[]

[tri_interwrapper_monolithic_staggered]
  type = Exodiff
  input = tri_interwrapper_monolithic_staggered.i
  exodiff = tri_interwrapper_monolithic_staggered_out.e
  custom_cmp = 'test.exodiff'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system will examine the monolithic staggered grif interwrapper solver for triangular assemblies'
[]

[psbt_regression_test_explicit]
  type = CSVDiff
  input = psbt_explicit.i
  csvdiff = psbt_explicit_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using an explicit algorithm.'
[]

[psbt_regression_test_explicit_friction2]
  type = CSVDiff
  input = psbt_explicit_v2.i
  csvdiff = psbt_explicit_v2_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using the non-default friction/mixing model.'
[]

[psbt_regression_test_explicit_staggered]
  type = CSVDiff
  input = psbt_explicit_staggered.i
  csvdiff = psbt_explicit_staggered_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using an explicit algorithm and a staggered pressure formulation.'
[]

[psbt_regression_test_implicit]
  type = CSVDiff
  input = psbt_implicit.i
  csvdiff = psbt_implicit_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using an implicit algorithm.'
[]

[psbt_regression_test_monolithic]
  type = CSVDiff
  input = psbt_monolithic.i
  csvdiff = psbt_monolithic_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using a monolithic algorithm.'
[]

[psbt_regression_test_full_monolithic]
  type = CSVDiff
  input = psbt_full_monolithic.i
  csvdiff = psbt_full_monolithic_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using a full-monolithic algorithm.'
[]

[psbt_regression_test_full_monolithic_staggered]
  type = CSVDiff
  input = psbt_full_monolithic_staggered.i
  csvdiff = psbt_full_monolithic_staggered_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  rel_err = 6e-5
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using a full-monolithic algorithm and a staggered pressure formulation.'
[]

[steady]
  type = Exodiff
  input = steady.i
  exodiff = steady_out.e
  recover = false
  abs_zero = 1e-7
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to restart after a steady state solution'
[]

[transient]
  type = Exodiff
  input = transient.i
  exodiff = transient_out.e
  prereq = 'steady'
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  requirement = 'The system shall be able to restart after a transient solution'
[]

[quad]
  type = Exodiff
  input = 'quad.i'
  exodiff = 'quad_out.e quad_out_viz0.e'
  recover = false
  max_parallel = 1
  requirement = 'The system shall be able to transger the subchannel/pin solution onto a detailed mesh'
[]

[quad_iw]
  type = Exodiff
  input = quad_iw.i
  exodiff = 'quad_iw_out.e quad_iw_out_viz0.e'
  custom_cmp = 'test.exodiff'
  recover = false
  abs_zero = 5e-6
  max_threads = 1
  max_parallel = 1
  cli_args = '-ksp_rtol 1e-3'
  requirement = 'The system shall be able to transger the interwrapper solution onto a detailed mesh'
[]