Heat Conduction Software Design Description

Introduction

The Software Design Description (SDD) for Heat Conduction is a description of the software product and serves as guideline that describes the architecture of the system and all of its parts.

Dependencies

The Heat Conduction application is developed using MOOSE and is based on various modules, as such the SDD for Heat Conduction is dependent upon the following documents.

• Framework System Design Description

Requirements Cross Reference

The following is a list of all design documents and the associated requirements for Heat Conduction.

• heat_conduction: ADHeatConduction
• F5.1.1AD heat conduction and the Jacobian shall be beautiful

  Specification: ad_heat_conduction:jacobian_test

  Design: ADHeatConduction

  Issue(s): #5658#12633

  16

• heat_conduction: ConvectiveFluxFunction
• F5.2.1The system shall allow prescribing a convective flux boundary condition using a constant heat transfer coefficient.

  Specification: convective_flux_function:constant

  Design: ConvectiveFluxFunction

  Issue(s): #14418

  20
• F5.2.2The system shall allow prescribing a convective flux boundary condition using a heat transfer coefficient that is a function of position and time.

  Specification: convective_flux_function:time_dependent

  Design: ConvectiveFluxFunction

  Issue(s): #14418

  20
• F5.2.3The system shall allow prescribing a convective flux boundary condition using a heat transfer coefficient that is a function of temperature.

  Specification: convective_flux_function:temperature_dependent

  Design: ConvectiveFluxFunction

  Issue(s): #14418

  20

• heat_conduction: PatchSidesetGenerator
• F5.3.1The system shall be able to divide a sideset into patches for more accurate radiative transfer modeling.

  Specification: generate_radiation_patch:generate_radiation_patch

  Design: PatchSidesetGenerator

  Issue(s): #14000

  20
• F5.3.2The system shall be able to use linear partitioner for subdividing sidesets into patches.

  Specification: generate_radiation_patch:generate_radiation_patch_linear

  Design: PatchSidesetGenerator

  Issue(s): #14000

  20
• F5.3.3The system shall be able to use centroid partitioner for subdividing sidesets into patches.

  Specification: generate_radiation_patch:generate_radiation_patch_centroid

  Design: PatchSidesetGenerator

  Issue(s): #14000

  20
• F5.3.4The system shall error when centroid partitioner is used but centroid_partitioner_direction is not provided.

  Specification: generate_radiation_patch:generate_radiation_patch_centroid_error

  Design: PatchSidesetGenerator

  Issue(s): #14000

  20

• heat_conduction: ConstantViewFactorSurfaceRadiation
• F5.4.1The system shall check consistency of boundary and emissivity entries

  Specification: gray_lambert_radiator:inconsistent_bnd_eps

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.2The system shall check consistency of boundary and view factor entries

  Specification: gray_lambert_radiator:inconsistent_bnd_view_factors

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.3The system shall check consistency of fixed_boundary_temperatures and fixed_temperature_boundary entries

  Specification: gray_lambert_radiator:inconsistent_iso_temperature

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.4The system shall check consistency of boundary and fixed_temperature_boundary entries

  Specification: gray_lambert_radiator:inconsistent_bnd_iso_bnd

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.5The system shall check consistency of boundary and adiabatic_boundary entries

  Specification: gray_lambert_radiator:inconsistent_bnd_adiabatic_bnd

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.6The system shall check consistency of the view_factors entry

  Specification: gray_lambert_radiator:incorrect_view_factor_shape

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.7The system shall check consistency of the view_factors entry

  Specification: gray_lambert_radiator:bad_rowsum

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  5
• F5.4.8The system shall compute radiative transfer between gray Lambert surfaces

  Specification: gray_lambert_radiator:gray_lambert_cavity

  Design: ConstantViewFactorSurfaceRadiationSurfaceRadiationVectorPostprocessorViewfactorVectorPostprocessor

  Issue(s): #13918

  20
• F5.4.9The system shall allow coupling radiative transfer between gray Lambert surfaces to solving heat conduction

  Specification: gray_lambert_radiator:coupled_heat_conduction

  Design: ConstantViewFactorSurfaceRadiation

  Issue(s): #13918

  20

• heat_conduction: SurfaceRadiationVectorPostprocessor
• F5.4.8The system shall compute radiative transfer between gray Lambert surfaces

  Specification: gray_lambert_radiator:gray_lambert_cavity

  Design: ConstantViewFactorSurfaceRadiationSurfaceRadiationVectorPostprocessorViewfactorVectorPostprocessor

  Issue(s): #13918

  20

• heat_conduction: ViewfactorVectorPostprocessor
• F5.4.8The system shall compute radiative transfer between gray Lambert surfaces

  Specification: gray_lambert_radiator:gray_lambert_cavity

  Design: ConstantViewFactorSurfaceRadiationSurfaceRadiationVectorPostprocessorViewfactorVectorPostprocessor

  Issue(s): #13918

  20

• heat_conduction: GrayLambertNeumannBC
• F5.4.10The system shall allow reconstructing the spatial distribution of the emission component on a radiation boundary via the T4 law

  Specification: gray_lambert_radiator:coupled_heat_conduction_emission_reconstruction

  Design: GrayLambertNeumannBC

  Issue(s): #13918

  20

• heat_conduction: ViewFactorObjectSurfaceRadiation
• F5.4.11The system shall compute radiative transfer between gray Lambert surfaces when the view factors are provided by a userobject

  Specification: gray_lambert_radiator:gray_lambert_cavity_automatic_vf

  Design: ViewFactorObjectSurfaceRadiation

  Issue(s): #13918

  20
• F5.4.12The system shall compute radiative transfer between gray Lambert surfaces when the view factors are provided by a userobject

  Specification: gray_lambert_radiator:gray_lambert_cavity_automatic_vf_3D

  Design: ViewFactorObjectSurfaceRadiation

  Issue(s): #13918

  20

• heat_conduction: GapConductance
• F5.5.1MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:perfect

  Design: GapConductance

  Issue(s): #6750

  20
• F5.5.2MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:perfectQ8

  Design: GapConductance

  Issue(s): #6750

  20
• F5.5.3MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:perfectQ9

  Design: GapConductance

  Issue(s): #6750

  20
• F5.5.17Optionally a constant attenuation shall be applied to compute the gap conductance below a gap length threshold.

  Specification: heat_conduction/min_gap:min_gap_order_zero

  Design: GapConductanceGapHeatTransfer

  Issue(s): #13221

  20
• F5.5.18Optionally 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.

  Specification: heat_conduction/min_gap:min_gap_order_one

  Design: GapConductanceGapHeatTransfer

  Issue(s): #13221

  20

• heat_conduction: Thermal Contact Action
• F5.5.4MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:nonmatching

  Design: Thermal Contact Action

  Issue(s): #6750

  20
• F5.5.5MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:second_order

  Design: Thermal Contact Action

  Issue(s): #6750

  20
• F5.5.6MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:moving

  Design: Thermal Contact Action

  Issue(s): #6750

  20
• F5.5.7MOOSE shall compute the heat transfer across small gaps for supported FEM orders and quadratures, moving interfaces, and non-matching meshes.

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:gap_conductivity_property

  Design: Thermal Contact Action

  Issue(s): #6750

  20
• F5.5.8MOOSE shall throw an error if the gap conductance model is used with uniform mesh refinement

  Specification: heat_conduction/2d_quadrature_gap_heat_transfer:gap_conductivity_property_r1_error

  Design: Thermal Contact Action

  Issue(s): #13043

  20
• F5.5.9The system shall support thermal contact with linear 3d hexahedral elements

  Specification: heat_conduction/3d_quadrature_gap_heat_transfer:nonmatching

  Design: Thermal Contact Action

  Issue(s): #6750

  20
• F5.5.10The system shall support thermal contact with second-order 3d hexahedral elements

  Specification: heat_conduction/3d_quadrature_gap_heat_transfer:second

  Design: Thermal Contact Action

  Issue(s): #6750

  20
• F5.5.11The system shall support thermal contact with 3d hexahedral elements where the surfaces move relative to one another

  Specification: heat_conduction/3d_quadrature_gap_heat_transfer:moving

  Design: Thermal Contact Action

  Issue(s): #6750

  20

• heat_conduction: CoupledConvectiveHeatFluxBC
• F5.5.12The system shall provide convective heat flux boundary condition where far-field temperature and convective heat transfer coefficient are given as constant variables

  Specification: heat_conduction/coupled_convective_heat_flux:const_hw

  Design: CoupledConvectiveHeatFluxBC

  Issue(s): #11631

  20
• F5.5.13The system shall provide convective heat flux boundary condition where far-field temperature and convective heat transfer coefficient are given as spatially varying variables

  Specification: heat_conduction/coupled_convective_heat_flux:coupled_convective_heat_flux

  Design: CoupledConvectiveHeatFluxBC

  Issue(s): #11631

  20
• F5.5.14The 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

  Specification: heat_conduction/coupled_convective_heat_flux:coupled_convective_heat_flux_two_phase

  Design: CoupledConvectiveHeatFluxBC

  Issue(s): #11631

  20
• F5.5.15The system shall report an error if the number of alpha components does not match the number of T_infinity components.

  Specification: heat_conduction/coupled_convective_heat_flux:not_enough_alpha

  Design: CoupledConvectiveHeatFluxBC

  Issue(s): #11631

  20
• F5.5.16The system shall report an error if the number of htc components does not match the number of T_infinity components.

  Specification: heat_conduction/coupled_convective_heat_flux:not_enough_htc

  Design: CoupledConvectiveHeatFluxBC

  Issue(s): #11631

  20

• heat_conduction: GapHeatTransfer
• F5.5.17Optionally a constant attenuation shall be applied to compute the gap conductance below a gap length threshold.

  Specification: heat_conduction/min_gap:min_gap_order_zero

  Design: GapConductanceGapHeatTransfer

  Issue(s): #13221

  20
• F5.5.18Optionally 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.

  Specification: heat_conduction/min_gap:min_gap_order_one

  Design: GapConductanceGapHeatTransfer

  Issue(s): #13221

  20
• F5.11.1The ThermalContact system shall enforce heat transfer across a meshed gap in a 2D plane geometry.

  Specification: meshed_gap_thermal_contact:test

  Design: ThermalContact SystemGapHeatTransfer

  Issue(s): #716

  20
• F5.11.4The ThermalContact system shall enforce heat transfer across a meshed circular annulus in a 2D plane geometry.

  Specification: meshed_gap_thermal_contact:annulus

  Design: ThermalContact SystemGapHeatTransfer

  Issue(s): #716

  20

• heat_conduction: AnisoHeatConductionMaterial
• F5.6.1The system shall allow the use of an anisotropic heat conduction material set by postprocessors.

  Specification: heat_conduction_ortho:test

  Design: AnisoHeatConductionMaterial

  Issue(s): #2674

  20

• heat_conduction: HeatConduction
• F5.7.1The system shall compute a tri-linear temperature field

  Specification: heat_conduction_patch:test

  Design: HeatConduction

  Issue(s): #6750

  20
• F5.7.2The system shall compute a bi-linear temperature field for an axisymmetric problem with quad8 elements

  Specification: heat_conduction_patch:test_rz_quad8

  Design: HeatConduction

  Issue(s): #6750

  20
• F5.7.3The system shall compute a bi-linear temperature field for an axisymmetric problem

  Specification: heat_conduction_patch:test_rz

  Design: HeatConduction

  Issue(s): #6750

  20
• F5.7.4The system shall compute a tri-linear temperature field with hex20 elements

  Specification: heat_conduction_patch:test_hex20

  Design: HeatConduction

  Issue(s): #6750

  20
• F5.7.5The system shall compute a tri-linear temperature field with hex20 elements using an anisotropic thermal conductivity model with isotropic thermal conductivities supplied

  Specification: heat_conduction_patch:test_hex20_aniso

  Design: HeatConduction

  Issue(s): #6750

  20
• F5.20.1Heat conduction shall match the answer from an analytical solution

  Specification: verify_against_analytical:1D_transient

  Design: HeatConduction

  Issue(s): #5975

  20
• F5.20.2Heat conduction from an AD kernel shall get the same answer as a traditional kernel

  Specification: verify_against_analytical:ad_1D_transient

  Design: HeatConduction

  Issue(s): #5658#12633

  20
• F5.20.3AD heat conduction and the Jacobian shall be beautiful

  Specification: verify_against_analytical:ad_1D_transient_jacobian

  Design: HeatConduction

  Issue(s): #5658#12633

  12
• F5.20.4Heat conduction shall match the answer from an analytical solution

  Specification: verify_against_analytical:2D_steady_state

  Design: HeatConduction

  Issue(s): #8194

  20
• F5.20.5Heat conduction from an AD kernel shall get the same answer as a traditional kernel

  Specification: verify_against_analytical:ad_2D_steady_state

  Design: HeatConduction

  Issue(s): #5658#12633

  20
• F5.20.6AD heat conduction and the Jacobian shall be beautiful

  Specification: verify_against_analytical:ad_2D_steady_state_jacobian

  Design: HeatConduction

  Issue(s): #5658#12633

  12

• heat_conduction: HeatSource
• F5.8.1MOOSE shall reproduce an analytical solution of a heat source in a 1D ceramic bar

  Specification: heat_source_bar:heat_source_bar

  Design: HeatSource

  Issue(s): #2582

  20

• heat_conduction: ADMatHeatSource
• F5.8.2MOOSE shall reproduce an analytical solution of a heat source in a 1D ceramic bar using AD kernels

  Specification: heat_source_bar:ad_heat_source_bar

  Design: ADMatHeatSource

  Issue(s): #12633

  20
• F5.8.3The AD heat conduction and heat source Jacobian shall be beautiful

  Specification: heat_source_bar:ad_heat_source_bar_jacobian

  Design: ADMatHeatSource

  Issue(s): #5658#12633

  16

• heat_conduction: HomogenizedHeatConduction
• F5.9.1The system shall compute homogenized thermal conductivity using the asymptotic expansion homogenization approach

  Specification: homogenization:heatConduction_test

  Design: HomogenizedHeatConduction

  Issue(s): #6750

  20
• F5.13.1The system shall compute total heat flux from heat transfer coefficient and temperature difference

  Specification: postprocessors:convective_ht_side_integral

  Design: HomogenizedHeatConduction

  Issue(s): #14390

  20

• heat_conduction: JouleHeatingSource
• F5.10.1The system shall compute Joule heating

  Specification: joule_heating:test

  Design: JouleHeatingSource

  Issue(s): #8220

  20

• heat_conduction: ThermalContact System
• F5.11.1The ThermalContact system shall enforce heat transfer across a meshed gap in a 2D plane geometry.

  Specification: meshed_gap_thermal_contact:test

  Design: ThermalContact SystemGapHeatTransfer

  Issue(s): #716

  20
• F5.11.2The ThermalContact system shall correctly enforce heat transfer across a meshed gap in a 2D plane geometry using a prescribed constant conductance.

  Specification: meshed_gap_thermal_contact:constant_conductance

  Design: ThermalContact SystemGapConductanceConstant

  Issue(s): #13061

  20
• F5.11.3The 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

  Specification: meshed_gap_thermal_contact:constant_conductance_quadrature

  Design: ThermalContact SystemGapConductanceConstant

  Issue(s): #13061

  20
• F5.11.4The ThermalContact system shall enforce heat transfer across a meshed circular annulus in a 2D plane geometry.

  Specification: meshed_gap_thermal_contact:annulus

  Design: ThermalContact SystemGapHeatTransfer

  Issue(s): #716

  20

• heat_conduction: GapConductanceConstant
• F5.11.2The ThermalContact system shall correctly enforce heat transfer across a meshed gap in a 2D plane geometry using a prescribed constant conductance.

  Specification: meshed_gap_thermal_contact:constant_conductance

  Design: ThermalContact SystemGapConductanceConstant

  Issue(s): #13061

  20
• F5.11.3The 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

  Specification: meshed_gap_thermal_contact:constant_conductance_quadrature

  Design: ThermalContact SystemGapConductanceConstant

  Issue(s): #13061

  20

• heat_conduction: ElementIntegralVariablePostprocessor
• F5.12.1The system shall computed an integrated value on elements in parallel

  Specification: parallel_element_pps_test:test

  Design: ElementIntegralVariablePostprocessor

  Issue(s): #861

  2

• heat_conduction: Radiation Transfer Action
• F5.14.1The system shall provide an action to set up radiative heat transfer problems using the net radiation method

  Specification: radiation_transfer_action:radiative_transfer_action

  Design: Radiation Transfer Action

  Issue(s): #13918

  16
• F5.14.2The system shall provide an action to set up radiative heat transfer problems using the net radiation method

  Specification: radiation_transfer_action:radiative_transfer_no_action

  Design: Radiation Transfer Action

  Issue(s): #13918

  16

• heat_conduction: InfiniteCylinderRadiativeBC
• F5.15.1Moose shall be able to model radiative transfer from a cylindrical surface as boundary condition.

  Specification: radiative_bcs:radiative_bc_cyl

  Design: InfiniteCylinderRadiativeBC

  Issue(s): #13053

  20

• heat_conduction: Heat Conduction Module
• F5.16.1MOOSE shall run a simulation with heat conduction, a heat source, thermal contact, and boundary conditions.

  Specification: recover:recover_1

  Design: Heat Conduction Module

  Issue(s): #10079

  20
• F5.16.2MOOSE shall run a short simulation with heat conduction, a heat source, thermal contact, and boundary conditions.

  Specification: recover:recover_2

  Design: Heat Conduction Module

  Issue(s): #10079

  20
• F5.16.3MOOSE 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.

  Specification: recover:recover_3

  Design: Heat Conduction Module

  Issue(s): #10079

  20
• F5.16.4MOOSE shall run a simulation with heat conduction, a heat source, thermal contact, and boundary conditions with automatic differentiation.

  Specification: recover:ad_recover_1

  Design: Heat Conduction Module

  Issue(s): #10079

  20
• F5.16.5MOOSE shall run a short simulation with heat conduction, a heat source, thermal contact, and boundary conditions with automatic differentiation.

  Specification: recover:ad_recover_2

  Design: Heat Conduction Module

  Issue(s): #10079

  20
• F5.16.6MOOSE 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.

  Specification: recover:ad_recover_3

  Design: Heat Conduction Module

  Issue(s): #10079

  20
• F5.16.7AD heat conduction and the Jacobian shall be beautiful

  Specification: recover:ad_recover_jacobian

  Design: Heat Conduction Module

  Issue(s): #5658#12633

• heat_conduction: SemiconductorLinearConductivity
• F5.17.1The system shall compute conductivity of semiconductors according to the Steinhart-Hart equation

  Specification: semiconductor_linear_conductivity:test

  Design: SemiconductorLinearConductivity

  Issue(s): #10278

  20

• heat_conduction: SideSetHeatTransferKernel
• F5.18.1The system shall solve the side set heat transfer model with:
  1. discontinuous finite elements,
  2. bulk gap temperature as an auxiliary variable,
  3. temperature dependent gap conductivity, and
  4. block restricted continuous finite element variables.

  Specification: sideset_heat_transfer:group

  Design: SideSetHeatTransferKernelSideSetHeatTransferMaterial

  Issue(s): #14519

  20202020
• F5.18.2MOOSE shall throw an error if the inputted boundary does not exist.

  Specification: sideset_heat_transfer:1D_gap_err

  Design: SideSetHeatTransferKernelSideSetHeatTransferMaterial

  Issue(s): #14519

  20

• heat_conduction: SideSetHeatTransferMaterial
• F5.18.1The system shall solve the side set heat transfer model with:
  1. discontinuous finite elements,
  2. bulk gap temperature as an auxiliary variable,
  3. temperature dependent gap conductivity, and
  4. block restricted continuous finite element variables.

  Specification: sideset_heat_transfer:group

  Design: SideSetHeatTransferKernelSideSetHeatTransferMaterial

  Issue(s): #14519

  20202020
• F5.18.2MOOSE shall throw an error if the inputted boundary does not exist.

  Specification: sideset_heat_transfer:1D_gap_err

  Design: SideSetHeatTransferKernelSideSetHeatTransferMaterial

  Issue(s): #14519

  20

• heat_conduction: SpecificHeatConductionTimeDerivative
• F5.19.1The system shall compute the time derivative term of the heat equation

  Specification: transient_heat:test

  Design: SpecificHeatConductionTimeDerivative

  Issue(s): #7759

  20

• heat_conduction: UnobstructedPlanarViewFactor
• F5.21.1The system shall compute view factors for unobstructed, planar surfaces in radiative exchange

  Specification: view_factors:view_factor_cube

  Design: UnobstructedPlanarViewFactor

  Issue(s): #13918

  20
• F5.21.2The system shall compute view factors for unobstructed, planar surfaces in radiative exchange and normalize view factors from each surface to sum to one

  Specification: view_factors:view_factor_cube_normalized

  Design: UnobstructedPlanarViewFactor

  Issue(s): #13918

  20

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

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

[./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'
  issues = "#14000"
[../]

[./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'
  issues = "#14000"
[../]

[./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'
  issues = "#14000"
[../]

[./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'
  issues = "#14000"
[../]

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

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

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

[./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 when the view factors are provided by a userobject'
  design = 'ViewFactorObjectSurfaceRadiation.md'
[../]

[./perfect]
  type = 'Exodiff'
  input = 'perfect.i'
  exodiff = 'perfect_out.e'
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  design = "/materials/GapConductance.md"
  issues = "#6750"
[../]

[./perfectQ8]
  type = 'Exodiff'
  input = 'perfectQ8.i'
  exodiff = 'perfectQ8_out.e'
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  design = "/materials/GapConductance.md"
  issues = "#6750"
[../]

[./perfectQ9]
  type = 'Exodiff'
  input = 'perfectQ9.i'
  exodiff = 'perfectQ9_out.e'
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  design = "/materials/GapConductance.md"
  issues = "#6750"
[../]

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

[./nonmatching]
  type = 'Exodiff'
  input = 'nonmatching.i'
  exodiff = 'nonmatching_out.e'
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./second_order]
  type = 'Exodiff'
  input = 'second_order.i'
  exodiff = 'second_order_out.e'
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./moving]
  type = 'Exodiff'
  input = 'moving.i'
  exodiff = 'moving_out.e'
  allow_warnings = true
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  design = "/actions/ThermalContactAction.md"
  issues = "#6750"
[../]

[./gap_conductivity_property]
  type = 'Exodiff'
  input = 'gap_conductivity_property.i'
  exodiff = 'gap_conductivity_property_out.e'
  requirement = "MOOSE shall compute the heat transfer across small gaps for supported
                   FEM orders and quadratures, moving interfaces, and non-matching meshes."
  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 = "MOOSE 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.'
[../]

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

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

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

[./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"
  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"
  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
  petsc_version = '>=3.9.4'
  method = 'OPT'
  cli_args = 'Executioner/num_steps=5'
  requirement = 'AD heat conduction and the Jacobian shall be beautiful'
  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"
  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"
  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
  petsc_version = '>=3.9.4'
  method = 'OPT'
  requirement = 'AD heat conduction and the Jacobian shall be beautiful'
  issues = "#5658 #12633"
[../]

[./heat_source_bar]
  type = 'Exodiff'
  input = 'heat_source_bar.i'
  exodiff = 'heat_source_bar_out.e'
  abs_zero = 1e-7
  requirement = 'MOOSE 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 = 'MOOSE 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
  petsc_version = '>=3.9.4'
  requirement = 'The AD heat conduction and heat source Jacobian shall be beautiful'
  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'
[../]

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

[./test]
  type = 'Exodiff'
  input = 'transient_jouleheating.i'
  exodiff = 'transient_jouleheating_out.e'
  requirement = 'The system shall compute Joule heating'
[../]

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

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

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

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

[./radiative_transfer_no_action]
  type = 'Exodiff'
  input = ' radiative_transfer_no_action.i'
  exodiff = 'radiative_transfer_no_action_out.e'
  requirement = 'The system shall provide an action to set up radiative heat transfer problems using the net radiation method'
  petsc_version = '>=3.9.4'
[../]

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

[./recover_1]
  type = Exodiff
  input = recover.i
  exodiff = recover_out.e
  superlu = true
  requirement = 'MOOSE 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 = 'MOOSE 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 = 'MOOSE 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 = 'MOOSE 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 = 'MOOSE 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 = 'MOOSE 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.'
[../]

[./ad_recover_jacobian]
  type = 'PetscJacobianTester'
  input = 'ad_recover.i'
  ratio_tol = 1e-7
  difference_tol = 1e-5
  run_sim = True
  petsc_version = '>=3.9'
  method = 'OPT'
  cli_args = 'Executioner/num_steps=2'
  skip = 'AD Thermal Contact not yet implemented'
  requirement = 'AD heat conduction and the Jacobian shall be beautiful'
  issues = "#5658 #12633"
[../]

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

[group]
  requirement = 'The system shall solve the side set heat transfer model with:'
  [./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 = 'MOOSE shall throw an error if the inputted boundary does not exist.'
[../]

[group]
  requirement = 'The system shall solve the side set heat transfer model with:'
  [./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 = 'MOOSE shall throw an error if the inputted boundary does not exist.'
[../]

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

[./view_factor_cube]
  type = CSVDiff
  input = 'view_factor_cube.i'
  csvdiff = 'view_factor_cube_out.csv'
  cli_args = 'UserObjects/view_factor/normalize_view_factor=false'
  design = 'UnobstructedPlanarViewFactor.md'
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces in radiative exchange'
[../]

[./view_factor_cube_normalized]
  type = CSVDiff
  input = 'view_factor_cube.i'
  csvdiff = 'view_factor_cube_normalized.csv'
  cli_args = 'Outputs/file_base=view_factor_cube_normalized'
  design = 'UnobstructedPlanarViewFactor.md'
  requirement = 'The system shall compute view factors for unobstructed, planar surfaces in radiative exchange and normalize view factors from each surface to sum to one'
[../]