Thermal Hydraulics Requirements Traceability Matrix

[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 = 'CSVDiff'
  input = '1phase.i'
  csvdiff = '1phase_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = '1phase.i'
  csvdiff = '1phase_out.csv'
[]

[1phase]
  type = 'CSVDiff'
  input = '1phase.i'
  csvdiff = '1phase_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = '1phase.i'
  csvdiff = '1phase_out.csv'
[]

[sum]
  type = 'CSVDiff'
  input = 'sum.i'
  csvdiff = 'sum_out.csv'
[]

[weighted_average]
  type = 'CSVDiff'
  input = 'weighted_average.i'
  csvdiff = 'weighted_average_out.csv'
[]

[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'
[]

[warning]
  type = 'RunApp'
  input = 'test.i'
  cli_args = "--color off"
  expect_out = "1 warning:\n\s+- warn: This is a warning"
  method = 'OPT'
  allow_test_objects = true
  allow_warnings = true
[]

[loop_identification]
  type = 'RunApp'
  input = 'loop_identification.i'
  cli_args = 'Debug/print_component_loops=true'
  expect_out = "Loop 1:
\s+corechannel
\s+hx:primary
\s+j1
\s+j2
\s+j3
\s+j4
\s+pipe1
\s+pipe2

\s+Loop 2:
\s+hx:secondary
\s+inlet
\s+outlet"
[]

[phy:]
  type = 'Exodiff'
  input = 'phy.test.i'
  exodiff = 'phy.test_out.e'
[]

[enumeration_option]
  type = RunException
  input = 'phy.test.i'
  cli_args = "GlobalParams/closures=none"
  allow_deprecated = False
  expect_err = 'The closures system now uses objects created in the input file'
[]

[err:missing_f_1phase]
  type = 'RunException'
  input = 'err.missing_f_1phase.i'
  expect_err = "pipe: When using simple closures, the parameter 'f' must be provided."
[]

[enumeration_option]
  type = RunException
  input = 'err.missing_f_1phase.i'
  cli_args = "GlobalParams/closures=simple Components/pipe/f=0"
  allow_deprecated = False
  expect_err = 'The closures system now uses objects created in the input file'
[]

[err:nonexisting_component]
  type = 'RunException'
  input = 'err.nonexisting_component.i'
  expect_err = "cmp: The component 'pipe' does not exist"
  allow_test_objects = true
[]

[err:setup_status]
  type = 'RunException'
  input = 'err.setup_status.i'
  expect_err = "pipe: The component setup status \(2\) is less than the required status \(3\)"
  allow_test_objects = true
[]

[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.'
[]

[phy:position]
  type = 'Exodiff'
  input = 'phy.position.i'
  exodiff = 'phy.position_out.e'
  custom_cmp = 'phy.position.exodiff'
[]

[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.'
[]

[phy:f_fn_3eqn]
  type = 'CSVDiff'
  input = 'phy.f_fn.3eqn.i'
  csvdiff = 'phy.f_fn.3eqn_out.csv'
  rel_err = 1e-10
[]

[phy:par_fn]
  type = 'CSVDiff'
  input = 'phy.par_fn.i'
  csvdiff = 'phy.par_fn_out.csv'
  rel_err = 1e-10
[]

[phy:sub_discretization]
  type = 'Exodiff'
  input = 'phy.sub_discretization.i'
  exodiff = 'phy.sub_discretization_out.e'
  custom_cmp = 'phy.sub_discretization.exodiff'
  # matching geometry at t=0
  recover = false
[]

[err:wrong_fp]
  type = 'RunException'
  input = 'err.wrong_fp.i'
  expect_err = "pipe: Supplied fluid properties must be for 1-phase fluids."
[]

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

[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.'
[]

[err:connecting_to_non_existent_component]
  type = 'RunException'
  input = 'err.connecting_to_non_existent_component.i'
  expect_err = "inlet_1p: Trying to connect to a component 'pipe', but there is no such component in the simulation. Please check your spelling."
[]

[err:old_connection_format]
  type = 'RunException'
  input = 'err.connection_format.i'
  expect_err = "You are using the old connection format 'comp_name\(end\)'\. Please update your input file to the new one 'comp_name:end'\."
  cli_args = "Components/left_wall/input='pipe(in)'"
[]

[err:incorrect_connection_format]
  type = 'RunException'
  input = 'err.connection_format.i'
  expect_err = "Incorrect connection specified 'pipe_in'\. Valid connection format is 'component_name:end_type'\."
  cli_args = "Components/left_wall/input='pipe_in'"
[]

[phy:form_loss_1phase]
  type = 'Exodiff'
  input = 'phy.form_loss_1phase.parent.i'
  exodiff = 'phy.form_loss_1phase.parent_out_slave0.e'
  rel_err = 1e-10
  recover = false
[]

[phy:form_loss_1phase]
  type = 'CSVDiff'
  input = 'phy.form_loss_1phase.i'
  csvdiff = 'phy.form_loss_1phase_out.csv'
  rel_err = 1e-3
[]

[phy:mass_1phase]
  type = 'CSVDiff'
  input = 'phy.conservation_free_boundary_1phase.i'
  csvdiff = 'phy.conservation_free_boundary_1phase_out_mass.csv'
  cli_args = 'Outputs/out_mass/type=CSV Outputs/out_mass/execute_postprocessors_on=final Outputs/out_mass/show=mass_conservation'
  abs_zero = 1e-10
  rel_err = 1e-15
[]

[phy:momentum_1phase]
  type = 'CSVDiff'
  input = 'phy.conservation_free_boundary_1phase.i'
  csvdiff = 'phy.conservation_free_boundary_1phase_out_momentum.csv'
  cli_args = 'Outputs/out_momentum/type=CSV Outputs/out_momentum/execute_postprocessors_on=final Outputs/out_momentum/show=momentum_conservation'
  abs_zero = 1e-6
  rel_err = 1e-15
  prereq = 'phy:mass_1phase'
[]

[phy:energy_1phase]
  type = 'CSVDiff'
  input = 'phy.conservation_free_boundary_1phase.i'
  csvdiff = 'phy.conservation_free_boundary_1phase_out_energy.csv'
  cli_args = 'Outputs/out_energy/type=CSV Outputs/out_energy/execute_postprocessors_on=final Outputs/out_energy/show=energy_conservation'
  abs_zero = 1e-2
  rel_err = 1e-15
  prereq = 'phy:momentum_1phase'
[]

[1phase_heavy]
  type = Exodiff
  input = 'gate_valve_1phase.i'
  exodiff = 'gate_valve_1phase_out.e'
  rel_err = 5e-5
  abs_zero = 1e-8
  heavy = true
[]

[1phase_light]
  type = Exodiff
  input = 'gate_valve_1phase.i'
  cli_args = '
      Executioner/num_steps=5
      Outputs/file_base=gate_valve_1phase_light'
  exodiff = 'gate_valve_1phase_light.e'
  rel_err = 1e-5
  abs_zero = 1e-8
[]

[1phase_light_flipped]
  type = Exodiff
  input = 'gate_valve_1phase.i'
  cli_args = "
      Components/pipe3_valve/connections='pipe3A:out pipe3B:out'
      Components/pipe3B/position='2 0 0'
      Components/pipe3B/orientation='-1 0 0'
      Components/pipe3_outlet/input='pipe3B:in'
      Executioner/num_steps=5
      Outputs/file_base=gate_valve_1phase_light_flipped"
  exodiff = 'gate_valve_1phase_light_flipped.e'
  rel_err = 1e-5
  abs_zero = 1e-8
[]

[err:not_2_connections]
  type = 'RunException'
  input = 'gate_valve_1phase.i'
  cli_args = "Components/pipe2_valve/connections='pipe1:in pipe2A:out pipe2B:in'"
  expect_err = "The number of connections \(3\) must equal 2"
[]

[err:not_parallel]
  type = 'RunException'
  input = 'gate_valve_1phase.i'
  cli_args = "Components/pipe3B/orientation='1 1 1'"
  expect_err = "The connected channels must be parallel at the junction"
[]

[warn:slope_reconstruction_used]
  type = 'RunException'
  input = 'gate_valve_1phase.i'
  cli_args = 'GlobalParams/rdg_slope_reconstruction=minmod'
  expect_err = "GateValve1Phase cannot perform slope reconstruction"
[]

[err:2nd_order]
  type = 'RunException'
  input = 'err.2nd_order.i'
  expect_err = "Cannot use TRAP quadrature rule with 2nd order elements.  Use SIMPSON or GAUSS instead."
[]

[phy:power_density_shape]
  type = 'Exodiff'
  input = 'phy.cylinder_power_shape_aux_var.i'
  exodiff = 'phy.cylinder_power_shape_aux_var_out.e'
  rel_err = 1e-4
  abs_zero = 1e-9

  requirement = "The system shall provide the heat source shape when power density is supplied"
[]

[phy:conservation_from_file_3d]
  type = 'CSVDiff'
  input = 'phy.conservation_from_file_3d.i'
  csvdiff = 'phy.conservation_from_file_3d_out.csv'
  rel_err = 1e-4
  abs_zero = 1e-9
  # because output has execute_on = 'final'
  recover = false
[]

[err:no_hs]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hgen/hs=hgen Components/hgen/regions='fuel'"
  expect_err = "hgen: The component 'hgen' is not of type 'HeatStructureInterface'"
[]

[err:wrong_block_name]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hgen/hs=hs Components/hgen/regions='asdf'"
  expect_err = "hgen: Region 'asdf' does not exist in heat structure 'hs'."
[]

[phy:power_component_shape_cylindrical]
  type = 'Exodiff'
  input = 'phy.cylinder_power_shape_fn.i'
  exodiff = 'phy.cylinder_power_shape_fn_out.e'
  rel_err = 1e-4
  abs_zero = 1e-9

  requirement = "The system shall provide the heat source shape for a cylindrical heat structure when power component and power shape function are supplied"
[]

[phy:power_component_shape_plate]
  type = 'Exodiff'
  input = 'phy.plate.i'
  exodiff = 'phy.plate_out.e'
  rel_err = 1e-4
  abs_zero = 1e-9

  requirement = "The system shall provide the heat source shape for a plate heat structure when power component and power shape function are supplied"
[]

[phy:conservation_plate]
  type = CSVDiff
  input = 'phy.conservation.i'
  cli_args = "
      Components/heat_structure/type=HeatStructurePlate
      Components/heat_structure/depth=2
      Postprocessors/E_tot/type=ADHeatStructureEnergy
      Postprocessors/E_tot/plate_depth=2
      Outputs/file_base=phy.conservation_plate"
  csvdiff = 'phy.conservation_plate.csv'
  abs_zero = 4e-7
  # because output has execute_on = 'initial final'
  recover = false
[]

[phy:conservation_cylinder]
  type = CSVDiff
  input = 'phy.conservation.i'
  cli_args = "
      Components/heat_structure/type=HeatStructureCylindrical
      Components/heat_structure/inner_radius=1.0
      Postprocessors/E_tot/type=ADHeatStructureEnergyRZ
      Postprocessors/E_tot/axis_point='0 1 0'
      Postprocessors/E_tot/axis_dir='1 0 0'
      Outputs/file_base=phy.conservation_cylinder"
  csvdiff = 'phy.conservation_cylinder.csv'
  abs_zero = 1e-7
  # because output has execute_on = 'initial final'
  recover = false
[]

[phy:conservation_from_file_3d]
  type = CSVDiff
  input = 'phy.conservation_from_file_3d.i'
  csvdiff = 'phy.conservation_from_file_3d_out.csv'
  abs_zero = 1e-7
  # because output has execute_on = 'final'
  recover = false
[]

[err:no_hs]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hgen/hs=reactor Components/hgen/regions='fuel' Components/hgen/power=reactor"
  expect_err = "hgen: The component 'reactor' is not of type 'HeatStructureInterface'"
[]

[err:wrong_block_name]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hgen/hs=hs Components/hgen/regions='asdf' Components/hgen/power=reactor"
  expect_err = "hgen: Region 'asdf' does not exist in heat structure 'hs'."
[]

[err:wrong_power_comp]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hgen/hs=hs Components/hgen/regions='fuel' Components/hgen/power=hs"
  expect_err = "hgen: The component 'hs' is not of type 'TotalPowerBase'"
[]

[phy:conservation_1phase]
  type = CSVDiff
  input = 'phy.conservation.1phase.i'
  csvdiff = 'phy.conservation.1phase_out.csv'
  abs_zero = 1e-9
  # because output has execute_on = 'final'
  recover = false
[]

[err:nonexistent_flow_channel]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hs/flow_channel=asdf"
  expect_err = "hs: The component 'asdf' does not exist"
[]

[err:not_a_flow_channel]
  type = 'RunException'
  input = 'err.base.i'
  cli_args = "Components/hs/flow_channel=total_power"
  expect_err = "hs: The component 'total_power' is not of type 'FlowChannel1Phase'"
[]

[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.'
[]

[axial_regions]
  type = Exodiff
  input = 'axial_regions.i'
  cli_args = "Components/hs/axial_region_names='region1 region2 region3'"
  exodiff = 'axial_regions_out.e'
  rel_err = 1e-5
[]

[phy:variable_init_t]
  type = Exodiff
  input = 'phy.variable_init_t.i'
  exodiff = 'phy.variable_init_t_out.e'
[]

[phy:sub_discretization]
  type = Exodiff
  input = 'phy.sub_discretization.i'
  exodiff = 'phy.sub_discretization_out.e'
[]

[inner_radial_boundary]
  type = CSVDiff
  input = 'inner_radial_boundary.i'
  csvdiff = 'inner_radial_boundary_out_test_vpp_0000.csv'
  recover = false
[]

[interior_axial_boundaries]
  type = CSVDiff
  input = 'interior_axial_boundaries.i'
  csvdiff = 'interior_axial_boundaries_out.csv'
  recover = false
[]

[2nd_order]
  type = 'RunApp'
  input = '2nd_order.i'
  expect_out = "Solve Converged"
[]

[err:missing_axial_regions_parameter]
  type = RunException
  input = 'axial_regions.i'
  expect_err = "If there is more than 1 axial region, then the parameter 'axial_region_names' must be specified."
[]

[err:incorrect_size_axial_regions_parameter]
  type = RunException
  input = 'axial_regions.i'
  cli_args = "Components/hs/axial_region_names='region1 region2'"
  expect_err = "The number of entries in parameter 'axial_region_names' \(2\) must equal the number of entries of parameter 'length' \(3\)"
[]

[err:no_2nd_order_with_trap]
  type = 'RunException'
  input = 'err.no_2nd_order_with_trap.i'
  expect_err = "Cannot use TRAP quadrature rule with 2nd order elements.  Use SIMPSON or GAUSS instead."
[]

[err:no_T_ic]
  type = 'RunException'
  input = 'err.no_T_ic.i'
  expect_err = "hs: Missing initial condition for temperature."
[]

[err:incorrect_size_of_n_part_elems]
  type = 'RunException'
  input = 'err.no_T_ic.i'
  cli_args = 'Components/hs/n_part_elems=1'
  expect_err = "hs: The number of entries in parameter 'names' \(3\) must equal the number of entries of parameter 'n_part_elems' \(1\)"
[]

[err:incorrect_size_of_widths]
  type = 'RunException'
  input = 'err.no_T_ic.i'
  cli_args = 'Components/hs/widths=1'
  expect_err = "hs: The number of entries in parameter 'names' \(3\) must equal the number of entries of parameter 'widths' \(1\)"
[]

[err:incorrect_size_of_materials]
  type = 'RunException'
  input = 'err.no_T_ic.i'
  cli_args = 'Components/hs/materials=a'
  expect_err = "hs: The number of entries in parameter 'names' \(3\) must equal the number of entries of parameter 'materials' \(1\)"
[]

[phy:rz]
  type = Exodiff
  input = 'phy.rz.i'
  exodiff = 'phy.rz_out.e'
  custom_cmp = 'phy.rz.exodiff'
[]

[phy:rz_mats]
  type = Exodiff
  input = 'phy.rz_mats.i'
  exodiff = 'phy.rz_out.e'
  custom_cmp = 'phy.rz.exodiff'
  cli_args = 'Outputs/file_base=phy.rz_out'
  prereq = phy:rz
[]

[phy:standalone]
  type = Exodiff
  input = 'phy.standalone.i'
  exodiff = 'phy.standalone_out.e'
[]

[err:2d_mesh]
  type = RunException
  input = 'phy.standalone.i'
  cli_args = 'Components/blk/file=square.e'
  expect_err = "blk: File '(.+)square\.e' does not contain a 3D mesh\."
[]

[phy:test]
  type = Exodiff
  input = 'phy.test.i'
  exodiff = 'phy.test_out.e'
[]

[phy:test_mats]
  type = Exodiff
  input = 'phy.test_mats.i'
  exodiff = 'phy.test_out.e'
  cli_args = 'Outputs/file_base=phy.test_out'
  prereq = 'phy:test'
[]

[err:mixed_heat_modes]
  type = 'RunException'
  input = 'err.mixed_heat_modes.i'
  expect_err = "Heat sources for a flow channel must be all of temperature type or all of heat flux type"
[]

[phy:q_wall_transfer_3eqn]
  type = Exodiff
  input = 'phy.q_wall_transfer_3eqn.parent.i'
  exodiff = 'phy.q_wall_transfer_3eqn.parent_out.e phy.q_wall_transfer_3eqn.parent_out_thm0.e'
[]

[phy:T_wall_transfer_3eqn]
  type = Exodiff
  input = 'phy.T_wall_transfer_3eqn.parent.i'
  exodiff = 'phy.T_wall_transfer_3eqn.parent_out.e phy.T_wall_transfer_3eqn.parent_out_thm0_out.e'
  cli_args = 'thm:Components/hxconn/initial_T_wall=300'
[]

[phy:T_wall_transfer_ext_3eqn]
  type = Exodiff
  input = 'phy.T_wall_transfer_3eqn.parent.i'
  exodiff = 'phy.T_wall_transfer_ext_3eqn.e phy.T_wall_transfer_ext_3eqn_thm0_out.e'
  cli_args = "Outputs/file_base=phy.T_wall_transfer_ext_3eqn Transfers/T_to_slave/variable=T_ext thm:AuxVariables/T_ext thm:AuxVariables/T_ext/initial_condition=300 thm:Components/hxconn/T_ext=T_ext thm:Outputs/out/show='T_ext'"
[]

[phy:T_wall_transfer_elem_3eqn]
  type = Exodiff
  input = 'phy.T_wall_transfer_elem_3eqn.parent.i'
  exodiff = 'phy.T_wall_transfer_elem_3eqn.parent_out.e phy.T_wall_transfer_elem_3eqn.parent_out_thm0_out.e'
[]

[phy:q_wall_multiple_3eqn]
  type = 'CSVDiff'
  input = 'phy.q_wall_multiple_3eqn.i'
  csvdiff = 'phy.q_wall_multiple_3eqn_out.csv'
[]

[phy.energy_heatflux_ss_1phase]
  type = 'CSVDiff'
  input = 'phy.energy_heatflux_ss_1phase.i'
  csvdiff = 'phy.energy_heatflux_ss_1phase_out.csv'
  abs_zero = 1e-8
[]

[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'
[]

[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.'
[]

[phy.energy_walltemperature_ss_1phase]
  type = 'CSVDiff'
  input = 'phy.energy_walltemperature_ss_1phase.i'
  csvdiff = 'phy.energy_walltemperature_ss_1phase_out.csv'
  abs_zero = 1e-9
[]

[clg:Hw]
  type = 'CSVDiff'
  input = 'clg.Hw.i'
  csvdiff = 'clg.Hw_out.csv'
[]

[clg:T_wall]
  type = 'CSVDiff'
  input = 'clg.T_wall.i'
  csvdiff = 'clg.T_wall_out.csv'
[]

[err:no_phf]
  type = RunException
  input = 'err.no_phf.i'
  expect_err = "hx\: Invalid side \'START\'\. This side does not have unit perimeter\. You probably want to use \'INNER\' or \'OUTER\' side instead\."
[]

[energy_conservation_cylindrical]
  type = 'CSVDiff'
  input = 'cylindrical.i'
  csvdiff = 'cylindrical_out.csv'
  abs_zero = 1e-6
[]

[energy_conservation_plate]
  type = 'CSVDiff'
  input = 'plate.i'
  csvdiff = 'plate_out.csv'
  abs_zero = 1e-6
[]

[energy_conservation_from_file_3d]
  type = 'CSVDiff'
  input = 'from_file_3d.i'
  csvdiff = 'from_file_3d_out.csv'
  abs_zero = 1e-6
  # because output has execute_on = 'final'
  recover = false
[]

[energy_conservation_plate]
  type = Exodiff
  input = 'plate.parent.i'
  exodiff = 'plate.parent_out.e plate.parent_out_thm0.e'
[]

[test]
  type = Exodiff
  input = 'phy.parent.i'
  exodiff = 'phy.parent_out.e phy.parent_out_thm0.e'
[]

[energy_conservation_cylindrical]
  type = 'CSVDiff'
  input = 'cylindrical.i'
  csvdiff = 'cylindrical_out.csv'
  abs_zero = 1e-6
[]

[energy_conservation_plate]
  type = 'CSVDiff'
  input = 'plate.i'
  csvdiff = 'plate_out.csv'
  abs_zero = 1e-6
[]

[energy_conservation_from_file_3d]
  type = 'CSVDiff'
  input = 'from_file_3d.i'
  csvdiff = 'from_file_3d_out.csv'
  abs_zero = 1e-6
  # because output has execute_on = 'final'
  recover = false
[]

[energy_conservation_cylindrical]
  type = 'CSVDiff'
  input = 'cylindrical.i'
  csvdiff = 'cylindrical_out.csv'
  abs_zero = 1e-6
[]

[energy_conservation_plate]
  type = 'CSVDiff'
  input = 'plate.i'
  csvdiff = 'plate_out.csv'
  abs_zero = 2e-6
[]

[energy_conservation_from_file_3d]
  type = 'CSVDiff'
  input = 'from_file_3d.i'
  csvdiff = 'from_file_3d_out.csv'
  abs_zero = 1e-6
  # because output has execute_on = 'final'
  recover = false
[]

[err:no_bnd]
  type = RunException
  input = 'err.no_bnd.i'
  expect_err = "hs_boundary\: The heat structure side of the heat structure \'hs\' corresponding to the boundary name \'hs:inner\' has a zero perimeter\. This can be caused by applying the boundary on the axis of symmetry of a cylindrical heat structure\."
[]

[phy:densityvelocity_3eqn]
  type = 'Exodiff'
  input = 'phy.densityvelocity_3eqn.i'
  exodiff = 'phy.densityvelocity_3eqn_out.e'
  recover = false
[]

[clg:densityvelocity_3eqn]
  type = 'CSVDiff'
  input = 'clg.densityvelocity_3eqn.i'
  csvdiff = 'clg.densityvelocity_3eqn_out.csv'
  rel_err = 1e-5
[]

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

[phy:massflowrate_3eqn]
  type = 'Exodiff'
  input = 'phy.massflowrate_3eqn.i'
  exodiff = 'phy.massflowrate_3eqn.e'
  group = 'inlet_mass_flow_rate outlet pipe 1phase'
  max_time = 1000
[]

[phy:reversed_flow]
  type = 'Exodiff'
  input = 'phy.reversed_flow.i'
  exodiff = 'phy.reversed_flow.e'
[]

[clg:ctrl_m_dot_3eqn_rdg]
  type = 'CSVDiff'
  input = 'clg.ctrl_m_dot_3eqn_rdg.i'
  csvdiff = 'clg.ctrl_m_dot_3eqn_rdg_out.csv'
[]

[clg:ctrl_T_3eqn_rdg]
  type = 'CSVDiff'
  input = 'clg.ctrl_T_3eqn_rdg.i'
  csvdiff = 'clg.ctrl_T_3eqn_rdg_out.csv'
[]

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

[jacobian-water97]
  type = PetscJacobianTester
  run_sim = True
  ratio_tol = 1e-5
  difference_tol = 1e6
  cli_args = 'Executioner/num_steps=3 Components/pipe/n_elems=3 -snes_test_err 1e-10'
  input = jac.massflowrate_3eqn_water97.i
[]

[phy:h_rhou_3eqn]
  type = 'RunException'
  input = 'phy.h_rhou_3eqn.i'
  expect_err = 'This component does not work with rDG, yet\.'
[]

[phy:stagnation_p_T_steady_3eqn]
  type = 'Exodiff'
  input = 'phy.stagnation_p_T_steady_3eqn.i'
  exodiff = 'phy.stagnation_p_T_steady_3eqn_out.e'
[]

[phy:stagnation_p_T_transient_3eqn]
  type = 'Exodiff'
  input = 'phy.stagnation_p_T_transient_3eqn.i'
  exodiff = 'phy.stagnation_p_T_transient_3eqn_out.e'
[]

[phy:p0T0_3eqn]
  type = 'Exodiff'
  input = 'phy.p0T0_3eqn.i'
  exodiff = 'phy.p0T0_3eqn.e'
  max_time = 500
  group = 'inlet_p0_T0 outlet pipe 1phase'
[]

[clg:ctrl_p0_3eqn]
  type = 'CSVDiff'
  input = 'clg.ctrl_p0_3eqn.i'
  csvdiff = 'clg.ctrl_p0_3eqn_out.csv'
[]

[clg:ctrl_T0_3eqn]
  type = 'CSVDiff'
  input = 'clg.ctrl_T0_3eqn.i'
  csvdiff = 'clg.ctrl_T0_3eqn_out.csv'
[]

[phy:velocity_t_3eqn]
  type = 'Exodiff'
  input = 'phy.velocity_t_3eqn.i'
  exodiff = 'phy.velocity_t_3eqn.e'
  custom_cmp = 'rdg_1phase.exodiff'
  group = 'outlet pipe 1phase rdg'
  max_time = 600
[]

[phy:reversed_flow]
  type = 'Exodiff'
  input = 'phy.reversed_flow.i'
  exodiff = 'phy.reversed_flow.e'
[]

[clg:velocity_t_3eqn]
  type = 'CSVDiff'
  input = 'clg.velocity_t_3eqn.i'
  csvdiff = 'clg.velocity_t_3eqn_out.csv'
[]

[phy:no_junction_equivalency_1phase_normal]
  type = 'CSVDiff'
  input = 'junction_one_to_one_1phase.i'
  cli_args = "
      Postprocessors/rhouA_right/type=SideAverageValue
      Postprocessors/rhouA_right/variable=rhouA
      Postprocessors/rhouA_right/boundary=right_channel:in
      Postprocessors/rhouA_right/execute_on='initial timestep_end'"
  csvdiff = 'junction_one_to_one_1phase_out.csv'
[]

[phy:no_junction_equivalency_1phase_flipped]
  type = 'CSVDiff'
  input = 'junction_one_to_one_1phase.i'
  cli_args = "
      Components/right_channel/position='1 0 0'
      Components/right_channel/orientation='-1 0 0'
      Components/junction/connections='left_channel:out right_channel:out'
      Components/right_boundary/input=right_channel:in
      Postprocessors/rhoA_right/boundary=right_channel:out
      Postprocessors/rhouA_right_negative/type=SideAverageValue
      Postprocessors/rhouA_right_negative/variable=rhouA
      Postprocessors/rhouA_right_negative/boundary=right_channel:out
      Postprocessors/rhouA_right_negative/execute_on='initial timestep_end'
      Postprocessors/rhouA_right/type=ScalePostprocessor
      Postprocessors/rhouA_right/scaling_factor=-1
      Postprocessors/rhouA_right/value=rhouA_right_negative
      Postprocessors/rhouA_right/execute_on='initial timestep_end'
      Postprocessors/rhoEA_right/boundary=right_channel:out"
  csvdiff = 'junction_one_to_one_1phase_out.csv'
  prereq = 'phy:no_junction_equivalency_1phase_normal'
[]

[phy:unequal_area_1phase]
  type = 'Exodiff'
  input = 'constriction_1phase.i'
  exodiff = 'constriction_1phase_out.e'
[]

[phy:no_junction_1phase_check]
  type = 'RunApp'
  input = 'no_junction_1phase.i'
  check_input = true
  recover = false
[]

[equal_area_x_direction]
  type = 'CSVDiff'
  input = 'equal_area_with_junction.i'
  # 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'
[]

[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
[]

[phy:shower_syntax]
  type = 'RunApp'
  input = 'phy.shower.i'
  check_input = true
  recover = false
  max_threads = 1
[]

[phy:shower]
  type = 'CSVDiff'
  input = 'phy.shower.i'
  csvdiff = 'phy.shower_out.csv'
  abs_zero = 8e-8
  rel_err = 0
  heavy = true
[]

[calorically_imperfect_gas]
  type = 'CSVDiff'
  input = 'junction_with_calorifically_imperfect_gas.i'
  csvdiff = 'junction_with_calorifically_imperfect_gas_out.csv'
  abs_zero = 4e-5
  rel_err = 1e-5
[]

[conservation_1phase]
  type = 'CSVDiff'
  input = 'conservation.i'
  csvdiff = 'conservation_out.csv'
  abs_zero = 2e-8
[]

[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"
[]

[err.non_parallel_channels]
  type = 'RunException'
  input = 'equal_area_with_junction.i'
  cli_args = "--error Components/junction/volume=1.0 Components/pipe2/orientation='0 1 0'"
  expect_err = "junction: Connected flow channels are not parallel, use VolumeJunction1Phase component instead."
[]

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

[phy:solidwall_outlet_3eqn]
  type = 'Exodiff'
  input = 'phy.solidwall_outlet_3eqn.i'
  exodiff = 'phy.solidwall_outlet_3eqn.e'
  custom_cmp = 'phy.solidwall_outlet_3eqn.exodiff'
  group = 'solid_wall outlet pipe 1phase'
  # see #530 why
  max_threads = 1
[]

[clg:ctrl_p_3eqn]
  type = 'CSVDiff'
  input = 'clg.ctrl_p_3eqn.i'
  csvdiff = 'clg.ctrl_p_3eqn_out.csv'
[]

[jacobian]
  type = 'PetscJacobianTester'
  input = 'jacobian.i'
  ratio_tol = 1e-6
  difference_tol = 15
  cli_args = 'Debug/check_jacobian=true'
  max_parallel = 1
[]

[phy:pump_mass_energy_conservation]
  type = 'CSVDiff'
  input = 'pump_mass_energy_conservation.i'
  csvdiff = 'pump_mass_energy_conservation_out.csv'
  group = 'pump1phase'
[]

[phy:pump_mass_energy_conservation_orientation]
  type = 'CSVDiff'
  input = 'pump_mass_energy_conservation.i'
  cli_args = "Components/pipe1/orientation='1 1 1'
                Outputs/file_base=pump_mass_energy_conservation_out_orientation
                Outputs/out/show='E_change'"
  # csvdiff made with Components/pipe1/orientation='1 0 0'
  csvdiff = 'pump_mass_energy_conservation_out_orientation.csv'
  group = 'pump1phase'
[]

[phy:pressure_check]
  type = 'CSVDiff'
  input = 'pump_pressure_check.i'
  csvdiff = 'pump_pressure_check_out.csv'
  abs_zero = 5e-7
  group = 'pump1phase'
[]

[phy:pressure_check_orientation]
  type = 'CSVDiff'
  input = 'pump_pressure_check.i'
  cli_args = "Components/pipe1/orientation='20 0 0'
                Outputs/file_base=pump_pressure_check_out_orientation"
  # csvdiff made with Components/pipe1/orientation='1 0 0'
  csvdiff = 'pump_pressure_check_out_orientation.csv'
  abs_zero = 5e-7
  group = 'pump1phase'
[]

[phy:pipe_friction_pump_head_balance_x_direction]
  type = 'Exodiff'
  input = 'pipe_friction_pump_head_balance.i'
  exodiff = 'pipe_friction_pump_head_balance_out_x.e'
  group = 'pump1phase'
[]

[phy:pipe_friction_pump_head_balance_y_direction]
  type = 'Exodiff'
  input = 'pipe_friction_pump_head_balance.i'
  cli_args = "Components/pipe1/orientation='0 1 0'
                Components/pump/initial_vel_x=0
                Components/pump/initial_vel_y=1
                Outputs/file_base=pipe_friction_pump_head_balance_out_y"
  exodiff = 'pipe_friction_pump_head_balance_out_y.e'
  group = 'pump1phase'
[]

[phy:pump_loop]
  type = 'Exodiff'
  input = 'pump_loop.i'
  exodiff = 'pump_loop_out.e'
  rel_err = 8e-5
  abs_zero = 1e-9
  group = 'pump1phase'
[]

[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'
[]

[err:no_connected_components]
  type = 'RunException'
  input = 'err.no_connected_components.i'
  expect_err = "shaft: No components are connected to the shaft."
[]

[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'
[]

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

[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'
[]

[phy:conservation]
  type = 'CSVDiff'
  input = 'phy.conservation.i'
  csvdiff = 'phy.conservation_csv.csv'
  override_columns = 'energy_diff mass_diff momentum_diff'
  override_abs_zero = '1e-6 2e-7 1e-4'
  override_rel_err = '0 0 0'
[]

[phy:test]
  type = 'Exodiff'
  input = 'phy.test.i'
  exodiff = 'phy.test_out.e'
[]

[clg:test]
  type = 'CSVDiff'
  input = 'clg.test.i'
  csvdiff = 'clg.test_csv.csv'
[]

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

[phy:3eqn]
  type = 'Exodiff'
  input = 'phy.3eqn.i'
  exodiff = 'phy.3eqn_out.e'
[]

[jacobian]
  type = 'PetscJacobianTester'
  input = 'jacobian.i'
  ratio_tol = 2e-5
  difference_tol = 2100
  cli_args = 'Debug/check_jacobian=true'
  max_parallel = 1
[]

[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'
[]

[phy:constant_power]
  type = 'CSVDiff'
  input = 'phy.constant_power.i'
  csvdiff = 'phy.constant_power_out.csv'
[]

[clg:power]
  type = 'CSVDiff'
  input = 'clg.power.i'
  csvdiff = 'clg.power_out.csv'
[]

[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'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[err.negative_tau]
  type = 'RunException'
  input = 'test.i'
  cli_args = 'ControlLogic/delay_ctrl/tau=-1'
  expect_err = "Negative values of 'tau' are not allowed\."
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[non_existent_control_data]
  type = 'RunException'
  input = 'non_existent_control_data.i'
  expect_err = "Control data 'wrong' was requested, but was not declared by any active control object."
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[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'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-4
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  allow_test_objects = true
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  allow_test_objects = true
  recover = false
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[no_error]
  type = 'RunApp'
  input = 'terminate.i'
  expect_out = 'Threshold exceeded'

  group = 'controls'
[]

[throw_error]
  type = 'RunException'
  input = 'terminate.i'
  cli_args = '
      ControlLogic/terminate/throw_error=true
      ControlLogic/terminate/termination_message="Custom error message"'
  expect_err = 'Custom error message'

  group = 'controls'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  allow_test_objects = True
[]

[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'
[]

[function_values]
  type = 'CSVDiff'
  input = 'piecewise_function.i'
  csvdiff = 'piecewise_function_out.csv'
[]

[error_inconsistent_number_of_points]
  type = 'RunException'
  input = 'piecewise_function.i'
  cli_args = "Functions/function_piecewise/axis_coordinates='3'"
  expect_err = "The number of entries in the parameter 'axis_coordinates' must equal the number of entries in the parameter 'functions' minus one"
[]

[error_points_nonascending_order]
  type = 'RunException'
  input = 'piecewise_function.i'
  cli_args = "Functions/function_piecewise/axis_coordinates='5 3'"
  expect_err = "The entries in the parameter 'axis_coordinates' must be in ascending order"
[]

[error_self_reference]
  type = 'RunException'
  input = 'piecewise_function.i'
  cli_args = "Functions/function_piecewise/functions='function_left function_piecewise function_end'"
  expect_err = "This function cannot use its own name in the 'functions' parameter"
[]

[space_weighted]
  type = 'Exodiff'
  input = 'space.i'
  exodiff = 'space_weighted.e'
[]

[time_weighted]
  type = 'CSVDiff'
  input = 'time.i'
  csvdiff = 'time_weighted.csv'
[]

[space_cubic]
  type = 'Exodiff'
  input = 'space.i'
  cli_args = '
      Functions/transition_fn/type=CubicTransitionFunction
      Functions/transition_fn/function1_derivative_end_point=0
      Functions/transition_fn/function2_derivative_end_point=0
      Outputs/file_base=space_cubic'
  exodiff = 'space_cubic.e'
[]

[time_cubic]
  type = 'CSVDiff'
  input = 'time.i'
  cli_args = '
      Functions/transition_fn/type=CubicTransitionFunction
      Functions/transition_fn/function1_derivative_end_point=0
      Functions/transition_fn/function2_derivative_end_point=0
      Outputs/file_base=time_cubic'
  csvdiff = 'time_cubic.csv'
[]

[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 = 'PetscJacobianTester'
  input = 'convection_heat_transfer_bc.i'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'convection_heat_transfer_rz_bc.i'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'external_app_convection_heat_transfer_bc.i'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'external_app_convection_heat_transfer_rz_bc.i'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'radiation_heat_flux_bc.i'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'radiation_heat_flux_rz_bc.i'
  ratio_tol = 1e-6
  difference_tol = 1e-4
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[mass_free_constraint]
  type = 'PetscJacobianTester'
  input = 'mass_free_constraint.i'
  ratio_tol = 5e-6
  difference_tol = 5e-6
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_momentum_flux]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_momentum_flux.i'
  ratio_tol = 6e-8
  difference_tol = 0.09
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_momentum_area_gradient]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_momentum_area_gradient.i'
  ratio_tol = 0.0005
  difference_tol = 3e5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_momentum_friction]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_momentum_friction.i'
  ratio_tol = 4e-8
  difference_tol = 1
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_energy_flux]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_energy_flux.i'
  ratio_tol = 2e-8
  difference_tol = 2e-4
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_energy_wall_heating]
  type = 'PetscJacobianTester'
  input = 'one_d_energy_wall_heating.i'
  ratio_tol = 2e-7
  difference_tol = 8e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_momentum_gravity]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_momentum_gravity.i'
  ratio_tol = 5e-8
  difference_tol = 5e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_energy_gravity]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_energy_gravity.i'
  ratio_tol = 5e-8
  difference_tol = 5e-2
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[one_d_3eqn_energy_heat_source]
  type = 'PetscJacobianTester'
  input = 'one_d_3eqn_energy_heat_source.i'
  ratio_tol = 1e-20
  difference_tol = 1e-20
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[ad_solid_material]
  type = 'PetscJacobianTester'
  input = 'ad_solid_material.i'
  ratio_tol = 1e-7
  difference_tol = 1e-6
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[constant_material]
  type = 'PetscJacobianTester'
  input = 'constant_material.i'
  ratio_tol = 1e-14
  difference_tol = 1e-14
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[density_3eqn_material]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=rho Kernels/test_kernel/coupled_variables='rhoA' JacobianTest1Phase/snes_test_err=1e-8"
  ratio_tol = 5e-9
  difference_tol = 5e-10
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[dynamic_viscosity]
  type = 'PetscJacobianTester'
  input = 'dynamic_viscosity.i'
  ratio_tol = 1e-7
  difference_tol = 1e-7
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[linear_test_material]
  type = 'PetscJacobianTester'
  input = 'linear_test_material.i'
  ratio_tol = 5e-14
  difference_tol = 5e-14
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[pressure_3eqn_material]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=p"
  ratio_tol = 1e-8
  difference_tol = 8e-4
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[specific_internal_energy_3eqn]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=e JacobianTest1Phase/snes_test_err=1e-8"
  ratio_tol = 5e-8
  difference_tol = 5e-7
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[specific_volume_3eqn]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=v JacobianTest1Phase/snes_test_err=1e-8"
  ratio_tol = 5e-8
  difference_tol = 5e-14
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[velocity_3eqn]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=v Kernels/test_kernel/coupled_variables='rhoA rhouA' JacobianTest1Phase/snes_test_err=1e-8"
  ratio_tol = 5e-8
  difference_tol = 5e-8
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[enthalpy_3eqn]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=h JacobianTest1Phase/snes_test_err=1e-8"
  ratio_tol = 1e-7
  difference_tol = 5e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[reynolds_number]
  type = 'PetscJacobianTester'
  input = 'reynolds_number.i'
  ratio_tol = 5e-8
  difference_tol = 1e-8
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[temperature_3eqn_material]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=T"
  ratio_tol = 2e-8
  difference_tol = 3e-10
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[total_enthalpy_3eqn]
  type = 'PetscJacobianTester'
  input = 'fluid_properties_3eqn.i'
  cli_args = "Kernels/test_kernel/material_property=H JacobianTest1Phase/snes_test_err=1e-8"
  ratio_tol = 5e-8
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[wall_friction_function]
  type = 'PetscJacobianTester'
  input = 'wall_friction_function.i'
  ratio_tol = 1e-14
  difference_tol = 1e-14
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[wall_friction_churchill]
  type = 'PetscJacobianTester'
  input = 'wall_friction_churchill.i'
  ratio_tol = 2e-10
  difference_tol = 2e-10
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[rhoV]
  type = 'PetscJacobianTester'
  input = '1phase.i'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[rhouV]
  type = 'PetscJacobianTester'
  input = '1phase.i'
  cli_args = '
      ScalarKernels/scalar_kernel/variable=rhouV
      ScalarKernels/scalar_kernel/equation_index=1'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[rhovV]
  type = 'PetscJacobianTester'
  input = '1phase.i'
  cli_args = '
      ScalarKernels/scalar_kernel/variable=rhovV
      ScalarKernels/scalar_kernel/equation_index=2'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[rhowV]
  type = 'PetscJacobianTester'
  input = '1phase.i'
  cli_args = '
      ScalarKernels/scalar_kernel/variable=rhowV
      ScalarKernels/scalar_kernel/equation_index=3'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[rhoEV]
  type = 'PetscJacobianTester'
  input = '1phase.i'
  cli_args = '
      ScalarKernels/scalar_kernel/variable=rhoEV
      ScalarKernels/scalar_kernel/equation_index=4'
  ratio_tol = 1e-6
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'shaft_connected_motor.i'
  ratio_tol = 1e-8
  difference_tol = 1e-7
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[test]
  type = 'PetscJacobianTester'
  input = 'shaft_connected_test_component.i'
  ratio_tol = 5e-8
  difference_tol = 1e-5
  max_threads = 1
  max_parallel = 1
  allow_test_objects = true
[]

[ad_average_wall_temperature_3eqn]
  type = 'CSVDiff'
  input = 'ad_average_wall_temperature_3eqn.i'
  csvdiff = 'ad_average_wall_temperature_3eqn_out.csv'
[]

[ad_alt_def]
  type = 'CSVDiff'
  input = 'ad_average_wall_temperature_3eqn.i'
  cli_args = "Outputs/file_base=alt_def_out Materials/const_materials/prop_values='0 0' AuxKernels/P_hf1_auxkernel/value=0.5e-5 AuxKernels/P_hf2_auxkernel/value=0.5e-5"
  csvdiff = 'alt_def_out.csv'
[]

[ad_err.n_hw_sources]
  type = 'RunException'
  input = 'ad_average_wall_temperature_3eqn.i'
  cli_args = "Materials/T_wall_avg_material/Hw_sources='Hw1' Materials/T_wall_avg_material/P_hf_sources='P_hf1 P_hf2'"
  expect_err = "T_wall_avg_material: The number of wall heat transfer coefficient values must equal the number of wall temperature values"
[]

[ad_err.n_phf_sources]
  type = 'RunException'
  input = 'ad_average_wall_temperature_3eqn.i'
  cli_args = "Materials/T_wall_avg_material/Hw_sources='Hw1 Hw2' Materials/T_wall_avg_material/P_hf_sources='P_hf1'"
  expect_err = "T_wall_avg_material: The number of heated perimeter values must equal the number of wall temperature values"
[]

[test]
  type = 'CSVDiff'
  input = 'ad_convection_heat_flux.i'
  csvdiff = 'ad_convection_heat_flux_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'ad_convection_heat_flux_hs.i'
  csvdiff = 'ad_convection_heat_flux_hs_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'ad_material_function_product.i'
  csvdiff = 'ad_material_function_product_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test_name]
  type = 'CSVDiff'
  cli_args = "Materials/Hw_material/Hw=htc AuxKernels/Hw_ak/property=htc"
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[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'
[]

[average_wall_temperature_3eqn]
  type = 'CSVDiff'
  input = 'average_wall_temperature_3eqn.i'
  csvdiff = 'average_wall_temperature_3eqn_out.csv'
[]

[alt_def]
  type = 'CSVDiff'
  input = 'average_wall_temperature_3eqn.i'
  cli_args = "Outputs/file_base=alt_def_out Materials/const_materials/prop_values='0 0' AuxKernels/P_hf1_auxkernel/value=0.5e-5 AuxKernels/P_hf2_auxkernel/value=0.5e-5"
  csvdiff = 'alt_def_out.csv'
[]

[err.n_hw_sources]
  type = 'RunException'
  input = 'average_wall_temperature_3eqn.i'
  cli_args = "Materials/T_wall_avg_material/Hw_sources='Hw1' Materials/T_wall_avg_material/P_hf_sources='P_hf1 P_hf2'"
  expect_err = "T_wall_avg_material: The number of wall heat transfer coefficient values must equal the number of wall temperature values"
[]

[err.n_phf_sources]
  type = 'RunException'
  input = 'average_wall_temperature_3eqn.i'
  cli_args = "Materials/T_wall_avg_material/Hw_sources='Hw1 Hw2' Materials/T_wall_avg_material/P_hf_sources='P_hf1'"
  expect_err = "T_wall_avg_material: The number of heated perimeter values must equal the number of wall temperature values"
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[churchill]
  type = 'CSVDiff'
  input = 'churchill.i'
  csvdiff = 'churchill_out.csv'
  recover = false
[]

[churchill_ad]
  type = 'CSVDiff'
  input = 'churchill_ad.i'
  csvdiff = 'churchill_ad_out.csv'
  recover = false
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test_name]
  type = 'CSVDiff'
  cli_args = "Materials/Hw_material/Hw=htc AuxKernels/Hw_ak/property=htc"
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[single_block]
  type = 'Exodiff'
  input = 'single_block.i'
  exodiff = 'single_block_out.e-s005'
  abs_zero = 1e-7
  recover = false
  mesh_mode = 'replicated'
[]

[count_iterations]
  type = 'CSVDiff'
  input = 'count_iterations.i'
  csvdiff = 'count_iterations_out.csv'
[]

[coupling_mD_nonoverlapping]
  type = Exodiff
  input = 'parent_non_overlapping.i'
  exodiff = 'parent_non_overlapping_out.e'
  recover = false
[]

[displaced_components]
  type = 'Exodiff'
  input = 'displaced_components.i'
  exodiff = 'displaced_components_out.e'
  custom_cmp = 'displaced_components.exodiff'
  recover = false
[]

[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'
[]

[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'
[]

[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 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'
[]

[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'
[]

[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'
[]

[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_in.e'
  cli_args = '--mesh-only'
  recover = False
[]

[part1]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'part1_out.e'
  cli_args = '
      GlobalParams/initial_T=507
      GlobalParams/initial_p=6.e6
      GlobalParams/initial_vel=0
      GlobalParams/initial_vel_x=0
      GlobalParams/initial_vel_y=0
      GlobalParams/initial_vel_z=0
      Components/hs/initial_T=Ts_init
      Components/hx3ext/initial_T_wall=520
      Executioner/num_steps=5
      Outputs/checkpoint=true
      Outputs/file_base=part1_out'
  recover = false
  abs_zero = 1e-7
[]

[part2]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'part2_out.e'
  cli_args = '
      Problem/restart_file_base=part1_out_cp/0005
      Executioner/num_steps=10
      Outputs/file_base=part2_out'
  abs_zero = 1e-7
  prereq = 'part1'
  recover = false
[]

[part2_ics]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'part2_ics_out.e'
  cli_args = '
      GlobalParams/initial_T=507
      GlobalParams/initial_p=6.e6
      GlobalParams/initial_vel=0.1
      GlobalParams/initial_vel_x=0.1
      GlobalParams/initial_vel_y=0.1
      GlobalParams/initial_vel_z=0.1
      Components/hs/initial_T=508
      Components/hx3ext/initial_T_wall=510
      Components/inlet/vel=0.1
      Components/temp_outside/T=508
      Problem/restart_file_base=part1_out_cp/0005
      Executioner/num_steps=10
      Outputs/file_base=part2_ics_out'
  prereq = 'part1'
  recover = false
  abs_zero = 1e-7
[]

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

[surrogate_power_profile]
  type = 'Exodiff'
  input = 'surrogate_power_profile.i'
  exodiff = 'surrogate_power_profile_out.e'
  prereq = 'power_profile'
  recover = false
[]

[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 = JSONDiff
  input = 'test.i'
  jsondiff = 'test_map.json'
[]

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

[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'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[test]
  type = 'CSVDiff'
  input = 'flow_junction_flux_1phase.i'
  csvdiff = 'flow_junction_flux_1phase_out.csv'
  recover = false
  abs_zero = 2e-8
[]

[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.'
[]

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

[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_convection.i'
  csvdiff = 'heat_rate_convection.csv'
  recover = false
[]

[test]
  type = 'CSVDiff'
  input = 'heat_rate_convection_1phase.i'
  csvdiff = 'heat_rate_convection_1phase_out.csv'
  recover = false
[]

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

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

[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 = 'heat_rate_radiation.i'
  csvdiff = 'heat_rate_radiation.csv'
  recover = false
[]

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

[heat_structure_energy_plate]
  type = 'CSVDiff'
  input = 'heat_structure_energy_plate.i'
  csvdiff = 'heat_structure_energy_plate.csv'
  recover = false
[]

[heat_structure_energy_plate_T_ref]
  type = 'CSVDiff'
  input = 'heat_structure_energy_plate.i'
  cli_args = 'Postprocessors/E_tot/T_ref=0.5 Outputs/file_base=heat_structure_energy_plate_T_ref'
  csvdiff = 'heat_structure_energy_plate_T_ref.csv'
  recover = false
[]

[heat_structure_energy_cylinder]
  type = 'CSVDiff'
  input = 'heat_structure_energy_cylinder.i'
  csvdiff = 'heat_structure_energy_cylinder.csv'
  rel_err = 1e-5
  recover = false
[]

[non_existent_par_name]
  type = 'RunException'
  input = 'non_existent_par_name.i'
  expect_err = "Either component 'inlet' does not exist or parameter 'p' does not exist in that component."
[]

[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.'
[]

[test]
  type = 'CSVDiff'
  input = 'side_flux_integral_rz.i'
  csvdiff = 'side_flux_integral_rz_out.csv'
  recover = false
[]

[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'
[]

[without_junction]
  type = 'Exodiff'
  input = 'area_constriction.i'
  exodiff = 'area_constriction_out.e'
  max_parallel = 1
[]

[with_junction]
  type = 'Exodiff'
  input = 'area_constriction_junction.i'
  exodiff = 'area_constriction_junction_out.e'
  max_parallel = 1
[]

[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'
[]

[1phase]
  type = 'Exodiff'
  input = '1phase.i'
  exodiff = '1phase.e'
  cli_args = 'Executioner/num_steps=5'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
[]

[test]
  type = CSVDiff
  input = 'freefall.i'
  csvdiff = 'freefall_out.csv'
[]

[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.'
[]

[1phase]
  type = 'CSVDiff'
  input = 'mms_1phase.i'
  csvdiff = 'mms_1phase_out.csv'
  allow_test_objects = true
  recover = false
  max_parallel = 1
[]

[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.'
[]

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

[sod_shock_tube_erk2_heavy]
  type = 'Exodiff'
  input = 'sod_shock_tube.i'
  exodiff = 'sod_shock_tube_erk2_heavy.e'
  cli_args = '
      Executioner/TimeIntegrator/order=2
      Outputs/file_base=sod_shock_tube_erk2_heavy
      Outputs/execute_on=final'
  abs_zero = 1e-8
  rel_err = 1e-5
  heavy = true
  recover = false
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'
[]

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

[non_x_direction_1phase]
  type = 'Exodiff'
  input = 'sod_shock_tube.i'
  exodiff = 'non_x_direction_1phase.e'
  cli_args = '
      Functions/p_ic_fn/axis=y
      Functions/T_ic_fn/axis=y
      Components/pipe/position="1 0 0"
      Components/pipe/orientation="0 1 0"
      Executioner/TimeIntegrator/order=2
      Executioner/num_steps=5
      Outputs/file_base=non_x_direction_1phase'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'
[]

[sod_shock_tube_be]
  type = 'Exodiff'
  input = 'sod_shock_tube.i'
  exodiff = 'sod_shock_tube_be.e'
  cli_args = '
      Preconditioning/pc/type=SMP
      Preconditioning/pc/full=true
      Executioner/TimeIntegrator/type=ImplicitEuler
      Executioner/solve_type=NEWTON
      Executioner/num_steps=5
      Outputs/file_base=sod_shock_tube_be'
  abs_zero = 1e-7
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg implicit pipe free_boundary'
[]

[square_wave_minmod]
  type = 'Exodiff'
  input = 'square_wave.i'
  exodiff = 'square_wave_minmod.e'
  cli_args = 'Executioner/num_steps=5 Outputs/file_base=square_wave_minmod'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'
[]

[square_wave_none]
  type = 'Exodiff'
  input = 'square_wave.i'
  exodiff = 'square_wave_none.e'
  cli_args = 'Executioner/num_steps=5 GlobalParams/rdg_slope_reconstruction=none Outputs/file_base=square_wave_none'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'
[]

[square_wave_mc]
  type = 'Exodiff'
  input = 'square_wave.i'
  exodiff = 'square_wave_mc.e'
  cli_args = 'Executioner/num_steps=5 GlobalParams/rdg_slope_reconstruction=mc Outputs/file_base=square_wave_mc'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'
[]

[square_wave_superbee]
  type = 'Exodiff'
  input = 'square_wave.i'
  exodiff = 'square_wave_superbee.e'
  cli_args = 'Executioner/num_steps=5 GlobalParams/rdg_slope_reconstruction=superbee Outputs/file_base=square_wave_superbee'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe free_boundary'
[]

[3eqn]
  type = 'Exodiff'
  input = '3eqn.i'
  exodiff = '3eqn_out.e'
  abs_zero = 1e-6
  max_parallel = 1
[]

[woodward_colella_blast_wave]
  type = 'Exodiff'
  input = 'woodward_colella_blast_wave.i'
  exodiff = 'woodward_colella_blast_wave.e'
  cli_args = 'Executioner/num_steps=5'
  abs_zero = 1e-10
  rel_err = 1e-5
  max_parallel = 1
  group = '1phase rdg explicit pipe solid_wall'
[]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
[]

[postprocessor_source]
  type = 'CSVDiff'
  input = 'postprocessor_source.i'
  csvdiff = 'postprocessor_source_out.csv'
[]

[test]
  type = CSVDiff
  input = 'function_element_loop_integral_uo.i'
  csvdiff = 'function_element_loop_integral_uo_out.csv'
  allow_test_objects = true
[]

[test]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
[]

[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."
[]

[cubic_nonad]
  type = 'CSVDiff'
  input = 'smooth_transition.i'
  csvdiff = 'cubic_nonad_test_vpp_0000.csv'
  allow_test_objects = True
  recover = False
[]

[weighted_nonad]
  type = 'CSVDiff'
  input = 'smooth_transition.i'
  cli_args = "Materials/test_mat/transition_type=weighted Outputs/file_base=weighted_nonad"
  csvdiff = 'weighted_nonad_test_vpp_0000.csv'
  allow_test_objects = True
  recover = False
[]

[cubic_derivative]
  type = 'PetscJacobianTester'
  input = 'smooth_transition_derivative.i'
  ratio_tol = 1e-9
  difference_tol = 1e-9
  max_threads = 1
  allow_test_objects = True
[]

[weighted_derivative]
  type = 'PetscJacobianTester'
  input = 'smooth_transition_derivative.i'
  cli_args = "Materials/test_mat/transition_type=weighted"
  ratio_tol = 1e-7
  difference_tol = 1e-7
  max_threads = 1
  allow_test_objects = True
[]

[ad_cubic]
  type = 'CSVDiff'
  input = 'ad_smooth_transition.i'
  csvdiff = 'ad_cubic_test_vpp_0000.csv'
  cli_args = "
      Materials/test_mat/transition_type=cubic
      Outputs/file_base=ad_cubic"
  allow_test_objects = True
  recover = False
[]

[ad_weighted]
  type = 'CSVDiff'
  input = 'ad_smooth_transition.i'
  csvdiff = 'ad_weighted_test_vpp_0000.csv'
  allow_test_objects = True
  recover = False
[]

[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.'
[]

[sampler_1d_vector]
  type = 'CSVDiff'
  input = 'sampler_1d_vector.i'
  csvdiff = 'out_test_property_vpp_0001.csv'
  allow_test_objects = true
[]