SubChannel Requirements Traceability Matrix

[SCMQuadPower_cells_aligned_with_pins]
  type = CSVDiff
  input = 'test_with_pins.i'
  csvdiff = 'test_with_pins_out_line_check_0001.csv test_with_pins_out.csv'
  recover = False
  abs_zero = 1e-7
  requirement = 'The system will populate a sinusoidal linear heat rate on fuel pins in a quadrilateral lattice using an initial condition and auxiliary kernel'
[]

[SCMQuadPower_cells_not_aligned_with_pins]
  type = CSVDiff
  input = 'test_with_pins.i'
  csvdiff = 'not_aligned_with_pins_line_check_0001.csv not_aligned_with_pins.csv'
  cli_args = 'num_cells=11 Outputs/file_base=not_aligned_with_pins'
  recover = False
  abs_zero = 1e-7
  requirement = 'The system will populate a sinusoidal linear heat rate on fuel pins in a quadrilateral lattice using an initial condition and auxiliary kernel in a mesh that has axial divisions that are not aligned with the heated region entrance and exit heights.'
[]

[SCMQuadPower_pin_dependent_axial_shape]
  type = CSVDiff
  input = 'two_pin_shapes.i'
  csvdiff = 'two_pin_shapes_out.csv two_pin_shapes_out_q_prime_profiles_0001.csv'
  recover = False
  abs_zero = 1e-7
  issues = '#34082'
  requirement = 'The system will assign a different axial linear heat rate shape to each fuel pin in a quadrilateral lattice from a single function of the pin position, while conserving the power assigned to each pin by the radial power distribution, using an initial condition and auxiliary kernel.'
[]

[SCMQuadPower_postprocessor_power]
  type = CSVDiff
  input = 'postprocessor_power.i'
  csvdiff = 'postprocessor_power_out.csv'
  recover = False
  abs_zero = 1e-7
  issues = '#34082'
  requirement = 'The system will assign the fuel pin linear heat rate in a quadrilateral lattice using an initial condition when the total power is supplied by a postprocessor, so that the linear heat rate integrated over the pins equals that power.'
[]

[SCMQuadPower_time_dependent_power]
  type = CSVDiff
  input = 'transient_power.i'
  csvdiff = 'transient_power_out.csv'
  recover = False
  abs_zero = 1e-6
  issues = '#34082'
  requirement = 'The system will assign the current value of a time-dependent total power postprocessor to the fuel pins in a quadrilateral lattice using an auxiliary kernel, without lagging it by one time step.'
[]

[SCMQuadPower_power_execute_on_warning]
  type = RunException
  input = 'transient_power.i'
  cli_args = 'Postprocessors/power/execute_on=timestep_end Executioner/num_steps=1'
  expect_err = 'The power postprocessor does not execute on every execute_on flag of this object, so q_prime may use a lagged power'
  issues = '#34082'
  requirement = 'The system will warn when the total power postprocessor of the quadrilateral lattice power auxiliary kernel does not execute on all of its execution flags.'
[]

[SCMQuadPower_power_force_postaux_warning]
  type = RunException
  input = 'transient_power.i'
  cli_args = 'Postprocessors/power/force_postaux=true Executioner/num_steps=1'
  expect_err = 'The power postprocessor sets force_postaux, so it executes after this object and q_prime uses a lagged power'
  issues = '#34082'
  requirement = 'The system will warn when the total power postprocessor of the quadrilateral lattice power auxiliary kernel is forced to execute after auxiliary kernels.'
[]

[SCMQuadPower_verbose]
  type = RunApp
  input = 'transient_power.i'
  cli_args = 'AuxKernels/q_prime_AUX/verbose=true Executioner/num_steps=1 Outputs/csv=false'
  expect_out = 'Aux Power correction vector'
  issues = '#34082'
  requirement = 'The system will report the per-pin power correction of the quadrilateral lattice power auxiliary kernel when requested.'
[]

[SCMQuadPower_receiver_power_no_warning]
  type = RunApp
  input = 'transient_power.i'
  # A Receiver with an execute_on and force_postaux that would otherwise trigger both warnings.
  # Its value is set by a transfer, so neither check applies and warnings remain errors here.
  cli_args = 'Postprocessors/received_power/type=Receiver Postprocessors/received_power/default=1e4 Postprocessors/received_power/execute_on=timestep_end Postprocessors/received_power/force_postaux=true AuxKernels/q_prime_AUX/power=received_power Executioner/num_steps=1 Outputs/csv=false'
  recover = False
  issues = '#34082'
  requirement = 'The system will not warn about the execution order of the total power postprocessor of the quadrilateral lattice power auxiliary kernel when that postprocessor is a Receiver filled by a transfer.'
[]

[SCMQuadPower_ic_power_not_initial_error]
  type = RunException
  input = 'postprocessor_power.i'
  cli_args = 'Postprocessors/power_pp/execute_on=timestep_end'
  expect_err = "The 'execute_on' parameter of the power postprocessor must include 'initial'"
  issues = '#34082'
  requirement = 'The system will report an error when the total power postprocessor of the quadrilateral lattice power initial condition does not execute on initial, since it would not have been computed when the initial condition is set.'
[]

[pin_power_on_initial]
  type = RunApp
  input = 'test_with_pins.i'
  cli_args = "Postprocessors/Total_power_SCMPinPowerPostprocessor/execute_on='initial timestep_end' Outputs/file_base=pin_power_on_initial"
  recover = False
  issues = '#34154'
  design = 'SCMPinPowerPostprocessor.md'
  requirement = 'The system shall compute the heat rate from the fuel pins to the coolant before the first subchannel solve.'
[]

[SCMTriPower_cells_aligned_with_pins]
  type = CSVDiff
  input = 'test_with_pins.i'
  capabilities = 'method!=dbg'
  csvdiff = 'test_with_pins_out.csv test_with_pins_out_line_check_0001.csv'
  recover = False
  abs_zero = 1e-7
  rel_err = 1e-5
  requirement = 'The system will populate a sinusoidal linear heat rate on pins in a triangular lattice using an initial condition and auxiliary kernel'
[]

[SCMTriPower_cells_not_aligned_with_pins]
  type = CSVDiff
  input = 'test_with_pins.i'
  capabilities = 'method!=dbg'
  csvdiff = 'not_aligned_with_pins.csv not_aligned_with_pins_line_check_0001.csv'
  cli_args = 'num_cells=31 SubChannel/T_relaxation=0.5 Outputs/file_base=not_aligned_with_pins'
  recover = False
  abs_zero = 1e-7
  rel_err = 2e-5
  requirement = 'The system will populate a sinusoidal linear heat rate on pins in a triangular lattice using an initial condition and auxiliary kernel in a mesh that has axial divisions that are not aligned with the heated region entrance and exit heights.'
[]

[SCMTriPower_pin_dependent_axial_shape]
  type = CSVDiff
  input = 'seven_pin_shapes.i'
  csvdiff = 'seven_pin_shapes_out.csv seven_pin_shapes_out_q_prime_profiles_0001.csv'
  recover = False
  abs_zero = 1e-7
  issues = '#34082'
  requirement = 'The system will assign a different axial linear heat rate shape to each pin in a triangular lattice from a single function of the pin position, while conserving the power assigned to each pin by the radial power distribution, using an initial condition and auxiliary kernel.'
[]

[SCMTriPower_postprocessor_power]
  type = CSVDiff
  input = 'postprocessor_power.i'
  csvdiff = 'postprocessor_power_out.csv'
  recover = False
  abs_zero = 1e-7
  issues = '#34082'
  requirement = 'The system will assign the pin linear heat rate in a triangular lattice using an initial condition when the total power is supplied by a postprocessor, so that the linear heat rate integrated over the pins equals that power.'
[]

[SCMTriPower_power_execute_on_warning]
  type = RunException
  input = 'postprocessor_power.i'
  # Adds the power aux kernel, executing on a flag the power postprocessor does not execute on
  cli_args = "AuxVariables/q_prime_aux/block=fuel_pins AuxKernels/q_prime_AUX/type=SCMTriPowerAux AuxKernels/q_prime_AUX/variable=q_prime_aux AuxKernels/q_prime_AUX/power=power_pp AuxKernels/q_prime_AUX/filename=seven_pin_power_profile.txt AuxKernels/q_prime_AUX/axial_heat_rate=axial_heat_rate Outputs/csv=false AuxKernels/q_prime_AUX/execute_on='initial timestep_end'"
  expect_err = 'The power postprocessor does not execute on every execute_on flag of this object, so q_prime may use a lagged power'
  issues = '#34082'
  requirement = 'The system will warn when the total power postprocessor of the triangular lattice power auxiliary kernel does not execute on all of its execution flags.'
[]

[SCMTriPower_power_force_postaux_warning]
  type = RunException
  input = 'postprocessor_power.i'
  cli_args = 'AuxVariables/q_prime_aux/block=fuel_pins AuxKernels/q_prime_AUX/type=SCMTriPowerAux AuxKernels/q_prime_AUX/variable=q_prime_aux AuxKernels/q_prime_AUX/power=power_pp AuxKernels/q_prime_AUX/filename=seven_pin_power_profile.txt AuxKernels/q_prime_AUX/axial_heat_rate=axial_heat_rate Outputs/csv=false AuxKernels/q_prime_AUX/execute_on=initial Postprocessors/power_pp/force_postaux=true'
  expect_err = 'The power postprocessor sets force_postaux, so it executes after this object and q_prime uses a lagged power'
  issues = '#34082'
  requirement = 'The system will warn when the total power postprocessor of the triangular lattice power auxiliary kernel is forced to execute after auxiliary kernels.'
[]

[SCMTriPower_verbose]
  type = RunApp
  input = 'postprocessor_power.i'
  cli_args = 'AuxVariables/q_prime_aux/block=fuel_pins AuxKernels/q_prime_AUX/type=SCMTriPowerAux AuxKernels/q_prime_AUX/variable=q_prime_aux AuxKernels/q_prime_AUX/power=power_pp AuxKernels/q_prime_AUX/filename=seven_pin_power_profile.txt AuxKernels/q_prime_AUX/axial_heat_rate=axial_heat_rate Outputs/csv=false AuxKernels/q_prime_AUX/execute_on=initial AuxKernels/q_prime_AUX/verbose=true'
  expect_out = 'Aux Power correction vector'
  issues = '#34082'
  requirement = 'The system will report the per-pin power correction of the triangular lattice power auxiliary kernel when requested.'
[]

[SCMTriPower_receiver_power_no_warning]
  type = RunApp
  input = 'postprocessor_power.i'
  # A Receiver with an execute_on and force_postaux that would otherwise trigger both warnings.
  # Its value is set by a transfer, so neither check applies and warnings remain errors here.
  cli_args = "AuxVariables/q_prime_aux/block=fuel_pins AuxKernels/q_prime_AUX/type=SCMTriPowerAux AuxKernels/q_prime_AUX/variable=q_prime_aux AuxKernels/q_prime_AUX/filename=seven_pin_power_profile.txt AuxKernels/q_prime_AUX/axial_heat_rate=axial_heat_rate Outputs/csv=false AuxKernels/q_prime_AUX/execute_on='initial timestep_end' Postprocessors/received_power/type=Receiver Postprocessors/received_power/default=7000 Postprocessors/received_power/execute_on=initial Postprocessors/received_power/force_postaux=true AuxKernels/q_prime_AUX/power=received_power"
  recover = False
  issues = '#34082'
  requirement = 'The system will not warn about the execution order of the total power postprocessor of the triangular lattice power auxiliary kernel when that postprocessor is a Receiver filled by a transfer.'
[]

[SCMTriPower_ic_power_not_initial_error]
  type = RunException
  input = 'postprocessor_power.i'
  cli_args = 'Postprocessors/power_pp/execute_on=timestep_end'
  expect_err = "The 'execute_on' parameter of the power postprocessor must include 'initial'"
  issues = '#34082'
  requirement = 'The system will report an error when the total power postprocessor of the triangular lattice power initial condition does not execute on initial, since it would not have been computed when the initial condition is set.'
[]

[blockedtest]
  type = CSVDiff
  input = 'blocked_test.i'
  csvdiff = 'blocked_test_out.csv'
  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'
[]

[Pressure_drop_counter_flow]
  type = CSVDiff
  input = 'one_pin_problem.i'
  csvdiff = 'counter_flow.csv'
  cli_args = 'Outputs/file_base=counter_flow'
  recover = False
  abs_zero = 1e-7
  requirement = 'The system will calculate the pressure drop across the pin with a counter_flow direction of gravity'
[]

[Pressure_drop_co_flow]
  type = CSVDiff
  input = 'one_pin_problem.i'
  csvdiff = 'co_flow.csv'
  cli_args = 'SubChannel/gravity=co_flow Outputs/file_base=co_flow'
  recover = False
  abs_zero = 1e-7
  requirement = 'The system will calculate the pressure drop across the pin with a co_flow direction of gravity'
[]

[Pressure_drop_no_gravity]
  type = CSVDiff
  input = 'one_pin_problem.i'
  csvdiff = 'no_gravity.csv'
  cli_args = 'SubChannel/gravity=none Outputs/file_base=no_gravity'
  recover = False
  abs_zero = 1e-7
  requirement = 'The system will calculate the pressure drop across the pin with no gravity'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system will calculate the geometric parameters for FCTF'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system shall calculate the flow area and wetted perimeter for a MARVEL-type microreactor core modeled in SCM'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system will apply initial massflow rate'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system will calculate the subchannel flow area for quadrilateral assemblies'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system will calculate the wetted perimeter of subchannels in a quadrilateral lattice'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system will calculate the subchannel flow area for triangular assemblies'
[]

[test]
  type = CSVDiff
  csvdiff = test_out.csv
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  abs_zero = 1e-6
  max_threads = 1
  input = 'test.i'
  requirement = 'The system will calculate the wetted perimeter of subchannels in a triangular lattice'
[]

[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 and fuel pins for quadrilateral assemblies (3x3)'
[]

[subchannel_number]
  type = 'RunException'
  input = 'coords.i'
  cli_args = 'Mesh/assembly/nx=1 Mesh/assembly/ny=1'
  expect_err = "The number of subchannels cannot be less than 2 in both directions\. Smallest assembly allowed is either 2X1 or 1X2\."
  requirement = 'The system shall check the minimum number of subchannels for detailed quadrilateral assemblies'
[]

[1x2]
  type = Exodiff
  input = coords.i
  exodiff = coords_1x2.e
  cli_args = "--mesh-only coords_1x2.e Mesh/assembly/nx=1 Mesh/assembly/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/assembly/nx=1 Mesh/assembly/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/assembly/nx=3 Mesh/assembly/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 subchannels and pins in triangular 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 and pins 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\. Smallest assembly allowed is either 2X1 or 1X2\."
  issues = '#29224'
  requirement = 'The system shall check the minimum number of subchannels'
[]

[test]
  type = Exodiff
  input = coords.i
  exodiff = coords_in.e
  cli_args = "--mesh-only"
  recover = false
  mesh_mode = 'replicated'
  requirement = 'The system will be able to create a rectangular duct mesh for square subchannel 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'
[]

[tricoords2]
  type = RunApp
  input = tricoords2.i
  cli_args = "--mesh-only"
  recover = false
  requirement = 'The system will create a 7-pin subchannel 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 be able to create a hexagonal duct mesh for triangular subchannel 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
  capabilities = 'method!=dbg'
[]

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

[pin_pos]
  type = JSONDiff
  input = 'pin_positions.i'
  jsondiff = 'pin_positions_out.json'
  requirement = 'The system shall be able to generate the list of positions for the subchannel assembly fuel pins.'
  max_threads = 1
  recover = false # steady solve
  # Distributed not supported
  mesh_mode = 'replicated'
[]

[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
  rel_err = 1e-5
  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
  requirement = 'The system will examine the tri subchannel solver for lead coolant'
[]

[pacio_bare_pins]
  type = RunException
  input = test_LBE-19pin.i
  cli_args = 'SCMClosures/Chen/friction_model=Pacio'
  expect_err = "The Pacio-Chen-Todreas friction correlation applies only to wire-wrapped triangular assemblies\."
  issues = '#33608'
  requirement = 'The system shall reject the Pacio-Chen-Todreas friction correlation for a bare pin triangular assembly.'
[]

[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
  rel_err = 3e-5
  # mdot-35 is the flow rate in a corner subchannel that is small relative to the assembly
  # flow and set by the iteratively-relaxed inter-channel crossflow (Wij), so it is more
  # sensitive to build/platform floating-point noise than the rest of the postprocessors
  override_columns = 'mdot-35'
  override_rel_err = '5e-5'
  override_abs_zero = '1e-6'
  max_threads = 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
  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
  # Tolerance added due to CIVET machine difference
  # See PR: Qprime #31380
  # and PR: various mods #31283
  rel_err = 3e-03
  allow_warnings = true
  requirement = 'The system will examine the subchannel implicit solver for a 19pin liquid sodium cooled assembly'
[]

[test_radial_conduction]
  type = CSVDiff
  input = test19_radial_conduction.i
  csvdiff = test19_radial_conduction_out.csv
  valgrind = NONE
  recover = false
  issues = '#33956'
  requirement = 'The system shall spread heat from a heated center pin to the outer subchannels of a 19-pin liquid sodium cooled assembly through radial heat conduction when turbulent mixing is disabled.'
[]

[invalid_segregated_relaxation]
  type = RunException
  input = test19_implicit.i
  cli_args = 'SubChannel/mass_flow_relaxation=0.5 Outputs/exodus=false Outputs/csv=false'
  expect_err = "The 'mass_flow_relaxation' parameter is only used by the implicit non-segregated solve\."
  requirement = 'The system shall reject monolithic relaxation controls when using a segregated solve.'
[]

[invalid_explicit_nonsegregated]
  type = RunException
  input = test19_explicit.i
  cli_args = 'SubChannel/segregated=false Outputs/exodus=false Outputs/csv=false'
  expect_err = "A non-segregated solve requires 'implicit = true'\."
  requirement = 'The system shall reject a non-segregated solve when implicit assembly is disabled.'
[]

[test_implicit_tiny_dt]
  type = RunApp
  input = test19_monolithic.i
  cli_args = 'Executioner/type=Transient Executioner/dt=1e-5 Executioner/end_time=1e-5 Executioner/abort_on_solve_fail=true SubChannel/interpolation_scheme=upwind SubChannel/enthalpy_subcycles=2 SubChannel/rtol=1e-10 SubChannel/atol=1e-12 SubChannel/mass_flow_relaxation=0.3 SubChannel/pressure_relaxation=0.8 Outputs/exodus=false Outputs/csv=false'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  max_threads = 1
  max_parallel = 1
  requirement = 'The coupled subchannel solver shall converge a transient sodium case with a 10 microsecond time step, upwind spatial interpolation, and enthalpy subcycling.'
[]

[test_saturated_sodium_tiny_dt]
  type = RunApp
  input = test19_monolithic.i
  cli_args = 'FluidProperties/sodium/type=SodiumSaturationFluidProperties Executioner/type=Transient Executioner/dt=1e-5 Executioner/end_time=1e-5 Executioner/abort_on_solve_fail=true SubChannel/P_tol=1e-5 SubChannel/interpolation_scheme=upwind SubChannel/enthalpy_subcycles=2 SubChannel/rtol=1e-10 SubChannel/atol=1e-12 SubChannel/mass_flow_relaxation=0.3 SubChannel/pressure_relaxation=0.8 Outputs/exodus=false Outputs/csv=false'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  allow_warnings = true
  max_threads = 1
  max_parallel = 1
  requirement = 'The coupled subchannel solver shall converge a transient saturated sodium case with a 10 microsecond time step, upwind spatial interpolation, and enthalpy subcycling.'
[]

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

[test_monolithic_pacio]
  type = CSVDiff
  input = test19_monolithic_pacio.i
  csvdiff = test19_monolithic_pacio_out.csv
  abs_zero = 1e-7
  rel_err = 7e-06
  valgrind = NONE
  recover = false
  design = 'SCMFrictionChenTodreas.md SCMMixingChenTodreas.md'
  issues = '#33608'
  requirement = 'The system will examine the relaxed subchannel monolithic solver for a 19pin liquid sodium cooled assembly with the Pacio-Chen-Todreas friction and mixing parameterizations'
[]

[case1_scm]
  type = CSVDiff
  input = case1_scm.i
  csvdiff = case1_scm_out.csv
  cli_args = '--error --compute-device=cpu'
  rel_err = 5e-5
  allow_warnings = true
  recover = false
  requirement = 'The system shall solve a 7-pin triangular subchannel assembly and preserve the center, edge, and corner outlet temperature and mass flow rate.'
[]

[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
  requirement = 'The system will examine the coupling using SCM as main app'
[]

[quad_pin_diameter_ic_recomputes_geometry]
  type = CSVDiff
  input = 'quad_pin_diameter.i'
  csvdiff = 'quad_pin_diameter_out.csv'
  recover = False
  requirement = 'The system shall recompute quadrilateral subchannel flow area and wetted perimeter from the Dpin initial condition when it differs from the mesh-generator pin diameter.'
[]

[tri_pin_diameter_ic_recomputes_geometry]
  type = CSVDiff
  input = 'tri_pin_diameter.i'
  capabilities = 'method!=dbg'
  csvdiff = 'tri_pin_diameter_out.csv'
  recover = False
  requirement = 'The system shall recompute triangular subchannel flow area and wetted perimeter from the Dpin initial condition when it differs from the mesh-generator pin diameter.'
[]

[dittus_boelter]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'dittus_boelter.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the default Dittus Boelter correlation for circular tubes.'
  cli_args = 'SubChannel/pin_HTC_closure=dittus-boelter SubChannel/duct_HTC_closure=dittus-boelter Outputs/file_base=dittus_boelter FluidProperties/sodium/type=Water97FluidProperties SCMClosures/dittus-boelter/correction_factor=none T_in=300.0'
[]

[dittus_boelter_Presser]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'dittus_boelter_Presser.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Dittus Boelter correlation and the Presser correction factor.'
  cli_args = 'SubChannel/pin_HTC_closure=dittus-boelter SubChannel/duct_HTC_closure=dittus-boelter Outputs/file_base=dittus_boelter_Presser FluidProperties/sodium/type=Water97FluidProperties SCMClosures/dittus-boelter/correction_factor=Presser T_in=300.0'
[]

[dittus_boelter_Weisman]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'dittus_boelter_Weisman.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Dittus Boelter correlation and the Weisman correction factor.'
  cli_args = 'SubChannel/pin_HTC_closure=dittus-boelter SubChannel/duct_HTC_closure=dittus-boelter Outputs/file_base=dittus_boelter_Weisman FluidProperties/sodium/type=Water97FluidProperties SCMClosures/dittus-boelter/correction_factor=Weisman T_in=300.0'
[]

[gnielinski]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'gnielinski.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Gnielinski correlation.'
  cli_args = 'SubChannel/pin_HTC_closure=gnielinski SubChannel/duct_HTC_closure=gnielinski Outputs/file_base=gnielinski'
[]

[invalid_wire_diameter_error]
  type = RunException
  allow_test_objects = true
  input = 'XX09_SCM_SS17.i'
  expect_err = 'The wire diameter must not exceed the pin-to-pin gap'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'reject wire diameters that exceed the pin-to-pin gap.'
  cli_args = 'SubChannel/pin_HTC_closure=gnielinski SubChannel/duct_HTC_closure=gnielinski wire_diameter=0.003 Outputs/file_base=invalid_wire_diameter_error'
[]

[negative_htc_error]
  type = RunException
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  expect_err = 'The heat transfer coefficient must be non-negative and finite\.'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'and reject non-physical heat transfer coefficients from HTC closures.'
  cli_args = 'SubChannel/pin_HTC_closure=negative-htc SubChannel/duct_HTC_closure=negative-htc Outputs/file_base=negative_htc_error'
[]

[kazimi_carelli]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'kazimi-carelli.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Kazimi-Carelli correlation.'
  cli_args = 'SubChannel/pin_HTC_closure=kazimi-carelli SubChannel/duct_HTC_closure=gnielinski Outputs/file_base=kazimi-carelli'
[]

[schad_modified]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'schad-modified.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Schad-Modified correlation.'
  cli_args = 'SubChannel/pin_HTC_closure=schad-modified SubChannel/duct_HTC_closure=gnielinski Outputs/file_base=schad-modified'
[]

[graber_rieger]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'graber-rieger.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Graber-Rieger correlation.'
  cli_args = 'SubChannel/pin_HTC_closure=graber-rieger SubChannel/duct_HTC_closure=gnielinski Outputs/file_base=graber-rieger'
[]

[borishanskii]
  type = CSVDiff
  input = 'XX09_SCM_SS17.i'
  allow_test_objects = true
  csvdiff = 'borishanskii.csv'
  capabilities = 'method!=dbg'
  abs_zero = 1e-5
  rel_err = 5e-5
  valgrind = NONE
  recover = false
  max_parallel = 1
  detail = 'using the Borishanskii correlation.'
  cli_args = 'SubChannel/pin_HTC_closure=borishanskii SubChannel/duct_HTC_closure=gnielinski Outputs/file_base=borishanskii'
[]

[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
  # The crossflow solve does not converge in every iteration; remove once #34306 is fixed
  allow_warnings = true
  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
  # The crossflow solve does not converge in every iteration; remove once #34306 is fixed
  allow_warnings = true
  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
  # The crossflow solve does not converge in every iteration; remove once #34306 is fixed
  allow_warnings = true
  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
  # The crossflow solve does not converge in every iteration; remove once #34306 is fixed
  allow_warnings = true
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using an implicit algorithm.'
[]

[psbt_pressure_iteration_exhaustion]
  type = RunException
  input = psbt_monolithic.i
  cli_args = 'SubChannel/P_maxit=1 Outputs/file_base=psbt_monolithic_p_maxit'
  expect_err = 'Reached maximum number of axial pressure iterations.*Solve Did NOT Converge!'
  capabilities = 'method!=dbg'
  valgrind = NONE
  recover = false
  max_threads = 1
  issues = '#34272'
  requirement = 'The system shall report nonconvergence and return a failing exit code when the pressure iteration limit is reached.'
[]

[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
  requirement = 'The system shall be able to solve a PSBT-type case with the subchannel solver using a monolithic algorithm.'
[]

[steady]
  type = Exodiff
  input = steady.i
  exodiff = steady_out.e
  recover = false
  abs_zero = 1e-7
  max_threads = 1
  # The crossflow solve does not converge in every iteration; remove once #34306 is fixed
  allow_warnings = true
  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
  # The crossflow solve does not converge in every iteration; remove once #34306 is fixed
  allow_warnings = true
  requirement = 'The system shall be able to restart after a transient solution'
[]

[transient_time_integrator_warning]
  type = RunApp
  input = transient.i
  prereq = 'steady'
  cli_args = 'Executioner/scheme=bdf2 Outputs/file_base=transient_time_integrator_warning'
  allow_warnings = true
  expect_out = "The subchannel solver always uses implicit (backward) Euler time integration; the requested 'BDF2' time integrator is ignored."
  match_literal = true
  max_threads = 1
  requirement = 'The system shall warn when a subchannel transient input requests a MOOSE time integrator other than implicit Euler.'
[]

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

[wrong_transfer_type]
  type = RunException
  input = 'quad.i'
  cli_args = "Transfers/T_transfer/transfer_type=subchannel Outputs/file_base=wrong_transfer_type"
  expect_err = "The variable 'T' does not have DOFs on subchannel nodes"
  requirement = 'The system shall report an error when a variable is transferred with the wrong subchannel transfer type'
[]

[legacy_pin_transfer]
  type = RunApp
  input = 'legacy_pin_transfer.i'
  recover = false
  requirement = 'The system shall preserve pin-transfer behavior when using the deprecated pin solution transfer name'
[]