Thermal Hydraulics System Design Description

[without_positions]
  type = Exodiff
  input = 'master.i'
  exodiff = 'master_out.e master_out_thm0_out.e'

  detail = 'Without sub-app positions provided.'
[]

[with_positions]
  type = Exodiff
  input = 'master.i'
  cli_args = "CoupledHeatTransfers/right/positions='0 0 0'"
  prereq = 'test/without_positions'
  exodiff = 'master_out.e master_out_thm0_out.e'

  detail = 'With sub-app positions provided.'
[]

[multiple_phases]
  type = Exodiff
  input = 'master.i'
  cli_args = "
        AuxVariables/kappa/family=MONOMIAL
        AuxVariables/kappa/order=CONSTANT
        AuxVariables/kappa/initial_condition=0
        CoupledHeatTransfers/right/T_fluid='T_fluid T_fluid'
        CoupledHeatTransfers/right/htc='htc htc'
        CoupledHeatTransfers/right/kappa='kappa kappa'
        CoupledHeatTransfers/right/T_fluid_user_objects='T_uo T_uo'
        CoupledHeatTransfers/right/htc_user_objects='Hw_uo Hw_uo'
        CoupledHeatTransfers/right/kappa_user_objects='kappa_uo kappa_uo'
        MultiApps/thm/input_files=sub_2phase.i
        Outputs/show='T T_fluid htc T_wall_avg T_fluid_avg htc_avg'"
  prereq = 'test/with_positions'
  exodiff = 'master_out.e master_out_thm0_out.e'

  detail = 'With multiple phases.'
[]

[misaligned]
  type = RunException
  input = 'master.i'
  cli_args = "CoupledHeatTransfers/right/orientation='0.01 1 0'"
  expect_err = "The direction given by the parameter 'orientation' must be aligned with the x, y, or z axis"

  requirement = 'The system shall report an error for the coupled heat transfer action if the mesh is not aligned with the x, y, or z axis.'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-10
  requirement = 'The system shall compute convective heat flux between fluid and wall temperature for 1-phase flow'
  design = 'ConvectiveHeatFlux1PhaseAux.md'
  issues = '#60'
[]

[test]
  type = 'RunApp'
  input = 'test.i'
  expect_out = "Loop 1:
\s+hx/primary
\s+pri_inlet
\s+pri_outlet

\s+Loop 2:
\s+hx/secondary
\s+sec_inlet
\s+sec_outlet"

  requirement = 'The system shall allow nesting components into groups in input files'
  design = 'ComponentGroup.md'
  issues = '#94'
[]

[free_boundary]
  type = 'RunException'
  input = 'free_boundary.i'
  expect_err = "Deprecated component. Use FreeBoundary1Phase or FreeBoundary2Phase instead."
  design = 'FreeBoundary.md'
  requirement = 'The system shall report an error if the FreeBoundary component is used.'
[]

[gate_valve]
  type = 'RunException'
  input = 'gate_valve.i'
  expect_err = "Deprecated component. Use GateValve1Phase or GateValve2Phase instead."
  design = 'GateValve.md'
  requirement = 'The system shall report an error if the GateValve component is used.'
[]

[solid_wall]
  type = 'RunException'
  input = 'solid_wall.i'
  expect_err = "Deprecated component. Use SolidWall1Phase or SolidWall2Phase instead."
  design = 'SolidWall.md'
  requirement = 'The system shall report an error if the SolidWall component is used.'
[]

[junction_one_to_one]
  type = 'RunException'
  input = 'junction_one_to_one.i'
  expect_err = "Deprecated component. Use JunctionOneToOne1Phase or JunctionOneToOne2Phase instead."
  design = 'JunctionOneToOne.md'
  requirement = 'The system shall report an error if the JunctionOneToOne component is used.'
[]

[heat_generation]
  type = 'RunException'
  input = 'heat_generation.i'
  expect_err = "'HeatGeneration' component is deprecated, use 'HeatSourceFromTotalPower' or 'HeatSourceFromPowerDensity' instead."
  design = 'HeatGeneration.md'
  requirement = 'The system shall report an error if the HeatGeneration component is used.'
[]

[heat_source_volumetric]
  type = 'RunException'
  input = 'heat_source_volumetric.i'
  expect_err = "Deprecated component. Use HeatSourceVolumetric1Phase or HeatSourceVolumetric2Phase instead."
  design = 'HeatSourceVolumetric.md'
  requirement = 'The system shall report an error if the HeatSourceVolumetric component is used.'
[]

[prescribed_reactor_power]
  type = 'RunException'
  input = 'prescribed_reactor_power.i'
  expect_err = "total_power: Deprecated component, use 'type = TotalPower' instead."
  design = 'PrescribedReactorPower.md'
  requirement = 'The system shall report an error if the PrescribedReactorPower component is used.'
[]

[mesh]
  type = RunApp
  input = 'mesh.i'
  cli_args = '--mesh-only'
  # Mesh-only runs should not test recovery
  recover = False

  requirement = 'The system shall generate a mesh for the FileMeshComponent test.'
[]

[test]
  type = Exodiff
  prereq = 'mesh'
  input = 'file_mesh_component.i'
  exodiff = 'file_mesh_component.e'
  # Prereq does not test recover
  recover = False

  requirement = 'The system shall provide a component that loads a mesh from an ExodusII file.'
[]

[file_not_readable]
  type = RunException
  input = 'file_mesh_component.i'
  cli_args = 'Components/hs_external/file=nonexistent_file.e'
  expect_err = "The file .* could not be opened."

  requirement = 'The system shall report an error for FileMeshComponent when the file is not readable.'
[]

[generate_mesh]
  type = Exodiff
  input = rectangle.i
  exodiff = rectangle.e
  recover = false # No transient, just mesh generation
  cli_args = "--mesh-only rectangle.e"
  requirement = 'The system shall generate a mesh for the FlowComponentNS tests.'
[]

[pincns]
  type = Exodiff
  input = flow_component_ns_pinc.i
  exodiff = flow_component_ns_pinc_out.e
  abs_zero = 5e-6
  prereq = 'generate_mesh'
  method = "!dbg"
  valgrind = HEAVY
  requirement = 'The system shall model the porous, incompressible Navier-Stokes equations with a finite volume discretization, using a component.'
[]

[cylindrical]
  type = Exodiff
  input = 'heat_structure_2d_coupler.i'
  cli_args = "
      Components/hs3/inner_radius=0.7
      Postprocessors/E_tot/axis_dir='1 0 0'
      Postprocessors/E_tot/axis_point='0 0 0'"
  exodiff = 'cylindrical.e'
  requirement = 'The system shall be able to couple two 2D cylindrical heat structures.'
[]

[plate]
  type = Exodiff
  input = 'heat_structure_2d_coupler.i'
  cli_args = "
      Components/hs1/type=HeatStructurePlate
      Components/hs2/type=HeatStructurePlate
      Components/hs3/type=HeatStructurePlate
      Components/hs1/depth=1.0
      Components/hs2/depth=1.0
      Components/hs3/depth=1.0
      Components/hs3/position='0.5 0.7 0'
      Postprocessors/E_tot/type=ADHeatStructureEnergy
      Postprocessors/E_tot/plate_depth=1.0
      Outputs/file_base=plate"
  exodiff = 'plate.e'
  requirement = 'The system shall be able to couple two 2D plate heat structures.'
[]

[separated]
  type = CSVDiff
  input = 'separated.i'
  csvdiff = 'separated_out.csv'
  issues = '#19851 #23853'
  requirement = 'The system shall be able to couple two 2D cylindrical heat structures on separated surfaces.'
[]

[missing_boundary]
  type = RunException
  input = 'heat_structure_2d_coupler.i'
  cli_args = "
        Components/hs3/inner_radius=0.7
        Postprocessors/E_tot/axis_dir='1 0 0'
        Postprocessors/E_tot/axis_point='0 0 0'
        Components/hs_coupling_1_2/secondary_boundary=hs3:inner"
  expect_err = "The heat structure 'hs1' does not have the boundary 'hs3:inner'"

  detail = 'the provided heat structure boundary does not exist.'
[]

[type_mismatch]
  type = RunException
  input = 'heat_structure_2d_coupler.i'
  cli_args = "
        Components/hs3/inner_radius=0.7
        Postprocessors/E_tot/axis_dir='1 0 0'
        Postprocessors/E_tot/axis_point='0 0 0'
        Components/hs1/type=HeatStructurePlate
        Components/hs1/depth=1.0"
  expect_err = 'The coupled heat structures must have the same type'

  detail = 'the types of the coupled heat structures do not match.'
[]

[invalid_hs_type]
  type = RunException
  input = 'heat_structure_2d_coupler.i'
  cli_args = "
        Components/hs3/inner_radius=0.7
        Postprocessors/E_tot/axis_dir='1 0 0'
        Postprocessors/E_tot/axis_point='0 0 0'
        Components/hs4/type=HeatStructureFromFile3D
        Components/hs4/file=../heat_structure_from_file_3d/box.e
        Components/hs4/position='0 0 0'
        Components/hs4/initial_T=300
        Materials/mat/type=ADGenericConstantMaterial
        Materials/mat/block='hs4:brick'
        Materials/mat/prop_names='density specific_heat thermal_conductivity'
        Materials/mat/prop_values='1 1 1'
        Components/hs_coupling_1_2/primary_heat_structure=hs4"
  expect_err = "The component 'hs4' is not of type 'HeatStructureBase'"

  detail = 'the types of either coupled heat structure is invalid.'
[]

[mesh_mismatch]
  type = RunException
  input = 'heat_structure_2d_coupler.i'
  cli_args = "
        Components/hs1/orientation='1 1 0'
        Components/hs3/inner_radius=0.8
        Postprocessors/E_tot/axis_dir='1 0 0'
        Postprocessors/E_tot/axis_point='0 0 0'"
  expect_err = "The primary and secondary boundaries must be aligned"

  detail = 'the boundary meshes are not aligned.'
[]

[physics]
  type = CSVDiff
  input = 'heat_structure_2d_radiation_coupler_rz.i'
  csvdiff = 'heat_structure_2d_radiation_coupler_rz.csv'
  abs_zero = 1e-5
  requirement = 'The system shall be able to couple two 2D cylindrical heat structures via radiation and conserve energy.'
[]

[missing_boundary]
  type = RunException
  input = 'heat_structure_2d_radiation_coupler_rz.i'
  cli_args = "Components/hs_coupler/primary_boundary=fake_boundary"
  expect_err = "The heat structure 'hs1' does not have the boundary 'fake_boundary'"

  detail = 'the provided heat structure boundary does not exist.'
[]

[invalid_hs_type]
  type = RunException
  input = 'heat_structure_2d_radiation_coupler_rz.i'
  cli_args = "
        Components/hs1/type=HeatStructurePlate
        Components/hs1/depth=1.0"
  expect_err = "The primary and secondary heat structures must be of a type inherited from 'HeatStructureCylindricalBase'"

  detail = 'the type of either coupled heat structure is invalid.'
[]

[mesh_not_aligned]
  type = RunException
  input = 'heat_structure_2d_radiation_coupler_rz.i'
  cli_args = "Components/hs1/position='0.1 0.1 0'"
  expect_err = "The primary and secondary boundaries must be aligned"

  detail = 'the boundary meshes are not aligned.'
[]

[phy:T_wall_transfer_3eqn_x]
  type = Exodiff
  input = 'phy.T_wall_transfer_3eqn_x.i'
  exodiff = 'phy.T_wall_transfer_3eqn_x_out.e'
  recover = false

  requirement = "The system shall be able to couple a flow channel and heat structure aligned with the x-axis."
  issues = '#19754'
[]

[phy:T_wall_transfer_3eqn_x_lengths]
  type = Exodiff
  input = 'phy.T_wall_transfer_3eqn_x.i'
  exodiff = 'phy.T_wall_transfer_3eqn_x_lengths.e'
  cli_args = "
      Components/pipe1/length='0.6 0.3 0.1'
      Components/pipe1/n_elems='36 18 6'
      Components/hs/length='0.1 0.3 0.6'
      Components/hs/n_elems='6 18 36'
      Components/hs/axial_region_names='1 3 6'
      Outputs/file_base=phy.T_wall_transfer_3eqn_x_lengths"
  recover = false

  requirement = "The system shall be able to couple a flow channel and heat structure with non-uniform meshes and opposite directions."
  issues = '#19754'
[]

[phy:T_wall_transfer_3eqn_y]
  type = Exodiff
  input = 'phy.T_wall_transfer_3eqn_y.i'
  exodiff = 'phy.T_wall_transfer_3eqn_y_out.e'
  recover = false

  requirement = "The system shall be able to couple a flow channel and heat structure aligned with the y-axis."
  issues = '#20747'
[]

[phy:T_wall_transfer_3eqn_z]
  type = Exodiff
  input = 'phy.T_wall_transfer_3eqn_z.i'
  exodiff = 'phy.T_wall_transfer_3eqn_z_out.e'
  recover = false

  requirement = "The system shall be able to couple a flow channel and heat structure aligned with the z-axis."
  issues = '#19754'
[]

[phy:conservation_1phase_plate]
  type = CSVDiff
  input = 'phy.conservation_1phase.i'
  csvdiff = 'phy.conservation_1phase_plate_rdg.csv'
  abs_zero = 1e-12
  cli_args = "
      Components/heat_transfer/P_hf=2
      Components/heat_structure/type=HeatStructurePlate
      Components/heat_structure/depth=2
      Postprocessors/E_heat_structure/type=ADHeatStructureEnergy
      Postprocessors/E_heat_structure/plate_depth=2
      Outputs/file_base=phy.conservation_1phase_plate_rdg"
  # because of output has execute_on = 'initial final'
  recover = false

  requirement = "The system shall conserve energy when a flow channel is coupled to a plate heat structure."
  issues = '#19754'
[]

[phy:conservation_1phase_cylinder]
  type = CSVDiff
  input = 'phy.conservation_1phase.i'
  csvdiff = 'phy.conservation_1phase_cylinder_rdg.csv'
  abs_zero = 1e-12
  cli_args = "
      Components/pipe/position='0 2 0'
      Components/heat_transfer/P_hf=1.5
      Components/heat_structure/type=HeatStructureCylindrical
      Components/heat_structure/inner_radius=1.0
      Components/heat_structure/offset_mesh_by_inner_radius=true
      Postprocessors/E_heat_structure/type=ADHeatStructureEnergyRZ
      Postprocessors/E_heat_structure/axis_point='0 2 0'
      Postprocessors/E_heat_structure/axis_dir='1 0 0'
      Outputs/file_base=phy.conservation_1phase_cylinder_rdg"
  # because of output has execute_on = 'initial final'
  recover = false

  requirement = "The system shall conserve energy when a flow channel is coupled to a cylindrical heat structure."
  issues = '#19754'
[]

[phy:heat_structure_multiple_3eqn]
  type = 'CSVDiff'
  input = 'phy.heat_structure_multiple_3eqn.i'
  csvdiff = 'phy.heat_structure_multiple_3eqn_rdg_out.csv'
  cli_args = "
      Outputs/file_base=phy.heat_structure_multiple_3eqn_rdg_out"
  abs_zero = 1e-5
  rel_err = 0
  max_time = 6000
  recover = false

  requirement = "The system shall conserve energy when a flow channel is coupled to several heat structures."
  issues = '#19754'

[]

[phy.energy_heatstructure_ss_1phase]
  type = 'CSVDiff'
  input = 'phy.energy_heatstructure_ss_1phase.i'
  csvdiff = 'phy.energy_heatstructure_ss_1phase_out.csv'
  abs_zero = 5e-7
  recover = false

  requirement = "The system shall conserve energy after reaching steady-state when a flow channel is coupled to a heat structure."
  issues = '#19754'
[]

[err:not_a_pipe]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/hx/flow_channel=inlet Components/hx/Hw=300"
  expect_err = "hx: The component 'inlet' is not of type 'FlowChannelBase'"

  requirement = "The system shall throw an error if the flow channel component is not of type 'FlowChannelBase'."
  issues = '#19754'
[]

[err:not_a_hs]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/hx/hs=inlet Components/hx/Hw=300 "
  expect_err = "hx: The component 'inlet' is not of type 'HeatStructureBase'"

  requirement = "The system shall throw an error if the heat structure component is not of type 'HeatStructureBase'."
  issues = '#19754'
[]

[err:wrong_hs_side]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/hx/hs_side=asdf"
  expect_err = "hx: The parameter 'hs_side' was given an invalid value"

  requirement = "The system shall throw an error if the provided heat structure side is invalid."
  issues = '#19754'
[]

[err:elems_mismatch]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/pipe/n_elems=1 Components/hx/Hw=300"
  expect_err = "hx: The number of elements in component 'pipe' is 1, but the number of axial elements in component 'hs' is 2. They must be the same."

  requirement = "The system shall throw an error if the heat structure and flow channel components don't have the same number of axial elements."
  issues = '#19754'
[]

[err:length_mismatch]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/pipe/length=3 Components/hx/Hw=300"
  expect_err = "hx: The length of component 'pipe' is 3, but the length of component 'hs' is 4. They must be the same."

  requirement = "The system shall throw an error if the heat structure and flow channel components don't have the same length."
  issues = '#19754'
[]

[err:wrong_position]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/pipe/position='0 0.2 0.5' Components/hx/Hw=300 Components/pipe/n_elems=10 Components/hs/n_elems=10"
  expect_err = "hx: The centers of the elements of flow channel 'pipe' do not align with the centers of the specified heat structure side."

  requirement = "The system shall throw an error if the center of the elements of the flow channel component don't align with the centers of the specified heat structure side."
  issues = '#19754'
[]

[err:wrong_orientation]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/pipe/orientation='0 1 0' Components/hx/Hw=300"
  expect_err = "hx: The centers of the elements of flow channel 'pipe' do not align with the centers of the specified heat structure side."

  requirement = "The system shall throw an error if the coupled flow channel and heat structure components don't have the same orientation."
  issues = '#19754'
[]

[err:zero_p_hs_radius]
  type = 'RunException'
  input = 'err.1phase.i'
  cli_args = "Components/pipe/position='0 0.1 0' Components/hx/hs_side=inner Components/hx/Hw=300"
  expect_err = "hx: 'hs_side' parameter is set to 'INNER', but this side of the heat structure 'hs' has radius of zero."

  requirement = "The system shall throw an error if the flow channel is coupled to the inner side of a heat structure that has a zero inner radius."
  issues = '#19754'
[]

[err:missing_hw]
  type = 'RunException'
  input = 'err.1phase.i'
  expect_err = "hx: The parameter 'Hw' must be provided when using simple closures."

  requirement = "The system shall throw an error if the heat transfer coefficient is not specified with simple closures."
  issues = '#19754'
[]

[jac:cylindrical_top_side_1phase]
  type = 'PetscJacobianTester'
  input = 'jac.1phase.i'
  ratio_tol = 1e-5
  difference_tol = 55
  cli_args = "Debug/check_jacobian=true Components/hs/type=HeatStructureCylindrical Components/hx/hs_side=outer"
  max_parallel = 1

  requirement = "The system shall compute jacobians when a flow channel is coupled to the outer side of a cylindrical heat structure."
  issues = '#19754'
[]

[jac:cylindrical_bottom_side_1phase]
  type = 'PetscJacobianTester'
  input = 'jac.1phase.i'
  ratio_tol = 1e-5
  difference_tol = 55
  cli_args = "Debug/check_jacobian=true Components/pipe/position='0 0 0' Components/hs/inner_radius=0.1 Components/hs/offset_mesh_by_inner_radius=true Components/hs/type=HeatStructureCylindrical Components/hx/hs_side=inner"
  max_parallel = 1

  requirement = "The system shall compute jacobians when a flow channel is coupled to the inner side of a cylindrical heat structure."
  issues = '#19754'
[]

[jac:plate_top_side_1phase]
  type = 'PetscJacobianTester'
  input = 'jac.1phase.i'
  ratio_tol = 1e-5
  difference_tol = 55
  cli_args = "Debug/check_jacobian=true Components/hs/type=HeatStructurePlate Components/hs/depth=1 Components/hx/hs_side=outer"
  max_parallel = 1

  requirement = "The system shall compute jacobians when a flow channel is coupled to the outer side of a plate heat structure."
  issues = '#19754'
[]

[jac:plate_bottom_side_1phase]
  type = 'PetscJacobianTester'
  input = 'jac.1phase.i'
  ratio_tol = 1e-5
  difference_tol = 55
  cli_args = "Debug/check_jacobian=true Components/pipe/position='0 0 0' Components/hs/depth=1 Components/hs/type=HeatStructurePlate Components/hx/hs_side=inner"
  max_parallel = 1

  requirement = "The system shall compute jacobians when a flow channel is coupled to the inner side of a plate heat structure."
  issues = '#19754'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in volume junctions'
  design = 'HeatTransferFromHeatStructure1Phase.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in volume junctions'
  design = 'HeatTransferFromHeatStructure1Phase.md SolutionIC.md'
  issues = '#20553'
[]

[phy:conservation]
  type = CSVDiff
  input = 'phy.conservation.i'
  csvdiff = 'phy.conservation.csv'
  abs_zero = 1e-9
  recover = false

  requirement = 'The system shall conserve energy when using HeatTransferFromHeatStructure3D1Phase.'
[]

[phy:conservation_ss]
  type = CSVDiff
  input = 'phy.conservation_ss.i'
  csvdiff = 'phy.conservation_ss.csv'
  abs_zero = 1e-6
  max_time = 6000
  recover = false

  requirement = 'The system shall allow to connect multiple flow channels to a single boundary in HeatTransferFromHeatStructure3D1Phase.'
  issues = '#19831 #20298'
[]

[phy:conservation_inv]
  type = CSVDiff
  input = 'phy.conservation.i'
  csvdiff = 'phy.conservation.csv'
  cli_args = "Components/fch1/position='0 0 1' Components/fch1/orientation='0 0 -1'"
  abs_zero = 1e-6
  max_time = 6000
  prereq = phy:conservation_ss

  requirement = 'The system shall allow to connect flow channels that have negative orientation to a HeatTransferFromHeatStructure3D1Phase component.'
  issues = '#19831 #19879'
[]

[err:not_a_pipe]
  type = 'RunException'
  input = 'phy.conservation.i'
  cli_args = "Components/ht/flow_channels='in1'"
  expect_err = "ht: The component 'in1' is not of type 'FlowChannel1Phase'."

  requirement = "The system shall throw an error if a flow channel connected to a HeatTransferFromHeatStructure3D1Phase component is not a FlowChannel1Phase."
[]

[err:fch_orientation]
  type = 'RunException'
  input = 'phy.conservation.i'
  cli_args = "Components/fch1/orientation='1 0 1'"
  expect_err = "ht: The flow channel 'fch1' must be aligned with the x-, y-, or z- axis."

  requirement = "The system shall throw an error if a flow channel connected to a HeatTransferFromHeatStructure3D1Phase component is not aligned with the x-, y-, or z- axis."
[]

[err:not_3d_hs]
  type = 'RunException'
  input = 'err.not_a_3d_hs.i'
  expect_err = "ht: The component 'blk' is not a HeatStructureFromFile3D component."

  requirement = "The system shall throw an error if the heat structure connected to a HeatTransferFromHeatStructure3D1Phase component is not a HeatStructureFromFile3D component."
[]

[err:non_existent_boundary]
  type = 'RunException'
  input = 'phy.conservation.i'
  cli_args = "Components/ht/boundary=asdf"
  expect_err = "ht: The boundary 'asdf' \(.+\) was not found."

  requirement = "The system shall throw an error if the heat structure boundary connected to a HeatTransferFromHeatStructure3D1Phase component doesn't exist."
[]

[err:differently_aligned_channels]
  type = 'RunException'
  input = 'phy.conservation.i'
  cli_args = "Components/fch1/orientation='0 0 1' Components/fch2/orientation='0 1 0'"
  expect_err = "ht: Flow channel 'fch2' has a different axis alignment \(1\). Make sure all flow channels are aligned with the same axis."

  requirement = "The system shall throw an error if the flow channels connected to a HeatTransferFromHeatStructure3D1Phase component are not aligned with the same axis."
[]

[err:different_n_elems]
  type = 'RunException'
  input = 'phy.conservation.i'
  cli_args = "Components/fch1/orientation='0 0 1' Components/fch2/orientation='0 0 1' Components/fch2/n_elems=7"
  expect_err = "ht: Flow channel 'fch2' has 7 elements which is inconsistent with the rest of the flow channels. Make sure all flow channels have the same number of elements."

  requirement = "The system shall throw an error if the flow channels connected to a HeatTransferFromHeatStructure3D1Phase component don't have the same lnumber of elements."
[]

[err:different_lengths]
  type = 'RunException'
  input = 'phy.conservation.i'
  cli_args = "Components/fch1/orientation='0 0 1' Components/fch2/orientation='0 0 1' Components/fch2/length=2"
  expect_err = "ht: Flow channel 'fch2' has length equal to 2 which is inconsistent with the rest of the flow channels. Make sure all flow channels have the length."
  requirement = "The system shall throw an error if the flow channels connected to a HeatTransferFromHeatStructure3D1Phase component don't have the same length."
[]

[jac]
  type = 'PetscJacobianTester'
  input = 'jac.1phase.i'
  ratio_tol = 1.3e-2
  difference_tol = 35000
  cli_args = "Debug/check_jacobian=true"
  max_parallel = 1
  requirement = 'The system shall correctly compute Jacobians for HeatTransferFromHeatStructure3D1Phase.'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false
  mesh_mode = 'replicated'

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in heat transfer from 3D heat structures'
  design = 'HeatTransferFromHeatStructure3D1Phase.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false
  mesh_mode = 'replicated'

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in heat transfer from 3D heat structures'
  design = 'HeatTransferFromHeatStructure3D1Phase.md SolutionIC.md'
  issues = '#20553'
[]

[clg:head]
  type = 'CSVDiff'
  input = 'clg.head.i'
  csvdiff = 'clg.head_out.csv'
  group = 'pump1phase'
  requirement = 'The system shall allow for controlling the pump head'
  design = 'Pump1Phase.md'
  issues = '#684'
[]

[jacobian]
  type = 'PetscJacobianTester'
  input = 'jacobian.i'
  ratio_tol = 2e-8
  difference_tol = 20
  cli_args = 'Debug/check_jacobian=true'
  max_parallel = 1
  requirement = 'The system shall allow for controlling the pump head'
  design = 'Pump1Phase.md'
  issues = '#684'
[]

[phy:mass_energy_conservation]
  type = 'CSVDiff'
  input = 'shaft_motor_compressor.i'
  csvdiff = 'shaft_motor_compressor_out.csv'
  cli_args = "Outputs/out/type=CSV Outputs/out/execute_on=final Outputs/out/show='mass_conservation energy_conservation'"
  abs_zero = 1e-5
  rel_err = 0

  requirement = 'The system shall conserve mass and energy when using ShaftConnectedCompressor1Phase.'
[]

[phy:loop]
  type = 'Exodiff'
  input = 'shaft_motor_compressor.i'
  exodiff = 'shaft_motor_compressor_out.e'
  rel_err = 4e-3
  abs_zero = 1e-9
  cli_args = 'Outputs/exodus=true'

  requirement = 'The system shall be able to model a compressor with ShaftConnectedCompressor1Phase.'
[]

[runs_with_zero_shaft_speed]
  type = 'RunApp'
  input = 'shaft_motor_compressor.i'
  cli_args = "
      Components/shaft/initial_speed=0
      GlobalParams/initial_vel=0
      GlobalParams/initial_vel_x=0
      Executioner/num_steps=1"
  recover = False

  requirement = 'The system shall allow ShaftConnectedCompressor1Phase to run with a zero shaft speed.'
  issues = '#19863 #20008 #20009'
[]

[jac:test]
  type = 'PetscJacobianTester'
  input = 'jac.test.i'
  ratio_tol = 4e-3
  difference_tol = 8
  cli_args = 'Debug/check_jacobian=true'
  allow_test_objects = true
  max_threads = 1
  max_parallel = 1

  requirement = 'The system shall correctly compute Jacobians for ShaftConnectedCompressor1Phase.'
[]

[err:not_connected_to_shaft]
  type = 'RunException'
  input = 'shaft_motor_compressor.i'
  cli_args = "Components/shaft/connected_components=''"
  expect_err = "This component must be connected to a shaft"

  requirement = 'The system shall throw an error if ShaftConnectedCompressor1Phase is not connected to a shaft component.'
  issues = '#19863 #19998'
[]

[err:no_initial_speed]
  type = 'RunException'
  input = 'test.i'
  cli_args = 'Components/hs/initial_T=300'
  expect_err = "shaft: The `initial_speed` parameter is missing."

  requirement = 'The system shall throw an error if the initial shaft speed is not provided and the application is not restarting.'
[]

[err:not_connected_to_shaft]
  type = 'RunException'
  input = 'test.i'
  cli_args = "
      Components/hs/initial_T=300
      Components/shaft/initial_speed=0
      Components/shaft/connected_components=''"
  expect_err = "This component must be connected to a shaft"

  requirement = 'The system shall throw an error if ShaftConnectedMotor is not connected to a shaft component.'
  issues = '#19833 #19998'
[]

[restart_part1]
  type = 'CSVDiff'
  csvdiff = part1_out.csv
  input = 'test.i'
  cli_args = '
      Components/hs/initial_T=300
      Components/shaft/initial_speed=0
      Outputs/checkpoint=true
      Outputs/file_base=part1_out'
  recover = false

  requirement = 'The system shall be able to model a motor connected to a shaft.'
[]

[restart_part2]
  type = 'CSVDiff'
  csvdiff = part2_out.csv
  input = 'test.i'
  cli_args = '
      Outputs/file_base=part2_out
      Problem/restart_file_base=part1_out_cp/LATEST
      Executioner/num_steps=2
      Outputs/file_base=part2_out'
  prereq = 'restart_part1'
  recover = false

  requirement = 'The system shall be able to execute a restart a simulation involving a shaft-connected motor.'
[]

[clg_test_torque]
  type = 'CSVDiff'
  csvdiff = clg.test.torque.csv
  input = 'clg.test.i'
  cli_args = 'ControlLogic/motor_ctrl/parameter=torque ControlLogic/motor_ctrl/function=torque_fn Postprocessors/test/parameter=torque Outputs/file_base=clg.test.torque'

  requirement = 'The system shall allow the torque of a shaft-connected motor to be controlled.'
[]

[clg_test_inertia]
  type = 'CSVDiff'
  csvdiff = clg.test.inertia.csv
  input = 'clg.test.i'
  cli_args = 'ControlLogic/motor_ctrl/parameter=inertia ControlLogic/motor_ctrl/function=inertia_fn Postprocessors/test/parameter=inertia Outputs/file_base=clg.test.inertia'

  requirement = 'The system shall allow the inertia of a shaft-connected motor to be controlled.'
[]

[phy:mass_energy_conservation]
  type = 'CSVDiff'
  input = 'shaft_motor_pump.i'
  csvdiff = 'shaft_motor_pump_out.csv'
  cli_args = "Outputs/out/type=CSV Outputs/out/execute_on=final Outputs/out/show='mass_conservation energy_conservation'"
  abs_zero = 1e-8
  rel_err = 0

  requirement = 'The system shall conserve mass and energy when using ShaftConnectedPump1Phase.'
[]

[phy:loop]
  type = 'Exodiff'
  input = 'shaft_motor_pump.i'
  exodiff = 'shaft_motor_pump_out.e'
  rel_err = 1e-3
  abs_zero = 1e-9
  cli_args = 'Outputs/exodus=true'

  requirement = 'The system shall be able to model a pump with ShaftConnectedPump1Phase.'
[]

[phy:coastdown]
  type = 'Exodiff'
  input = 'pump_coastdown.i'
  exodiff = 'pump_coastdown_out.e'
  rel_err = 1e-3
  abs_zero = 1e-9
  heavy = true

  requirement = 'The system shall be able to model a pump coastdown with ShaftConnectedPump1Phase.'
[]

[jacobian]
  type = 'PetscJacobianTester'
  input = 'jacobian.i'
  ratio_tol = 1e-7
  difference_tol = 50
  cli_args = 'Debug/check_jacobian=true'
  allow_test_objects = true
  max_threads = 1
  max_parallel = 1

  requirement = 'The system shall correctly compute Jacobians for ShaftConnectedPump1Phase.'
[]

[err:not_connected_to_shaft]
  type = 'RunException'
  input = 'shaft_motor_pump.i'
  cli_args = "Components/shaft/connected_components=''"
  expect_err = "This component must be connected to a shaft"

  requirement = 'The system shall throw an error if ShaftConnectedPump1Phase is not connected to a shaft component.'
  issues = '#19833 #19998'
[]

[phy:mass_energy_conservation]
  type = 'CSVDiff'
  input = 'shaft_motor_turbine.i'
  csvdiff = 'shaft_motor_turbine_out.csv'
  cli_args = "Outputs/out/type=CSV Outputs/out/execute_on=final Outputs/out/show='mass_conservation energy_conservation'"
  abs_zero = 1e-5
  rel_err = 0

  requirement = 'The system shall conserve mass and energy when using ShaftConnectedTurbine1Phase.'
[]

[phy:loop]
  type = 'Exodiff'
  input = 'shaft_motor_turbine.i'
  exodiff = 'shaft_motor_turbine_out.e'
  rel_err = 1.1
  abs_zero = 1e-7
  cli_args = 'Outputs/exodus=true'

  requirement = 'The system shall be able to model a turbine with ShaftConnectedTurbine1Phase.'
[]

[phy:startup]
  type = 'Exodiff'
  input = 'turbine_startup.i'
  exodiff = 'turbine_startup_out.e'
  rel_err = 4e-3
  abs_zero = 1e-9
  cli_args = 'Outputs/exodus=true Outputs/velocity_as_vector=false'
  heavy = true

  requirement = 'The system shall be able to model a turbine startup with ShaftConnectedTurbine1Phase.'
[]

[jac:test]
  type = 'PetscJacobianTester'
  input = 'jac.test.i'
  ratio_tol = 4e-2
  difference_tol = 7e3
  cli_args = 'Debug/check_jacobian=true'
  allow_test_objects = true
  max_threads = 1
  max_parallel = 1

  requirement = 'The system shall correctly compute Jacobians for ShaftConnectedTurbine1Phase.'
[]

[err:not_connected_to_shaft]
  type = 'RunException'
  input = 'shaft_motor_turbine.i'
  cli_args = "Components/shaft/connected_components=''"
  expect_err = "This component must be connected to a shaft"

  requirement = 'The system shall throw an error if ShaftConnectedTurbine1Phase is not connected to a shaft component.'
  issues = '#19876 #19998'
[]

[not_implemented]
  type = 'RunException'
  input = 'err.i'
  expect_err = "This component is deprecated."
  design = 'SupersonicInlet.md'
  requirement = 'The system shall report an error if the SupersonicInlet component is used.'
  issues = '#20383'
[]

[equal_area_x_direction]
  type = 'CSVDiff'
  input = 'equal_area_with_junction.i'
  cli_args = 'Components/junction/volume=1.0'
  # The gold file was actually produced with equal_area_no_junction.i and
  # then renamed to equal_area_with_junction_out.csv.
  csvdiff = 'equal_area_with_junction_out.csv'
  detail = '2 pipes of equal area in the x direction,'
[]

[equal_area_not_x_direction]
  type = 'CSVDiff'
  input = 'equal_area_with_junction.i'
  cli_args = '
        Components/junction/volume=1.0
        Functions/T0/axis=y
        Components/pipe1/orientation="0 1 0"
        Components/junction/position="0 1.02 0"
        Components/junction/initial_vel_x=0
        Components/junction/initial_vel_y=1
        Components/pipe2/position="0 1.04 0"
        Components/pipe2/orientation="0 1 0"
        Postprocessors/junction_rhou/variable=junction:rhovV'
  csvdiff = 'equal_area_with_junction_out.csv'
  prereq = 'test_phy/equal_area_x_direction'
  detail = '2 pipes of equal area not in the x direction,'
[]

[phy:unequal_area]
  type = 'CSVDiff'
  input = 'phy.unequal_area.i'
  csvdiff = 'phy.unequal_area_out.csv'
  abs_zero = 5e-5
  rel_err = 0
  recover = false
  detail = '2 pipes of unequal area,'
[]

[phy:deadend]
  type = 'CSVDiff'
  input = 'phy.deadend.i'
  csvdiff = 'phy.deadend_out.csv'
  abs_zero = 1e-5
  rel_err = 0
  detail = '3 pipes, 1 of which going to a dead-end,'
[]

[phy:shower_syntax]
  type = 'RunApp'
  input = 'phy.shower.i'
  check_input = true
  recover = false
  detail = '2 pipes with different temperatures mixing together into a third pipe with correct syntax,'
[]

[phy:shower]
  type = 'CSVDiff'
  input = 'phy.shower.i'
  csvdiff = 'phy.shower_out.csv'
  abs_zero = 1e-9
  rel_err = 0
  heavy = true
  detail = '2 pipes with different temperatures mixing together into a third pipe with correct results,'
[]

[calorically_imperfect_gas]
  type = 'CSVDiff'
  input = 'junction_with_calorifically_imperfect_gas.i'
  csvdiff = 'junction_with_calorifically_imperfect_gas_out.csv'
  abs_zero = 5e-5
  rel_err = 1e-4
  detail = 'pipes with the calorically imperfect gas fluid properties,'
[]

[phy.form_loss]
  type = 'CSVDiff'
  input = 'phy.form_loss.i'
  csvdiff = 'phy.form_loss.csv'
  cli_args = 'Components/junction/K=1 Outputs/file_base=phy.form_loss'
  rel_err = 2e-4
  recover = false
  detail = 'by specifying a constant loss coefficient and using the area of the first connected pipe, '
[]

[phy.form_loss_Aref]
  type = 'CSVDiff'
  input = 'phy.form_loss.i'
  csvdiff = 'phy.form_loss_Aref.csv'
  cli_args = 'Components/junction/K=1 Components/junction/A_ref=0.2 Outputs/file_base=phy.form_loss_Aref'
  rel_err = 2e-4
  recover = false
  detail = 'by specifying a constant loss coefficient and a reference flow area, and'
[]

[phy.form_loss_crtl]
  type = 'CSVDiff'
  input = 'phy.form_loss.i'
  csvdiff = 'phy.form_loss_crtl.csv'
  cli_args = 'Components/junction/K=1 ControlLogic/active=K_crtl Outputs/execute_on=TIMESTEP_END Outputs/file_base=phy.form_loss_crtl'
  rel_err = 2e-4
  recover = false
  detail = 'by specifying the loss coefficient through the control system.'
[]

[conservation_1phase]
  type = 'CSVDiff'
  input = 't_junction_1phase.i'
  csvdiff = 't_junction_1phase_out.csv'
  abs_zero = 1e-10
  requirement = 'The system shall conserve mass and energy when a VolumJunction1Phase component is used'
[]

[err.missing_ics]
  type = 'RunException'
  input = 'err.missing_ics.i'
  expect_err = "junction: The following initial condition parameters are missing: initial_p initial_T initial_vel_x initial_vel_y initial_vel_z"
  requirement = 'The system shall throw an error if initial conditions for the VolumeJunction1Phase component are missing.'
[]

[err.missing_K]
  type = 'RunException'
  input = 'phy.form_loss.i'
  cli_args = '--error Components/junction/A_ref=1'
  expect_err = "junction: Parameter 'A_ref' is specified, but 'K' is not specified, so the junction will behave as if there were no form loss."
  requirement = 'The system shall throw an error if the parameter "A_ref" is specifed and the paramter "K" is not specified.'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  requirement = 'The system shall provide a control that evaluates a parsed function'
  design = 'ParsedFunctionControl.md'
  issues = '#93'
[]

[no_latch]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'no_latch.csv'
  cli_args = 'Outputs/file_base=no_latch'
  group = 'controls'
  requirement = 'The system shall provide a unit trip component that report true if the trip condition was met and false otherwise.'
  design = 'UnitTripControl.md'
  issues = '#619'
[]

[latch]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'latch.csv'
  cli_args = 'ControlLogic/trip_ctrl/latch=true Outputs/file_base=latch'
  requirement = 'The system shall provide a unit trip component that stays in tripped state after the trip happened.'
  design = 'UnitTripControl.md'
  issues = '#619'
  group = 'controls'
[]

[err:not_boolean]
  type = 'RunException'
  input = 'test.i'
  cli_args = 'ControlLogic/trip_ctrl/condition=a'
  expect_err = 'The user-provided condition expression did not return a boolean value \(0 or 1\)'
  requirement = 'The system shall report an error when an unit trip condition does not evaluate as boolean value.'
  design = 'UnitTripControl.md'
  issues = '#619'
  group = 'controls'
[]

[compute_axial_coordinate]
  type = Exodiff
  input = 'discrete_line_segment_interface.i'
  exodiff = 'axial_coord.e'
  allow_test_objects = true
  # 0 time steps taken
  recover = false

  requirement = 'The system shall provide an interface to compute an axial coordinate from an arbitrary spatial point.'
[]

[compute_radial_coordinate]
  type = Exodiff
  input = 'discrete_line_segment_interface.i'
  cli_args = "
      AuxKernels/testvar_aux/test_type=radial_coord
      Outputs/file_base=radial_coord"
  exodiff = 'radial_coord.e'
  allow_test_objects = true
  # 0 time steps taken
  recover = false

  requirement = 'The system shall provide an interface to compute a radial coordinate from an arbitrary spatial point.'
[]

[get_axial_section_index]
  type = Exodiff
  input = 'discrete_line_segment_interface.i'
  cli_args = "
      AuxKernels/testvar_aux/test_type=axial_section_index
      Outputs/file_base=axial_section_index"
  exodiff = 'axial_section_index.e'
  allow_test_objects = true
  # 0 time steps taken
  recover = false

  requirement = 'The system shall provide an interface to get the axial section index for an arbitrary spatial point.'
[]

[get_axial_element_index]
  type = Exodiff
  input = 'discrete_line_segment_interface.i'
  cli_args = "
      AuxVariables/testvar/family=MONOMIAL
      AuxVariables/testvar/order=CONSTANT
      AuxKernels/testvar_aux/test_type=axial_element_index
      Executioner/Quadrature/type=MONOMIAL
      Executioner/Quadrature/order=CONSTANT
      Outputs/file_base=axial_element_index"
  exodiff = 'axial_element_index.e'
  allow_test_objects = true
  # 0 time steps taken
  recover = false

  requirement = 'The system shall provide an interface to get the axial element index for an arbitrary spatial point.'
[]

[compute_axial_coordinate_invalid_axial_coord]
  type = RunException
  input = 'discrete_line_segment_interface.i'
  cli_args = "AuxKernels/testvar_aux/length=5.0"
  expect_err = "testvar_aux: The point \(x,y,z\)=\(.*\) has an invalid axial coordinate \(.*\). Valid axial coordinates are in the range \(0,5\)."
  allow_test_objects = true

  detail = 'if an invalid axial coordinate is provided.'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  requirement = "The system shall be able to compute the convective heat transfer coefficient using the Wolf-McCarthy correlation."
  issues = '#23709'
  design = 'ADWallHeatTransferCoefficientWolfMcCarthyMaterial.md'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in flow channels'
  design = 'FlowChannel1Phase.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in flow channels'
  design = 'FlowChannel1Phase.md SolutionIC.md'
  issues = '#20553'
[]

[explicit]
  type = 'Exodiff'
  input = 'lax_shock_tube.i'
  exodiff = 'lax_shock_tube.e'
  cli_args = 'Executioner/num_steps=5 Executioner/TimeIntegrator/order=2'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'

  detail = 'using an explicit temporal discretization, and'
[]

[implicit]
  type = 'Exodiff'
  input = 'lax_shock_tube.i'
  exodiff = 'lax_shock_tube_implicit.e'
  cli_args = '
        Preconditioning/pc/type=SMP
        Preconditioning/pc/full=true
        Executioner/TimeIntegrator/type=BDF2
        Executioner/solve_type=NEWTON
        Executioner/dt=0.01
        Executioner/num_steps=5
        Outputs/file_base=lax_shock_tube_implicit'
  abs_zero = 1e-8
  rel_err = 1e-5
  max_parallel = 1
  method = opt

  detail = 'using an implicit temporal discretization.'
[]

[test]
  type = XMLDiff
  input = 'natural_circulation.i'
  xmldiff = 'natural_circulation_out.xml'
  rel_err = 1e-4 # parallel/threading tests require this
  heavy = true
  recover = false # run to steady-state; recover does not know when half transient ends
  requirement = 'The system shall simulate a natural circulation loop using flow channels and junctions.'
[]

[without_junction]
  type = CSVDiff
  input = 'pressure_drop.i'
  csvdiff = 'pressure_drop_csv_pressure_vpp_FINAL.csv'
  rel_err = 1e-6
  cli_args = 'GlobalParams/rdg_slope_reconstruction=minmod Executioner/num_steps=5'

  detail = 'without a junction, and'
[]

[with_junction]
  type = CSVDiff
  input = 'pressure_drop_with_junction.i'
  csvdiff = 'pressure_drop_with_junction_csv_pressure_vpp_FINAL.csv'
  rel_err = 1e-6
  cli_args = 'GlobalParams/rdg_slope_reconstruction=minmod Executioner/num_steps=5'

  detail = 'with a junction.'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in flow channels'
  design = 'FlowChannel1Phase.md SolutionIC.md'
  issues = '#20553'
[]

[non_existent_block]
  type = 'RunException'
  input = 'err.non_existent_block.i'
  prereq = steady_state
  recover = false
  expect_err = "Block 'asdf' does not exist in the file '.+'"

  requirement = 'The system shall report an error when a block is non found in the restart ExodusII file'
  design = 'SolutionIC.md'
  issues = '#20526'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in heat structures'
  design = 'heat_structure.md SolutionIC.md'
  issues = '#20465'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in 3D heat structures'
  design = 'HeatStructureFromFile3D.md SolutionIC.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in volume junctions'
  design = 'HeatTransferFromHeatStructure1Phase.md SolutionIC.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false
  mesh_mode = 'replicated'

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in heat transfer from 3D heat structures'
  design = 'HeatTransferFromHeatStructure3D1Phase.md SolutionIC.md'
  issues = '#20553'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in heat structures'
  design = 'heat_structure.md'
  issues = '#20465'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in heat structures'
  design = 'heat_structure.md SolutionIC.md'
  issues = '#20465'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in 3D heat structures'
  design = 'HeatStructureFromFile3D.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in 3D heat structures'
  design = 'HeatStructureFromFile3D.md SolutionIC.md'
  issues = '#20553'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in shaft'
  design = 'Shaft.md'
  issues = '#20553'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in shaft'
  design = 'Shaft.md ScalarSolutionIC.md'
  issues = '#20553'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in shaft'
  design = 'Shaft.md ScalarSolutionIC.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in volume junctions'
  design = 'flow_junction.md ScalarSolutionIC.md'
  issues = '#20553'
[]

[steady_state]
  type = 'Exodiff'
  input = 'steady_state.i'
  exodiff = 'steady_state_out.e'
  recover = false

  requirement = 'The system shall be able to produce an exodus file for setting initial conditions in volume junctions'
  design = 'flow_junction.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  prereq = steady_state
  recover = false

  requirement = 'The system shall be able to use an exodus file for setting initial conditions in volume junctions'
  design = 'flow_junction.md ScalarSolutionIC.md'
  issues = '#20553'
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  cli_args = '-r 1'
  recover = false

  issues = "#226"
  design = "THMCreateMeshAction.md"
  requirement = "The system shall uniform refine mesh when specifid on the command line"
[]

[disable]
  type = RunApp
  input = 'disable_scalars_in_console.i'
  absent_out = 'Scalar Variable Values:'
  requirement = 'The system shall be able to disable the output of scalar variables to the console.'
[]

[allow]
  type = RunApp
  input = 'disable_scalars_in_console.i'
  cli_args = 'Outputs/disable_scalars_in_console=false'
  expect_out = 'Scalar Variable Values:'
  requirement = 'The system shall be able to allow the output of scalar variables to the console.'
[]

[test]
  type = CSVDiff
  input = 'element_integral_material_property_rz.i'
  csvdiff = 'element_integral_material_property_rz.csv'
  rel_err = 0
  abs_zero = 1e-12

  requirement = 'This system shall compute an RZ integral of a material property.'
  design = 'ADElementIntegralMaterialPropertyRZ.md'
  issues = '#23420'
[]

[err:rz_domain]
  type = 'RunException'
  input = 'err.rz_domain.i'
  expect_err = 'el_int: This is a THM-specific object can be applied only on subdomains with Cartesian coordinate system\. If your domain has RZ cooridnate system associated with it, you need to use the object without the RZ suffix to obtain the desired result\.'
  recover = false
  issues = '#215'
  design = 'RZSymmetry.md'
  requirement = 'The system should report an error when users set subdomain-restricted RZ-symmtrical THM-specific objects on RZ-subdomains.'
[]

[err:rz_domain]
  type = 'RunException'
  input = 'err.rz_domain.i'
  expect_err = 'el_int: This is a THM-specific object can be applied only on subdomains with Cartesian coordinate system\. If your domain has RZ cooridnate system associated with it, you need to use the object without the RZ suffix to obtain the desired result\.'
  issues = '#215'
  design = 'RZSymmetry.md'
  requirement = 'The system should error out when users set boundary-restricted RZ-symmtrical THM-specific objects on RZ-subdomains.'
[]

[test]
  type = CSVDiff
  input = 'heat_rate_conduction_rz.i'
  csvdiff = 'heat_rate_conduction_rz.csv'
  rel_err = 0

  requirement = 'The system shall compute the heat conduction rate across an RZ boundary.'
[]

[test]
  type = 'CSVDiff'
  input = 'heat_rate_heat_flux.i'
  csvdiff = 'heat_rate_heat_flux_out.csv'
  recover = false # no time steps

  requirement = 'The system shall compute the heat rate for a user-provided heat flux function.'
  issues = '#24261'
[]

[test]
  type = 'CSVDiff'
  input = 'heat_rate_heat_flux_rz.i'
  csvdiff = 'heat_rate_heat_flux_rz_out.csv'
  recover = false # no time steps

  requirement = 'The system shall compute the heat rate for a user-provided heat flux function for a cylindrical boundary.'
  issues = '#24261'
[]

[test]
  type = 'CSVDiff'
  input = 'shaft_connected_component_postprocessor.i'
  csvdiff = 'shaft_connected_component_postprocessor_out.csv'
  recover = false

  requirement = 'The system shall provide a post-processor to retrieve the torque and moment of inertia from a shaft-connected component.'
[]

[Isp_1ph]
  type = 'CSVDiff'
  input = 'Isp_1ph.i'
  csvdiff = 'Isp_1ph_out.csv'

  issues = '#189'
  design = 'SpecificImpulse1Phase.md'
  requirement = 'The system shall compute specific impulse from conditions on a boundary'
[]

[open]
  requirement = 'The system shall be able to model an open Brayton cycle'
  type = CSVDiff
  input = 'open_brayton_cycle.i'
  csvdiff = 'open_brayton_cycle.csv'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  cli_args = 'motor_ramp_up_duration=1 motor_ramp_down_duration=1 post_motor_time=1'
[]

[closed]
  requirement = 'The system shall be able to model a closed Brayton cycle'
  type = CSVDiff
  input = 'closed_brayton_cycle.i'
  csvdiff = 'closed_brayton_cycle.csv'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  cli_args = 'motor_ramp_up_duration=1 motor_ramp_down_duration=1 post_motor_time=1'
[]

[recuperated]
  type = RunApp
  input = 'recuperated_brayton_cycle.i'
  check_input = true
  max_parallel = 1
  recover = false
  requirement = 'The system shall be able to model an open recuperated Brayton cycle'
[]

[test]
  type = XMLDiff
  input = 'natural_circulation.i'
  xmldiff = 'natural_circulation_out.xml'
  rel_err = 1e-4 # parallel/threading tests require this
  heavy = true
  recover = false # run to steady-state; recover does not know when half transient ends
  requirement = 'The system shall simulate a natural circulation loop using flow channels and junctions.'
[]

[layered_area_change]
  type = Exodiff
  input = layered_flow_area_2D.i
  exodiff = layered_flow_area_2D_out.e
  design = 'userobjects/LayeredFlowAreaChange.md'
  requirement = "The system shall allow computing changes in channel flow areas from deformation."
[]

[non_ad]
  type = 'CSVDiff'
  input = 'sampler_1d_real.i'
  csvdiff = 'out_test_property_vpp_0001.csv'

  requirement = 'The system shall provide a vector post-processor to sample regular material properties in one or more blocks.'
[]

[ad]
  type = 'CSVDiff'
  input = 'ad_sampler_1d_real.i'
  csvdiff = 'ad_sampler_1d_real_out_test_vpp_0000.csv'
  recover = False

  requirement = 'The system shall provide a vector post-processor to sample AD material properties in one or more blocks.'
[]

[error_on_nonexistent_matprop]
  type = RunException
  input = 'sampler_1d_real.i'
  cli_args = "VectorPostprocessors/test_property_vpp/property=nonexistent_matprop"
  expect_err = "The material property 'nonexistent_matprop' does not exist."

  requirement = 'The system shall report an error if a non-existent material property is requested for the block material property sampler vector post-processor.'
[]