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

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

[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'
  ad_indexing_type = 'global'
[]

[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
  ad_indexing_type = 'global'
  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
  ad_indexing_type = 'global'
  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
  ad_indexing_type = 'global'
  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
  ad_indexing_type = 'global'
  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"
[]

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

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

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

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

[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'
  ad_indexing_type = 'global'
[]

[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'
  ad_indexing_type = 'global'
[]

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

[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'
  ad_indexing_type = 'global'
[]

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

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

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

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

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

[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'
  rel_err = 1.0e-3
  requirement = 'We shall be able to generate mortar numerical results for conduction and radiationthrough cylinders that are close to those generated by the node-to-segment formulation.'
  ad_indexing_type = 'global'
[]

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

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

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

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

[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'
  ad_indexing_type = 'global'
[]

[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'
  ad_indexing_type = 'global'
[]

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

[large_gap_heat_transfer_test_cylinder_mortar_error]
  type = RunException
  input = 'large_gap_heat_transfer_test_cylinder_mortar_error.i'
  expect_err = 'gap_geometry_type = CYLINDER'
  requirement = 'We shall be able to generate a meaningful error message if the user does not define the cylinder points when a cylindrical geometry has been chosen.'
  ad_indexing_type = 'global'
[]

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

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

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

[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
  min_ad_size = 100
  mesh_mode = 'REPLICATED'
  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.'
  ad_indexing_type = 'global'
[]

[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
  min_ad_size = 100
  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.'
  ad_indexing_type = 'global'
[]

[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
  min_ad_size = 100
  ad_indexing_type = 'global'
  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
  min_ad_size = 100
  ad_indexing_type = 'global'
  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'
  ad_indexing_type = 'global'
  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.'
[]

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

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

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

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

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

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

[./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 heat flux across a gap using the ThermalContact methods'
  design = 'source/bcs/GapHeatTransfer.md'
  issues = '#1609'
[../]

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

[./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 = 1e-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'
[../]

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

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

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

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

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

[./test]
  type = 'Exodiff'
  input = 'transient_heat.i'
  exodiff = 'transient_heat_out.e'
  requirement = 'The system shall compute the time derivative term of the heat equation'
[../]

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

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

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