Heat Transfer Requirements Traceability Matrix

[./hex8]
  type = 'CSVDiff'
  input = 'nafems_t3_hex_template.i'
  cli_args = 'Outputs/file_base=nafems_t3_hex8_coarse_out'
  csvdiff = 'nafems_t3_hex8_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'HEX8 elements'
[../]

[./hex20]
  type = 'CSVDiff'
  input = 'nafems_t3_hex_template.i'
  cli_args = 'Mesh/elem_type=HEX20 Postprocessors/target_temp/nodeid=47 Outputs/file_base=nafems_t3_hex20_coarse_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_hex20_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'HEX20 elements'
[../]

[./hex27]
  type = 'CSVDiff'
  input = 'nafems_t3_hex_template.i'
  cli_args = 'Mesh/elem_type=HEX27 Postprocessors/target_temp/nodeid=66 Outputs/file_base=nafems_t3_hex27_coarse_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_hex27_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'HEX27 elements'
[../]

[./edge2]
  type = 'CSVDiff'
  input = 'nafems_t3_edge_template.i'
  cli_args = 'Outputs/file_base=nafems_t3_edge2_coarse_out'
  csvdiff = 'nafems_t3_edge2_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'EDGE2 elements'
[../]

[./edge3]
  type = 'CSVDiff'
  input = 'nafems_t3_edge_template.i'
  cli_args = 'Mesh/elem_type=EDGE3 Postprocessors/target_temp/nodeid=7 Outputs/file_base=nafems_t3_edge3_coarse_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_edge3_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'EDGE3 elements'
[../]

[./quad4]
  type = 'CSVDiff'
  input = 'nafems_t3_quad_template.i'
  cli_args = 'Outputs/file_base=nafems_t3_quad4_coarse_out'
  csvdiff = 'nafems_t3_quad4_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'QUAD4 elements'
[../]

[./quad8]
  type = 'CSVDiff'
  input = 'nafems_t3_quad_template.i'
  cli_args = 'Mesh/elem_type=QUAD8 Postprocessors/target_temp/nodeid=19 Outputs/file_base=nafems_t3_quad8_coarse_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_quad8_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'QUAD8 elements'
[../]

[./quad9]
  type = 'CSVDiff'
  input = 'nafems_t3_quad_template.i'
  cli_args = 'Mesh/elem_type=QUAD9 Postprocessors/target_temp/nodeid=22 Outputs/file_base=nafems_t3_quad9_coarse_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_quad9_coarse_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'QUAD9 elements'
[../]

[./hex8]
  type = 'CSVDiff'
  input = 'nafems_t3_hex_template.i'
  cli_args = 'Mesh/nx=10 Postprocessors/target_temp/nodeid=35 Outputs/file_base=nafems_t3_hex8_fine_out'
  csvdiff = 'nafems_t3_hex8_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'HEX8 mesh'
[../]

[./hex20]
  type = 'CSVDiff'
  input = 'nafems_t3_hex_template.i'
  cli_args = 'Mesh/nx=10 Mesh/elem_type=HEX20 Postprocessors/target_temp/nodeid=95 Outputs/file_base=nafems_t3_hex20_fine_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_hex20_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'HEX20 mesh'
[../]

[./hex27]
  type = 'CSVDiff'
  input = 'nafems_t3_hex_template.i'
  cli_args = 'Mesh/nx=10 Mesh/elem_type=HEX27 Postprocessors/target_temp/nodeid=138 Outputs/file_base=nafems_t3_hex27_fine_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_hex27_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'HEX27 mesh'
[../]

[./edge2]
  type = 'CSVDiff'
  input = 'nafems_t3_edge_template.i'
  cli_args = 'Mesh/nx=10 Postprocessors/target_temp/nodeid=8 Outputs/file_base=nafems_t3_edge2_fine_out'
  csvdiff = 'nafems_t3_edge2_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'EDGE2 mesh'
[../]

[./edge3]
  type = 'CSVDiff'
  input = 'nafems_t3_edge_template.i'
  cli_args = 'Mesh/nx=10 Mesh/elem_type=EDGE3 Postprocessors/target_temp/nodeid=15 Outputs/file_base=nafems_t3_edge3_fine_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_edge3_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'EDGE3 mesh'
[../]

[./quad4]
  type = 'CSVDiff'
  input = 'nafems_t3_quad_template.i'
  cli_args = 'Mesh/nx=10 Postprocessors/target_temp/nodeid=17 Outputs/file_base=nafems_t3_quad4_fine_out'
  csvdiff = 'nafems_t3_quad4_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'QUAD4 mesh'
[../]

[./quad8]
  type = 'CSVDiff'
  input = 'nafems_t3_quad_template.i'
  cli_args = 'Mesh/nx=10 Mesh/elem_type=QUAD8 Postprocessors/target_temp/nodeid=39 Outputs/file_base=nafems_t3_quad8_fine_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_quad8_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'QUAD8 mesh'
[../]

[./quad9]
  type = 'CSVDiff'
  input = 'nafems_t3_quad_template.i'
  cli_args = 'Mesh/nx=10 Mesh/elem_type=QUAD9 Postprocessors/target_temp/nodeid=46 Outputs/file_base=nafems_t3_quad9_fine_out Variables/temp/order=SECOND'
  csvdiff = 'nafems_t3_quad9_fine_out.csv'
  verification = 'nafems_t3_verif.md'
  mesh_mode = REPLICATED
  detail = 'QUAD9 mesh'
[../]

[./flux]
  type = 'CSVDiff'
  input = 'flux.i'
  csvdiff = 'flux_out.csv'
  detail = 'and match hand calculations for flux through a boundary.'
[../]

[./equilibrium]
  type = 'CSVDiff'
  input = 'equilibrium.i'
  csvdiff = 'equilibrium_out.csv'
  detail = 'and approach a constant far-field temperature value over time as heat flux decreases.'
[../]

[./coupled]
  type = 'CSVDiff'
  input = 'coupled.i'
  csvdiff = 'coupled_out.csv'
  detail = 'and couple a temperature dependent far-field temperature and heat transfer coefficient.'
  allow_test_objects = true
[../]

[./jacobian_test]
  type = 'PetscJacobianTester'
  input = 'test.i'
  ratio_tol = 1e-7
  difference_tol = 1e-5
  run_sim = True
  requirement = 'AD heat conduction and the Jacobian shall be beautiful'
  design = "ADHeatConduction.md"
  issues = "#5658 #12633"
  allow_test_objects = true
[../]

[spatial_csv]
  type = CSVDiff
  command = mes_spatial.py
  recover = false
  csvdiff = 'cartesian_test_no1_first_order.csv'
  requirement = 'The MOOSE solutions shall converge to the analytic solutions with an expected order of accuracy (two for linear, three for quadratic) where a standard set of heat conduction problems is used for code verification.'
  valgrind = none
  method = opt
  heavy = true
  required_python_packages = 'sympy'
[]

[./convection]
  type = Exodiff
  input = conjugate_heat_transfer.i
  exodiff = conjugate_heat_transfer_out.e
  abs_zero = 1e-8
  requirement = "The system shall correctly model convection heat transfer across internal sidesets aka conjugate heat transfer."
  issues = '#15114'
  design = 'ConjugateHeatTransfer.md'
[../]

[./constant]
  type = 'CSVDiff'
  input = 'convective_flux_function.i'
  csvdiff = 'convective_flux_function_out.csv'
  requirement = 'The system shall allow prescribing a convective flux boundary condition using a constant heat transfer coefficient.'
[../]

[./time_dependent]
  prereq = constant
  type = 'CSVDiff'
  input = 'convective_flux_function.i'
  cli_args = "BCs/right/coefficient='t*10.0'"
  csvdiff = 'convective_flux_function_out.csv'
  requirement = 'The system shall allow prescribing a convective flux boundary condition using a heat transfer coefficient that is a function of position and time.'
[../]

[./temperature_dependent]
  prereq = time_dependent
  type = 'CSVDiff'
  input = 'convective_flux_function.i'
  cli_args = "BCs/right/coefficient='t/15.0' BCs/right/coefficient_function_type=TEMPERATURE"
  csvdiff = 'convective_flux_function_out.csv'
  requirement = 'The system shall allow prescribing a convective flux boundary condition using a heat transfer coefficient that is a function of temperature.'
[../]

[./flux]
  type = 'CSVDiff'
  input = 'flux.i'
  csvdiff = 'flux_out.csv'
  detail = 'and match hand calculations for flux through a boundary.'
[../]

[./equilibrium]
  type = 'CSVDiff'
  input = 'equilibrium.i'
  csvdiff = 'equilibrium_out.csv'
  detail = 'and approach a constant far-field temperature value over time as heat flux decreases.'
[../]

[./coupled]
  type = 'CSVDiff'
  input = 'coupled.i'
  csvdiff = 'coupled_out.csv'
  detail = 'and couple a temperature dependent far-field temperature and heat transfer coefficient.'
[../]

[2d]
  type = Exodiff
  input = 2d.i
  exodiff = 2d_out.e
  requirement = 'The system shall be able to apply a directed flux to the surface of a 2d geometry with and without self shadowing'
[]

[3d]
  type = Exodiff
  input = 3d.i
  exodiff = 3d_out.e
  installation_type = in_tree
  requirement = 'The system shall be able to apply a directed flux to the surface of a 3d geometry with and without self shadowing'
  recover = false # outputs only on FINAL; recover is fully tested through the 2d case
[]

[quad4]
  type = CSVDiff
  input = 2d_elem.i
  csvdiff = 2d_quad4.csv
  cli_args = 'GlobalParams/elem_type=quad4 Outputs/file_base=2d_quad4'
  detail = 'QUAD4'
[]

[quad8]
  type = CSVDiff
  input = 2d_elem.i
  csvdiff = 2d_quad8.csv
  cli_args = 'GlobalParams/elem_type=quad8 GlobalParams/order=second Outputs/file_base=2d_quad8'
  detail = 'QUAD8'
[]

[quad9]
  type = CSVDiff
  input = 2d_elem.i
  csvdiff = 2d_quad9.csv
  cli_args = 'GlobalParams/elem_type=quad9 GlobalParams/order=second Outputs/file_base=2d_quad9'
  detail = 'QUAD9'
[]

[tri3]
  type = CSVDiff
  input = 2d_elem.i
  csvdiff = 2d_tri3.csv
  cli_args = 'GlobalParams/elem_type=tri3 Outputs/file_base=2d_tri3'
  detail = 'TRI3'
[]

[tri6]
  type = CSVDiff
  input = 2d_elem.i
  csvdiff = 2d_tri6.csv
  cli_args = 'GlobalParams/elem_type=tri6 GlobalParams/order=second Outputs/file_base=2d_tri6'
  detail = 'TRI6'
[]

[tri7]
  type = CSVDiff
  input = 2d_elem.i
  csvdiff = 2d_tri7.csv
  cli_args = 'GlobalParams/elem_type=tri7 GlobalParams/order=third Outputs/file_base=2d_tri7'
  detail = 'TRI7'
[]

[tet4]
  type = CSVDiff
  input = 3d_elem.i
  csvdiff = 3d_tet4.csv
  cli_args = 'GlobalParams/elem_type=tet4 Outputs/file_base=3d_tet4'
  detail = 'TET4'
[]

[tet10]
  type = CSVDiff
  input = 3d_elem.i
  csvdiff = 3d_tet10.csv
  cli_args = 'GlobalParams/elem_type=tet10 GlobalParams/order=second Outputs/file_base=3d_tet10'
  detail = 'TET10'
  heavy = true
[]

[tet14]
  type = CSVDiff
  input = 3d_elem.i
  csvdiff = 3d_tet14.csv
  cli_args = 'GlobalParams/elem_type=tet14 GlobalParams/order=third Outputs/file_base=3d_tet14'
  detail = 'TET14'
  heavy = true
[]

[hex8]
  type = CSVDiff
  input = 3d_elem.i
  csvdiff = 3d_hex8.csv
  cli_args = 'GlobalParams/elem_type=hex8 Outputs/file_base=3d_hex8'
  detail = 'HEX8'
[]

[hex20]
  type = CSVDiff
  input = 3d_elem.i
  csvdiff = 3d_hex20.csv
  cli_args = 'GlobalParams/elem_type=hex20 GlobalParams/order=second Outputs/file_base=3d_hex20'
  detail = 'HEX20'
[]

[hex27]
  type = CSVDiff
  input = 3d_elem.i
  csvdiff = 3d_hex27.csv
  cli_args = 'GlobalParams/elem_type=hex27 GlobalParams/order=second Outputs/file_base=3d_hex27'
  detail = 'HEX27'
[]

[prism6]
  type = RunApp
  input = 3d_elem.i
  cli_args = 'GlobalParams/elem_type=prism6 Outputs/file_base=3d_prism6'
  detail = 'PRISM6'
[]

[prism15]
  type = RunApp
  input = 3d_elem.i
  cli_args = 'GlobalParams/elem_type=prism15 GlobalParams/order=second  Outputs/file_base=3d_prism15'
  detail = 'PRISM15'
[]

[prism18]
  type = RunApp
  input = 3d_elem.i
  cli_args = 'GlobalParams/elem_type=prism18 GlobalParams/order=second  Outputs/file_base=3d_prism18'
  detail = 'PRISM18'
[]

[pyramid5]
  type = RunApp
  input = 3d_elem.i
  cli_args = 'GlobalParams/elem_type=pyramid5 Outputs/file_base=3d_pyramid5'
  detail = 'PYRAMID5'
[]

[pyramid13]
  type = RunApp
  input = 3d_elem.i
  cli_args = 'GlobalParams/elem_type=pyramid13 GlobalParams/order=second  Outputs/file_base=3d_pyramid13'
  detail = 'PYRAMID13'
  heavy = true
[]

[pyramid14]
  type = RunApp
  input = 3d_elem.i
  cli_args = 'GlobalParams/elem_type=pyramid14 GlobalParams/order=second  Outputs/file_base=3d_pyramid14'
  detail = 'PYRAMID14'
  heavy = true
[]

[./test]
  type = 'CSVDiff'
  input = 'function_heat_source.i'
  csvdiff = 'function_heat_source_out.csv'
  requirement = 'The system shall produce a moving heat source where its path is function dependent'
[../]

[P1_radiation_material_test]
  type = CSVDiff
  input = 'P1_radiation_material_test.i'
  csvdiff = 'P1_radiation_material_test_out.csv'
  requirement = 'The system shall be able to compute a P1 diffusion coefficient from opacity and effective scattering.'
[]

[test]
  type = CSVDiff
  input = 'convection_heat_flux.i'
  csvdiff = 'convection_heat_flux_out.csv'
  requirement = 'The system shall be able to compute a convection heat flux functor material property.'
[]

[test]
  type = Exodiff
  input = 'cylindrical_gap_heat_flux_functor_material.i'
  exodiff = 'cylindrical_gap_heat_flux_functor_material_out.e'
  requirement = 'The system shall compute functor material properties for heat flux across a cylindrical gap.'
  design = 'source/functormaterials/CylindricalGapHeatFluxFunctorMaterial.md'
  issues = '#25113'
[]

[test]
  type = CSVDiff
  input = 'fin_efficiency.i'
  csvdiff = 'fin_efficiency_out.csv'
  requirement = 'The system shall be able to compute a fin efficiency functor material property.'
[]

[test]
  type = CSVDiff
  input = 'fin_enhancement_factor.i'
  csvdiff = 'fin_enhancement_factor_out.csv'
  requirement = 'The system shall be able to compute a fin heat transfer enhancement factor functor material property.'
[]

[test]
  type = Exodiff
  input = fv_functor_convective_heat_flux.i
  exodiff = 'fv_functor_convective_heat_flux_out.e'
  requirement = 'The system shall provide a convective heat flux boundary condition which uses functors as heat transfer coefficients and far-field temperature values'
[]

[1D]
  type = CSVDiff
  input = 'rad_istothermal_medium_1d.i'
  csvdiff = 'rad_istothermal_medium_1d_out.csv'
  abs_zero = 1e-10
  detail = " in 1D."
[]

[2D]
  type = CSVDiff
  input = 'rad_istothermal_medium_2d.i'
  csvdiff = 'rad_istothermal_medium_2d_out.csv'
  abs_zero = 1e-10
  detail = " in 2D."
[]

[infinite_cylinder_radiation]
  issues = '#18626'
  design = 'FVInfiniteCylinderRadiativeBC.md'
  requirement = 'The system shall be able to solve a heat conduction problem with boundary conditions representing radiation to an infinite cylinder.'
  type = Exodiff
  input = test.i
  exodiff = 'test_out.e'
[]

[matprop]
  type = Exodiff
  input = test.i
  exodiff = 'test_out.e'
  allow_warnings = true # thermal resistance bc has an inner solve that may not always converge enough
  detail = 'using regular material properties.'
[]

[matprop_rz]
  type = Exodiff
  input = test.i
  exodiff = 'test_rz_out.e'
  allow_warnings = true # thermal resistance bc has an inner solve that may not always converge enough
  detail = 'using regular material properties in RZ geometry.'
  # We remove the thermal resistance BC from the top, since the conduction formulation there
  # is incorrect. We also move the thermal resistance BC to the outer cylinder, makes more sense than at r=0
  cli_args = "Problem/coord_type=RZ FVBCs/left/boundary=right FVBCs/left/geometry=cylindrical FVBCs/left/inner_radius=2 FVBCs/other/boundary='left bottom' Outputs/file_base=test_rz_out"
[]

[functor_matprop]
  type = Exodiff
  input = test_functor.i
  exodiff = 'test_functor_out.e'
  allow_warnings = true # thermal resistance bc has an inner solve that may not always converge enough
  detail = 'using functor material properties.'
[]

[functor_matprop_rz]
  type = Exodiff
  input = test_functor.i
  exodiff = 'test_functor_rz_out.e'
  allow_warnings = true # thermal resistance bc has an inner solve that may not always converge enough
  detail = 'using functor material properties in RZ geometry.'
  # We remove the thermal resistance BC from the top, since the conduction formulation there
  # is incorrect. We also move the thermal resistance BC to the outer cylinder, makes more sense than at r=0
  cli_args = "Problem/coord_type=RZ FVBCs/left/boundary=right FVBCs/left/geometry=cylindrical FVBCs/left/inner_radius=2 FVBCs/other/boundary='left bottom' Outputs/file_base=test_functor_rz_out"
[]

[1d_P1_rad_heat_trasfer]
  type = CSVDiff
  input = 'radiation_istothermal_medium_1d.i'
  csvdiff = 'radiation_istothermal_medium_1d_out.csv'
  abs_zero = 1e-4
  requirement = "P1 Radiation heat transfer should match an analytical solution in 1D."
[]

[large_gap_heat_transfer_test_sphere]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_sphere.i'
  csvdiff = 'large_gap_heat_transfer_test_sphere_out.csv'
  abs_zero = 1e-6
  rel_err = 1e-6
  requirement = "Energy balance must be fulfilled for the heat transfer of concentric spheres involving radiation, when the gap distance is not negligible with respect to the body main dimensions."
[]

[large_gap_heat_transfer_test_rz_cylinder]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_rz_cylinder.i'
  csvdiff = 'large_gap_heat_transfer_test_rz_cylinder_out.csv'
  abs_zero = 1e-6
  rel_err = 1e-6
  requirement = "Energy balance must be fulfilled for the heat transfer of concentric cylinders involving radiation in two-dimensions, when the gap distance is not negligible with respect to the body main dimensions."
[]

[large_gap_heat_transfer_test_cylinder]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_cylinder.i'
  csvdiff = 'large_gap_heat_transfer_test_cylinder_out.csv'
  abs_zero = 1e-6
  rel_err = 1e-6
  requirement = "Energy balance must be fulfilled for the heat transfer of concentric cylinders involving radiation in two-dimensions with axisymmetry, when the gap distance is not negligible with respect to the body main dimensions."
[]

[1D]
  type = Exodiff
  input = 1D.i
  exodiff = 1D_out.e
  requirement = "The system shall compute thermal contact in 1D with axisymmetric coordinates."
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#27216'
[]

[3D]
  type = Exodiff
  input = 'gap_heat_transfer_htonly_test.i'
  exodiff = 'gap_heat_transfer_htonly_test_out.e'
  abs_zero = 1e-6
  requirement = "Thermal contact shall solve plate heat transfer for a constant conductivity gap in 3D"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#1609'
[]

[syntax]
  type = Exodiff
  input = 'gap_heat_transfer_htonly_syntax.i'
  exodiff = 'gap_heat_transfer_htonly_syntax_out.e'
  abs_zero = 1e-6
  requirement = "Thermal contact shall solve plate heat transfer for a constant conductivity gap in 3D using the Modules/HeatConduction/Thermal contact syntax"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#1609'
[]

[3D_Iters]
  type = Exodiff
  input = 'gap_heat_transfer_htonly_it_plot_test.i'
  exodiff = 'out_it_plot.e'
  custom_cmp = 'gap_heat_transfer_htonly_it_plot_test.cmp'
  requirement = "Thermal contact shall solve plate heat transfer for a constant conductivity gap in 3D at each iteration"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#1609'
[]

[RZ]
  type = Exodiff
  input = 'gap_heat_transfer_htonly_rz_test.i'
  exodiff = 'gap_heat_transfer_htonly_rz_test_out.e'
  abs_zero = 1e-9
  rel_err = 5e-3
  compiler = 'CLANG GCC'
  requirement = "Thermal contact shall solve cylindrical and plate heat transfer for a constant conductivity gap in 2D axisymmetric coordinates"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#5104'
[]

[ZR]
  type = Exodiff
  input = 'gap_heat_transfer_htonly_rz_test.i'
  exodiff = 'gap_heat_transfer_htonly_zr_test_out.e'
  abs_zero = 1e-9
  rel_err = 5e-3
  compiler = 'CLANG GCC'
  requirement = "Thermal contact shall solve cylindrical heat transfer for a constant conductivity gap in 2D axisymmetric coordinates where the axial axis is along the x-direction"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#12071'
  cli_args = "Problem/rz_coord_axis=X Mesh/active='file rotate' Outputs/file_base=gap_heat_transfer_htonly_zr_test_out"
[]

[RSpherical]
  type = Exodiff
  input = 'gap_heat_transfer_htonly_rspherical.i'
  exodiff = 'gap_heat_transfer_htonly_rspherical_out.e'
  custom_cmp = 'gap_heat_transfer_htonly_rspherical.exodiff'
  compiler = 'CLANG GCC'
  requirement = "Thermal contact shall solve spherical heat transfer for a constant conductivity gap in 1D spherically symmetric coordinates"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#1609 #5104'
[]

[cyl3D]
  type = Exodiff
  input = 'cyl3D.i'
  exodiff = 'cyl3D_out.e'
  superlu = true
  requirement = "Thermal contact shall solve cylindrical heat transfer for a constant conductivity gap in 3D"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#6161'
[]

[cyl2D]
  type = Exodiff
  input = 'cyl2D.i'
  exodiff = 'cyl2D_out.e'
  superlu = true
  requirement = "Thermal contact shall solve cylindrical heat transfer for a constant conductivity gap in the x-y plane"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#6161'
[]

[sphere3D]
  type = Exodiff
  input = 'sphere3D.i'
  exodiff = 'sphere3D_out.e'
  superlu = true
  requirement = "Thermal contact shall solve spherical heat transfer for a constant conductivity gap in 3D"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#6161'
[]

[sphere2DRZ]
  type = Exodiff
  input = 'sphere2DRZ.i'
  exodiff = 'sphere2DRZ_out.e'
  superlu = true
  requirement = "Thermal contact shall solve spherical heat transfer for a constant conductivity gap in 2D axisymmetric coordinates"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#6161'
[]

[cyl2D_xz]
  type = Exodiff
  input = 'cyl2D_xz.i'
  exodiff = 'cyl2D_xz_out.e'
  superlu = true
  requirement = "Thermal contact shall solve cylindrical heat transfer for a constant conductivity gap in the x-z plane"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#11913'
[]

[cyl2D_yz]
  type = Exodiff
  input = 'cyl2D_yz.i'
  exodiff = 'cyl2D_yz_out.e'
  superlu = true
  requirement = "Thermal contact shall solve cylindrical heat transfer for a constant conductivity gap in the y-z plane"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#11913'
[]

[planar_xy]
  type = Exodiff
  input = 'planar_xy.i'
  exodiff = 'planar_xy_out.e'
  custom_cmp = 'planar.cmp'
  superlu = true
  requirement = "Thermal contact shall solve plate heat transfer for a constant conductivity gap in the x-y plane"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#11913'
[]

[planar_xz]
  type = Exodiff
  input = 'planar_xz.i'
  exodiff = 'planar_xz_out.e'
  custom_cmp = 'planar.cmp'
  superlu = true
  requirement = "Thermal contact shall solve plate heat transfer for a constant conductivity gap in the x-z plane"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#11913'
[]

[planar_yz]
  type = Exodiff
  input = 'planar_yz.i'
  exodiff = 'planar_yz_out.e'
  custom_cmp = 'planar.cmp'
  superlu = true
  requirement = "Thermal contact shall solve plate heat transfer for a constant conductivity gap in the y-z plane"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#11913'
[]

[corner_wrap]
  type = Exodiff
  input = 'corner_wrap.i'
  exodiff = 'corner_wrap_out.e'
  requirement = "When the 'check_boundary_restricted' option is set to false, the thermal contact system shall solve problems in which multiple faces of an element are in one of the contact surfaces, but provide an information message that contact variables may have issues in those areas."
  expect_out = "This parameter is set to 'false'. Although thermal contact will be correctly enforced, the contact-related output may have issues in cases where where more than one face of an element belongs to a contact surface because the values from only one of the faces will be reported."
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#23058'
[]

[corner_wrap_err_check_true]
  prereq = corner_wrap
  type = RunException
  input = 'corner_wrap.i'
  cli_args = ThermalContact/thermal_contact/check_boundary_restricted=true
  requirement = "When the 'check_boundary_restricted' option is set to true, the AuxKernels set up by the ThermalContact system shall generate an error when multiple faces of an element are in one of the contact surfaces."
  expect_err = "Boundary restricted auxiliary kernel 'gap_value_thermal_contact_0' has element "
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#23058'
[]

[corner_wrap_err_quadrature]
  prereq = corner_wrap_err_check_true
  type = RunException
  input = 'corner_wrap.i'
  cli_args = ThermalContact/thermal_contact/quadrature=false
  requirement = "If 'quadrature=false' the thermal contact system shall generate an error if the user also sets 'check_boundary_restricted=true'."
  expect_err = "This parameter cannot be 'false' when 'quadrature=false'"
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#23058'
[]

[test]
  type = Exodiff
  input = 'gap_heat_transfer_mortar.i'
  exodiff = 'gap_heat_transfer_mortar_out.e'
  custom_cmp = 'one-variable-temp.cmp'
  partial = True
  map = False
  requirement = "We shall be able to produce the expected result for a gap conductance test case using the mortar method."
[]

[modular]
  type = Exodiff
  input = 'modular_gap_heat_transfer_mortar.i'
  exodiff = 'modular_gap_heat_transfer_mortar_out.e'
  custom_cmp = 'one-variable-temp.cmp'
  partial = True
  map = False
  requirement = "We shall be able to produce the expected result for a gap conductance test case using the mortar method using the modular gap flux system."
[]

[modular_multiple]
  type = 'Exodiff'
  input = 'gap_heat_transfer_radiation_test.i'
  exodiff = 'gap_heat_transfer_radiation_test_out.e'
  abs_zero = 1e-06
  requirement = "We shall be able to produce the expected result for a combined gap conductance and radiative heat transfer test case using the mortar method using the modular gap flux system"
[]

[displaced]
  type = Exodiff
  input = 'gap_heat_transfer_mortar_displaced.i'
  exodiff = 'gap_heat_transfer_mortar_displaced_out.e'
  custom_cmp = 'one-variable-temp.cmp'
  partial = True
  map = False
  requirement = "We shall be able to run the mortar method on a displaced mesh, supplying the displacements with **constant** auxiliary variables"
[]

[modular_displaced]
  type = Exodiff
  input = 'modular_gap_heat_transfer_mortar_displaced.i'
  exodiff = 'modular_gap_heat_transfer_mortar_displaced_out.e'
  custom_cmp = 'one-variable-temp.cmp'
  partial = True
  map = False
  requirement = "We shall be able to produce the expected result for a gap conductance test case using the mortar method using the modular gap flux system with a displaced mesh."
[]

[displaced_rz]
  type = Exodiff
  input = 'gap_heat_transfer_mortar_displaced.i'
  exodiff = 'gap_heat_transfer_mortar_displaced_rz_out.e'
  cli_args = 'Outputs/file_base=gap_heat_transfer_mortar_displaced_rz_out Problem/coord_type=RZ'
  requirement = "The system shall accurately calculate axisymmetric coordinates on mortar finite element segments"
[]

[bc_gap_heat_transfer_displaced_radiation]
  type = CSVDiff
  input = 'bc_gap_heat_transfer_displaced_radiation.i'
  csvdiff = 'bc_gap_heat_transfer_displaced_radiation_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to generate node-to-segment numerical results for radiation through plates and use it as a reference for mortar-based constraints."
[]

[modular_gap_heat_transfer_mortar_displaced_radiation]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_out_NodalTemperature_0001.csv'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for radiation through plates that are close to those generated by the node-to-segment formulation."
[]

[bc_gap_heat_transfer_displaced_conduction]
  type = CSVDiff
  input = 'bc_gap_heat_transfer_displaced_conduction.i'
  csvdiff = 'bc_gap_heat_transfer_displaced_conduction_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to generate node-to-segment numerical results for conduction through plates and use it as a reference for mortar-based constraints."
[]

[modular_gap_heat_transfer_mortar_displaced_conduction]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_conduction.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_conduction_out_NodalTemperature_0001.csv'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction through plates that are close to those generated by the node-to-segment formulation."
[]

[large_gap_heat_transfer_test_cylinder]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_cylinder.i'
  csvdiff = 'large_gap_heat_transfer_test_cylinder_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to generate node-to-segment numerical results for conduction and radiation through cylinders and use it as a reference for mortar-based constraints."
[]

[large_gap_heat_transfer_test_cylinder_mortar]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_cylinder_mortar.i'
  csvdiff = 'large_gap_heat_transfer_test_cylinder_mortar_out_NodalTemperature_0001.csv'
  cli_args = 'Constraints/ced/cylinder_axis_point_1="0 0 0" Constraints/ced/cylinder_axis_point_2="0 0 5"'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation through cylinders that are close to those generated by the node-to-segment formulation, where the cylinder axis is deduced automatically."
[]

[large_gap_heat_transfer_test_cylinder_mortar_auto]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_cylinder_mortar.i'
  csvdiff = 'large_gap_heat_transfer_test_cylinder_mortar_out_NodalTemperature_0001.csv'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation through cylinders that are close to those generated by the node-to-segment formulation."
[]

[large_gap_heat_transfer_test_2d_sphere_mortar]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_cylinder_mortar.i'
  csvdiff = 'large_gap_heat_transfer_test_cylinder_mortar_out_NodalTemperature_0001.csv'
  cli_args = 'Constraints/ced/sphere_origin="0 0 0" Constraints/ced/gap_geometry_type=SPHERE'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation through spheres that are close to those generated by the node-to-segment formulation."
[]

[large_gap_heat_transfer_test_2d_sphere_mortar_auto]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_cylinder_mortar.i'
  csvdiff = 'large_gap_heat_transfer_test_cylinder_mortar_out_NodalTemperature_0001.csv'
  cli_args = 'Constraints/ced/gap_geometry_type=SPHERE'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation through spheres that are close to those generated by the node-to-segment formulation, where the sphere origin is deduced automatically."
[]

[large_gap_heat_transfer_test_rz_cylinder]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_rz_cylinder.i'
  csvdiff = 'large_gap_heat_transfer_test_rz_cylinder_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to generate node-to-segment numerical results for conduction and radiation through cylinders with axisymmetry and use it as a reference for mortar-based constraints."
[]

[large_gap_heat_transfer_test_rz_cylinder_mortar]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_rz_cylinder_mortar.i'
  csvdiff = 'large_gap_heat_transfer_test_rz_cylinder_mortar_out_NodalTemperature_0001.csv'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction and radiationthrough cylinders with axisymmetry that are close to those generated by the node-to-segment formulation."
[]

[large_gap_heat_transfer_test_sphere]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_sphere.i'
  csvdiff = 'large_gap_heat_transfer_test_sphere_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to generate node-to-segment numerical results for conduction and radiation through concentric spheres with axisymmetry and use it as a reference for mortar-based constraints."
[]

[large_gap_heat_transfer_test_sphere_mortar]
  type = CSVDiff
  input = 'large_gap_heat_transfer_test_sphere_mortar.i'
  csvdiff = 'large_gap_heat_transfer_test_sphere_mortar_out_NodalTemperature_0001.csv'
  rel_err = 0.06
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation through concentric spheres with axisymmetry that are close to those generated by the node-to-segment formulation."
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation in two dimensions. This test is used as a reference for computing separate gap physics, i.e. for the use of multiple heat flux Lagrange multipliers"
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_separate]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_separate.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_separate_out_NodalTemperature_0001.csv'
  requirement = "We shall be able to reproduce heat transfer mortar results when the gap physics (i.e. radiation and conduction) are separated in two constraint classes with independent Lagrange multipliers"
[]

[large_gap_heat_transfer_test_sphere_mortar_geometry_error]
  type = RunException
  input = 'large_gap_heat_transfer_test_sphere_mortar_error.i'
  expect_err = 'gap_geometry_type = SPHERE'
  requirement = "We shall be able to generate a meaningful error message if the user does not define the sphere origin when a spherical geometry has been chosen."
[]

[modular_gap_heat_transfer_mortar_displaced_conduction_function]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_conduction_function.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_conduction_function_out_NodalTemperature_0001.csv'
  rel_err = 1.0e-3
  requirement = "We shall be able to generate mortar numerical results for conduction through plates that are close to those generated by the node-to-segment formulation and using the function feature to enrich the evolution of the gap conductance."
[]

[gap_heat_transfer_3d]
  input = 'gap_heat_transfer_3D.i'
  type = CSVDiff
  csvdiff = 'gap_heat_transfer_3D_out_NodalTemperature_0001.csv'
  rel_err = 1e-3
  abs_zero = 1e-8
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to generate node-to-segment numerical results for conduction and radiation between two blocks in 3D and use it as a reference for mortar-based constraints."
[]

[gap_heat_transfer_3d_hex20]
  input = 'gap_heat_transfer_3D.i'
  type = CSVDiff
  csvdiff = 'gap_heat_transfer_3D_hex20_out_NodalTemperature_0001.csv'
  cli_args = 'Mesh/left_block/elem_type=HEX20 Mesh/right_block/elem_type=HEX20 Outputs/file_base=gap_heat_transfer_3D_hex20_out'
  rel_err = 1e-3
  abs_zero = 1e-8
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to generate node-to-segment numerical results for conduction and radiation between two blocks in 3D using HEX20 elements and use it as a reference for mortar-based constraints."
[]

[gap_heat_transfer_3d_mortar]
  input = 'gap_heat_transfer_3D_mortar.i'
  type = CSVDiff
  csvdiff = 'gap_heat_transfer_3D_mortar_out_NodalTemperature_0001.csv'
  rel_err = 4e-2
  abs_zero = 1e-8
  mesh_mode = 'REPLICATED'
  min_ad_size = 100
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation between two blocks in 3D and match reasonably well with the node-to-segment approach with a refined mesh."
  method = '!dbg'
  valgrind = 'none'
[]

[gap_heat_transfer_3d_mortar_hex20]
  input = 'gap_heat_transfer_3D_mortar.i'
  type = CSVDiff
  csvdiff = 'gap_heat_transfer_3D_mortar_hex20_out_NodalTemperature_0001.csv'
  cli_args = 'Mesh/left_block/elem_type=HEX20 Mesh/right_block/elem_type=HEX20 Outputs/file_base=gap_heat_transfer_3D_mortar_hex20_out'
  rel_err = 4e-2
  abs_zero = 1e-8
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to generate mortar numerical results for conduction and radiation between two blocks in 3D using HEX20 elements and match reasonably well with the node-to-segment approach with a refined mesh."
  valgrind = 'heavy'
  method = '!dbg'
[]

[gap_heat_transfer_sphere3d]
  input = 'gap_heat_transfer_sphere3D.i'
  type = CSVDiff
  csvdiff = 'gap_heat_transfer_sphere3D_out_NodalTemperature_0001.csv'
  rel_err = 1e-3
  abs_zero = 1e-8
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to generate node-to-segment numerical results using thermal contact that resolves spherical heat transfer for a constant conductivity gap in 3D using HEX20 elements"
[]

[gap_heat_transfer_sphere3d_mortar]
  input = 'gap_heat_transfer_sphere3D_mortar.i'
  type = CSVDiff
  csvdiff = 'gap_heat_transfer_sphere3D_mortar_out_NodalTemperature_0001.csv'
  rel_err = 4e-2
  abs_zero = 1e-8
  mesh_mode = 'REPLICATED'
  requirement = "We shall be able to generate mortar numerical results using thermal contact that resolves spherical heat transfer for a constant conductivity gap in 3D using HEX20 elements and match reasonably well with the node-to-segment approach with a refined mesh."
[]

[closed_gap_prescribed_pressure]
  type = CSVDiff
  input = 'closed_gap_prescribed_pressure.i'
  csvdiff = 'closed_gap_prescribed_pressure_out.csv'
  cli_args = "Mesh/allow_renumbering=false"
  detail = "and the material thermal conductivities and hardness values consistent with an analytical solution for the temperatures at the interface."
[]

[xyz]
  input = 'fv_modular_gap_heat_transfer_mortar_radiation_conduction.i'
  type = Exodiff
  exodiff = 'fv_modular_gap_heat_transfer_mortar_radiation_conduction_out.e'
  detail = 'in a Cartesian coordinate system'
[]

[rz]
  input = 'fv_modular_gap_heat_transfer_mortar_radiation_conduction.i'
  type = Exodiff
  exodiff = 'fv_modular_gap_heat_transfer_mortar_radiation_conduction_rz_out.e'
  cli_args = "Outputs/file_base=fv_modular_gap_heat_transfer_mortar_radiation_conduction_rz_out Mesh/coord_type=RZ Mesh/inactive=''"
  detail = 'in an axisymmetric coordinate system'
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_verbose]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_verbose.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_verbose_out_NodalTemperature_0001.csv'
  requirement = 'We shall be able to leverage mortar constraint and user objects to describe gap heat transfer physics by spelling out those objects in the input file.'
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_out_NodalTemperature_0001.csv'
  requirement = 'We shall be able to leverage mortar constraints and user objects to describe gap heat transfer physics by using the mortar thermal action in MOOSE.'
[]

[modular_gap_heat_transfer_mortar_displaced_conduction_UOs_function]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_conduction_UOs_function.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_conduction_UOs_function_out_NodalTemperature_0001.csv'
  requirement = 'We shall be able to leverage mortar constraints and user objects to describe gap heat transfer physics by using the mortar thermal action in MOOSE and describe the gap conductance with a function of temperature.'
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_lowerd_exists]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_lowerd_exists.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_lowerd_exists_out_NodalTemperature_0001.csv'
  requirement = 'We shall be able to leverage mortar constraints and user objects to describe gap heat transfer physics by using the mortar thermal action in MOOSE when the lower-dimensional domains have already been appended to the mesh.'
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_existing_UOs]
  type = CSVDiff
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_existing_UOs.i'
  csvdiff = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_existing_UOs_out_NodalTemperature_0001.csv'
  requirement = 'We shall be able to leverage mortar constraints and user objects to describe gap heat transfer physics by using the mortar thermal action in MOOSE when the user objects are manually built by the user in the input file.'
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_existing_UOs_error]
  type = RunException
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_existing_UOs.i'
  cli_args = "MortarGapHeatTransfer/mortar_heat_transfer/primary_emissivity=1.0"
  expect_err = 'The mortar gap heat transfer action requires that the input file defines user objects'
  requirement = 'We shall be able to inform the user that he or she provided physics parameters for two ways of building gap heat transfer options and error out, to avoid having misleading input files.'
[]

[modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action_too_many_uos]
  type = RunException
  input = 'modular_gap_heat_transfer_mortar_displaced_radiation_conduction_action.i'
  cli_args = "MortarGapHeatTransfer/mortar_heat_transfer/gap_flux_options='CONDUCTION RADIATION CONDUCTION'"
  expect_err = 'You cannot choose to have more than one conduction or more than one radiation user'
  requirement = 'We shall be able to inform the user that the mortar gap heat transfer action cannot internally build more than one conduction or more than one radiation user object.'
[]

[test]
  type = 'Exodiff'
  input = 'gap_heat_transfer_radiation_test.i'
  exodiff = 'gap_heat_transfer_radiation_test_out.e'
  abs_zero = 1e-06
  requirement = 'The system shall be able to compute radiative heat flux across a gap using the ThermalContact methods.'
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#1609'
[]

[cylinder]
  type = CSVDiff
  input = cylinder.i
  csvdiff = cylinder_out.csv
  override_columns = 'error_1 error_2'
  override_abs_zero = '1e-9 1e-9'
  override_rel_err = '5e-3 5e-3'
  requirement = 'The system shall be able to compute radiative heat flux across a cylindrical gap using the ThermalContact methods.'
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#26627'
[]

[sphere]
  type = CSVDiff
  input = sphere.i
  csvdiff = sphere_out.csv
  override_columns = 'error_1 error_2'
  override_abs_zero = '1e-9 1e-9'
  override_rel_err = '5e-3 5e-3'
  requirement = 'The system shall be able to compute radiative heat flux across a spherical gap using the ThermalContact methods.'
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#26627'
[]

[test]
  type = CSVDiff
  input = perfect_transfer_gap.i
  csvdiff = perfect_transfer_gap_out.csv
  detail = 'the use of a penalty parameter, and'
[]

[error_check]
  type = RunException
  input = perfect_transfer_gap.i
  cli_args = "ThermalContact/thermal_contact_1/type=GapHeatTransfer"
  expect_err = "This parameter should only be set by the user when 'type=GapPerfectConductance'."
  detail = 'error if the penalty is set for other heat transfer types.'
[]

[./generate_radiation_patch]
  type = 'Exodiff'
  input = 'generate_radiation_patch.i'
  exodiff = 'generate_radiation_patch_in.e'
  requirement = "The system shall be able to divide a sideset into patches for more accurate radiative transfer modeling."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = '--mesh-only'
  mesh_mode = REPLICATED
  issues = "#14000"
  recover = false
[../]

[./generate_radiation_patch_linear]
  type = 'Exodiff'
  input = 'generate_radiation_patch.i'
  exodiff = 'generate_radiation_patch_linear.e'
  requirement = "The system shall be able to use linear partitioner for subdividing sidesets into patches."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = 'Mesh/patch/partitioner=linear --mesh-only generate_radiation_patch_linear.e'
  mesh_mode = REPLICATED
  issues = "#14000"
  recover = false
[../]

[./generate_radiation_patch_centroid]
  type = 'Exodiff'
  input = 'generate_radiation_patch.i'
  exodiff = 'generate_radiation_patch_centroid.e'
  requirement = "The system shall be able to use centroid partitioner for subdividing sidesets into patches."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = 'Mesh/patch/partitioner=centroid Mesh/patch/centroid_partitioner_direction=x --mesh-only generate_radiation_patch_centroid.e'
  mesh_mode = REPLICATED
  issues = "#14000"
  recover = false
[../]

[./generate_radiation_patch_centroid_error]
  type = RunException
  input = 'generate_radiation_patch.i'
  expect_err = "If using the centroid partitioner you _must_ specify centroid_partitioner_direction!"
  requirement = "The system shall report an error when centroid partitioner is used but centroid_partitioner_direction is not provided."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = 'Mesh/patch/partitioner=centroid --mesh-only generate_radiation_patch_centroid.e'
  mesh_mode = REPLICATED
  issues = "#14000"
  recover = false
[../]

[./generate_radiation_patch_grid]
  type = 'Exodiff'
  input = 'generate_radiation_patch.i'
  exodiff = 'generate_radiation_patch_grid.e'
  requirement = "The system shall be able to use a uniform grid for subdividing sidesets into patches."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = 'Mesh/patch/partitioner=grid --mesh-only generate_radiation_patch_grid.e'
  mesh_mode = REPLICATED
  issues = "#15829"
  recover = false
[../]

[./generate_radiation_patch_grid_2D]
  type = 'Exodiff'
  input = 'generate_radiation_patch_grid_2D.i'
  exodiff = 'generate_radiation_patch_grid_2D_in.e'
  requirement = "The system shall be able to use a uniform grid for subdividing 1D sidesets into patches."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = 'Mesh/patch/partitioner=grid --mesh-only'
  mesh_mode = REPLICATED
  issues = "#15829"
  recover = false
[../]

[./generate_radiation_patch_grid_2D_overpart]
  type = 'Exodiff'
  input = 'generate_radiation_patch_grid_2D.i'
  exodiff = 'generate_radiation_patch_grid_2D_overpart.e'
  requirement = "The system shall be able to adjust the number of patches of partitions that end up empty."
  design = 'source/meshgenerators/PatchSidesetGenerator.md'
  cli_args = 'Mesh/patch/n_patches=61 Mesh/patch/partitioner=grid --mesh-only generate_radiation_patch_grid_2D_overpart.e'
  mesh_mode = REPLICATED
  issues = "#15829"
  recover = false
[../]

[./inconsistent_bnd_eps]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/emissivity="1"'
  expect_err = 'The number of entries must match the number of boundary entries.'
  requirement = 'The system shall check consistency of boundary and emissivity entries.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./inconsistent_bnd_view_factors]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/boundary="top bottom"
                UserObjects/gray_lambert/adiabatic_boundary="bottom"
                UserObjects/gray_lambert/fixed_temperature_boundary="top"
                UserObjects/gray_lambert/fixed_boundary_temperatures="550"
                UserObjects/gray_lambert/emissivity="1 1"'
  expect_err = 'Leading dimension of view_factors must be equal to number of side sets.'
  requirement = 'The system shall check consistency of boundary and view factor entries.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./inconsistent_iso_temperature]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/boundary="top bottom"
                UserObjects/gray_lambert/fixed_temperature_boundary="top"
                UserObjects/gray_lambert/emissivity="1 1"
                UserObjects/gray_lambert/view_factors="0.5 0.5 ; 0.5 0.5"'
  expect_err = 'fixed_boundary_temperatures and fixed_temperature_boundary must have the same length.'
  requirement = 'The system shall check consistency of fixed_boundary_temperatures and fixed_temperature_boundary entries.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./inconsistent_bnd_iso_bnd]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/boundary="top bottom"
                UserObjects/gray_lambert/fixed_temperature_boundary="right"
                UserObjects/gray_lambert/fixed_boundary_temperatures="300"
                UserObjects/gray_lambert/emissivity="1 1"
                UserObjects/gray_lambert/view_factors="0.5 0.5 ; 0.5 0.5"'
  expect_err = 'fixed_temperature_boundary must be a subset of boundary.'
  requirement = 'The system shall check consistency of boundary and fixed_temperature_boundary entries.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./inconsistent_bnd_adiabatic_bnd]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/boundary="top bottom"
                UserObjects/gray_lambert/fixed_temperature_boundary="top"
                UserObjects/gray_lambert/fixed_boundary_temperatures="300"
                UserObjects/gray_lambert/emissivity="1 1"
                UserObjects/gray_lambert/view_factors="0.5 0.5 ; 0.5 0.5"'
  expect_err = 'adiabatic_boundary must be a subset of boundary.'
  requirement = 'The system shall check consistency of boundary and adiabatic_boundary entries.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./incorrect_view_factor_shape]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/view_factors="0.25 0.25 0.25 0.25; 0.5 0.5 ; 1 1 1 1; 1 1 1 1"'
  expect_err = 'view_factors must be provided as square array. Row 1 has 2 entries.'
  requirement = 'The system shall check consistency of the view_factors entry shape.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./bad_rowsum]
  type = 'RunException'
  input = 'gray_lambert_cavity.i'
  cli_args = 'UserObjects/gray_lambert/view_factors="0.25 0.25 0.25 0.25; 0.2 0.25 0.25 0.25 ; 0.25 0.25 0.25 0.25; 0.25 0.25 0.25 0.25"'
  expect_err = 'view_factors row 1 sums to 0.95.'
  compiler = '!GCC'
  requirement = 'The system shall check consistency of the view_factors entry norm.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./gray_lambert_cavity]
  type = CSVDiff
  input = 'gray_lambert_cavity.i'
  csvdiff = 'gray_lambert_cavity_out.csv gray_lambert_cavity_out_lambert_vpp_0001.csv gray_lambert_cavity_out_view_factors_0001.csv'
  cli_args = 'Outputs/csv=true'
  requirement = 'The system shall compute radiative transfer between gray Lambert surfaces.'
  design = 'ConstantViewFactorSurfaceRadiation.md SurfaceRadiationVectorPostprocessor.md ViewfactorVectorPostprocessor.md'
[../]

[./coupled_heat_conduction]
  type = Exodiff
  input = 'coupled_heat_conduction.i'
  exodiff = 'coupled_heat_conduction_out.e'
  requirement = 'The system shall allow coupling radiative transfer between gray Lambert surfaces to solving heat conduction.'
  design = 'ConstantViewFactorSurfaceRadiation.md'
[../]

[./coupled_heat_conduction_emission_reconstruction]
  type = Exodiff
  input = 'coupled_heat_conduction.i'
  exodiff = 'coupled_heat_conduction_emission_er.e'
  cli_args = 'BCs/radiation/reconstruct_emission=true Outputs/file_base=coupled_heat_conduction_emission_er'
  requirement = 'The system shall allow reconstructing the spatial distribution of the emission component on a radiation boundary via the T4 law.'
  design = 'GrayLambertNeumannBC.md'
[../]

[./gray_lambert_cavity_automatic_vf]
  type = CSVDiff
  input = 'gray_lambert_cavity_automatic_vf.i'
  csvdiff = 'gray_lambert_cavity_automatic_vf_out.csv'
  requirement = 'The system shall compute radiative transfer between gray Lambert surfaces when the view factors are provided by a userobject.'
  design = 'ViewFactorObjectSurfaceRadiation.md'
  mesh_mode = REPLICATED
[../]

[./gray_lambert_cavity_automatic_vf_3D]
  type = CSVDiff
  input = 'gray_lambert_cavity_automatic_vf_3D.i'
  csvdiff = 'gray_lambert_cavity_automatic_vf_3D_out.csv'
  requirement = 'The system shall compute radiative transfer between gray Lambert surfaces in 3D when the view factors are provided by a userobject.'
  design = 'ViewFactorObjectSurfaceRadiation.md'
  mesh_mode = REPLICATED
[../]

[./perfect]
  type = 'Exodiff'
  input = 'perfect.i'
  exodiff = 'perfect_out.e'
  requirement = "The system shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures (QUAD4)."
  design = "/materials/GapConductance.md"
  issues = "#6750"
[../]

[./perfectQ8]
  type = 'Exodiff'
  input = 'perfectQ8.i'
  exodiff = 'perfectQ8_out.e'
  requirement = "The system shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures (QUAD8)"
  design = "/materials/GapConductance.md"
  issues = "#6750"
[../]

[./perfectQ9]
  type = 'Exodiff'
  input = 'perfectQ9.i'
  exodiff = 'perfectQ9_out.e'
  requirement = "The system shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures (QUAD9)"
  design = "/materials/GapConductance.md"
  issues = "#6750"
[../]

[./nonmatching]
  type = 'Exodiff'
  input = 'nonmatching.i'
  exodiff = 'nonmatching_out.e'
  requirement = "The system shall compute the heat transfer across small gaps for non-matching meshes."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./second_order]
  type = 'Exodiff'
  input = 'second_order.i'
  exodiff = 'second_order_out.e'
  requirement = "The system shall compute the heat transfer across small gaps for second order FEM bases."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./moving]
  type = 'Exodiff'
  input = 'moving.i'
  exodiff = 'moving_out.e'
  allow_warnings = true
  requirement = "The system shall compute the heat transfer across small gaps for moving interfaces."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./gap_conductivity_property]
  type = 'Exodiff'
  input = 'gap_conductivity_property.i'
  exodiff = 'gap_conductivity_property_out.e'
  requirement = "The system shall compute the heat transfer across small gaps with a specified gap conductivity."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./gap_conductivity_property_r1_error]
  type = RunException
  input = 'gap_conductivity_property.i'
  cli_args = 'Mesh/uniform_refine=1'
  expect_err = 'GapConductance does not work with uniform mesh refinement.'
  requirement = "The system shall throw an error if the gap conductance model is used
                   with uniform mesh refinement"
  design = "/actions/ThermalContactAction.md"
  issues = "#13043"
[../]

[perfect_prereq]
  type = 'CheckFiles'
  prereq = perfect
  input = 'perfect.i'
  cli_args = '--split-mesh=2'
  mesh_mode = 'REPLICATED'
  check_files = 'perfect.cpa.gz/2/header.gz perfect.cpa.gz/2/split-2-0.gz perfect.cpa.gz/2/split-2-1.gz'
  recover = false
  requirement = "The system shall generate a parallel mesh split across 2 processes."
  design = "Mesh/splitting.md"
  issues = "#27203"
[]

[perfect_split_mesh]
  type = 'Exodiff'
  prereq = perfect_prereq
  input = 'perfect_split.i'
  exodiff = 'perfect_out.e'
  cli_args = 'Outputs/file_base=perfect_out'
  requirement = "The system shall read in parallel a mesh split across 2 processes."
  design = "Mesh/splitting.md"
  min_parallel = 2
  max_parallel = 2
  issues = "#27203"
[]

[./nonmatching]
  type = 'Exodiff'
  input = 'nonmatching.i'
  exodiff = 'nonmatching_out.e'
  requirement = 'The system shall support thermal contact with linear 3d hexahedral elements'
[../]

[./second]
  type = 'Exodiff'
  input = 'second.i'
  exodiff = 'second_out.e'
  requirement = 'The system shall support thermal contact with second-order 3d hexahedral elements'
[../]

[./moving]
  type = 'Exodiff'
  input = 'moving.i'
  exodiff = 'moving_out.e'
  valgrind = 'HEAVY'
  allow_warnings = true
  requirement = 'The system shall support thermal contact with 3d hexahedral elements where the surfaces move relative to one another'
[../]

[./const_hw]
  type = 'Exodiff'
  input = 'const_hw.i'
  exodiff = 'const_hw_out.e'
  requirement = 'The system shall provide convective heat flux boundary condition where far-field temperature and convective heat transfer coefficient are given as constant variables'
[../]

[./coupled_convective_heat_flux]
  type = 'Exodiff'
  input = 'coupled_convective_heat_flux.i'
  exodiff = 'coupled_convective_heat_flux_out.e'
  requirement = 'The system shall provide convective heat flux boundary condition where far-field temperature and convective heat transfer coefficient are given as spatially varying variables'
[../]

[./coupled_convective_heat_flux_two_phase]
  type = 'Exodiff'
  input = 'coupled_convective_heat_flux_two_phase.i'
  exodiff = 'coupled_convective_heat_flux_two_phase_out.e'
  requirement = 'The system shall provide convective heat flux boundary condition for multi-phase fluids where far-field temperatures and convective heat transfer coefficients are given as spatially varying variables'
[../]

[./not_enough_alpha]
  type = RunException
  input = 'coupled_convective_heat_flux_two_phase.i'
  cli_args = "BCs/right/alpha='alpha_liquid'"
  expect_err = '\(BCs/right/alpha\).*?The number of coupled components does not match the number of `T_infinity` components.'
  requirement = 'The system shall report an error if the number of `alpha` components does not match the number of `T_infinity` components.'
[../]

[./not_enough_htc]
  type = RunException
  input = 'coupled_convective_heat_flux_two_phase.i'
  cli_args = "BCs/right/htc='Hw_liquid'"
  expect_err = '\(BCs/right/htc\).*?The number of coupled components does not match the number of `T_infinity` components.'
  requirement = 'The system shall report an error if the number of `htc` components does not match the number of `T_infinity` components.'
[../]

[./on_off]
  type = 'Exodiff'
  input = 'on_off.i'
  exodiff = 'on_off_out.e'
  requirement = 'The system shall enable scaling of the total heat flux of the convective heat flux boundary condition'
  issues = '#15421'
[../]

[./min_gap_order_zero]
  type = CSVDiff
  csvdiff = min_gap_order_zero.csv
  input = min_gap.i
  cli_args = 'ThermalContact/left_to_right/min_gap_order=0 Outputs/file_base=min_gap_order_zero'
  issues = '#13221'
  design = 'GapConductance.md GapHeatTransfer.md'
  mesh_mode = REPLICATED
  requirement = 'Optionally a constant attenuation shall be applied to compute the gap conductance below a gap length threshold.'
[../]

[./min_gap_order_one]
  type = CSVDiff
  csvdiff = min_gap_order_one.csv
  input = min_gap.i
  cli_args = 'ThermalContact/left_to_right/min_gap_order=1 Outputs/file_base=min_gap_order_one'
  issues = '#13221'
  mesh_mode = REPLICATED
  design = 'GapConductance.md GapHeatTransfer.md'
  requirement = 'Optionally a linear Taylor expansion of the inverse gap length shall be applied as the attenuation to compute the gap conductance below a gap length threshold.'
[../]

[./test]
  type = 'Exodiff'
  input = 'heat_conduction_ortho.i'
  exodiff = 'heat_conduction_ortho_out.e'
  max_parallel = 1
  design = 'AnisoHeatConductionMaterial.md'
  issues = '#2674'
  requirement = 'The system shall allow the use of an anisotropic heat conduction material set by postprocessors.'
[../]

[./test]
  type = 'Exodiff'
  input = 'heat_conduction_patch.i'
  exodiff = 'heat_conduction_patch_out.e'
  max_parallel = 1
  requirement = 'The system shall compute a tri-linear temperature field'
[../]

[./test_rz_quad8]
  type = 'Exodiff'
  input = 'heat_conduction_patch_rz_quad8.i'
  exodiff = 'heat_conduction_patch_rz_quad8_out.e'
  max_parallel = 1
  requirement = 'The system shall compute a bi-linear temperature field for an axisymmetric problem with quad8 elements'
[../]

[./test_rz]
  type = 'Exodiff'
  input = 'heat_conduction_patch_rz.i'
  exodiff = 'heat_conduction_patch_rz_out.e'
  max_parallel = 1
  requirement = 'The system shall compute a bi-linear temperature field for an axisymmetric problem'
[../]

[./test_hex20]
  type = 'Exodiff'
  input = 'heat_conduction_patch_hex20.i'
  exodiff = 'heat_conduction_patch_hex20_out.e'
  max_parallel = 1
  requirement = 'The system shall compute a tri-linear temperature field with hex20 elements'
[../]

[./test_hex20_aniso]
  type = 'Exodiff'
  input = 'heat_conduction_patch_hex20_aniso.i'
  exodiff = 'heat_conduction_patch_hex20_out.e'
  max_parallel = 1
  prereq = test_hex20
  requirement = 'The system shall compute a tri-linear temperature field with hex20 elements using an anisotropic thermal conductivity model with isotropic thermal conductivities supplied'
[../]

[./heat_source_bar]
  type = 'Exodiff'
  input = 'heat_source_bar.i'
  exodiff = 'heat_source_bar_out.e'
  abs_zero = 1e-7
  requirement = 'The system shall reproduce an analytical solution of a heat source in a 1D ceramic bar.'
  design = "HeatSource.md"
  issues = "#2582"
[../]

[./ad_heat_source_bar]
  type = 'Exodiff'
  input = 'ad_heat_source_bar.i'
  exodiff = 'ad_heat_source_bar_out.e'
  abs_zero = 1e-7
  requirement = 'The system shall reproduce an analytical solution of a heat source in a 1D ceramic bar using AD kernels.'
  design = "ADMatHeatSource.md"
  issues = "#12633"
  custom_cmp = 'ad_heat_source_bar.cmp'
[../]

[./ad_heat_source_bar_jacobian]
  type = 'PetscJacobianTester'
  input = 'ad_heat_source_bar.i'
  ratio_tol = 2e-7
  difference_tol = 2e-2
  run_sim = True
  requirement = 'The system shall produce correct Jacobians for the AD heat conduction and heat source kernel objects.'
  design = "ADMatHeatSource.md"
  issues = "#5658 #12633"
  prereq = 'ad_heat_source_bar' # same input file
[../]

[heatConduction_test]
  type = 'Exodiff'
  input = 'heatConduction2D.i'
  exodiff = 'heatConduction2D_out.e'
  max_parallel = 1
  design = 'HomogenizedHeatConduction.md'
  requirement = 'The system shall compute homogenized thermal conductivity using the asymptotic expansion homogenization approach'
  issues = '#6750'
[]

[homogenize_tc_hexagonal]
  type = 'CSVDiff'
  input = 'homogenize_tc_hex.i'
  csvdiff = 'homogenize_tc_hex_outp_csv.csv'
  requirement = 'The system shall compute homogenized thermal conductivity using the asymptotic expansion homogenization approach for hexagonal geometries.'
  design = 'HomogenizedThermalConductivity.md'
  issues = '#22919'
[]

[homogenize_tensor_tc]
  type = 'Exodiff'
  input = 'heatConduction2D_tensor_tc.i'
  exodiff = 'heatConduction2D_tensor_tc_out.e'
  design = 'AnisoHomogenizedHeatConduction.md'
  requirement = 'The system shall compute homogenized thermal conductivity using the asymptotic expansion homogenization approach when the heterogeneous thermal conductivities are given as tensor'
  issues = '#22919'
[]

[homogenize_aniso_tensor_tc]
  type = 'Exodiff'
  input = 'heatConduction2D_tensor_tc.i'
  exodiff = 'homogenize_aniso_tensor_tc.e'
  design = 'AnisoHomogenizedHeatConduction.md'
  cli_args = 'Outputs/file_base=homogenize_aniso_tensor_tc Materials/heat2/thermal_conductivity="100 1 0 0 80 8 0 0 0"'
  requirement = 'The system shall compute homogenized thermal conductivity using the asymptotic expansion homogenization approach when the heterogeneous thermal conductivities are given as anisotropic tensor'
  issues = '#22919'
[]

[constraint_joule_heating_single_material_insulated]
  type = CSVDiff
  input = "constraint_joule_heating_single_material_insulated.i"
  csvdiff = "constraint_joule_heating_single_material_insulated_out.csv"
  cli_args = "Mesh/allow_renumbering=false"
  detail = " between two separate blocks of the same material constistent with the analytical solution (850.595K) in the case of thermal insulation across the interface."
[]

[constraint_joule_heating_single_material]
  type = CSVDiff
  input = "constraint_joule_heating_single_material.i"
  csvdiff = "constraint_joule_heating_single_material_out.csv"
  cli_args = "Mesh/allow_renumbering=false"
  detail = " between two separate blocks of the same material as equivalent to the steady state case when thermal contact across the interface is not considered."
[]

[constraint_joule_heating_offset_single_material_insulated]
  type = CSVDiff
  input = "constraint_joule_heating_offset_single_material_insulated.i"
  csvdiff = "constraint_joule_heating_offset_single_material_insulated_out.csv"
  cli_args = "Mesh/allow_renumbering=false"
  detail = "as non-zero only along the portion of the interface at which the gap between the two blocks is closed and as zero along the open-gap portion of the interface."
[]

[constraint_joule_heating_dual_material_insulated]
  type = CSVDiff
  input = "constraint_joule_heating_dual_material_insulated.i"
  csvdiff = "constraint_joule_heating_dual_material_insulated_out.csv"
  cli_args = "Mesh/allow_renumbering=false"
  detail = " that, when the blocks at the interface have different material properties, is different in each material while the Joule heating source into each block is equivalent from the common electric potential drop across the interface."
[]

[constraint_joule_heating_dual_material]
  type = CSVDiff
  input = "constraint_joule_heating_dual_material.i"
  csvdiff = "constraint_joule_heating_dual_material_out.csv"
  cli_args = "Mesh/allow_renumbering=false"
  detail = "as a function of both the Joule heating and thermal contact at the interface, which results in a smaller temperature difference between the two blocks at the interface."
[]

[transient_joule_heating_constraint]
  type = CSVDiff
  input = "transient_joule_heating_constraint.i"
  csvdiff = "transient_joule_heating_constraint_out.csv"
  cli_args = "Mesh/allow_renumbering=false"
  requirement = "We shall be able to converge towards the steady-state interface temperature profile due to Joule heating and thermal contact in a transient simulation, given sufficient simulation time."
[]

[./joule_heating]
  type = 'Exodiff'
  input = 'transient_jouleheating.i'
  exodiff = 'transient_jouleheating_out.e'
  requirement = 'The system shall compute Joule heating'
  issues = '#8220'
  design = 'JouleHeatingSource.md'
[../]

[./ad_joule_heating]
  type = 'Exodiff'
  input = 'transient_ad_jouleheating.i'
  exodiff = 'transient_ad_jouleheating_out.e'
  requirement = 'The system shall compute Joule heating using automatic differentiation'
  issues = '#15536'
  design = 'ADJouleHeatingSource.md'
[../]

[./ad_joule_heating_jacobian]
  type = 'PetscJacobianTester'
  input = 'transient_ad_jouleheating.i'
  cli_args = 'Executioner/end_time=1'
  ratio_tol = 3.2e-6
  difference_tol = 0.002
  requirement = 'The system shall compute a perfect jacobian for Joule heating using automatic differentiation'
  issues = '#15536'
  design = 'ADJouleHeatingSource.md'
[../]

[test]
  type = Exodiff
  input = test.i
  exodiff = test_out.e
  requirement = 'The system shall be able to model an equilibrium between an incoming heat flux from a focused beam (e.g. laser), which is described by a Gaussian shape, and outgoing heat flux due to radiative losses.'
[]

[./test]
  type = 'Exodiff'
  input = 'meshed_gap_thermal_contact.i'
  exodiff = 'meshed_gap_thermal_contact_out.e'
  requirement = "The ThermalContact system shall enforce heat transfer
    across a meshed gap in a 2D plane geometry."
[../]

[./constant_conductance]
  type = 'Exodiff'
  input = 'meshed_gap_thermal_contact_constant_conductance.i'
  exodiff = 'meshed_gap_thermal_contact_constant_conductance_out.e'
  issues = '#13061'
  design = 'syntax/ThermalContact/index.md source/materials/GapConductanceConstant.md'
  requirement = "The ThermalContact system shall correctly enforce heat transfer
    across a meshed gap in a 2D plane geometry using a prescribed constant conductance."
[../]

[./constant_conductance_quadrature]
  prereq = constant_conductance
  type = 'Exodiff'
  input = 'meshed_gap_thermal_contact_constant_conductance.i'
  cli_args = 'ThermalContact/gap_conductance/quadrature=true Outputs/file_base=meshed_gap_thermal_contact_constant_conductance_quadrature_out'
  exodiff = 'meshed_gap_thermal_contact_constant_conductance_quadrature_out.e'
  issues = '#13061'
  design = 'syntax/ThermalContact/index.md source/materials/GapConductanceConstant.md'
  requirement = "The ThermalContact system shall correctly enforce heat transfer
    across a meshed gap in a 2D plane geometry using a prescribed constant conductance with the quadrature option"
[../]

[./annulus]
  type = 'Exodiff'
  input = 'meshed_annulus_thermal_contact.i'
  exodiff = 'meshed_annulus_thermal_contact_out.e'
  requirement = "The ThermalContact system shall enforce heat transfer
    across a meshed circular annulus in a 2D plane geometry."
[../]

[multiple_contact_pairs]
  type = CSVDiff
  input = 'multiple_contact_pairs.i'
  csvdiff = 'multiple_contact_pairs_out.csv'
  abs_zero = 1.0e-6
  rel_err = 1.0e-5
  installation_type = in_tree
  requirement = 'Heat transfer module action shall allow for providing multiple contact pairs.'
[]

[multiple_radiation_cavities]
  type = 'Exodiff'
  issues = '#16954'
  design = 'RadiationTransferAction.md'
  requirement = 'The system shall support the the modeling of radiative heat transfer with multiple radiation cavities.'
  input = 'multiple_radiation_cavities.i'
  exodiff = 'multiple_radiation_cavities_out.e'
  mesh_mode = REPLICATED
[]

[./test]
  type = 'Exodiff'
  input = 'parallel_element_pps_test.i'
  exodiff = 'out.e'
  min_parallel = 2
  platform = 'DARWIN'
  requirement = 'The system shall computed an integrated value on elements in parallel'
[../]

[default]
  type = CSVDiff
  input = test_fv.i
  csvdiff = 'default.csv'
  cli_args = 'Outputs/file_base=default'
  detail = 'using the default initial condition,'
[]

[user_ics]
  type = CSVDiff
  input = test_fv.i
  csvdiff = 'user_ics.csv'
  cli_args = 'Physics/HeatConduction/FiniteVolume/h1/initial_temperature=100 Outputs/file_base=user_ics'
  detail = 'with a user-defined initial condition,'
[]

[restart_with_user_ics]
  type = CSVDiff
  input = test_fv.i
  csvdiff = 'restart_user_ics.csv'
  cli_args = "Physics/HeatConduction/FiniteVolume/h1/initial_temperature=100
                  Problem/restart_file_base=default_cp/LATEST Problem/allow_initial_conditions_with_restart=true
                  Outputs/file_base=restart_user_ics"
  detail = 'when performing a regular checkpoint restart, but still obeying the user-defined initial condition,'
[]

[restart_from_file]
  type = CSVDiff
  input = test_fv.i
  csvdiff = 'from_file.csv'
  cli_args = "Mesh/active='fmg_restart' Physics/HeatConduction/FiniteVolume/h1/initialize_variables_from_mesh_file=true Outputs/file_base=from_file"
  detail = 'when performing manual restart from a mesh file, ignoring the default initial condition.'
[]

[error]
  type = RunException
  input = test_fv.i
  cli_args = 'Physics/HeatConduction/FiniteVolume/h1/initial_temperature=100 Physics/HeatConduction/FiniteVolume/h1/initialize_variables_from_mesh_file=true'
  expect_err = 'Initial temperature should not be set if the variables should be initialized from the mesh file'
  requirement = 'The system shall error if the user specifies initial conditions while also requesting variables be loaded from a mesh file.'
[]

[base]
  type = Exodiff
  input = 'test_cg.i'
  exodiff = 'test_cg_out.e'
  detail = 'with Dirichlet and Neumann boundary conditions and a heat source defined using a variable,'
[]

[functor_heat]
  type = Exodiff
  input = 'test_cg.i'
  exodiff = 'test_cg_functor_hs_out.e'
  cli_args = "AuxVariables/Q_0/initial_condition=0
                  Physics/HeatConduction/FiniteElement/h1/heat_source_var='Q_0'
                  Physics/HeatConduction/FiniteElement/h1/heat_source_functor='100'
                  Outputs/hide=Q_0
                  Outputs/file_base=test_cg_functor_hs_out"
  detail = 'with a heat source defined using a functor,'
[]

[convective_bc]
  type = Exodiff
  input = 'test_cg_convective_bc.i'
  exodiff = 'test_cg_convective_bc_out.e'
  detail = 'with convective heat flux boundary conditions.'
[]

[base]
  type = Exodiff
  input = 'test_fv.i'
  exodiff = 'test_fv_out.e'
  detail = 'with Dirichlet and Neumann boundary conditions and a heat source defined using a variable,'
[]

[functor_heat]
  type = Exodiff
  input = 'test_fv.i'
  exodiff = 'test_fv_functor_hs_out.e'
  cli_args = "Physics/HeatConduction/FiniteVolume/h1/heat_source_var='0'
                  Physics/HeatConduction/FiniteVolume/h1/heat_source_functor='100'
                  Outputs/file_base=test_fv_functor_hs_out"
  detail = 'with a heat source defined using a functor,'
[]

[convective_bc]
  type = Exodiff
  input = 'test_fv_convective_bc.i'
  exodiff = 'test_fv_convective_bc_out.e'
  detail = 'with convective heat flux boundary conditions.'
[]

[./convective_ht_side_integral]
  type = 'CSVDiff'
  input = 'convective_ht_side_integral.i'
  csvdiff = 'convective_ht_side_integral_out.csv'
  design = 'ConvectiveHeatTransferSideIntegral.md'
  requirement = 'The system shall compute total heat flux from heat transfer coefficient and temperature difference'
  issues = '#14390'
[../]

[./ad_convective_ht_side_integral]
  type = 'CSVDiff'
  input = 'ad_convective_ht_side_integral.i'
  csvdiff = 'ad_convective_ht_side_integral_out.csv'
  design = 'ConvectiveHeatTransferSideIntegral.md'
  requirement = 'The system shall compute total heat flux from heat transfer coefficient and temperature difference for AD variables'
  issues = '#14390 #17020'
[../]

[radiative_transfer_action_analytical]
  type = 'Exodiff'
  input = 'radiative_transfer_action.i'
  exodiff = 'radiative_transfer_action_out.e'
  requirement = 'The system shall provide an action to set up radiative heat transfer problems using the net radiation method for cavities with unobstructed, planar faces.'
  petsc_version = '>=3.9.4'
  mesh_mode = REPLICATED
  cli_args = 'GrayDiffuseRadiation/cavity/view_factor_calculator=analytical'
[]

[radiative_transfer_action_raytracing]
  type = 'Exodiff'
  input = 'radiative_transfer_action.i'
  exodiff = 'radiative_transfer_action_raytracing_out.e'
  requirement = 'The system shall provide an action to set up radiative heat transfer problems using the net radiation method and allow computing view factors using raytracing.'
  petsc_version = '>=3.9.4'
  mesh_mode = REPLICATED
  cli_args = 'Outputs/file_base=radiative_transfer_action_raytracing_out GrayDiffuseRadiation/cavity/view_factor_calculator=ray_tracing'
[]

[bnd_names]
  type = 'Exodiff'
  input = 'radiative_transfer_action.i'
  exodiff = 'radiative_transfer_action_out_bnd_names.e'
  cli_args = 'GrayDiffuseRadiation/cavity/boundary="inner_bottom 5 inner_right inner_top"
                GrayDiffuseRadiation/cavity/adiabatic_boundary="inner_top"
                GrayDiffuseRadiation/cavity/fixed_temperature_boundary="inner_bottom"
                Outputs/file_base=radiative_transfer_action_out_bnd_names'
  requirement = 'The system shall allow the specification of boundary names and ids in the modeling of radiative heat transfer.'
  petsc_version = '>=3.9.4'
  mesh_mode = REPLICATED
[]

[no_action]
  type = 'Exodiff'
  input = 'radiative_transfer_no_action.i'
  exodiff = 'radiative_transfer_no_action_out.e'
  requirement = 'The system shall ensure that results between manually created radiative transfer inputs and inputs that use the radiation transfer action are identical.'
  petsc_version = '>=3.9.4'
  mesh_mode = REPLICATED
[]

[external_boundary_analytical]
  type = 'Exodiff'
  input = 'radiative_transfer_action_external_boundary.i'
  exodiff = 'radiative_transfer_action_external_boundary_out.e'
  petsc_version = '>=3.9.4'
  mesh_mode = REPLICATED
  requirement = 'The system shall provide an action to set up radiative heat transfer problems where sidesets participating in the radiative exchange are external faces of the domain, with view factors computed by simple quadrature rules for cavities with unobstructed, planar faces.'
[]

[external_boundary_ray_tracing]
  type = 'Exodiff'
  input = 'radiative_transfer_action_external_boundary_ray_tracing.i'
  exodiff = 'radiative_transfer_action_external_boundary_ray_tracing_out.e'
  petsc_version = '>=3.9.4'
  mesh_mode = REPLICATED
  requirement = 'The system shall provide an action to set up radiative heat transfer problems where sidesets participating in the radiative exchange are external faces of the domain, with view factors computed by ray tracing.'
[]

[unfolded]
  type = 'Exodiff'
  input = 'cavity_with_pillars.i'
  exodiff = 'cavity_with_pillars_out.e'
  mesh_mode = REPLICATED
  petsc_version = '>=3.9.4'

  detail = 'unfolding the problem at the symmetry boundary and'
[]

[symmetry_bc]
  type = 'Exodiff'
  input = 'cavity_with_pillars_symmetry_bc.i'
  exodiff = 'cavity_with_pillars_symmetry_bc_out.e'
  mesh_mode = REPLICATED
  petsc_version = '>=3.9.4'

  detail = 'by using a symmetry boundary condition.'
[]

[./radiative_bc_cyl]
  type = CSVDiff
  input = 'radiative_bc_cyl.i'
  csvdiff = 'radiative_bc_cyl_out.csv'
  requirement = "The system shall be able to model radiative transfer from a cylindrical surface as a boundary condition."
  design = 'source/bcs/InfiniteCylinderRadiativeBC.md'
  issues = "#13053"
  # see #21185
  allow_warnings = true
[../]

[ad_radiative_bc_cyl]
  type = CSVDiff
  input = 'ad_radiative_bc_cyl.i'
  csvdiff = 'radiative_bc_cyl_out.csv'
  requirement = "The system shall be able to model radiative transfer from a cylindrical surface as boundary condition with automated differentiation."
  design = 'source/bcs/ADInfiniteCylinderRadiativeBC.md'
  issues = "#13053"
  cli_args = "Outputs/file_base=radiative_bc_cyl_out"
  # see #21185
  allow_warnings = true
[]

[ad_radiative_bc_cyl_jacobian]
  type = PetscJacobianTester
  input = 'ad_radiative_bc_cyl.i'
  ratio_tol = 1e-7
  difference_tol = 2
  requirement = "The system shall be able to model radiative transfer from a cylindrical surface as boundary condition with automated differentiation and provide exact Jacobian."
  design = 'source/bcs/ADInfiniteCylinderRadiativeBC.md'
  issues = "#13053"
  # see #21185
  allow_warnings = true
[]

[./function_radiative_bc]
  type = Exodiff
  input = 'function_radiative_bc.i'
  exodiff = 'function_radiative_bc_out.e'
  requirement = "The system shall be able to model radiative heat transfer using a user-specified emissivity function."
  design = 'source/bcs/FunctionRadiativeBC.md'
  issues = "#13053"
  mesh_mode = REPLICATED
[../]

[ad_function_radiative_bc]
  type = Exodiff
  input = 'ad_function_radiative_bc.i'
  exodiff = 'function_radiative_bc_out.e'
  requirement = "The system shall be able to model radiative heat transfer using a user-specified emissivity function with automated differentiation."
  design = 'source/bcs/ADFunctionRadiativeBC.md'
  issues = "#13053"
  mesh_mode = REPLICATED
  cli_args = "Outputs/file_base=function_radiative_bc_out"
[]

[ad_function_radiative_bc_jacobian]
  type = PetscJacobianTester
  input = 'ad_function_radiative_bc.i'
  ratio_tol = 1e-7
  difference_tol = 5e-5
  cli_args = "Executioner/num_steps=1"
  requirement = "The system shall be able to model radiative heat transfer using a user-specified emissivity function with automated differentiation and provide exact Jacobian."
  design = 'source/bcs/ADFunctionRadiativeBC.md'
  issues = "#13053"
  mesh_mode = REPLICATED
[]

[recover_1]
  type = Exodiff
  input = recover.i
  exodiff = recover_out.e
  superlu = true
  requirement = 'The system shall run a simulation with heat conduction, a heat source, thermal contact, and boundary conditions.'
[]

[recover_2]
  type = RunApp
  input = recover.i
  prereq = recover_1
  cli_args = 'Outputs/checkpoint=true Executioner/num_steps=5'
  superlu = true
  requirement = 'The system shall run a short simulation with heat conduction, a heat source, thermal contact, and boundary conditions.'
[]

[recover_3]
  type = Exodiff
  input = recover.i
  exodiff = recover_out.e
  prereq = recover_2
  delete_output_before_running = false
  cli_args = '--recover recover_out_cp/0005'
  superlu = true
  requirement = 'The system shall be able to recover from a short simulation and reproduce a the full time scale simulation with heat conduction, a heat source, thermal contact, and boundary conditions.'
[]

[ad_recover_1]
  type = Exodiff
  input = ad_recover.i
  exodiff = ad_recover_out.e
  superlu = true
  requirement = 'The system shall run a simulation with heat conduction, a heat source, thermal contact, and boundary conditions with automatic differentiation.'
[]

[ad_recover_2]
  type = RunApp
  input = ad_recover.i
  prereq = ad_recover_1
  cli_args = 'Outputs/checkpoint=true Executioner/num_steps=5'
  superlu = true
  requirement = 'The system shall run a short simulation with heat conduction, a heat source, thermal contact, and boundary conditions with automatic differentiation.'
[]

[ad_recover_3]
  type = Exodiff
  input = ad_recover.i
  exodiff = ad_recover_out.e
  prereq = ad_recover_2
  delete_output_before_running = false
  cli_args = '--recover ad_recover_out_cp/0005'
  superlu = true
  requirement = 'The system shall be able to recover from a short simulation and reproduce a the full time scale simulation with heat conduction, a heat source, thermal contact, and boundary conditions with automatic differentiation.'
[]

[./test]
  type = 'Exodiff'
  input = 'steinhart-hart_linear.i'
  exodiff = 'steinhart-hart_linear_out.e'
  requirement = 'The system shall compute conductivity of semiconductors according to the Steinhart-Hart equation'
[../]

[./1D_gap]
  type = 'Exodiff'
  input = 'gap_thermal_1D.i'
  exodiff = 'gap_thermal_1D_out.e'
  abs_zero = 1e-6
  detail = 'discontinuous finite elements, '
[../]

[./1D_gap_Tbulk_var]
  type = 'Exodiff'
  input = 'gap_thermal_1D.i'
  exodiff = 'gap_thermal_1D_Tbulk_var.e'
  cli_args = 'InterfaceKernels/active="gap_var"
                  Outputs/file_base=gap_thermal_1D_Tbulk_var'
  abs_zero = 1e-6
  detail = 'bulk gap temperature as an auxiliary variable, '
[../]

[./1D_gap_ktemp]
  type = 'Exodiff'
  input = 'gap_thermal_ktemp_1D.i'
  exodiff = 'gap_thermal_ktemp_1D_out.e'
  abs_zero = 1e-6
  detail = 'temperature dependent gap conductivity, and'
[../]

[./CFEM_gap]
  type = 'Exodiff'
  input = 'cfem_gap.i'
  exodiff = 'cfem_gap_out.e'
  abs_zero = 1e-6
  detail = 'block restricted continuous finite element variables.'
[../]

[./1D_gap_err]
  type = RunException
  input = 'gap_thermal_1D.i'
  cli_args = 'DGKernels/dg_diff/exclude_boundary="fake_interface another_fake last_one"'
  expect_err = "DGKernel 'dg_diff' contains the following boundary names that do not exist on the mesh: fake_interface, another_fake, last_one"
  requirement = 'The system shall throw an error if the specified boundary does not exist when using the DGDiffusion DGKernel.'
[../]

[2d]
  type = Exodiff
  input = 2d.i
  exodiff = 2d_out.e
  requirement = 'The system shall be able to compute thermal compliance and its sensitivity in the heat transfer module.'
[]

[steady_2d]
  type = 'Exodiff'
  input = 'steady_2d.i'
  exodiff = 'steady_2d_out.e'
  requirement = 'The system shall model steady state heat transfer with an interface between two domains in 2D.'
  issues = '#21988'
  design = 'ThinLayerHeatTransfer.md'
  mesh_mode = 'REPLICATED'
  recover = false
[]

[steady_3d]
  type = 'Exodiff'
  input = 'steady_3d.i'
  exodiff = 'steady_3d_out.e'
  requirement = 'The system shall model steady state heat transfer with an interface between two domains in 3D.'
  issues = '#21988'
  design = 'ThinLayerHeatTransfer.md'
  mesh_mode = 'REPLICATED'
  recover = false
[]

[transient_2d]
  type = 'Exodiff'
  input = 'transient_2d.i'
  exodiff = 'transient_2d_out.e'
  requirement = 'The system shall model transient heat transfer with an interface between two domains in 2D.'
  issues = '#21988'
  design = 'ThinLayerHeatTransfer.md'
  mesh_mode = 'REPLICATED'
  recover = false
[]

[transient_3d]
  type = 'Exodiff'
  input = 'transient_3d.i'
  exodiff = 'transient_3d_out.e'
  requirement = 'The system shall model transient heat transfer with an interface between two domains in 3D.'
  issues = '#21988'
  design = 'ThinLayerHeatTransfer.md'
  mesh_mode = 'REPLICATED'
  recover = false
[]

[jacobian]
  type = 'PetscJacobianTester'
  input = 'transient_2d.i'
  run_sim = 'True'
  petsc_version = '>=3.9'
  ratio_tol = 2e-7
  issues = '#21988'
  design = 'ThinLayerHeatTransfer.md'
  difference_tol = 1e-6
  mesh_mode = 'REPLICATED'
  recover = false
  requirement = 'The Jacobian for the ThinLayerHeatTransfer calculations shall provide perfect jacobians.'
[]

[transient_heat]
  type = 'Exodiff'
  input = 'transient_heat.i'
  exodiff = 'transient_heat_out.e'
  requirement = 'The system shall compute the time derivative term of the heat equation'
  issues = '#7759'
  design = 'SpecificHeatConductionTimeDerivative.md'
[]

[transient_heat_derivatives]
  type = 'CSVDiff'
  input = 'transient_heat_derivatives.i'
  csvdiff = 'transient_heat_derivatives_out.csv'
  requirement = 'The system shall compute the time derivative term of the heat equation and utilize the derivative of the specific heat, thermal conductivity, and density to solve a heat conduction problem'
  issues = '#26647'
  design = 'HeatConductionTimeDerivative.md'
[]

[transient_heat_t_limit]
  type = 'RunApp'
  input = 'transient_heat_derivatives.i'
  expect_out = '(is below min_T|specified minimum)'
  cli_args = 'Materials/constant/min_T=10'
  requirement = 'The system shall use a lower temperature limit to compute the specific heat'
  issues = '#26647'
  design = 'HeatConductionTimeDerivative.md'
[]

[line]
  type = 'Exodiff'
  input = 'line.i'
  exodiff = 'line_out.e'
  rel_err = 1e-6
  abs_zero = 1e-10
  detail = 'using 1D truss elements.'
[]

[strip]
  type = 'Exodiff'
  input = 'strip.i'
  exodiff = 'strip_out.e'
  rel_err = 1e-6
  abs_zero = 1e-10
  detail = 'using 2D continuum elements in a 2D medium.'
[]

[rectangle_with_strip]
  type = 'Exodiff'
  input = 'rectangle_w_strip.i'
  exodiff = 'rectangle_w_strip_out.e'
  rel_err = 1e-6
  abs_zero = 1e-10
  detail = 'using 2D continuum elements for the bar contiguously meshed with a 2D medium.'
[]

[rectangle_with_line]
  type = 'Exodiff'
  input = 'rectangle_w_line.i'
  exodiff = 'rectangle_w_line_out.e'
  rel_err = 1e-6
  abs_zero = 1e-10
  detail = 'using 1D truss elements embedded in a 2D medium and connected using constraints.'
[]

[block_with_bar]
  type = 'Exodiff'
  input = 'block_w_bar.i'
  exodiff = 'block_w_bar_out.e'
  rel_err = 1e-6
  abs_zero = 1e-10
  detail = 'using 3D continuum elements for the bar contiguously meshed with a 3D medium.'
[]

[block_with_line]
  type = 'Exodiff'
  input = 'block_w_line.i'
  exodiff = 'block_w_line_out.e'
  rel_err = 1e-6
  abs_zero = 1e-10
  detail = 'using 1D truss elements embedded in a 3D medium and connected using constraints.'
[]

[plotting]
  requirement = 'The system shall generate comparison plots of the thermal solutions for a bar embedded in a continuum represented various ways.'
  type = RunCommand
  method = opt
  command = "python3 temp_plots.py"
  required_python_packages = 'matplotlib'
[]

[./1D_transient]
  type = 'Exodiff'
  input = '1D_transient.i'
  exodiff = '1D_transient_out.e'
  abs_zero = 1e-8
  requirement = "Heat conduction shall match the answer from an analytical solution in 1D"
  issues = "#5975"
[../]

[./ad_1D_transient]
  type = 'Exodiff'
  input = 'ad_1D_transient.i'
  exodiff = 'ad_1D_transient_out.e'
  abs_zero = 1e-8
  requirement = "Heat conduction from an AD kernel shall get the same answer as a traditional kernel in 1D"
  issues = "#5658 #12633"
[../]

[./ad_1D_transient_jacobian]
  type = 'PetscJacobianTester'
  input = 'ad_1D_transient.i'
  ratio_tol = 1e-7
  difference_tol = 1e-5
  run_sim = True
  cli_args = 'Executioner/num_steps=5'
  requirement = 'AD heat conduction and the Jacobian shall be beautiful in 1D'
  issues = "#5658 #12633"
[../]

[./2D_steady_state]
  type = 'Exodiff'
  input = '2d_steady_state.i'
  exodiff = '2d_steady_state_out.e'
  abs_zero = 1e-8
  requirement = "Heat conduction shall match the answer from an analytical solution in 2D"
  issues = "#8194"
[../]

[./ad_2D_steady_state]
  type = 'Exodiff'
  input = 'ad_2d_steady_state.i'
  exodiff = 'ad_2d_steady_state_out.e'
  abs_zero = 1e-8
  requirement = "Heat conduction from an AD kernel shall get the same answer as a traditional kernel in 2D"
  issues = "#5658 #12633"
[../]

[./ad_2D_steady_state_jacobian]
  type = 'PetscJacobianTester'
  input = 'ad_2d_steady_state.i'
  ratio_tol = 1e-7
  difference_tol = 1e-5
  run_sim = True
  requirement = 'AD heat conduction and the Jacobian shall be beautiful in 2D'
  issues = "#5658 #12633"
[../]

[unnormalized]
  type = CSVDiff
  input = 'view_factor_3d.i'
  csvdiff = 'view_factor_3d_unnormalized_out.csv'
  cli_args = 'Outputs/file_base=view_factor_3d_unnormalized_out
                UserObjects/unobstructed_vf/normalize_view_factor=false'
  design = 'UnobstructedPlanarViewFactor.md'
  mesh_mode = REPLICATED
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces without normalization.'
[]

[obstructed]
  type = CSVDiff
  input = 'view_factor_obstructed.i'
  csvdiff = 'view_factor_obstructed_out.csv'
  design = 'RayTracingViewFactor.md'
  requirement = 'The system shall compute view factors for cavities with obstruction using ray tracing.'
[]

[analytical2D]
  type = CSVDiff
  input = 'view_factor_2d.i'
  csvdiff = 'view_factor_2d_out.csv'
  design = 'UnobstructedPlanarViewFactor.md'
  mesh_mode = REPLICATED
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces in two-dimensional meshes using simple quadrature rules.'
[]

[ray2D]
  type = CSVDiff
  input = 'view_factor_2d.i'
  csvdiff = 'view_factor_2d_ray_out.csv'
  cli_args = 'Outputs/file_base=view_factor_2d_ray_out
                UserObjects/active=\'vf_study rt_vf\'
                RayBCs/vf/type=ViewFactorRayBC
                RayBCs/vf/boundary=\'left right bottom top\'
                GlobalParams/view_factor_object_name=rt_vf'
  design = 'RayTracingViewFactor.md'
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces in two-dimensional meshes using ray tracing.'
[]

[analytical3D]
  type = CSVDiff
  input = 'view_factor_3d.i'
  csvdiff = 'view_factor_3d_out.csv'
  design = 'UnobstructedPlanarViewFactor.md'
  mesh_mode = REPLICATED
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces in three-dimensional meshes using simple quadrature rules.'
[]

[ray3D]
  type = CSVDiff
  input = 'view_factor_3d.i'
  csvdiff = 'view_factor_3d_ray_out.csv'
  cli_args = 'Outputs/file_base=view_factor_3d_ray_out
                UserObjects/active=\'vf_study rt_vf\'
                RayBCs/vf/type=ViewFactorRayBC
                RayBCs/vf/boundary=\'left right front back bottom top\'
                GlobalParams/view_factor_object_name=rt_vf'
  design = 'RayTracingViewFactor.md'
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces in three-dimensional meshes using ray tracing.'
[]

[ray3D_nonplanar]
  type = RunApp
  input = 'view_factor_3d_non_planar_face.i'
  expect_out = '732 rays were skipped as they exited the mesh at their starting point through non-planar sides.'
  design = 'RayTracingViewFactor.md'
  requirement = 'The system shall be able to skip rays that exit the mesh when starting from non-planar faces in three-dimensional problems.'
  recover = false
  mesh_mode = REPLICATED
[]

[cavity_with_pillars]
  type = CSVDiff
  input = 'cavity_with_pillars.i'
  csvdiff = 'cavity_with_pillars_out.csv'
  requirement = 'The system shall support ensure that symmetry boundary conditions provide exactly the same answer as unfolding the problem about its axis of symmetry.'
  #see 21185
  allow_warnings = true
[]

[cavity_with_pillars_symmetry_bc]
  type = CSVDiff
  input = 'cavity_with_pillars_symmetry_bc.i'
  csvdiff = 'cavity_with_pillars_symmetry_bc_out.csv'
  requirement = 'The system shall support symmetry boundary conditions for view factor calculations.'
  #see 21185
  allow_warnings = true
[]