Navier Stokes Requirement Traceability Matrix

[./advection_bc]
  type = 'Exodiff'
  input = 'advection_bc.i'
  exodiff = 'advection_bc_out.e'
  requirement = "The system shall compute inflow and outflow boundary conditions for advected variables"
[../]

[check_too_few_components]
  type = 'RunException'
  expect_err = 'Number of components of velocity_vector must be at least equal to the mesh dimension'
  input = '2d_advection_bc.i'
  cli_args = 'Outputs/exodus=false'
  requirement = 'We shall error if the user provides less velocity components than the mesh dimension'
[]

[check_too_many_components]
  type = RunException
  expect_err = 'You cannot supply more than 3 components for the velocity vector'
  input = 2d_advection_bc.i
  cli_args = "Outputs/exodus=false BCs/outflow_term/velocity_vector='1 0 0 0'"
  requirement = 'We shall error if the user provides more than 3 velocity components'
[]

[check_more_components_than_mesh_dim]
  type = Exodiff
  input = 2d_advection_bc.i
  exodiff = '2d_advection_bc_out.e'
  cli_args = "BCs/outflow_term/velocity_vector='1 0 0'"
  requirement = 'We shall allow the user to supply more velocity components than the mesh dimension (up to 3 components)'
[]

[./bump]
  type = 'Exodiff'
  input = 'bump.i'
  exodiff = 'bump_out.e'
  # This test has very small rhov and vel_y values that may report
  # relative differences when this test is run in parallel, so we
  # use a custom comparison file to switch those components to
  # absolute comparisons.
  custom_cmp = 'bump_out.cmp'
  threading = '!pthreads'
  requirement = 'The system shall be able to solve the Euler equations for subsonic flow with a bump in the channel.'
[../]

[./hllc_sod_shocktube_1D_benchmark]
  type = Exodiff
  input = 'hllc_sod_shocktube.i'
  exodiff = 'hllc_sod_shocktube_out.e'
  design = 'CNSFVHLLCBase.md'
  issues = '#16758'
  requirement = 'The system shall be able to solve the 1D Sod shock-tube benchmark problem using an HLLC scheme to compute convective fluxes.'
  cli_args = 'Executioner/num_steps=10'
[../]

[kt-van-leer-primitive]
  type = Exodiff
  input = transient-porous-kt-primitive.i
  exodiff = transient-porous-kt-primitive_exo.e
  ad_indexing_type = global
  requirement = 'The system shall be able to solve the steady Euler equations in a heated channel using Kurganov-Tadmor with linearly reconstructed data with Van-Leer limiting for the convection term and a primitive variable set and show a flat momentum profile'
  mumps = true
  recover = false # See https://github.com/idaholab/moose/issues/17906
[]

[./fv_implicit_bcs]
  type = Exodiff
  input = hllc_sod_shocktube.i
  exodiff = hllc_sod_shocktube_out.e
  abs_zero = 1e-9
  requirement = "The system shall be able to impose boundary advective fluxes for HLLC discretizations that use implicit/interior cell information."
[../]

[1d-free-flow-hllc]
  type = PythonUnitTest
  input = test.py
  test_case = Test1DFreeFlowHLLC
  requirement = 'The system shall exhibit first order convergence for all variables for the free-flow Euler equations using a HLLC discretization scheme for the advection flux and with specified temperature and momentum at one boundary and specified pressure at another boundary.'
  required_python_packages = 'pandas matplotlib'
  method = "!dbg"
  design = 'CNSFVHLLCBase.md'
[]

[constant]
  type = PythonUnitTest
  input = test.py
  test_case = Test1DPorousHLLC
  required_python_packages = 'pandas matplotlib'
  method = "!dbg"
  detail = 'constant porosity'
[]

[varying]
  type = PythonUnitTest
  input = test.py
  test_case = Test1DVaryingEpsPorousHLLC
  required_python_packages = 'pandas matplotlib'
  method = "!dbg"
  detail = 'varying porosity'
[]

[central_difference]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DPorousKTPrimitiveCD
  required_python_packages = 'pandas matplotlib'
  detail = 'when using central differencing to interpolate cell center values to faces and display second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[upwind]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DPorousKTPrimitiveUpwind
  required_python_packages = 'pandas matplotlib'
  detail = 'when using directional upwinding to interpolate cell center values to faces and display first order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[vanLeer]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DPorousKTPrimitiveVanLeer
  required_python_packages = 'pandas matplotlib'
  detail = 'when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[central_difference]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DPorousKTConservedCD
  required_python_packages = 'pandas matplotlib'
  detail = 'when using central differencing to interpolate cell center values to faces and display second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[upwind]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DPorousKTConservedUpwind
  required_python_packages = 'pandas matplotlib'
  detail = 'when using directional upwinding to interpolate cell center values to faces and display first order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[vanLeer]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DPorousKTConservedVanLeer
  required_python_packages = 'pandas matplotlib'
  detail = 'when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[central_difference]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKTPrimitiveCD
  required_python_packages = 'pandas matplotlib'
  detail = 'when using central differencing to interpolate cell center values to faces and display second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[upwind]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKTPrimitiveUpwind
  required_python_packages = 'pandas matplotlib'
  detail = 'when using directional upwinding to interpolate cell center values to faces and display first order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[vanLeer]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKTPrimitiveVanLeer
  required_python_packages = 'pandas matplotlib'
  detail = 'when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[central_difference]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKNPPrimitiveCD
  required_python_packages = 'pandas matplotlib'
  detail = 'when using central differencing to interpolate cell center values to faces and display second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[upwind]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKNPPrimitiveUpwind
  required_python_packages = 'pandas matplotlib'
  detail = 'when using directional upwinding to interpolate cell center values to faces and display first order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[vanLeer]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKNPPrimitiveVanLeer
  required_python_packages = 'pandas matplotlib'
  detail = 'when using linear interpolation of cell center values to faces with Van-Leer limiting and display at least second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[central_difference]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKTMixedCD
  required_python_packages = 'pandas matplotlib'
  detail = 'when using central differencing to interpolate cell center values to faces and display second order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[upwind]
  type = PythonUnitTest
  input = test.py
  test_case = TestBasic1DVaryingPorousKTMixedUpwind
  required_python_packages = 'pandas matplotlib'
  detail = 'when using directional upwinding to interpolate cell center values to faces and display first order convergence'
  ad_indexing_type = 'global'
  method = "!dbg"
[]

[hllc_with_volume_source]
  type = Exodiff
  input = 'straight-channel-hllc.i'
  exodiff = straight-channel-hllc_out.e
  detail = 'when the source has a cell-centered volumetric discretization'
  method = "!dbg"
[]

[conserved]
  type = Exodiff
  input = 'regular-straight-channel.i'
  exodiff = 'regular-straight-channel_out.e'
  detail = 'conservation of mass when no sources are present'
  cli_args = "FVKernels/inactive='drag'"
  method = "!dbg"
[]

[non_conserved]
  input = 'regular-straight-channel.i'
  exodiff = 'regular-straight-channel-sources_out.e'
  detail = 'violation of conservation of mass when sources are present'
  type = Exodiff
  cli_args = 'Outputs/file_base=regular-straight-channel-sources_out'
  method = "!dbg"
[]

[non_conserved_finer]
  input = 'regular-straight-channel.i'
  exodiff = 'regular-straight-channel-sources-nx-150_out.e'
  detail = 'lesser violation of conservation of mass when sources are present and the mesh is refined'
  type = Exodiff
  cli_args = 'Outputs/file_base=regular-straight-channel-sources-nx-150_out Mesh/gen_mesh/nx=150'
  method = "!dbg"
[]

[1d-hllc]
  type = PythonUnitTest
  input = test.py
  test_case = Test1DHLLC
  requirement = 'The system shall be able to solve compressible fluid flow kernels for mass, momentum, and energy with the addition of diffusion and display first order convergence for all variables when using a HLLC scheme for the convection terms.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  method = "!dbg"
[]

[ray_1pt7e4]
  type = RunApp
  input = average-boussinesq.i
  detail = '1.7e4'
  strumpack = true
  ad_indexing_type = 'global'
[]

[ray_1pt4e5]
  type = RunApp
  input = average-boussinesq.i
  cli_args = 'Mesh/gen/xmax=2e-2 Mesh/gen/ymax=2e-2 Outputs/file_base=ra-1e5'
  detail = '1.4e5'
  strumpack = true
  ad_indexing_type = 'global'
[]

[fv_specified_pressure_out]
  type = Exodiff
  input = subsonic_nozzle_fixed_inflow_hllc.i
  exodiff = subsonic_nozzle_fixed_inflow_hllc_out.e
  abs_zero = 1e-9
  requirement = "The system shall be able to model subsonic nozzle flow using an HLLC discretization with a specified outlet pressure."
  recover = false # See https://github.com/idaholab/moose/issues/17906
[]

[scalar_advection]
  type = Exodiff
  input = mass-frac-advection.i
  exodiff = mass-frac-advection_out.e
  requirement = 'The system shall be able to advect a scalar using density and velocity computed through solution of the Euler equations.'
  ad_indexing_type = 'global'
  abs_zero = 3e-9
[]

[./hllc_sod_shocktube_2D]
  type = RunApp
  input = 'hllc_sod_shocktube_2D.i'
  check_input = True
  requirement = "The system shall be able to run a two-dimensional version of Sod's shocktube problem."
[../]

[./supersonic_nozzle_hllc]
  type = Exodiff
  input = supersonic_nozzle_hllc.i
  exodiff = supersonic_nozzle_hllc_out.e
  cli_args = 'Executioner/end_time=0.1'
  abs_zero = 1e-9
  requirement = "The system shall be able to model supersonic nozzle flow using an HLLC advective flux discretization and with inlet boundary conditions based on stagnation temperature and stagnation pressure."
[../]

[continuous_eps_kt]
  type = Exodiff
  input = implicit-euler-basic-kt-primitive.i
  exodiff = implicit-euler-basic-kt-primitive_out.e
  detail = 'the Kurganov-Tadmor scheme'
  ad_indexing_type = 'global'
  recover = false # See https://github.com/idaholab/moose/issues/17906
[]

[continuous_eps_hllc]
  type = Exodiff
  input = hllc.i
  exodiff = hllc_out.e
  detail = 'the HLLC scheme'
  ad_indexing_type = 'global'
[]

[twod_y_channel_upwind_frictional_porosity_function]
  type = Exodiff
  input = rotated-2d-bkt-function-porosity.i
  exodiff = rotated-2d-bkt-function-porosity_out.e
  requirement = 'The system shall be able to solve a two-dimensional y-channel problem with frictional drag and a series of porosity jumps smoothed into a continuous porosity function, using the Kurganov-Tadmor scheme for computing intercell convective fluxes with upwind limiting interpolation (e.g. the + cell centroid value is used as the + side value at the face).'
  ad_indexing_type = 'global'
  mumps = true
  recover = false # See https://github.com/idaholab/moose/issues/17906
  method = "!dbg"
  valgrind = 'none'
  design = 'PCNSFVKT.md PNSFVPGradEpsilon.md PorousPrimitiveVarMaterial.md PNSFVMomentumFriction.md'
[]

[twod_y_channel_upwind_frictional_porosity_function_mixed]
  type = Exodiff
  input = rotated-2d-bkt-function-porosity-mixed.i
  exodiff = rotated-2d-bkt-function-porosity-mixed_out.e
  requirement = 'The system shall be able to solve a two-dimensional y-channel problem using a mixed variable set with frictional drag and a series of porosity jumps smoothed into a continuous porosity function, using the Kurganov-Tadmor scheme for computing intercell convective fluxes with upwind limiting interpolation (e.g. the + cell centroid value is used as the + side value at the face).'
  ad_indexing_type = 'global'
  mumps = true
  recover = false # See https://github.com/idaholab/moose/issues/17906
  method = "!dbg"
  valgrind = 'none'
  design = 'PCNSFVKT.md PNSFVPGradEpsilon.md PorousMixedVarMaterial.md PNSFVMomentumFriction.md'
[]

[deferred_correction]
  type = Exodiff
  input = dc.i
  exodiff = dc_out.e
  ad_indexing_type = 'global'
  requirement = 'The system shall support the deferred correction algorithm for transitioning from low-order to high-order representations of the convective flux during a transient simulation.'
  abs_zero = 1e-8
  design = 'PCNSFVKTDC.md'
[]

[./2D_symmetry_hllc]
  type = RunApp
  input = 2D_symmetry.i
  check_input = True
  requirement = 'The system shall be able to run a two-dimensional symmetric flow problem with an HLLC discretization for advection.'
[../]

[./HLLC_wave_speeds_1D]
  type = CSVDiff
  input = 'hllc_uo_1D.i'
  csvdiff = 'hllc_uo_1D_out_wave_speeds_0001.csv'
  requirement = "The system shall be able to compute wave speeds for HLLC Riemann solvers."
[../]

[./HLLC_wave_speeds_2D]
  type = CSVDiff
  input = 'hllc_uo_2D_tri.i'
  csvdiff = 'hllc_uo_2D_tri_out_wave_speeds_0001.csv'
  requirement = "The system shall be able to compute wave speeds for HLLC Riemann solvers in multiple dimensions."
[../]

[./local_porosity]
  type = Exodiff
  input = wall_heat_transfer.i
  exodiff = local_porosity_out.e
  abs_zero = 1e-9
  cli_args = "Outputs/file_base=local_porosity_out GlobalParams/splitting='porosity' GlobalParams/locality='local'"
  requirement = "The system shall provide a boundary condition to split a constant heat flux according to local values of porosity."
[../]

[./global_porosity]
  type = Exodiff
  input = wall_heat_transfer.i
  exodiff = global_porosity_out.e
  abs_zero = 1e-9
  cli_args = "Outputs/file_base=global_porosity_out GlobalParams/splitting='porosity' GlobalParams/locality='global'"
  requirement = "The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of porosity."
[../]

[./local_k]
  type = Exodiff
  input = wall_heat_transfer.i
  exodiff = local_k_out.e
  abs_zero = 1e-9
  cli_args = "Outputs/file_base=local_k_out GlobalParams/splitting='thermal_conductivity' GlobalParams/locality='local'"
  requirement = "The system shall provide a boundary condition to split a constant heat flux according to local values of thermal conductivity."
[../]

[./global_k]
  type = Exodiff
  input = wall_heat_transfer.i
  exodiff = global_k_out.e
  abs_zero = 1e-9
  cli_args = "Outputs/file_base=global_k_out GlobalParams/splitting='thermal_conductivity' GlobalParams/locality='global'"
  requirement = "The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of thermal conductivity."
[../]

[./local_kappa]
  type = Exodiff
  input = wall_heat_transfer.i
  exodiff = local_kappa_out.e
  abs_zero = 1e-9
  cli_args = "Outputs/file_base=local_kappa_out GlobalParams/splitting='effective_thermal_conductivity' GlobalParams/locality='local'"
  requirement = "The system shall provide a boundary condition to split a constant heat flux according to local values of effective thermal conductivity."
[../]

[./global_kappa]
  type = Exodiff
  input = wall_heat_transfer.i
  exodiff = global_kappa_out.e
  abs_zero = 1e-9
  cli_args = "Outputs/file_base=global_kappa_out GlobalParams/splitting='effective_thermal_conductivity' GlobalParams/locality='global'"
  requirement = "The system shall provide a boundary condition to split a constant heat flux according to domain-averaged values of effective thermal conductivity."
[../]

[wall-function-bc]
  type = 'Exodiff'
  input = Re_t395.i
  exodiff = Re_t395_out.e
  method = "!dbg"
  requirement = 'The system shall be able to impose a wall shear stress at the wall according to the algebraic wall function.'
  ad_indexing_type = 'global'
  issues = '#18273'
  design = 'rans_theory.md INSFVWallFunctionBC.md'
  mesh_mode = REPLICATED
  recover = false
[]

[cavity]
  type = Exodiff
  input = convection_cavity.i
  exodiff = 'convection_cavity_out.e'
  ad_indexing_type = 'global'
  detail = 'for a cavity problem'
  max_parallel = 1
[]

[channel]
  type = Exodiff
  input = convection_channel.i
  exodiff = 'convection_channel_out.e'
  ad_indexing_type = 'global'
  detail = 'and for a channel problem.'
[]

[block_restricted_variables]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = 2d-rc_out.e
  requirement = 'The system shall be able to block-restrict all variables in a heated channel simulation with passive scalar advection.'
  ad_indexing_type = 'global'
[]

[1e3]
  type = Exodiff
  input = boussinesq.i
  exodiff = 1e3.e
  cli_args = 'rayleigh=1e3 Outputs/file_base=1e3'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e3 using a finite volume discretization.'
  valgrind = 'none'
  ad_indexing_type = 'global'
[]

[1e4]
  type = Exodiff
  input = boussinesq.i
  exodiff = 1e4.e
  cli_args = 'rayleigh=1e4 Outputs/file_base=1e4'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e4 using a finite volume discretization.'
  valgrind = 'none'
  ad_indexing_type = 'global'
[]

[1e5]
  type = Exodiff
  input = boussinesq.i
  exodiff = 1e5.e
  cli_args = 'rayleigh=1e5 Outputs/file_base=1e5'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e5 using a finite volume discretization.'
  valgrind = 'none'
  ad_indexing_type = 'global'
[]

[1e6]
  type = Exodiff
  input = boussinesq.i
  exodiff = 1e6.e
  cli_args = 'rayleigh=1e6 Outputs/file_base=1e6'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e6 using a finite volume discretization.'
  abs_zero = 1e-9
  valgrind = 'none'
  ad_indexing_type = 'global'
[]

[average-no-slip]
  type = 'Exodiff'
  input = 2d-average-no-slip.i
  exodiff = 2d-average-no-slip_out.e
  requirement = 'The system shall be able to solve incompressible Navier-Stokes channel flow with no-slip boundary conditions on the wall in an axisymmetric coordinate system using an average interpolation scheme for the velocity.'
  ad_indexing_type = 'global'
[]

[rc-rz-no-slip-mass-conservation]
  type = 'Exodiff'
  input = 2d-average-no-slip.i
  cli_args = "velocity_interp_method='rc' Outputs/file_base=2d-rc-no-slip"
  exodiff = 2d-rc-no-slip.e
  requirement = 'The system shall be able to solve incompressible Navier-Stokes channel flow with no-slip boundary conditions on the wall in an axisymmetric coordinate system using a Rhie-Chow interpolation scheme for the velocity.'
  ad_indexing_type = 'global'
[]

[rc-free-slip]
  type = 'Exodiff'
  input = 2d-rc-slip.i
  exodiff = 2d-rc-slip_out.e
  requirement = 'The system shall be able to solve incompressible Navier-Stokes channel flow with free-slip boundary conditions on the wall in an axisymmetric coordinate system using a Rhie-Chow interpolation scheme for the velocity.'
  ad_indexing_type = 'global'
[]

[rz-diverging-no-slip]
  type = 'Exodiff'
  input = diverging.i
  exodiff = 'diverging-no-slip.e'
  cli_args = "FVBCs/active='inlet-u inlet-v no-slip-wall-u no-slip-wall-v outlet-p axis-u axis-v' Outputs/file_base='diverging-no-slip' -pc_type lu -pc_factor_shift_type NONZERO"
  requirement = 'The system shall be able to solve a diverging channel problem in cylindrical coordinates with no slip boundary conditions.'
  ad_indexing_type = 'global'
[]

[rz-diverging-free-slip]
  type = 'Exodiff'
  input = diverging.i
  exodiff = 'diverging-free-slip.e'
  cli_args = "FVBCs/active='inlet-u inlet-v free-slip-wall-u free-slip-wall-v outlet-p axis-u axis-v' Outputs/file_base='diverging-free-slip' -pc_type lu -pc_factor_shift_type NONZERO"
  requirement = 'The system shall be able to solve a diverging channel problem in cylindrical coordinates with free slip boundary conditions.'
  ad_indexing_type = 'global'
  abs_zero = 1e-9
[]

[rc-xyz-no-slip-mass-conservation]
  type = 'Exodiff'
  input = '2d-average-no-slip.i'
  exodiff = '2d-rc-no-slip-xyz-mass-conservation.e'
  cli_args = "velocity_interp_method='rc' Outputs/file_base=2d-rc-no-slip-xyz-mass-conservation Problem/coord_type='XYZ'"
  requirement = 'The system shall conserve mass when solving a Cartesian channel flow problem with one symmetry boundary condition and one no-slip wall boundary condition.'
  ad_indexing_type = 'global'
[]

[1d-rc]
  type = 'Exodiff'
  input = 1d-rc.i
  exodiff = 1d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to model free-slip conditions in a 1D channel; specifically the tangential velocity shall have a uniform value of unity and the pressure shall not change.'
  ad_indexing_type = 'global'
[]

[free-slip]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = 2d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to model free-slip conditions in a channel; specifically the tangential velocity shall have a uniform value of unity, the normal velocity shall have a uniform value of zero, and the pressure shall not change.'
  ad_indexing_type = 'global'
[]

[no-slip]
  type = 'Exodiff'
  input = 2d-rc-no-slip.i
  exodiff = 2d-rc-no-slip_out.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a channel; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = '-pc_type lu -pc_factor_shift_type NONZERO'
[]

[scalar-transport]
  type = 'Exodiff'
  input = 2d-scalar-transport.i
  exodiff = 2d-scalar-transport_out.e
  method = "!dbg"
  requirement = 'The system shall be able to transport arbitrary scalar field variables in a fluid flow field.'
  ad_indexing_type = 'global'
  design = 'INSFVScalarFieldAdvection.md'
  issues = '#16732'
[]

[momentum-outflow-bcs]
  type = 'Exodiff'
  input = 2d-rc-no-slip-outflow-bcs.i
  exodiff = 2d-rc-no-slip-outflow-bcs_out.e
  requirement = 'The system shall be able to use flux boundary conditions for the momentum and match results produced by using flux kernels.'
  design = 'INSFVMomentumAdvectionOutflowBC.md'
  issues = '#16854'
  ad_indexing_type = 'global'
  abs_zero = 2e-9
[]

[extrapolated-outlet-pressure]
  type = 'Exodiff'
  input = 2d-rc-no-slip-extrapolated-outlet-pressure.i
  exodiff = 2d-rc-no-slip-extrapolated-outlet-pressure_out.e
  requirement = 'The system shall be able to extrapolate a pressure value at a fully developed outflow boundary and use a mean pressure approach to eliminate the nullspace for the pressure.'
  design = 'INSFVMassAdvectionOutflowBC.md'
  issues = '#16854'
  ad_indexing_type = 'global'
  abs_zero = 6e-9
  cli_args = '-pc_type lu -sub_pc_factor_shift_type NONZERO'
[]

[mixing-length]
  type = 'Exodiff'
  input = 2d-mixing-length.i
  exodiff = 2d-mixing-length_out.e
  requirement = 'The system shall be able to model the effect of Reynolds-averaged parameters on the momentum and passive scalar advection equations using a mixing length model'
  ad_indexing_type = 'global'
  design = 'navier_stokes/rans_theory.md'
  issues = '#16794 #16937'
  mesh_mode = REPLICATED
[]

[linear-friction]
  type = 'Exodiff'
  input = 2d-rc-friction.i
  exodiff = 2d-rc-friction_out.e
  requirement = 'The system shall be able to model linear volumetric friction in a channel.'
  issues = '#16872'
  ad_indexing_type = 'global'
  design = 'NSFVMomentumFriction.md'
[]

[quadratic-friction]
  type = 'Exodiff'
  input = 2d-rc-friction.i
  exodiff = 2d-rc-friction-quad_out.e
  requirement = 'The system shall be able to model quadratic volumetric friction in a channel.'
  issues = '#16872'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='2d-rc-friction-quad_out' FVKernels/inactive='u_friction_linear v_friction_linear'"
  design = 'NSFVMomentumFriction.md'
[]

[ambient-convection]
  type = 'Exodiff'
  input = 2d-rc-ambient-convection.i
  exodiff = 2d-rc-ambient-convection_out.e
  requirement = 'The system shall be able to model ambient volumetric convection in a channel.'
  issues = '#16948'
  ad_indexing_type = 'global'
  design = 'NSFVEnergyAmbientConvection.md'
[]

[tris]
  type = 'Exodiff'
  input = 2d-rc-no-slip.i
  exodiff = 2d-rc-no-slip-tris.e
  cli_args = 'Mesh/gen/elem_type=TRI3 GlobalParams/two_term_boundary_expansion=false Outputs/file_base=2d-rc-no-slip-tris'
  ad_indexing_type = 'global'
  detail = 'two-dimensional triangular elements'
[]

[tets]
  type = 'Exodiff'
  input = 3d-rc-no-slip.i
  exodiff = 3d-rc-no-slip_out.e
  ad_indexing_type = 'global'
  detail = 'three-dimensional tetrahedral elements'
  mumps = true
  heavy = true
[]

[3d-rc]
  type = 'Exodiff'
  input = 3d-rc.i
  exodiff = 3d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to model free-slip conditions in a 3D square channel; specifically the tangential velocity shall have a uniform value of unity and the pressure shall not change.'
  ad_indexing_type = 'global'
[]

[exo]
  type = 'Exodiff'
  input = 'lid-driven.i'
  exodiff = 'lid-driven_out.e'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity using the finite volume method.'
  ad_indexing_type = 'global'
[]

[nonsingular]
  type = 'RunApp'
  input = 'lid-driven.i'
  cli_args = 'Mesh/gen/nx=4 Mesh/gen/ny=4 Outputs/exodus=false -pc_type svd -pc_svd_monitor'
  expect_out = '0 of 49 singular values'
  absent_out = '\s+[1-9]+[0-9]* of 49 singular values'
  requirement = 'The system shall be able to solve an incompressible Navier-Stokes problem with dirichlet boundary conditions for all the normal components of velocity, using the finite volume method, and have a nonsingular system matrix.'
  ad_indexing_type = 'global'
[]

[jacobian]
  type = 'PetscJacobianTester'
  input = 'lid-driven.i'
  cli_args = 'Mesh/gen/nx=2 Mesh/gen/ny=2'
  requirement = 'The system shall be able to compute a perfect Jacobian when solving a lid-driven incompressible Navier-Stokes problem with the finite volume method.'
  run_sim = True
  ad_indexing_type = 'global'
  ratio_tol = 1e-7
  difference_tol = 1e-7
[]

[with-temp]
  type = 'Exodiff'
  input = 'lid-driven-with-energy.i'
  exodiff = 'lid-driven-with-energy_out.e'
  requirement = 'The system shall be able to transport scalar quantities using the simultaneously calculated velocity field from the incompressible Navier Stokes equations.'
  method = '!dbg'
  ad_indexing_type = 'global'
  cli_args = '-pc_type lu -pc_factor_shift_type NONZERO'
[]

[transient-with-temp]
  type = 'Exodiff'
  input = 'transient-lid-driven-with-energy.i'
  exodiff = 'transient-lid-driven-with-energy_out.e'
  requirement = 'The system shall be able to transport scalar quantities using the simultaneously calculated velocity field from the transient incompressible Navier Stokes equations.'
  method = '!dbg'
  ad_indexing_type = 'global'
  cli_args = '-pc_type lu -pc_factor_shift_type NONZERO'
[]

[capped_mixing_length]
  type = 'Exodiff'
  input = mixing_length_eddy_viscosity.i
  exodiff = mixing_length_eddy_viscosity_out.e
  method = "!dbg"
  requirement = 'The system shall be able to compute the turbulent viscosity based on the capped mixing length model and store it in a variable.'
  ad_indexing_type = 'global'
  issues = '#18666'
  design = 'rans_theory.md INSFVMixingLengthTurbulentViscosityAux.md'
  mesh_mode = REPLICATED
  recover = false
[]

[capped_mixing_length]
  type = 'Exodiff'
  input = mixing_length_total_viscosity.i
  exodiff = mixing_length_total_viscosity_out.e
  method = "!dbg"
  requirement = 'The system shall be able to calculate the material property comprising the total turbulent viscosity, based on the capped mixing length model.'
  ad_indexing_type = 'global'
  issues = '#18666'
  design = 'rans_theory.md MixingLengthTurbulentViscosityMaterial.md'
  mesh_mode = REPLICATED
  recover = false
[]

[average]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DAverage
  method = '!dbg'
  requirement = 'The system shall be able to solve a problem with channel-flow like boundary conditions in the coordinate system with an average interpolation for the velocity and demonstrate second order convergence in the velocity variables and first order convergence in the pressure variable.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[rc]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DRC
  method = '!dbg'
  min_parallel = 16
  requirement = 'The system shall be able to solve a problem with channel-flow like boundary conditions in the coordinate system with a Rhie-Chow interpolation for the velocity and demonstrate second order convergence in the velocity and pressure variables.'
  heavy = true
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  valgrind = 'none'
[]

[average-with-temp]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DAverageTemp
  method = '!dbg'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system, including energy, using an average interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable other than the pressure which shall demonstrate first order convergence.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[rc-with-temp]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DRCTemp
  method = '!dbg'
  min_parallel = 16
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system, including energy, using a RC interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable.'
  heavy = true
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  valgrind = 'none'
  max_time = 600
[]

[plane-poiseuille-average]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleAverage
  method = '!dbg'
  requirement = 'The system shall demonstrate global second order convergence for all variables on a rotated mesh when using an average interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 8
  heavy = true
[]

[plane-poiseuille-rc]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleRC
  method = '!dbg'
  requirement = 'The system shall demonstrate global second order convergence for all variables on a rotated mesh when using an RC interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 16
  heavy = true
[]

[plane-poiseuille-average-first]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleAverageFirst
  method = '!dbg'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 16
  heavy = true
  requirement = 'The system shall demonstrate global second order convergence for velocity variables and first order convergence for the pressure variable on a rotated mesh when using an average interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.'
[]

[plane-poiseuille-rc-first]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleRCFirst
  method = '!dbg'
  requirement = 'The system shall demonstrate global second order convergence for all variables on a rotated mesh when using an RC interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 16
  heavy = true
[]

[1d-average]
  type = PythonUnitTest
  input = test.py
  test_case = Test1DAverage
  method = '!dbg'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in one dimension with prescribed inlet velocity and outlet pressure and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using an average interpolation scheme for the velocity.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[2d-average]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DAverage
  method = '!dbg'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in two dimensions with prescribed inlet velocity and outlet pressure, free slip along the walls, and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using an average interpolation scheme for the velocity.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[2d-rc]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DRC
  method = '!dbg'
  min_parallel = 16
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in two dimensions with prescribed inlet velocity and outlet pressure, free slip along the walls, and implicit zero gradient boundary conditions elsewhere, and demonstrate second order convergence in both velocity and pressure when using a Rhie-Chow interpolation scheme for the velocity.'
  heavy = 'true'
  valgrind = 'none'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  max_time = 600
[]

[plane-poiseuille-average]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleAverage
  method = '!dbg'
  requirement = 'The system shall demonstrate global second order convergence for all variables when using an average interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 8
[]

[plane-poiseuille-rc]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleRC
  method = '!dbg'
  requirement = 'The system shall demonstrate global second order convergence for all variables when using an RC interpolation for the velocity and a two term Taylor series expansion for face values on non-Dirichlet boundaries.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 16
  heavy = true
[]

[plane-poiseuille-average-first]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleAverageFirst
  method = '!dbg'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 16
  heavy = true
  requirement = 'The system shall demonstrate global second order convergence for all variables when using an average interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.'
[]

[plane-poiseuille-rc-first]
  type = PythonUnitTest
  input = test.py
  test_case = TestPlanePoiseuilleRCFirst
  method = '!dbg'
  requirement = 'The system shall demonstrate global second order convergence for all variables when using an RC interpolation for the velocity and a one term Taylor series expansion for face values on non-Dirichlet boundaries.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  min_parallel = 16
  heavy = true
[]

[2d-average-with-temp]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DAverageTemp
  method = '!dbg'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations, including energy, using an average interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable.'
  heavy = 'true'
  valgrind = 'none'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[2d-rc-with-temp]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DRCTemp
  method = '!dbg'
  min_parallel = 16
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations, including energy, using a Rhie-Chow interpolation for the velocity, with a mix of Dirichlet and zero-gradient boundary conditions for each variable, and demonstrate second order convergence for each variable.'
  heavy = 'true'
  valgrind = 'none'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[2d-rc-diri]
  type = PythonUnitTest
  input = test.py
  test_case = Test2DRC
  method = '!dbg'
  min_parallel = 2
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in 2D cylindrical coordinates, using a Rhie-Chow scheme, dirichlet boundary conditions for both variables, and demonstrate second order convergence for the velocity and pressure.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  heavy = true
[]

[rc]
  type = PythonUnitTest
  input = test.py
  test_case = TestRC
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations using a Rhie-Chow interpolation scheme and produce second order convergence for all variables.'
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[capped_mixing_length]
  type = 'Exodiff'
  input = capped_mixing_length.i
  exodiff = capped_mixing_length_out.e
  method = "!dbg"
  requirement = 'The system shall be able to compute the turbulent viscosity based on the capped mixing length model.'
  ad_indexing_type = 'global'
  issues = '#18666'
  design = 'rans_theory.md WallDistanceMixingLengthAux.md'
  mesh_mode = REPLICATED
  recover = false
[]

[fluid_only]
  type = 'Exodiff'
  input = 2d-rc-heated.i
  exodiff = rc-heated.e
  method = "!dbg"
  abs_zero = 1e-9
  detail = 'with a Cartesian geometry, only modeling the fluid phase,'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-heated'"
[]

[rz]
  type = 'Exodiff'
  input = 2d-rc-heated.i
  exodiff = rc-heated-rz.e
  method = "!dbg"
  abs_zero = 1e-9
  detail = 'in rz geometry,'
  ad_indexing_type = 'global'
  cli_args = "Problem/coord_type='rz' Problem/rz_coord_axis='x' Outputs/file_base='rc-heated-rz'"
[]

[kappa]
  type = 'Exodiff'
  input = 2d-rc-heated-effective.i
  exodiff = rc-heated-effective.e
  method = "!dbg"
  abs_zero = 1e-9
  detail = 'with an effective diffusion coefficient,'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-heated-effective'"
[]

[solid-fluid]
  type = 'Exodiff'
  input = 2d-rc-heated.i
  exodiff = rc-heated-solid-fluid.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#16756'
  requirement = 'The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature.'
  ad_indexing_type = 'global'
  design = 'PINSFVEnergyAmbientConvection.md'
  cli_args = "Variables/inactive='' AuxVariables/inactive='temp_solid' FVKernels/inactive='' FVBCs/inactive='' Outputs/file_base='rc-heated-solid-fluid'"
[]

[boussinesq]
  type = 'Exodiff'
  input = 2d-rc-heated-boussinesq.i
  exodiff = rc-heated-boussinesq.e
  method = "!dbg"
  issues = '#16756'
  requirement = 'The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel with a Boussinesq approximation for the influence of temperature on density.'
  ad_indexing_type = 'global'
  design = 'PINSFVMomentumBoussinesq.md'
  cli_args = "Outputs/file_base='rc-heated-boussinesq'"
[]

[smooth-jump]
  type = 'Exodiff'
  input = 2d-rc-epsjump.i
  exodiff = rc-smooth-epsjump.e
  method = "!dbg"
  requirement = 'The system shall be able to model a smooth porosity gradient in a 2D channel.'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-smooth-epsjump' ICs/inactive='porosity_1 porosity_2' FVKernels/inactive='' GlobalParams/smooth_porosity=true"
[]

[1d-discontinuous-jump-average-average]
  type = 'Exodiff'
  input = 1d-rc-epsjump.i
  exodiff = rc-discontinuous-1d-epsjump-average-average.e
  requirement = 'The system shall be able to model a discontinuous porosity jump in a 1D channel with average interpolation of velocity and advected quantity.'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-discontinuous-1d-epsjump-average-average' velocity_interp_method='average' advected_interp_method='average'"
[]

[1d-discontinuous-jump-average-upwind]
  type = 'Exodiff'
  input = 1d-rc-epsjump.i
  exodiff = rc-discontinuous-1d-epsjump-average-upwind.e
  requirement = 'The system shall be able to model a discontinuous porosity jump in a 1D channel with average interpolation of velocity and upwinding of the advected quantity.'
  ad_indexing_type = 'global'
  max_parallel = 1
  cli_args = "Outputs/file_base='rc-discontinuous-1d-epsjump-average-upwind' velocity_interp_method='average' advected_interp_method='upwind' Executioner/nl_abs_tol=1e-13"
[]

[1d-discontinuous-jump-rc-average]
  type = 'Exodiff'
  input = 1d-rc-epsjump.i
  exodiff = rc-discontinuous-1d-epsjump-rc-average.e
  requirement = 'The system shall be able to model a discontinuous porosity jump in a 1D channel with Rhie Chow interpolation of velocity and averaging of the advected quantity.'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-discontinuous-1d-epsjump-rc-average' velocity_interp_method='rc' advected_interp_method='average' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[1d-discontinuous-jump-rc-upwind]
  type = 'Exodiff'
  input = 1d-rc-epsjump.i
  exodiff = rc-discontinuous-1d-epsjump-rc-upwind.e
  requirement = 'The system shall be able to model a discontinuous porosity jump in a 1D channel with Rhie Chow interpolation of velocity and upwinding of the advected quantity.'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-discontinuous-1d-epsjump-rc-upwind' velocity_interp_method='rc' advected_interp_method='upwind' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[discontinuous-jump]
  type = 'Exodiff'
  input = 2d-rc-epsjump.i
  exodiff = rc-discontinuous-epsjump.e
  method = "!dbg"
  requirement = 'The system shall be able to model a discontinuous porosity jump in a 2D channel.'
  recover = false
  max_parallel = 1
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-discontinuous-epsjump'"
[]

[free-slip]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc.e
  method = "!dbg"
  requirement = 'The system shall be able to model free-slip conditions in a porous media channel; specifically the tangential velocity shall have a uniform value of unity, the normal velocity shall have a uniform value of zero, and the pressure shall not change.'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc'"
[]

[free-slip-rz]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-rz.e
  method = "!dbg"
  requirement = 'The system shall be able to model free-slip conditions in a porous media cylindrical channel; specifically the tangential velocity shall have a uniform value of unity, the normal velocity shall have a uniform value of zero, and the pressure shall not change.'
  ad_indexing_type = 'global'
  cli_args = "Problem/coord_type='rz' Problem/rz_coord_axis='x' Outputs/file_base='rc-rz'"
[]

[no-slip]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "FVBCs/active='inlet-u inlet-v outlet-p no-slip-u no-slip-v' Outputs/file_base='rc-no-slip' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[no-slip-pressure-driven]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip-delta-pressure.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel with flow driven by a pressure differential; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "FVBCs/active='inlet-p outlet-p no-slip-u no-slip-v' Outputs/file_base='rc-no-slip-delta-pressure'"
[]

[no-slip-pressure-average]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip-mean-pressure.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel with a set mean pressure; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "Variables/inactive='' FVKernels/inactive='' FVBCs/active='inlet-u inlet-v outlet-p-novalue outlet-u outlet-v no-slip-u no-slip-v' Outputs/file_base='rc-no-slip-mean-pressure'"
[]

[no-slip-average]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip-average-velocity.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel using an average interpolation for velocity; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "velocity_interp_method='average' FVBCs/active='inlet-u inlet-v outlet-p no-slip-u no-slip-v' Outputs/file_base='rc-no-slip-average-velocity' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[no-slip-match]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip-match-insfv.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel with a porosity of 1; specifically, it should match a regular INSFV simulation results.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "AuxVariables/porosity/initial_condition=1 FVBCs/active='inlet-u inlet-v outlet-p no-slip-u no-slip-v' Outputs/file_base='rc-no-slip-match-insfv' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[symmetry]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip-symmetry.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel with reflective boundary conditions on one side; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "FVBCs/no-slip-u/boundary='top' FVBCs/no-slip-v/boundary='top' FVBCs/active='inlet-u inlet-v outlet-p no-slip-u no-slip-v symmetry-u symmetry-v symmetry-p' Outputs/file_base='rc-no-slip-symmetry'"
[]

[symmetry-rz]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = rc-no-slip-rz-symmetry.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a cylindrical porous media channel with reflective boundary conditions on one side; specifically, moving down the channel, the tangential velocity shall develop a parabolic profile.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "Problem/coord_type='rz' Problem/rz_coord_axis='x' FVBCs/no-slip-u/boundary='top' FVBCs/no-slip-v/boundary='top' FVBCs/active='inlet-u inlet-v outlet-p no-slip-u no-slip-v symmetry-u symmetry-v symmetry-p' Outputs/file_base='rc-no-slip-rz-symmetry' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[friction]
  type = 'Exodiff'
  input = 2d-rc-friction.i
  exodiff = rc-friction.e
  method = "!dbg"
  requirement = 'The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "FVBCs/active='inlet-u inlet-v outlet-p no-slip-u no-slip-v symmetry-u symmetry-v symmetry-p' Outputs/file_base='rc-friction' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[1D_continuous_porosity]
  type = PythonUnitTest
  input = test.py
  test_case = TestRC
  requirement = 'The system shall be able to solve the incompressible porous flow Navier-Stokes equations using a Rhie-Chow interpolation scheme in a 1D channel with a continuously varying porosity and produce second order convergence for all variables.'
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[2D_continuous_porosity]
  type = PythonUnitTest
  input = test.py
  test_case = TestRC_2D
  requirement = 'The system shall be able to solve the incompressible porous flow Navier-Stokes equations using a Rhie-Chow interpolation scheme in a 2D channel with a continuously varying porosity and produce second order convergence for all variables.'
  method = '!dbg'
  min_parallel = 16
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[rc]
  type = PythonUnitTest
  input = test.py
  test_case = TestRC
  requirement = 'The system shall be able to solve the incompressible porous flow Navier-Stokes equations in a 1D channel using a Rhie-Chow interpolation scheme and produce second order convergence for all variables.'
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[rc-2d]
  type = PythonUnitTest
  input = test.py
  test_case = TestRC_2D
  requirement = 'The system shall be able to solve the incompressible porous flow Navier-Stokes equations in a 2D channel using a Rhie-Chow interpolation scheme and produce second order convergence for all variables.'
  method = '!dbg'
  min_parallel = 16
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[ns_ics]
  type = 'Exodiff'
  input = 'test.i'
  exodiff = 'test_out.e'
  requirement = 'The system shall be able to set initial conditions for fluid flow variables.'
  design = 'NSInitialCondition.md'
[]

[cns_ics]
  type = 'Exodiff'
  input = 'pns_test.i'
  exodiff = 'pns_test_out.e'
  requirement = 'The system shall be able to set intial conditions for porous flow variables.'
  design = 'PNSInitialCondition.md'
[]

[./RZ_cone_no_parts]
  type = 'Exodiff'
  input = 'RZ_cone_no_parts.i'
  exodiff = 'RZ_cone_no_parts_out.e'
  issues = '#7651'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system while not integrating the pressure term by parts.'
[../]

[./RZ_cone_by_parts]
  type = 'Exodiff'
  input = 'RZ_cone_by_parts.i'
  exodiff = 'RZ_cone_by_parts_out.e'
  issues = '#7651'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in an RZ coordinate system while integrating the pressure term by parts.'
[../]

[./high_re]
  type = 'Exodiff'
  input = 'RZ_cone_high_reynolds.i'
  exodiff = 'steady_linear.e'
  cli_args = 'Outputs/out/file_base=steady_linear'
  issues = '#7651'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations for a high Reynolds number in an RZ coordinate system.'
[../]

[./jac]
  type = 'PetscJacobianTester'
  input = 'RZ_cone_stab_jac_test.i'
  ratio_tol = 1e-7
  difference_tol = 1e-5
  issues = '#7651'
  requirement = 'The system shall be able to compute an accurate Jacobian for the incompressible Navier-Stokes equations in an RZ coordinate system.'
[../]

[ad_rz_cone_by_parts]
  type = 'Exodiff'
  input = 'ad_rz_cone_by_parts.i'
  exodiff = 'ad_rz_cone_by_parts_out.e'
  requirement = 'The system shall be able to solve the transient incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and reproduce the results of a hand-coded Jacobian implementation.'
  issues = '#14901'
  custom_cmp = 'rz-cone.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[ad_rz_cone_no_parts]
  type = 'Exodiff'
  input = 'ad_rz_cone_no_parts.i'
  exodiff = 'ad_rz_cone_no_parts_out.e'
  requirement = 'The system shall be able to solve the transient incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts, using a traction form for the viscous term, and using a no-bc boundary condition, and reproduce the results of a hand-coded Jacobian implementation.'
  issues = '#14901'
  custom_cmp = 'rz-cone.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[ad_rz_cone_no_parts_steady]
  type = 'Exodiff'
  input = 'ad_rz_cone_no_parts_steady.i'
  exodiff = 'ad_rz_cone_no_parts_steady_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.'
  issues = '#14901'
[]

[ad_rz_cone_by_parts_steady]
  type = 'Exodiff'
  input = 'ad_rz_cone_by_parts_steady.i'
  exodiff = 'ad_rz_cone_by_parts_steady_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition'
  issues = '#14901'
[]

[ad_rz_cone_no_parts_steady_nobcbc]
  type = 'Exodiff'
  input = 'ad_rz_cone_no_parts_steady_nobcbc.i'
  exodiff = 'ad_rz_cone_no_parts_steady_nobcbc_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a NoBC outflow boundary condition.'
  issues = '#14901'
[]

[ad_rz_cone_by_parts_steady_nobcbc]
  type = 'Exodiff'
  input = 'ad_rz_cone_by_parts_steady_nobcbc.i'
  exodiff = 'ad_rz_cone_by_parts_steady_nobcbc_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a NoBC outflow boundary condition'
  issues = '#14901'
[]

[rz_cone_no_parts_steady]
  type = 'Exodiff'
  input = 'rz_cone_no_parts_steady.i'
  exodiff = 'rz_cone_no_parts_steady_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  issues = '#14901'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[rz_cone_by_parts_steady]
  type = 'Exodiff'
  input = 'rz_cone_by_parts_steady.i'
  exodiff = 'rz_cone_by_parts_steady_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  issues = '#14901'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[rz_cone_no_parts_steady_nobcbc]
  type = 'Exodiff'
  input = 'rz_cone_no_parts_steady_nobcbc.i'
  exodiff = 'rz_cone_no_parts_steady_nobcbc_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a NoBC outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  issues = '#14901'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[rz_cone_by_parts_steady_nobcbc]
  type = 'Exodiff'
  input = 'rz_cone_by_parts_steady_nobcbc.i'
  exodiff = 'rz_cone_by_parts_steady_nobcbc_out.e'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equations using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a NoBC outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  issues = '#14901'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[ad_rz_cone_no_parts_steady_supg_pspg]
  type = 'Exodiff'
  input = 'ad_rz_cone_no_parts_steady_stabilized.i'
  exodiff = 'ad_rz_cone_no_parts_steady_stabilized_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.'
[]

[rz_cone_no_parts_steady_supg_pspg]
  type = 'Exodiff'
  input = 'rz_cone_no_parts_steady_stabilized.i'
  exodiff = 'rz_cone_no_parts_steady_stabilized_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[ad_rz_cone_by_parts_steady_supg_pspg]
  type = 'Exodiff'
  input = 'ad_rz_cone_by_parts_steady_stabilized.i'
  exodiff = 'ad_rz_cone_by_parts_steady_stabilized_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition.'
[]

[rz_cone_by_parts_steady_supg_pspg]
  type = 'Exodiff'
  input = 'rz_cone_by_parts_steady_stabilized.i'
  exodiff = 'rz_cone_by_parts_steady_stabilized_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[ad_rz_cone_no_parts_steady_supg_pspg_second_order]
  # TODO: update the gold file once second order derivatives are implemented in libMesh for TypeNTensors
  type = 'Exodiff'
  input = 'ad_rz_cone_no_parts_steady_stabilized_second_order.i'
  exodiff = 'ad_rz_cone_no_parts_steady_stabilized_second_order_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a second order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.'
[]

[rz_cone_no_parts_steady_supg_pspg_second_order]
  type = 'Exodiff'
  input = 'rz_cone_no_parts_steady_stabilized_second_order.i'
  exodiff = 'rz_cone_no_parts_steady_stabilized_second_order_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a second order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while not integrating the pressure term by parts and applying a natural outflow boundary condition.'
[]

[ad_rz_cone_by_parts_steady_supg_pspg_second_order]
  # TODO: update the gold file once second order derivatives are implemented in libMesh for TypeNTensors
  type = 'Exodiff'
  input = 'ad_rz_cone_by_parts_steady_stabilized_second_order.i'
  exodiff = 'ad_rz_cone_by_parts_steady_stabilized_second_order_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a second order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition.'
[]

[rz_cone_by_parts_steady_supg_pspg_second_order]
  type = 'Exodiff'
  input = 'rz_cone_by_parts_steady_stabilized_second_order.i'
  exodiff = 'rz_cone_by_parts_steady_stabilized_second_order_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a second order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts and applying a natural outflow boundary condition.'
[]

[ad_jac]
  type = 'PetscJacobianTester'
  input = 'ad_rz_cone_stab_jac_test.i'
  ratio_tol = 1e-7
  difference_tol = 1e-5
  run_sim = True
  issues = '#14901'
  requirement = 'The system shall compute an accurate Jacobian using automatic differentiation when solving the incompressible Navier Stokes equations in an axisymmetric coordinate system with SUPG and PSPG stabilization'
[../]

[ad_rz_cone_by_parts_traction_steady_supg_pspg]
  type = 'Exodiff'
  input = 'ad_rz_cone_by_parts_traction_steady_stabilized.i'
  exodiff = 'ad_rz_cone_by_parts_traction_steady_stabilized_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts, using a traction form for the viscous term, and applying a natural outflow boundary condition.'
[]

[rz_cone_by_parts_traction_steady_supg_pspg]
  type = 'Exodiff'
  input = 'rz_cone_by_parts_traction_steady_stabilized.i'
  exodiff = 'rz_cone_by_parts_traction_steady_stabilized_out.e'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using a hand-coded Jacobian, standard variable implementation in an RZ coordinate system while integrating the pressure term by parts, using a traction form for the viscous term, and applying a natural outflow boundary condition and reproduce the results of the AD, vector variable implementation.'
  custom_cmp = 'rz-cone-with-globals.cmp' # To ignore the vector variable components automatically generated by libMesh/ParaView
[]

[ad_rz_cone_by_parts_traction_steady_supg_pspg_jac]
  type = 'PetscJacobianTester'
  difference_tol = 1e-6
  run_sim = True
  input = 'ad_rz_cone_by_parts_traction_steady_stabilized.i'
  issues = '#14901'
  requirement = 'The system shall be able to solve the steady incompressible Navier-Stokes equationswith SUPG and PSPG stabilization and a first order velocity basis using an automatic differentiation, vector variable implementation in an RZ coordinate system while integrating the pressure term by parts, using a traction form for the viscous term, and applying a natural outflow boundary condition and obtain a perfect Jacobian.'
[]

[two-mats-two-eqn-sets]
  type = Exodiff
  input = two-mats-two-eqn-sets.i
  exodiff = two-mats-two-eqn-sets_out.e
  requirement = 'The system shall be able to solve two different kernel sets with two different material domains.'
  ad_mode = 'sparse'
[]

[one-mat-two-eqn-sets]
  type = Exodiff
  input = one-mat-two-eqn-sets.i
  exodiff = one-mat-two-eqn-sets_out.e
  requirement = 'The system shall be able to solve two different kernel sets within one material domain.'
[]

[two-mats-one-eqn-set]
  type = Exodiff
  input = two-mats-one-eqn-set.i
  exodiff = two-mats-one-eqn-set_out.e
  requirement = 'The system shall be able to solve one kernel set with two different material domains.'
[]

[1e3]
  type = Exodiff
  input = benchmark.i
  exodiff = 1e3.e
  cli_args = 'rayleigh=1e3 Outputs/file_base=1e3'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e3.'
  valgrind = 'none'
  heavy = True
[]

[1e4]
  type = Exodiff
  input = benchmark.i
  exodiff = 1e4.e
  cli_args = 'rayleigh=1e4 Outputs/file_base=1e4'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e4.'
  valgrind = 'none'
  heavy = True
[]

[1e5]
  type = Exodiff
  input = benchmark.i
  exodiff = 1e5.e
  cli_args = 'rayleigh=1e5 Outputs/file_base=1e5'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e5.'
  valgrind = 'none'
  heavy = True
[]

[1e6]
  type = Exodiff
  input = benchmark.i
  exodiff = 1e6.e
  cli_args = 'rayleigh=1e6 Outputs/file_base=1e6'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e6.'
  abs_zero = 1e-9
  valgrind = 'none'
  heavy = True
[]

[exo]
  type = Exodiff
  input = boussinesq_square.i
  exodiff = boussinesq_square_out.e
  requirement = 'The system shall be able to simulate natural convection by adding the Boussinesq approximation to the incompressible Navier-Stokes equations.'
[]

[threaded_exo]
  type = Exodiff
  input = boussinesq_square.i
  exodiff = boussinesq_square_out.e
  requirement = 'The system shall be able to solve mass, momentum, and energy incompressible Navier-Stokes equations with multiple threads.'
  min_threads = 2
  issues = '#15713'
  recover = false # Already testing recover in 'exo' and if we do this we'll get errors about recovering with a different number of threads. Also unsure if we support threaded recovery
[]

[jac]
  type = PetscJacobianTester
  input = boussinesq_square.i
  run_sim = True
  difference_tol = 1e-6
  ratio_tol = 1e-7
  cli_args = "Mesh/gen/nx=4 Mesh/gen/ny=4 Outputs/active=''"
  requirement = 'The system shall have an accurate Jacobian provided by automatic differentiation when computing the Boussinesq approximation.'
[]

[exo_stab]
  type = Exodiff
  input = boussinesq_stabilized.i
  exodiff = boussinesq_stabilized_out.e
  requirement = 'The system shall be able to support SUPG and PSPG stabilization of the incompressible Navier Stokes equations including the Boussinesq approximation.'
[]

[threaded_exo_stab]
  type = Exodiff
  input = boussinesq_stabilized.i
  exodiff = boussinesq_stabilized_out.e
  requirement = 'The system shall be able to solve stablized mass, momentum, and energy incompressible Navier-Stokes equations with multiple threads.'
  min_threads = 2
  issues = '#15713'
  recover = false # Already testing recover in 'exo_stab' and if we do this we'll get errors about recovering with a different number of threads. Also unsure if we support threaded recovery
[]

[jac_stab]
  type = PetscJacobianTester
  input = boussinesq_stabilized.i
  run_sim = True
  difference_tol = 1e-5
  ratio_tol = 3e-7
  cli_args = "Mesh/gen/nx=8 Mesh/gen/ny=8 Outputs/active='' -snes_test_err 1e-9"
  requirement = 'The system shall have an accurate Jacobian provided by automatic differentiation when computing the Boussinesq approximation with SUPG and PSPG stabilization.'
  method = "!dbg" # this problem is highly nonlinear and takes ~21 nonlinear iterations to solve in serial. Hence this is a relatively long test
[]

[exo_stab_action]
  type = 'Exodiff'
  input = 'boussinesq_stabilized_action.i'
  exodiff = 'boussinesq_stabilized_out.e'
  prereq = 'exo_stab'
  issues = '#15159'
  requirement = 'The system shall be able to reproduce results of incompressible Navier-Stokes with Boussinesq approximation using a customized and condensed action syntax.'
  recover = false # We use the same output as exo_stab
[]

[threaded_exo_stab_action]
  type = 'Exodiff'
  input = 'boussinesq_stabilized_action.i'
  exodiff = 'boussinesq_stabilized_out.e'
  requirement = 'The system shall be able to solve mass, momentum, and energy incompressible Navier-Stokes equations with a custom action syntax using multiple threads.'
  min_threads = 2
  issues = '#15713'
  recover = false # We use the same output as exo_stab and also unsure if we support threaded recovery
[]

[steady]
  type = Exodiff
  input = steady.i
  exodiff = steady_out.e
  abs_zero = 1e-8
  requirement = 'The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a coupled variable.'
[]

[steady-jac]
  type = PetscJacobianTester
  input = steady.i
  cli_args = 'Mesh/gen/nx=5 Mesh/gen/ny=5'
  run_sim = True
  ratio_tol = 1e-7
  difference_tol = 1e-6
  requirement = 'The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a coupled variable.'
[]

[steady-function]
  type = Exodiff
  input = steady-function.i
  exodiff = steady-function_out.e
  abs_zero = 1e-8
  requirement = 'The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a vector function.'
[]

[steady-function-jac]
  type = PetscJacobianTester
  input = steady-function.i
  cli_args = 'Mesh/gen/nx=5 Mesh/gen/ny=5'
  run_sim = True
  ratio_tol = 1e-7
  difference_tol = 1e-6
  requirement = 'The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a vector function.'
[]

[steady-action]
  type = Exodiff
  input = steady-action.i
  exodiff = steady-action_out.e
  abs_zero = 1e-8
  requirement = 'The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a coupled variable, with the problem setup through automatic action syntax.'
[]

[steady-action-jac]
  type = PetscJacobianTester
  input = steady-action.i
  cli_args = 'Mesh/gen/nx=5 Mesh/gen/ny=5'
  run_sim = True
  ratio_tol = 1e-7
  difference_tol = 1e-6
  requirement = 'The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a coupled variable, with the problem setup through automatic action syntax.'
[]

[steady-action-function]
  type = Exodiff
  input = steady-action-function.i
  exodiff = steady-action-function_out.e
  abs_zero = 1e-8
  requirement = 'The system shall be able to apply an external force to the incompressible Navier-Stokes momentum equation through a vector function, with the problem setup through automatic action syntax.'
[]

[steady-action-function-jac]
  type = PetscJacobianTester
  input = steady-action-function.i
  cli_args = 'Mesh/gen/nx=5 Mesh/gen/ny=5'
  run_sim = True
  ratio_tol = 1e-7
  difference_tol = 1e-6
  requirement = 'The system shall be able to compute an accurate Jacobian when applying an external force to the incompressible Navier-Stokes momentum equation through a vector function, with the problem setup through automatic action syntax.'
[]

[gravity-object]
  type = Exodiff
  input = gravity-object.i
  exodiff = 'gravity-object_out.e'
  detail = 'provided through a dedicated object,'
[]

[var-and-func]
  type = Exodiff
  input = 'gravity-through-coupled-force.i'
  exodiff = 'gravity-through-coupled-force_out.e'
  detail = 'or through a generic object that can simultaneously add multiple forces through both a coupled variable and a function.'
[]

[two-vars]
  type = Exodiff
  input = 'gravity-through-coupled-force.i'
  cli_args = "Outputs/out/file_base=two-vars Kernels/inactive='momentum_coupled_forces_var_and_func momentum_coupled_forces_two_funcs'"
  exodiff = 'two-vars.e'
  detail = 'The generic object shall also be able to compute the forces solely through multiple coupled variables,'
[]

[two-funcs]
  type = Exodiff
  input = 'gravity-through-coupled-force.i'
  cli_args = "Outputs/out/file_base=two-funcs Kernels/inactive='momentum_coupled_forces_var_and_func momentum_coupled_forces_two_vars'"
  exodiff = 'two-funcs.e'
  detail = 'or solely through multiple vector functions.'
[]

[var-and-func-action]
  type = Exodiff
  input = 'gravity-through-coupled-force-action.i'
  cli_args = "Modules/IncompressibleNavierStokes/coupled_force_var=u Modules/IncompressibleNavierStokes/coupled_force_vector_function='vector_gravity_func'"
  exodiff = 'gravity-through-coupled-force-action_out.e'
  detail = 'The system shall be able to add the generic object through an automatic action syntax and provide two forces either through a coupled variable and a function,'
[]

[two-vars-action]
  type = Exodiff
  input = 'gravity-through-coupled-force-action.i'
  cli_args = "Outputs/out/file_base=two-vars-action Modules/IncompressibleNavierStokes/coupled_force_var='u gravity'"
  exodiff = 'two-vars-action.e'
  detail = 'two coupled variables,'
[]

[two-funcs-action]
  type = Exodiff
  input = 'gravity-through-coupled-force-action.i'
  cli_args = "Outputs/out/file_base=two-funcs-action Modules/IncompressibleNavierStokes/coupled_force_vector_function='vector_func vector_gravity_func'"
  exodiff = 'two-funcs-action.e'
  detail = 'or two functions.'
[]

[steady]
  type = 'Exodiff'
  input = 'steady.i'
  exodiff = 'steady_out.e'
  requirement = 'The system shall be able to model a volumetric heat source and included it in stabilization terms.'
  design = 'INSADEnergySource.md'
  issues = '#15500'
[]

[steady-action]
  type = 'Exodiff'
  input = 'steady-action.i'
  exodiff = 'steady-action_out.e'
  requirement = 'The system shall be able to build a volumetric heat source model using an automatic action syntax.'
  design = 'INSADEnergySource.md'
  issues = '#15500'
[]

[steady-var]
  type = 'Exodiff'
  input = 'steady-var.i'
  exodiff = 'steady-var_out.e'
  requirement = 'The system shall be able to model a volumetric heat source with a coupled variable and included it in stabilization terms.'
  design = 'INSADEnergySource.md'
  issues = '#15500'
[]

[steady-var-action]
  type = 'Exodiff'
  input = 'steady-var-action.i'
  exodiff = 'steady-var-action_out.e'
  requirement = 'The system shall be able to build a volumetric heat source model, provided through a coupled variable, using an automatic action syntax.'
  design = 'INSADEnergySource.md'
  issues = '#15500'
[]

[./gravity]
  type = 'Exodiff'
  input = 'gravity.i'
  exodiff = 'gravity_out.e'
  requirement = 'The system shall be able to model the effect of gravity on incompressible flow using a finite element discretization.'
[../]

[./jacobian_test]
  type = AnalyzeJacobian
  input = jacobian_test.i
  expect_out = '\nNo errors detected. :-\)\n'
  recover = false
  petsc_version = '<3.9.0 || >=3.9.4'
  mesh_mode = REPLICATED
  issues = '#13025'
  requirement = 'The system shall compute accurate Jacobians for the incompressible Navier-Stokes equation.'
[../]

[./jacobian_stabilized_test]
  type = AnalyzeJacobian
  input = jacobian_stabilized_test.i
  expect_out = '\nNo errors detected. :-\)\n'
  recover = false
  petsc_version = '<3.9.0 || >=3.9.4'
  mesh_mode = REPLICATED
  issues = '#13025'
  requirement = 'The system shall compute accurate Jacobians for the incompressible Navier-Stokes equation with stabilization.'
[../]

[./jacobian_traction_stabilized_test]
  type = AnalyzeJacobian
  input = jacobian_traction_stabilized.i
  expect_out = '\nNo errors detected. :-\)\n'
  recover = false
  petsc_version = '<3.9.0 || >=3.9.4'
  mesh_mode = REPLICATED
  issues = '#13025'
  requirement = 'The system shall compute accurate Jacobians for the incompressible Navier-Stokes
                   equation with stabilization with a traction boundary condition.'
[../]

[./wedge_dirichlet]
  type = 'Exodiff'
  input = 'wedge_dirichlet.i'
  exodiff = 'wedge_dirichlet_out.e'
  # This test pins the pressure to zero, so it can experience some
  # spurious diffs around p=0 if we are use relative tolerances.
  custom_cmp = 'wedge_dirichlet.cmp'
  allow_test_objects = true
  detail = 'with pressure Dirichlet boundary conditions'
[../]

[./wedge_natural]
  type = 'Exodiff'
  input = 'wedge_natural.i'
  exodiff = 'wedge_natural_out.e'
  allow_test_objects = true
  detail = 'and with natural advection boundary conditions.'
[../]

[mixed]
  # Note that the best test would be a MMS test which I usually setup with a PythonUnitTest but
  # those are expensive and take some developer hours to get tuned
  type = 'Exodiff'
  input = mixed.i
  exodiff = mixed_out.e
  requirement = 'The system shall support solving a steady energy equation and transient momentum equations and apply the correct stabilization.'
  rel_err = 7e-6
[]

[jac]
  type = PetscJacobianTester
  run_sim = True
  input = mixed.i
  difference_tol = 1e-6
  ratio_tol = 1e-7
  requirement = 'The system shall support solving a steady energy equation and transient momentum equations with correct stabilization and compute a perfect Jacobian.'
  cli_args = 'Mesh/gen/nx=2 Mesh/gen/ny=2'
[]

[./lid_driven]
  type = 'Exodiff'
  input = 'lid_driven.i'
  exodiff = 'lid_driven_out.e'
  custom_cmp = 'lid_driven.cmp'
  issues = '#000'
  requirement = 'We shall be able to solve a canonical lid-driven problem without stabilization, using mixed order finite elements for velocity and pressure.'
[../]

[./ad_lid_driven]
  type = 'Exodiff'
  input = 'ad_lid_driven.i'
  exodiff = 'lid_driven_out.e'
  custom_cmp = 'lid_driven.cmp'
  prereq = 'lid_driven'
  cli_args = 'Outputs/file_base=lid_driven_out'
  method = 'OPT'
  issues = '#13025'
  requirement = 'We shall be able to reproduce the results from the hand-coded lid-driven simulation using automatic differentiation objects.'
  recover = false # Same output as lid_driven
[../]

[./ad_lid_driven_mean_zero_pressure]
  type = 'Exodiff'
  input = 'ad_lid_driven_mean_zero_pressure.i'
  exodiff = 'ad_lid_driven_mean_zero_pressure_out.e'
  method = 'OPT'
  issues = '#15549'
  requirement = 'We shall be able to run lid-dirven simulation using a global mean-zero pressure constraint approach.'
[../]

[./ad_lid_driven_jacobian]
  type = 'PetscJacobianTester'
  input = 'ad_lid_driven.i'
  cli_args = 'Outputs/exodus=false Mesh/gen/nx=3 Mesh/gen/ny=3'
  run_sim = True
  difference_tol = 1e-7
  issues = '#13025'
  requirement = 'The Jacobian for the mixed-order INS problem shall be perfect when provided through automatic differentiation.'
[../]

[./lid_driven_stabilized]
  type = 'Exodiff'
  input = 'lid_driven_stabilized.i'
  exodiff = 'lid_driven_stabilized_out.e'
  custom_cmp = 'lid_driven_stabilized.cmp'
  issues = '#9687'
  requirement = 'We shall be able to solve the lid-driven problem using equal order shape functions with pressure-stabilized petrov-galerkin stabilization. We shall also demonstrate SUPG stabilization.'
[../]

[./ad_lid_driven_stabilized]
  type = 'Exodiff'
  input = 'ad_lid_driven_stabilized.i'
  exodiff = 'lid_driven_stabilized_out.e'
  custom_cmp = 'lid_driven_stabilized.cmp'
  prereq = 'lid_driven_stabilized'
  issues = '#13025'
  requirement = 'We shall be able to reproduce the hand-coded stabilized results with automatic differentiation objects.'
  recover = false # Same output as lid_driven_stabilized
[../]

[./ad_lid_driven_stabilized_jacobian]
  type = 'PetscJacobianTester'
  input = 'ad_lid_driven_stabilized.i'
  run_sim = True
  cli_args = 'Outputs/exodus=false Mesh/gen/nx=3 Mesh/gen/ny=3'
  difference_tol = 1e-7
  issues = '#13025'
  requirement = 'The Jacobian for the automatic differentiation stabilized lid-driven problem shall be perfect.'
[../]

[./still_unstable]
  type = 'RunApp'
  input = 'lid_driven_stabilized.i'
  expect_out = 'Aborting as solve did not converge'
  cli_args = 'GlobalParams/alpha=0 Outputs/exodus=false'
  method = 'OPT'
  issues = '#9687'
  requirement = 'Simulation with equal-order shape functions without pressure stabilization shall be unstable.'
[../]

[./lid_driven_chorin]
  type = 'Exodiff'
  input = 'lid_driven_chorin.i'
  exodiff = 'lid_driven_chorin_out.e'
  issues = '#000'
  requirement = 'We shall be able to solve the INS equations using the classical Chorin splitting algorithm.'
[../]

[./lid_driven_action]
  type = 'Exodiff'
  input = 'lid_driven_action.i'
  exodiff = 'lid_driven_out.e'
  custom_cmp = 'lid_driven.cmp'
  prereq = 'lid_driven'
  issues = '#15159'
  requirement = 'The system shall be able to reproduce unstabilized incompressible Navier-Stokes results with hand-coded Jacobian using a customized and condensed action syntax.'
  recover = false # Same output as lid_driven
[../]

[./lid_driven_stabilized_action]
  type = 'Exodiff'
  input = 'lid_driven_stabilized_action.i'
  exodiff = 'lid_driven_stabilized_out.e'
  custom_cmp = 'lid_driven_stabilized.cmp'
  prereq = 'lid_driven_stabilized'
  issues = '#15159'
  requirement = 'The system shall be able to reproduce stabilized incompressible Navier-Stokes results with hand-coded Jacobian using a customized and condensed action syntax.'
  recover = false # Same output as lid_driven_stabilized
[../]

[./ad_lid_driven_action]
  type = 'Exodiff'
  input = 'ad_lid_driven_action.i'
  exodiff = 'lid_driven_out.e'
  custom_cmp = 'lid_driven.cmp'
  method = 'OPT'
  prereq = 'ad_lid_driven'
  issues = '#15159'
  requirement = 'The system shall be able to reproduce unstabilized incompressible Navier-Stokes results with auto-differentiation using a customized and condensed action syntax.'
  recover = false # Same output as lid_driven
[../]

[./ad_lid_driven_stabilized_action]
  type = 'Exodiff'
  input = 'ad_lid_driven_stabilized_action.i'
  exodiff = 'lid_driven_stabilized_out.e'
  custom_cmp = 'lid_driven_stabilized.cmp'
  prereq = 'ad_lid_driven_stabilized'
  issues = '#15159'
  requirement = 'The system shall be able to reproduce stabilized incompressible Navier-Stokes results with auto-differentiation using a customized and condensed action syntax.'
  recover = false # Same output as lid_driven_stabilized
[../]

[ad_stabilized_energy_steady]
  type = 'Exodiff'
  input = 'ad_lid_driven_stabilized_with_temp.i'
  exodiff = 'ad_lid_driven_stabilized_with_temp_out.e'
  issues = '#15500'
  requirement = 'The system shall be able to solve a steady stabilized mass/momentum/energy incompressible Navier-Stokes formulation.'
[]

[ad_stabilized_energy_transient]
  type = 'Exodiff'
  input = 'ad_lid_driven_stabilized_with_temp_transient.i'
  exodiff = 'ad_lid_driven_stabilized_with_temp_transient_out.e'
  issues = '#15500'
  requirement = 'The system shall be able to solve a transient stabilized mass/momentum/energy incompressible Navier-Stokes formulation.'
[]

[ad_stabilized_energy_steady_action]
  type = 'Exodiff'
  input = 'ad_lid_driven_action_stabilized_steady.i'
  exodiff = 'ad_lid_driven_action_stabilized_steady_out.e'
  issues = '#15500'
  requirement = 'The system shall be able to solve a steady stabilized mass/momentum/energy incompressible Navier-Stokes formulation with action syntax.'
[]

[ad_stabilized_energy_transient_action]
  type = 'Exodiff'
  input = 'ad_lid_driven_action_stabilized_transient.i'
  exodiff = 'ad_lid_driven_action_stabilized_transient_out.e'
  issues = '#15500'
  requirement = 'The system shall be able to solve a transient stabilized mass/momentum/energy incompressible Navier-Stokes formulation with action syntax.'
[]

[ad_stabilized_transient_les]
  type = 'Exodiff'
  input = 'ad_lid_driven_les.i'
  exodiff = 'ad_lid_driven_les_out.e'
  issues = '#15757'
  requirement = 'The system shall be able to solve a transient incompressible Navier-Stokes with nonlinear Smagorinsky eddy viscosity.'
[]

[./4x4]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '4x4_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=4x4_alpha_1e-6 Mesh/gen/nx=4 Mesh/gen/ny=4 GlobalParams/alpha=1e-6'
  heavy = true
  detail = '4x4,'
[../]

[./8x8]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '8x8_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=8x8_alpha_1e-6 Mesh/gen/nx=8 Mesh/gen/ny=8 GlobalParams/alpha=1e-6'
  heavy = true
  detail = '8x8,'
[../]

[./16x16]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '16x16_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=16x16_alpha_1e-6 Mesh/gen/nx=16 Mesh/gen/ny=16 GlobalParams/alpha=1e-6'
  heavy = true
  detail = '16x16,'
[../]

[./32x32]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '32x32_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=32x32_alpha_1e-6 Mesh/gen/nx=32 Mesh/gen/ny=32 GlobalParams/alpha=1e-6'
  heavy = true
  detail = 'and 32x32 mesh.'
[../]

[./4x4]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '4x4_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=4x4_alpha_1e-3 Mesh/gen/nx=4 Mesh/gen/ny=4 GlobalParams/alpha=1e-3'
  heavy = true
  detail = '4x4,'
[../]

[./8x8]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '8x8_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=8x8_alpha_1e-3 Mesh/gen/nx=8 Mesh/gen/ny=8 GlobalParams/alpha=1e-3'
  heavy = true
  detail = '8x8,'
[../]

[./16x16]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '16x16_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=16x16_alpha_1e-3 Mesh/gen/nx=16 Mesh/gen/ny=16 GlobalParams/alpha=1e-3'
  heavy = true
  detail = '16x16,'
[../]

[./32x32]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '32x32_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=32x32_alpha_1e-3 Mesh/gen/nx=32 Mesh/gen/ny=32 GlobalParams/alpha=1e-3'
  heavy = true
  detail = 'and 32x32 mesh.'
[../]

[./4x4]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '4x4_alpha_1e0.e'
  cli_args = 'Outputs/file_base=4x4_alpha_1e0 Mesh/gen/nx=4 Mesh/gen/ny=4 GlobalParams/alpha=1e0'
  detail = '4x4,'
[../]

[./8x8]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '8x8_alpha_1e0.e'
  cli_args = 'Outputs/file_base=8x8_alpha_1e0 Mesh/gen/nx=8 Mesh/gen/ny=8 GlobalParams/alpha=1e0'
  detail = '8x8,'
[../]

[./16x16]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '16x16_alpha_1e0.e'
  cli_args = 'Outputs/file_base=16x16_alpha_1e0 Mesh/gen/nx=16 Mesh/gen/ny=16 GlobalParams/alpha=1e0'
  detail = '16x16,'
[../]

[./32x32]
  type = 'Exodiff'
  input = 'pspg_mms_test.i'
  exodiff = '32x32_alpha_1e0.e'
  cli_args = 'Outputs/file_base=32x32_alpha_1e0 Mesh/gen/nx=32 Mesh/gen/ny=32 GlobalParams/alpha=1e0'
  detail = 'and 32x32 mesh.'
[../]

[./4x4]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '4x4_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=4x4_alpha_1e-6 Mesh/gen/nx=4 Mesh/gen/ny=4 GlobalParams/alpha=1e-6'
  heavy = true
  detail = '4x4,'
[../]

[./8x8]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '8x8_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=8x8_alpha_1e-6 Mesh/gen/nx=8 Mesh/gen/ny=8 GlobalParams/alpha=1e-6'
  heavy = true
  detail = '8x8,'
[../]

[./16x16]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '16x16_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=16x16_alpha_1e-6 Mesh/gen/nx=16 Mesh/gen/ny=16 GlobalParams/alpha=1e-6'
  heavy = true
  detail = '16x16,'
[../]

[./32x32]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '32x32_alpha_1e-6.e'
  cli_args = 'Outputs/file_base=32x32_alpha_1e-6 Mesh/gen/nx=32 Mesh/gen/ny=32 GlobalParams/alpha=1e-6'
  heavy = true
  detail = 'and 32x32 mesh.'
  valgrind = 'none'
[../]

[./4x4]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '4x4_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=4x4_alpha_1e-3 Mesh/gen/nx=4 Mesh/gen/ny=4 GlobalParams/alpha=1e-3'
  heavy = true
  detail = '4x4,'
[../]

[./8x8]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '8x8_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=8x8_alpha_1e-3 Mesh/gen/nx=8 Mesh/gen/ny=8 GlobalParams/alpha=1e-3'
  heavy = true
  detail = '8x8,'
[../]

[./16x16]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '16x16_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=16x16_alpha_1e-3 Mesh/gen/nx=16 Mesh/gen/ny=16 GlobalParams/alpha=1e-3'
  heavy = true
  detail = '16x16,'
[../]

[./32x32]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '32x32_alpha_1e-3.e'
  cli_args = 'Outputs/file_base=32x32_alpha_1e-3 Mesh/gen/nx=32 Mesh/gen/ny=32 GlobalParams/alpha=1e-3'
  heavy = true
  detail = 'and 32x32 mesh.'
  valgrind = 'none'
[../]

[./4x4]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '4x4_alpha_1e0.e'
  cli_args = 'Outputs/file_base=4x4_alpha_1e0 Mesh/gen/nx=4 Mesh/gen/ny=4 GlobalParams/alpha=1e0'
  detail = '4x4,'
[../]

[./8x8]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '8x8_alpha_1e0.e'
  cli_args = 'Outputs/file_base=8x8_alpha_1e0 Mesh/gen/nx=8 Mesh/gen/ny=8 GlobalParams/alpha=1e0'
  detail = '8x8,'
[../]

[./16x16]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '16x16_alpha_1e0.e'
  cli_args = 'Outputs/file_base=16x16_alpha_1e0 Mesh/gen/nx=16 Mesh/gen/ny=16 GlobalParams/alpha=1e0'
  detail = '16x16,'
[../]

[./32x32]
  type = 'Exodiff'
  input = 'supg_mms_test.i'
  exodiff = '32x32_alpha_1e0.e'
  cli_args = 'Outputs/file_base=32x32_alpha_1e0 Mesh/gen/nx=32 Mesh/gen/ny=32 GlobalParams/alpha=1e0'
  detail = 'and 32x32 mesh.'
  valgrind = 'none'
[../]

[./adv_dominated_supg_stabilized]
  type = 'Exodiff'
  input = 'supg_adv_dominated_mms.i'
  exodiff = 'adv_dominated_stabilized.e'
  cli_args = 'GlobalParams/supg=true Outputs/file_base=adv_dominated_stabilized'
  requirement = 'The system shall be able to solve high Reynolds number incompressible flow problems through use of streamline upwind Petrov-Galerkin stabilization and with a Q2Q1 discretization'
  recover = false # We only execute output on final
[../]

[./adv_dominated_supg_pspg_stabilized]
  type = 'Exodiff'
  input = 'supg_pspg_adv_dominated_mms.i'
  exodiff = 'supg_pspg_adv_dominated_mms_exodus.e'
  # This problem uses SuperLU.
  superlu = true
  requirement = 'The system shall be able to solve high Reynolds number incompressible flow problems through use of streamline upwind and pressure stabilized Petrov-Galerkin and with a Q1Q1 discretization'
  recover = false # We only execute output on final
[../]

[malloc]
  type = RunException
  input = test.i
  expect_err = 'New nonzero at.*caused a malloc'
  requirement = 'The system shall allow MOOSE applications to specify nonzero malloc behavior; for the Navier-Stokes application, new nonzero allocations shall be errors.'
  allow_test_objects = True
  executable_pattern = 'navier_stokes'
[]

[kernels]
  type = 'Exodiff'
  input = 'open_bc_pressure_BC.i'
  exodiff = 'open_bc_out_pressure_BC.e'
  custom_cmp = open_bc_out_pressure_BC.cmp
  detail = 'using the kernel formulation,'
[../]

[action]
  type = 'Exodiff'
  input = 'open_bc_pressure_BC_action.i'
  exodiff = 'open_bc_out_pressure_BC.e'
  custom_cmp = open_bc_out_pressure_BC.cmp
  prereq = open_bc_pressure_BC/kernels
  detail = 'using the action formulation'
[../]

[fieldSplit]
  type = 'Exodiff'
  input = 'open_bc_pressure_BC_fieldSplit.i'
  exodiff = 'open_bc_out_pressure_BC_fieldSplit.e'
  custom_cmp = open_bc_out_pressure_BC_fieldSplit.cmp
  # This test uses the 'selfp' option for the fieldsplit
  # preconditioner, which requires at least PETSc 3.5.0.
  petsc_version = '>=3.5.0'
  detail = 'and using a field split preconditioning.'
[../]

[./2D]
  type = 'Exodiff'
  input = 'stagnation.i'
  cli_args = 'Outputs/file_base=2D_out'
  exodiff = '2D_out.e'
  detail = 'in 2D'
[../]

[./axisymmetric]
  type = 'Exodiff'
  input = 'stagnation.i'
  # Note: We need to specify a slightly different vertical velocity
  # boundary condition in the axisymmetric case in order to have
  # the "far field" flow conditions satisfy the divergence-free
  # constraint.
  cli_args = 'Outputs/file_base=axisymmetric_out Problem/coord_type=RZ Kernels/mass/type=INSMassRZ Kernels/x_momentum_space/type=INSMomentumLaplaceFormRZ Kernels/y_momentum_space/type=INSMomentumLaplaceFormRZ Functions/vel_y_inlet/value=-2*k*y'
  exodiff = 'axisymmetric_out.e'
  detail = 'and in 2D RZ axisymmetric simulations.'
[../]

[no_parts]
  type = 'Exodiff'
  input = 'velocity_inletBC_no_parts.i'
  exodiff = 'velocity_inletBC_no_parts_out.e'
  custom_cmp = velocity_inletBC_no_parts_out.cmp
  detail = 'with the regular volumetric integration of the momentum pressure term'
[]

[by_parts]
  type = 'Exodiff'
  input = 'velocity_inletBC_by_parts.i'
  exodiff = 'velocity_inletBC_by_parts_out.e'
  custom_cmp = velocity_inletBC_by_parts_out.cmp
  detail = 'and with the momentum pressure term integrated by parts.'
[]

[steady]
  type = 'Exodiff'
  input = 'steady.i'
  exodiff = 'steady_out.e'
  requirement = 'The system shall be able to model heat transfer from ambient surroundings using a volumetric approximation.'
  design = 'INSADEnergyAmbientConvection.md'
  issues = '#15500'
[]

[steady-action]
  type = 'Exodiff'
  input = 'steady-action.i'
  exodiff = 'steady-action_out.e'
  requirement = 'The system shall be able to build a simulation, modeling heat transfer from ambient surroundings, using an automated action syntax.'
  design = 'INSADEnergyAmbientConvection.md'
  issues = '#15500'
[]

[steady-action-no-temp-var]
  type = 'Exodiff'
  input = 'steady-action.i'
  exodiff = 'steady-action-no-temp-var.e'
  cli_args = "Variables/active='temperature' Outputs/file_base=steady-action-no-temp-var"
  requirement = 'The system shall be able to add a incompressible Navier-Stokes energy/temperature equation using an action, but use a temperature variable already added in the input file.'
  design = 'INSAction.md'
  issues = '#15607'
[]

[fe]
  type = CSVDiff
  input = conservation_INSFE.i
  csvdiff = straight_INSFE.csv
  abs_zero = 1e-8
  requirement = 'The system shall be able to compute mass and momentum flow rates at internal and external boundaries of a straight channel with a finite element incompressible Navier Stokes model.'
  cli_args = "Mesh/inactive='diverging_mesh' Outputs/file_base='straight_INSFE'"
[]

[fe_diverging]
  type = CSVDiff
  input = conservation_INSFE.i
  csvdiff = diverging_INSFE.csv
  abs_zero = 1e-8
  requirement = 'The system shall be able to compute mass and momentum flow rates at internal and external boundaries of a diverging channel with a finite element incompressible Navier Stokes model.'
  cli_args = "Outputs/file_base='diverging_INSFE'"
[]

[insfv_straight]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = straight_INSFV.csv
  abs_zero = 1e-8
  requirement = 'The system shall be able to compute flow rates and prove mass, momentum and energy conservation at internal and external boundaries of a frictionless heated straight channel with a finite volume incompressible Navier Stokes model.'
  cli_args = "Mesh/inactive='diverging_mesh' Outputs/file_base='straight_INSFV' -pc_type lu -pc_factor_shift_type NONZERO"
  ad_indexing_type = 'global'
[]

[insfv_quad_xy]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = diverging_INSFV_quad_xy.csv
  abs_zero = 1e-8
  detail = 'with a quadrilateral mesh in XY geometry, with mass flow measured using either a variable or material property,'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='diverging_INSFV_quad_xy'"
[]

[insfv_quad_rz]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = diverging_INSFV_quad_rz.csv
  abs_zero = 1e-8
  detail = 'with a quadrilateral mesh in RZ geometry,'
  cli_args = "Problem/coord_type=RZ Outputs/file_base='diverging_INSFV_quad_rz'"
  ad_indexing_type = 'global'
[]

[insfv_tri_xy]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = diverging_INSFV_tri_xy.csv
  abs_zero = 1e-8
  detail = 'with a triangular mesh in XY geometry,'
  cli_args = "Mesh/diverging_mesh/file='expansion_tri.e' Outputs/file_base='diverging_INSFV_tri_xy' GlobalParams/two_term_boundary_expansion=false"
  ad_indexing_type = 'global'
[]

[insfv_quad_xy_rc]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = diverging_INSFV_quad_xy_rc.csv
  abs_zero = 1e-8
  detail = 'with Rhie Chow velocity interpolation,'
  cli_args = "velocity_interp_method='rc' Outputs/file_base='diverging_INSFV_quad_xy_rc'"
  ad_indexing_type = 'global'
[]

[insfv_quad_xy_upwind]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = diverging_INSFV_quad_xy_upwind.csv
  abs_zero = 1e-8
  detail = 'with upwind interpolation of advected quantities,'
  cli_args = "advected_interp_method='upwind' Outputs/file_base='diverging_INSFV_quad_xy_upwind'"
  ad_indexing_type = 'global'
[]

[insfv_quad_xy_noslip]
  type = CSVDiff
  input = conservation_INSFV.i
  csvdiff = diverging_INSFV_quad_xy_noslip.csv
  abs_zero = 1e-8
  detail = 'and with no-slip boundary conditions, for which momentum and energy will be dissipated at the wall.'
  cli_args = "FVBCs/inactive='free-slip-u free-slip-v' Outputs/file_base='diverging_INSFV_quad_xy_noslip'"
  ad_indexing_type = 'global'
[]

[pinsfv_straight]
  type = CSVDiff
  input = conservation_PINSFV.i
  csvdiff = straight_PINSFV.csv
  abs_zero = 1e-8
  requirement = 'The system shall be able to compute flow rates and prove mass, momentum and energy conservation at internal and external boundaries of a frictionless heated straight channel with a finite volume porous media incompressible Navier Stokes model.'
  cli_args = "Mesh/inactive='diverging_mesh' Outputs/file_base='straight_PINSFV'"
  ad_indexing_type = 'global'
[]

[pinsfv_quad_xy]
  type = CSVDiff
  input = conservation_PINSFV.i
  csvdiff = diverging_PINSFV_quad_xy.csv
  abs_zero = 1e-8
  detail = 'with a quadrilateral mesh in XY geometry, with mass flow measured using either a variable or material property,'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='diverging_PINSFV_quad_xy'"
[]

[pinsfv_quad_rz]
  type = CSVDiff
  input = conservation_PINSFV.i
  csvdiff = diverging_PINSFV_quad_rz.csv
  abs_zero = 1e-8
  detail = 'with a quadrilateral mesh in RZ geometry,'
  cli_args = "Problem/coord_type=RZ Outputs/file_base='diverging_PINSFV_quad_rz'"
  ad_indexing_type = 'global'
[]

[pinsfv_quad_xy_rc]
  type = CSVDiff
  input = conservation_PINSFV.i
  csvdiff = diverging_PINSFV_quad_xy_rc.csv
  abs_zero = 1e-8
  detail = 'with Rhie Chow velocity interpolation,'
  cli_args = "velocity_interp_method='rc' Outputs/file_base='diverging_PINSFV_quad_xy_rc'"
  ad_indexing_type = 'global'
[]

[pinsfv_quad_xy_upwind]
  type = CSVDiff
  input = conservation_PINSFV.i
  csvdiff = diverging_PINSFV_quad_xy_upwind.csv
  abs_zero = 1e-8
  detail = 'with upwind interpolation of advected quantities,'
  cli_args = "advected_interp_method='upwind' Outputs/file_base='diverging_PINSFV_quad_xy_upwind'"
  ad_indexing_type = 'global'
[]

[pinsfv_quad_xy_noslip]
  type = CSVDiff
  input = conservation_PINSFV.i
  csvdiff = diverging_PINSFV_quad_xy_noslip.csv
  abs_zero = 1e-8
  detail = 'and with no-slip boundary conditions, for which momentum and energy will be dissipated at the wall.'
  cli_args = "FVBCs/inactive='free-slip-u free-slip-v' Outputs/file_base='diverging_PINSFV_quad_xy_noslip'"
  ad_indexing_type = 'global'
[]

[./tauOpt]
  type = 'Exodiff'
  input = 'tauOpt.i'
  exodiff = 'tauOpt_out.e'
  allow_test_objects = true
  detail = 'using the optimal tau stabilization,'
[../]

[./tauMod]
  type = 'Exodiff'
  input = 'tauOpt.i'
  exodiff = 'tauMod_out.e'
  cli_args = "GlobalParams/tau_type=mod Outputs/file_base=tauMod_out"
  allow_test_objects = true
  detail = 'using the modified tau stabilization,'
[../]

[./1d_error_test_supg]
  type = 'Exodiff'
  input = 'advection_error_testing.i'
  exodiff = 'advection_error_testing_exodus.e'
  allow_test_objects = true
  detail = 'and still satisfy MMS testing in 1D'
[../]

[./2d_error_test_supg]
  type = 'Exodiff'
  exodiff = '2d_advection_error_testing_exodus.e'
  input = '2d_advection_error_testing.i'
  allow_test_objects = true
  # This problem uses SuperLU.
  superlu = true
  detail = 'and in 2D.'
[../]