Navier Stokes Requirements Traceability Matrix

[fv]
  type = Exodiff
  input = fv-thermal.i
  exodiff = fv-thermal_out.e
  detail = 'finite volume discretizations'
  ad_indexing_type = 'global'
[]

[fe]
  type = Exodiff
  input = fe-thermal.i
  exodiff = fe-thermal_out.e
  detail = 'finite element discretizations'
  ad_indexing_type = 'global'
[]

[fv]
  type = Exodiff
  input = fv.i
  exodiff = fv_out.e
  detail = 'finite volume discretizations'
  ad_indexing_type = 'global'
[]

[fe]
  type = Exodiff
  input = fe.i
  exodiff = fe_out.e
  detail = 'finite element discretizations'
  ad_indexing_type = 'global'
[]

[./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'
  requirement = 'The system shall be able to solve the Euler equations for subsonic flow with a bump in the channel.'
[../]

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

[traction]
  type = 'Exodiff'
  input = 'ad-rz-displacements.i'
  detail = 'traction form'
  exodiff = 'ad-rz-displacements-traction.e'
  cli_args = "Outputs/file_base=ad-rz-displacements-traction GlobalParams/viscous_form='traction'"
[]

[laplace]
  type = 'Exodiff'
  input = 'ad-rz-displacements.i'
  detail = 'laplace form'
  exodiff = 'ad-rz-displacements-laplace.e'
  cli_args = "Outputs/file_base=ad-rz-displacements-laplace GlobalParams/viscous_form='laplace'"
[]

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

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

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

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

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

[exo]
  type = 'Exodiff'
  input = 'pipe-flow.i'
  requirement = 'The system shall be able to solve an axisymmetric pipe flow problem using a finite element discretization in which the axis of symmetry is the x-axis.'
  exodiff = 'pipe-flow_out.e'
[]

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

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

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

[pwcnsfv]
  type = PythonUnitTest
  input = test.py
  test_case = PWCNSFV
  required_python_packages = 'pandas matplotlib'
  ad_indexing_type = 'global'
  method = "!dbg"
  requirement = 'The system shall demonstrate first order convergence rates for pressure and superficial velocity when using an upwind interpolation for advected quantities in a weakly compressible formulation of the mass and momentum Euler equations.'
  issues = '#18215'
  design = 'PINSFVMomentumAdvection.md'
[]

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

[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
  valgrind = HEAVY
[]

[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
  valgrind = HEAVY
[]

[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 PCNSFVMomentumFriction.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 PCNSFVMomentumFriction.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
  valgrind = HEAVY
  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, using functor material properties."
  ad_indexing_type = 'global'
[]

[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, using functor material properties."
  ad_indexing_type = 'global'
[]

[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, using functor material properties."
  ad_indexing_type = 'global'
[]

[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, using functor material properties."
  ad_indexing_type = 'global'
[]

[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, using functor material properties."
  ad_indexing_type = 'global'
[]

[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, using functor material properties."
  ad_indexing_type = 'global'
[]

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

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

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

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

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

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

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

[friction_flow_diode]
  type = Exodiff
  issues = '#20695'
  design = 'INSFVFrictionFlowDiode.md'
  requirement = 'The system shall be able to model a simple friction-based volumetric flow diode'
  input = 'friction.i'
  exodiff = friction_out.e
  ad_indexing_type = 'global'
[]

[momentum-inlet-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of momentum inlet types does not match the number of inlet boundaries.'
  cli_args = "Modules/NavierStokesFV/momentum_inlet_types=''"
  expect_err = "Size is not the same as the number of inlet boundaries in 'inlet_boundaries'"
  ad_indexing_type = 'global'
[]

[momentum-inlet-function-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of momentum inlet functions does not match the number of inlet boundaries times the dimension of the mesh.'
  cli_args = "Modules/NavierStokesFV/momentum_inlet_function='1'"
  expect_err = "Size is not the same as the number of direction dependent boundaries in 'inlet_boundaries' times the mesh dimension"
  ad_indexing_type = 'global'
[]

[momentum-outlet-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of momentum outlet types does not match the number of outlet boundaries.'
  cli_args = "Modules/NavierStokesFV/momentum_outlet_types=''"
  expect_err = "Size is not the same as the number of outlet boundaries in 'outlet_boundaries'"
  ad_indexing_type = 'global'
[]

[pressure-outlet-function-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of pressure outlet functions is not the same the pressure outlet boundaries.'
  cli_args = "Modules/NavierStokesFV/pressure_function=''"
  expect_err = "Size is not the same as the number of pressure outlet boundaries!"
  ad_indexing_type = 'global'
[]

[momentum-wall-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of momentum wall types does not match the number of wall boundaries.'
  cli_args = "Modules/NavierStokesFV/momentum_wall_types=''"
  expect_err = "Size is not the same as the number of wall boundaries in 'wall_boundaries'"
  ad_indexing_type = 'global'
[]

[energy-inlet-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of energy inlet types does not match the number of inlet boundaries.'
  cli_args = "Modules/NavierStokesFV/energy_inlet_types=''"
  expect_err = "Size is not the same as the number of inlet boundaries in 'inlet_boundaries'"
  ad_indexing_type = 'global'
[]

[passive-scalar-inlet-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of passive scalar inlet types does not match the number of inlet boundaries.'
  cli_args = "Modules/NavierStokesFV/passive_scalar_inlet_types=''"
  expect_err = "Size is not the same as the number of inlet boundaries in 'inlet_boundaries'"
  ad_indexing_type = 'global'
[]

[passive-scalar-inlet-function-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of passive scalar inlet functions does not match the number of inlet boundaries.'
  cli_args = "Modules/NavierStokesFV/passive_scalar_inlet_function=''"
  expect_err = "Size is not the same as the number of boundaries in 'passive_scalar_inlet_types'"
  ad_indexing_type = 'global'
[]

[energy-wall-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of energy wall types does not match the number of wall boundaries.'
  cli_args = "Modules/NavierStokesFV/energy_wall_types=''"
  expect_err = "Size is not the same as the number of wall boundaries in 'wall_boundaries'"
  ad_indexing_type = 'global'
[]

[energy-wall-function-action-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of energy wall functions does not match the number of energy wall types.'
  cli_args = "Modules/NavierStokesFV/energy_wall_function=''"
  expect_err = "Size 0 is not the same as the number of Dirichlet/Neumann conditions in 'energy_wall_types' \(2\)"
  ad_indexing_type = 'global'
[]

[scalar-ic-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of defined initial conditions is different than the number of created scalar variables.'
  cli_args = "Modules/NavierStokesFV/passive_scalar_names='scalar1 scalar2' Modules/NavierStokesFV/initial_scalar_variables='0'"
  expect_err = "The number of initial conditions \(1\) is not equal to the number of self-generated variables \(2\) !"
  ad_indexing_type = 'global'
[]

[velocity-ic-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the number of components for the initial velocity is not 3.'
  cli_args = "Modules/NavierStokesFV/initial_velocity='0 0'"
  expect_err = "The number of velocity components in the NSFVAction initial condition is not 3!"
  ad_indexing_type = 'global'
[]

[porosity-correction-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if friction correction is requested with no porous medium treatment.'
  cli_args = "Modules/NavierStokesFV/use_friction_correction='true'"
  expect_err = "This parameter should not be defined if the porous medium treatment is disabled!"
  ad_indexing_type = 'global'
[]

[porosity-scaling-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if consistent scaling is defined without using friction correction'
  cli_args = "Modules/NavierStokesFV/consistent_scaling='true'"
  expect_err = "Consistent scaling should not be defined if friction correction is disabled!"
  ad_indexing_type = 'global'
[]

[porosity-smoothing-layer-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if porosity smoothing is requested without porous medium treatment.'
  cli_args = "Modules/NavierStokesFV/porosity_smoothing_layers=4"
  expect_err = "This parameter should not be defined if the porous medium treatment is disabled!"
  ad_indexing_type = 'global'
[]

[passive-scalar-diffusivity-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user defines an inappropriate number of passive scalar diffusivities'
  cli_args = "Modules/NavierStokesFV/passive_scalar_diffusivity='0.5 0.6'"
  expect_err = "The number of diffusivities defined is not equal to the number of passive scalar fields!"
  ad_indexing_type = 'global'
[]

[passive-scalar-source-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user defines an inappropriate number of passive scalar source functions'
  cli_args = "Modules/NavierStokesFV/passive_scalar_source='0.5 0.6'"
  expect_err = "The number of external sources defined is not equal to the number of passive scalar fields!"
  ad_indexing_type = 'global'
[]

[passive-scalar-coupled-source-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user defines an inappropriate number of passive scalar coupled source functions'
  cli_args = "Modules/NavierStokesFV/passive_scalar_coupled_source='v w'"
  expect_err = "The number of coupled sources defined is not equal to the number of passive scalar fields!"
  ad_indexing_type = 'global'
[]

[passive-scalar-coupled-source-coeff-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user defines an inappropriate number of passive scalar coupled source coefficients'
  cli_args = "Modules/NavierStokesFV/passive_scalar_coupled_source='v' Modules/NavierStokesFV/passive_scalar_coupled_source_coeff='0.1 0.2'"
  expect_err = "The number of coupled sources coefficients defined is not equal to the number of passive scalar fields!"
  ad_indexing_type = 'global'
[]

[passive-scalar-schmidt-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user defines an inappropriate number of passive scalar Schmidt numbers'
  cli_args = "Modules/NavierStokesFV/passive_scalar_schmidt_number='0.5 0.6'"
  expect_err = "The number of Schmidt numbers defined is not equal to the number of passive scalar fields!"
  ad_indexing_type = 'global'
[]

[velocity-variable-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user supplies a velocity variable which does not exist'
  cli_args = "Modules/NavierStokesFV/velocity_variable='vel_x vely' Variables/inactive='pressure T_fluid scalar'"
  expect_err = "Variable \(vely\) supplied to the NavierStokesFV action does not exist!"
  ad_indexing_type = 'global'
[]

[velocity-component-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user supplies an inappropriate number of externally created velocity components'
  cli_args = "Modules/NavierStokesFV/velocity_variable='vel_x vel_y vel_z vel_q' Variables/inactive='pressure T_fluid scalar'"
  expect_err = "The number of velocity variable names supplied to the NSFVAction is not 2 \(mesh dimension\) or 3!"
  ad_indexing_type = 'global'
[]

[pressure-variable-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user supplies a pressure variable which does not exist'
  cli_args = "Modules/NavierStokesFV/pressure_variable='pessure' Variables/inactive='vel_x vel_y T_fluid scalar'"
  expect_err = "Variable \(pessure\) supplied to the NavierStokesFV action does not exist!"
  ad_indexing_type = 'global'
[]

[fluid-temperature-variable-error]
  type = 'RunException'
  input = 2d-rc-error-action.i
  requirement = 'The system shall throw an error if the user supplies a fluid temperature variable which does not exist'
  cli_args = "Modules/NavierStokesFV/pressure_variable='T_fuid' Variables/inactive='vel_x vel_y pressure scalar'"
  expect_err = "Variable \(T_fuid\) supplied to the NavierStokesFV action does not exist!"
  ad_indexing_type = 'global'
[]

[sou]
  type = Exodiff
  input = 'test.i'
  exodiff = sou.e
  cli_args = 'GlobalParams/advected_interp_method=sou Outputs/file_base=sou'
  ad_indexing_type = 'global'
  detail = 'second-order upwind'
[]

[vanLeer]
  type = Exodiff
  input = 'test.i'
  exodiff = vanLeer.e
  cli_args = 'GlobalParams/advected_interp_method=vanLeer Outputs/file_base=vanLeer'
  ad_indexing_type = 'global'
  detail = 'van Leer'
[]

[min_mod]
  type = Exodiff
  input = 'test.i'
  exodiff = min_mod.e
  cli_args = 'GlobalParams/advected_interp_method=min_mod Outputs/file_base=min_mod'
  ad_indexing_type = 'global'
  detail = 'min-mod'
[]

[quick]
  type = Exodiff
  input = 'test.i'
  exodiff = quick.e
  cli_args = 'GlobalParams/advected_interp_method=quick Outputs/file_base=quick'
  ad_indexing_type = 'global'
  detail = 'QUICK'
[]

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

[1e6-action]
  type = Exodiff
  input = boussinesq-action.i
  exodiff = 1e6-action.e
  cli_args = 'rayleigh=1e6 Outputs/file_base=1e6-action'
  requirement = 'The system shall be able to reproduce benchmark results for a Rayleigh number of 1e6 using the INSFV actions.'
  valgrind = 'none'
  ad_indexing_type = 'global'
  issues = '#19742'
[]

[rho-error]
  type = RunException
  input = boussinesq-action.i
  cli_args = "Modules/NavierStokesFV/density='x'"
  expect_err = 'rho: The density in the boussinesq term is not constant!'
  requirement = 'The system should throw an error if the density is not a constant functor in case of Boussinesq treatment.'
  valgrind = 'none'
  ad_indexing_type = 'global'
  issues = '#19742'
[]

[wcnsfv]
  type = Exodiff
  input = wcnsfv.i
  exodiff = wcnsfv_out.e
  requirement = 'The system shall be able to model natural convection using a weakly compressible implementation.'
  ad_indexing_type = 'global'
  method = '!dbg'
  valgrind = 'none'
[]

[transient_wcnsfv]
  type = Exodiff
  input = transient-wcnsfv.i
  exodiff = transient-wcnsfv_out.e
  requirement = 'The system shall be able to model transient natural convection with a low Rayleigh number using a weakly compressible implementation.'
  ad_indexing_type = 'global'
  method = '!dbg'
  cli_args = "Outputs/out/execute_on='final'"
  recover = false # See https://github.com/idaholab/moose/issues/17906
  valgrind = 'none'
[]

[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 = 'CSVDiff'
  input = 2d-average-no-slip.i
  cli_args = "velocity_interp_method='rc' Outputs/file_base=2d-rc-no-slip"
  csvdiff = 2d-rc-no-slip.csv
  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 = 'CSVDiff'
  input = 2d-rc-slip.i
  csvdiff = 2d-rc-slip_out.csv
  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 = 'CSVDiff'
  input = diverging.i
  csvdiff = 'diverging-no-slip.csv'
  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-no-slip-tris]
  type = 'Exodiff'
  input = no-slip-tris.i
  exodiff = 'no-slip-tris_out.e'
  cli_args = "FVBCs/active='inlet-u inlet-v no-slip-wall-u no-slip-wall-v outlet-p axis-u axis-v' -pc_type lu -pc_factor_shift_type NONZERO"
  requirement = 'The system shall be able to solve a straight channel problem in cylindrical coordinates using triangular elements with no slip boundary conditions.'
  ad_indexing_type = 'global'
[]

[rz-no-slip-tris-action]
  type = 'Exodiff'
  input = no-slip-tris-action.i
  exodiff = 'no-slip-tris-action_out.e'
  cli_args = "-pc_type lu -pc_factor_shift_type NONZERO"
  requirement = 'The system shall be able to solve a straight channel problem in cylindrical coordinates using triangular elements with no slip boundary conditions and the NSFV action syntax.'
  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'
  valgrind = HEAVY
[]

[free-slip-action]
  type = 'Exodiff'
  input = 2d-rc-action.i
  exodiff = 2d-rc-action_out.e
  method = "!dbg"
  requirement = 'The system shall be able to model free-slip conditions in a channel using the NSFV action.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
[]

[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'
  valgrind = HEAVY
[]

[no-slip-action]
  type = 'Exodiff'
  input = 2d-rc-no-slip-action.i
  exodiff = 2d-rc-no-slip-action_out.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a channel using incompressible Navier Stokes action.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  issues = '#19472'
  cli_args = '-pc_type lu -pc_factor_shift_type NONZERO'
  valgrind = HEAVY
[]

[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'
  valgrind = HEAVY
[]

[scalar-transport-action]
  type = 'Exodiff'
  input = 2d-scalar-transport-action.i
  exodiff = 2d-scalar-transport-action_out.e
  method = "!dbg"
  requirement = 'The system shall be able to create and transport scalar field variables in a fluid flow field using the NSFVAction syntax.'
  ad_indexing_type = 'global'
  design = 'NSFVAction.md'
  issues = '#19472'
  valgrind = HEAVY
[]

[external-scalar-transport-action]
  type = 'Exodiff'
  input = 2d-scalar-transport-action.i
  exodiff = 2d-scalar-transport-action_out.e
  method = "!dbg"
  requirement = 'The system shall be able to transport an externally-generated scalar field variables in a fluid flow field using the NSFVAction syntax.'
  ad_indexing_type = 'global'
  cli_args = "Variables/inactive='' Modules/NavierStokesFV/velocity_variable='vel_x vel_y' Modules/NavierStokesFV/pressure_variable='pressure' Modules/NavierStokesFV/fluid_temperature_variable='T_fluid'"
  design = 'NSFVAction.md'
  issues = '#19472'
  valgrind = HEAVY
[]

[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
  valgrind = HEAVY
[]

[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
  valgrind = HEAVY
  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
[]

[mixing-length-action]
  type = 'Exodiff'
  input = 2d-mixing-length-action.i
  exodiff = 2d-mixing-length-action_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 and the NSFVAction syntax'
  ad_indexing_type = 'global'
  issues = '#19472'
  mesh_mode = REPLICATED
[]

[mixing-length-jac]
  type = 'PetscJacobianTester'
  input = 2d-mixing-length.i
  run_sim = True
  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 and show a perfect Jacobian'
  ad_indexing_type = 'global'
  design = 'navier_stokes/rans_theory.md'
  issues = '#16794 #16937'
  mesh_mode = REPLICATED
  ratio_tol = 3e-8
  difference_tol = 1e-6
  cli_args = "Mesh/gen/nx=2 Mesh/gen/ny=2"
  only_final_jacobian = True
[]

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

[linear-friction-action]
  type = 'Exodiff'
  input = 2d-rc-friction-action.i
  exodiff = 2d-rc-friction-action_out.e
  requirement = 'The system shall be able to model linear volumetric friction in a channel using NSFVAction syntax.'
  issues = '#19472'
  ad_indexing_type = 'global'
[]

[quadratic-friction-action]
  type = 'Exodiff'
  input = 2d-rc-friction-action.i
  exodiff = 2d-rc-friction-action-quad_out.e
  requirement = 'The system shall be able to model quadratic volumetric friction in a channel using NSFVAction syntax.'
  issues = '#19472'
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/friction_types='forchheimer' Outputs/file_base='2d-rc-friction-action-quad_out'"
[]

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

[ambient-convection-action]
  type = 'CSVDiff'
  input = 2d-rc-ambient-convection-action.i
  csvdiff = 2d-rc-ambient-convection-action_out.csv
  requirement = 'The system shall be able to model ambient volumetric convection in a channel using the Navier-Stokes action.'
  issues = '#19472'
  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'
  valgrind = HEAVY
  detail = 'two-dimensional triangular elements'
  abs_zero = 1e-9
[]

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

[singular_two_term_expansion]
  type = 'Exodiff'
  input = one-elem-wide-channel.i
  ad_indexing_type = 'global'
  exodiff = one-elem-wide-channel_out.e
  issues = '#16822'
  requirement = 'The system shall be able to compute gradients, when extrapolated boundary face values make the gradient computation singular, by catching the singularity error and re-running without doing boundary face value extrapolation.'
  # currently throws a warning which is non-deterministic with threads, see #20326
  max_threads = 1
[]

[discontinuous]
  type = 'Exodiff'
  ad_indexing_type = 'global'
  input = discontinuous-body-forces.i
  exodiff = discontinuous-body-forces_out.e
  design = INSFVRhieChowInterpolator.md
  requirement = 'The system shall show a monotone pressure profile in the presence of discontinuous body forces, in this case a transient from free flow to linear drag.'
[]

[transient]
  type = 'Exodiff'
  input = 2d-rc-transient.i
  exodiff = 2d-rc-transient_out.e
  method = "!dbg"
  requirement = 'The system shall be able to simulate a transient flow in a channel with enthalpy as an advected quantity.'
  ad_indexing_type = 'global'
  abs_zero = 6e-9
  valgrind = HEAVY
[]

[transient-action]
  type = 'Exodiff'
  input = 2d-rc-transient-action.i
  exodiff = 2d-rc-transient-action_out.e
  method = "!dbg"
  requirement = 'The system shall be able to simulate a transient flow in a channel with enthalpy as an advected quantity using the NSFV action syntax.'
  ad_indexing_type = 'global'
  abs_zero = 6e-9
  valgrind = HEAVY
[]

[bad_ro]
  type = RunException
  input = 'bad-ro.i'
  expect_err = "is not a INSFVMomentumResidualObject"
  requirement = 'The system shall error if a user specifies a residual object for the momentum equation that is not a incompressible Navier-Stokes finite volume momentum residual object.'
  ad_indexing_type = 'global'
[]

[bad_restriction]
  type = RunException
  input = bad-restriction.i
  expect_err = "Block restriction of interpolator user object.* doesn't match the block restriction of variable"
  requirement = 'The system shall error if the interpolation object has block restriction different from the nonlinear flow variables.'
  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'
[]

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

[point-pressure]
  type = 'Exodiff'
  input = 'lid-driven.i'
  exodiff = 'lid-driven-point-pressure.e'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity by fixing the point value of the pressure at a certain coordinate.'
  ad_indexing_type = 'global'
  cli_args = "FVKernels/mean_zero_pressure/type='FVPointValueConstraint' FVKernels/mean_zero_pressure/point='0.001 0.099 0.0' Outputs/file_base='lid-driven-point-pressure'"
[]

[point-pressure-action]
  type = 'Exodiff'
  input = 'lid-driven-action.i'
  exodiff = 'lid-driven-point-pressure-action.e'
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations in a lid-driven cavity by fixing the point value of the pressure at a certain coordinate using the NSFV action syntax.'
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/pinned_pressure_type=point-value Modules/NavierStokesFV/pinned_pressure_point='0.001 0.099 0.0' Outputs/file_base='lid-driven-point-pressure-action'"
[]

[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 = 'transient-lid-driven-with-energy.i'
  cli_args = 'Mesh/gen/nx=2 Mesh/gen/ny=2 Executioner/num_steps=1'
  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'
  valgrind = HEAVY
  cli_args = '-pc_type lu -pc_factor_shift_type NONZERO'
[]

[with-temp-action]
  type = 'Exodiff'
  input = 'lid-driven-with-energy-action.i'
  exodiff = 'lid-driven-with-energy-action_out.e'
  requirement = 'The system shall be able to get the same result as the enthalpy transport example using the NSFVAction to set up the run.'
  method = '!dbg'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  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'
  valgrind = 'none'
  cli_args = '-pc_type lu -pc_factor_shift_type NONZERO'
[]

[quiescent]
  type = 'Exodiff'
  input = 'rz-gravity-quiescent-fluid.i'
  exodiff = 'rz-gravity-quiescent-fluid_out.e'
  ad_indexing_type = 'global'
  requirement = 'The system shall yield a quiescent fluid in an axisymmetric coordinate system with a gravitational force applied and Rhie-Chow interpolation used for the velocity field.'
[]

[quiescent_jac]
  type = PetscJacobianTester
  run_sim = True
  input = rz-gravity-quiescent-fluid.i
  prereq = quiescent
  cli_args = 'Mesh/gen/nx=3 Mesh/gen/ny=3 Executioner/automatic_scaling=true Executioner/off_diagonals_in_auto_scaling=true'
  ad_indexing_type = 'global'
  requirement = 'The system shall compute an accurate Jacobian when a scaling factor is applied to a scalar variable.'
  difference_tol = 1e-7
[]

[transient]
  type = 'Exodiff'
  input = mixing_length_eddy_viscosity.i
  exodiff = mixing_length_eddy_viscosity_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  valgrind = 'none'
  recover = false
  abs_zero = 3e-10
  heavy = true
  detail = 'transient simulations'
[]

[steady]
  type = 'Exodiff'
  input = steady.i
  exodiff = steady_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  valgrind = 'none'
  recover = false
  abs_zero = 6e-10
  detail = 'steady simulations'
[]

[transient]
  type = 'Exodiff'
  input = mixing_length_total_viscosity.i
  exodiff = mixing_length_total_viscosity_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  valgrind = 'none'
  recover = false
  heavy = true
  detail = 'transient simulations'
  abs_zero = 5e-10
[]

[steady]
  type = 'Exodiff'
  input = steady.i
  exodiff = steady_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  valgrind = 'none'
  recover = false
  detail = 'steady simulations'
  abs_zero = 5e-10
[]

[steady_action]
  type = 'Exodiff'
  input = steady-action.i
  exodiff = steady-action_out.e
  method = "!dbg"
  detail = 'steady simulations with action syntax'
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  valgrind = 'none'
  recover = false
  abs_zero = 5e-10
[]

[two_term]
  type = PythonUnitTest
  input = 'test.py'
  test_case = Test2DRCWalls
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  requirement = 'The system shall exhibit second order convergence for all variables in a Cartesian, no-slip, channel-flow problem using a Rhie-Chow interpolation, including body forces, and two term boundary expansion for gradient and extrapolated boundary face value computation.'
[]

[two_term_symmetry]
  type = PythonUnitTest
  input = 'test.py'
  test_case = Test2DRCSymmetry
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  requirement = 'The system shall exhibit second order convergence for all variables in a Cartesian, no-slip, channel-flow problem with a symmetry axis using a Rhie-Chow interpolation, including body forces, and two term boundary expansion for gradient and extrapolated boundary face value computation.'
[]

[two_term_symmetry_rz]
  type = PythonUnitTest
  input = 'test.py'
  test_case = Test2DRZRCSymmetry
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  requirement = 'The system shall exhibit second order convergence for all variables in an axisymmetric, no-slip, channel-flow problem using a Rhie-Chow interpolation, including body forces, and two term boundary expansion for gradient and extrapolated boundary face value computation.'
[]

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

[vortex-skewness-corrected]
  type = PythonUnitTest
  input = test.py
  test_case = TestVortexSkewCorrected
  method = '!dbg'
  min_parallel = 2
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes, using a Rhie-Chow scheme and skewness-correction, Dirichlet boundary conditions for the velocity, and demonstrate second order convergence for the velocity and first order convergence for pressure.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  heavy = true
[]

[vortex-skewness-corrected-action]
  type = PythonUnitTest
  input = test.py
  test_case = TestVortexSkewCorrectedAction
  method = '!dbg'
  min_parallel = 2
  requirement = 'The system shall be able to solve the incompressible Navier-Stokes equations on triangular meshes using skewness-correction and the NSFV action syntax.'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
  heavy = true
[]

[run]
  type = RunApp
  input = skewed-vortex.i
  ad_indexing_type = 'global'
  requirement = 'The system shall be able to automatically expand its ghosting pattern when skew-corrected face gradients are involved in the simulation.'
[]

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

[multiapp]
  type = Exodiff
  input = 'fluid-flow.i'
  exodiff = 'fluid-flow_out.e fluid-flow_out_scalar0.e'
  requirement = 'The system shall be able to transport passive scalar quantities in an application different from the fluid flow.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  issues = '#16585'
[]

[average_yet_a_provided]
  type = RunException
  input = 'scalar-transport.i'
  cli_args = "UserObjects/rc/velocity_interp_method='average'"
  expect_err = 'Rhie Chow coefficients may not be specified for average velocity interpolation'
  requirement = 'The system shall error out if both average and Rhie Chow velocity interpolation parameters are provided.'
  issues = '#20294'
  ad_indexing_type = 'global'
[]

[RC_yet_a_not_provided]
  type = RunException
  input = 'scalar-transport.i'
  cli_args = "UserObjects/active='rc_bad' GlobalParams/rhie_chow_user_object='rc_bad'"
  expect_err = 'No INSFVKernels detected for the velocity variables*'
  requirement = 'The system shall error out if Rhie Chow velocity interpolation is requested but the a coefficients are not provided and cannot be computed from the momemtum equation.'
  issues = '#20294'
  ad_indexing_type = 'global'
[]

[vu]
  type = RunException
  input = 'scalar-transport.i'
  cli_args = "UserObjects/active='rc_bad' GlobalParams/rhie_chow_user_object='rc_bad' UserObjects/rc_bad/a_v=ay"
  expect_err = 'If the a_v coefficients are provided, then a_u must be provided'
  detail = 'for example with the Y coefficient but not the X one'
  ad_indexing_type = 'global'
[]

[bad_a_components_wu]
  type = RunException
  input = 'scalar-transport.i'
  cli_args = "UserObjects/active='rc_bad' GlobalParams/rhie_chow_user_object='rc_bad' UserObjects/rc_bad/a_w=ay"
  expect_err = 'If the a_w coefficients are provided, then a_u must be provided'
  detail = 'for example with the Z coefficient but not the X one'
  ad_indexing_type = 'global'
[]

[3d-no-face-caching]
  type = 'Exodiff'
  input = 3d-rc.i
  exodiff = 3d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to solve for flow in a 3D channel while not caching face variable values.'
  ad_indexing_type = 'global'
  cli_args = 'pressure_face_gradient_caching=false
                velocity_face_gradient_caching=false
                pressure_face_value_caching=false
                velocity_face_value_caching=false
                GlobalParams/two_term_boundary_expansion=false'
[]

[3d-no-gradient-caching]
  type = 'Exodiff'
  input = 3d-rc.i
  exodiff = 3d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to solve for flow in a 3D channel while not caching face variable gradients.'
  ad_indexing_type = 'global'
  cli_args = 'pressure_face_gradient_caching=false
                velocity_face_gradient_caching=false'
[]

[3d-pressure-gradient-caching]
  type = 'Exodiff'
  input = 3d-rc.i
  exodiff = 3d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to solve for flow in a 3D channel while only caching face pressure gradients.'
  ad_indexing_type = 'global'
  cli_args = 'velocity_face_gradient_caching=false
                pressure_face_value_caching=false
                velocity_face_value_caching=false
                GlobalParams/two_term_boundary_expansion=false'
[]

[3d-all-caching]
  type = 'Exodiff'
  input = 3d-rc.i
  exodiff = 3d-rc_out.e
  method = "!dbg"
  requirement = 'The system shall be able to solve for flow in a 3D channel while caching all variable face values and gradients.'
  ad_indexing_type = 'global'
[]

[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
  valgrind = 'none'
  recover = false
  abs_zero = 3e-10
  heavy = true
[]

[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'
  valgrind = HEAVY
  cli_args = "Outputs/file_base='rc-heated'"
[]

[rz]
  type = 'Exodiff'
  input = 2d-rc-heated.i
  exodiff = rc-heated-rz.e
  method = "!dbg"
  abs_zero = 2e-9
  detail = 'in rz geometry,'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  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'
  valgrind = HEAVY
  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'
  valgrind = HEAVY
  design = 'PINSFVEnergyAmbientConvection.md'
  cli_args = "Variables/inactive='' AuxVariables/inactive='T_solid' FVKernels/inactive='' FVBCs/inactive='' Outputs/file_base='rc-heated-solid-fluid'"
[]

[solid-fluid-action]
  type = 'Exodiff'
  input = 2d-rc-heated-action.i
  exodiff = rc-heated-solid-fluid-action.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#19472'
  design = 'NSFVAction.md'
  requirement = 'The system shall be able to solve for fluid energy diffusion, advection and convection with the solid phase in a 2D channel using the NSFV action syntax.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  cli_args = "Outputs/file_base='rc-heated-solid-fluid-action'"
[]

[solid-fluid-block]
  type = 'Exodiff'
  input = 2d-rc-heated.i
  exodiff = rc-heated-solid-fluid-block.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#16756'
  requirement = 'The system shall be able to solve for fluid energy diffusion, advection and block-restricted convection.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  design = 'PINSFVEnergyAmbientConvection.md'
  cli_args = "Variables/inactive='' AuxVariables/inactive='T_solid' FVKernels/inactive='' FVKernels/energy_convection/block='2' FVKernels/solid_energy_convection/block='2' FVKernels/energy_advection/advected_interp_method='upwind' FVBCs/inactive='' Outputs/file_base='rc-heated-solid-fluid-block'"
[]

[solid-fluid-block-action]
  type = 'Exodiff'
  input = 2d-rc-heated-action.i
  exodiff = rc-heated-solid-fluid-block-action.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#19472'
  design = 'NSFVAction.md'
  requirement = 'The system shall be able to solve for fluid energy diffusion, advection and block-restricted convection with the solid phase in a 2D channel using the NSFV action syntax.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  cli_args = "Outputs/file_base='rc-heated-solid-fluid-block-action' Modules/NavierStokesFV/ambient_convection_blocks='2' Modules/NavierStokesFV/energy_advection_interpolation='upwind' FVKernels/solid_energy_convection/block='2'"
[]

[transient]
  type = 'Exodiff'
  input = 2d-transient.i
  exodiff = 2d-transient_out.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#16756'
  requirement = 'The system shall be able to solve transient relaxations with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature.'
  ad_indexing_type = 'global'
  valgrind = 'none'
  design = 'PINSFVEnergyTimeDerivative.md INSFVMomentumTimeDerivative.md'
[]

[transient-action]
  type = 'Exodiff'
  input = 2d-transient-action.i
  exodiff = 2d-transient-action_out.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#19472'
  design = 'NSFVAction.md'
  requirement = 'The system shall be able to run transient simulations with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature using the NSFV action syntax.'
  ad_indexing_type = 'global'
  valgrind = 'none'
[]

[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'
  valgrind = HEAVY
  design = 'PINSFVMomentumBoussinesq.md'
  cli_args = "Outputs/file_base='rc-heated-boussinesq'"
[]

[boussinesq-action]
  type = 'Exodiff'
  input = 2d-rc-heated-boussinesq-action.i
  exodiff = rc-heated-boussinesq-action.e
  method = "!dbg"
  issues = '#19472'
  design = 'NSFVAction.md'
  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 using the NSFV action syntax.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  cli_args = "Outputs/file_base='rc-heated-boussinesq-action'"
[]

[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' Materials/inactive='jump smooth'"
[]

[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' Materials/inactive='jump smooth'"
[]

[reconstruct-porosity]
  type = 'Exodiff'
  input = 2d-rc-epsjump.i
  exodiff = reconstructed.e
  cli_args = "Outputs/file_base=reconstructed UserObjects/rc/smoothing_layers=20 AuxVariables/active='' ICs/inactive='porosity_1 porosity_2 porosity_continuous'"
  ad_indexing_type = 'global'
  requirement = 'The system shall exhibit a monotonic pressure profile when a number of reconstructions (essentially smoothing) are applied to the porosity.'
[]

[reconstruct-error]
  type = RunException
  input = 2d-rc-epsjump.i
  cli_args = 'UserObjects/rc/smoothing_layers=1'
  expect_err = 'If we are reconstructing porosity, then the input porosity to this user object cannot be a Moose variable'
  ad_indexing_type = 'global'
  requirement = 'The system shall error if the user attempts to create a reconstructed field from a variable.'
[]

[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-action]
  type = 'Exodiff'
  input = 2d-rc-action.i
  exodiff = rc-action.e
  method = "!dbg"
  issues = "#19472"
  requirement = 'The system should be able to run a porous medium channel computation using the NSFV action.'
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-action'"
[]

[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-driven-action]
  type = 'Exodiff'
  input = 2d-rc-action.i
  exodiff = rc-no-slip-delta-pressure-action.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 using the NSFVAction syntax.'
  abs_zero = 1e-9
  issues = "#19472"
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/inlet_boundaries='' Modules/NavierStokesFV/momentum_inlet_types='' Modules/NavierStokesFV/momentum_inlet_function='' Modules/NavierStokesFV/outlet_boundaries='left right' Modules/NavierStokesFV/momentum_outlet_types='fixed-pressure fixed-pressure' Modules/NavierStokesFV/momentum_wall_types='noslip noslip'  Modules/NavierStokesFV/pinned_pressure_value=0.01  Outputs/file_base='rc-no-slip-delta-pressure-action' Modules/NavierStokesFV/pressure_function='1 0'"
[]

[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-pressure-average-action]
  type = 'Exodiff'
  input = 2d-rc-action.i
  exodiff = rc-no-slip-mean-pressure-action.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a porous media channel using the NSFV action and a set mean pressure.'
  abs_zero = 1e-9
  issues = "#19472"
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/momentum_wall_types='noslip noslip' Modules/NavierStokesFV/momentum_outlet_types='zero-gradient'  Modules/NavierStokesFV/pin_pressure=true  Modules/NavierStokesFV/pinned_pressure_type=average  Modules/NavierStokesFV/pinned_pressure_value=0.01  Outputs/file_base='rc-no-slip-mean-pressure-action' Modules/NavierStokesFV/pressure_function=''"
[]

[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 = "GlobalParams/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-action]
  type = 'Exodiff'
  input = 2d-rc-action.i
  exodiff = rc-no-slip-symmetry-action.e
  method = "!dbg"
  requirement = 'The system shall be able to model no-slip conditions in a symmetric porous media channel using the NSFV action.'
  abs_zero = 1e-9
  issues = "#19472"
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/momentum_wall_types='noslip symmetry' Outputs/file_base='rc-no-slip-symmetry-action'"
[]

[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 with no slip boundary conditions on the channel walls.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-friction' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[friction-action]
  type = 'Exodiff'
  input = 2d-rc-friction-action.i
  exodiff = rc-friction-action.e
  method = "!dbg"
  requirement = 'The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models and the NSFV action input syntax.'
  abs_zero = 1e-9
  issues = "#19472"
  ad_indexing_type = 'global'
  cli_args = "Outputs/file_base='rc-friction-action' -pc_type lu -pc_factor_shift_type NONZERO"
[]

[friction-free-slip]
  type = 'Exodiff'
  input = 2d-rc-friction.i
  exodiff = rc-friction-free-slip.e
  method = "!dbg"
  requirement = 'The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models with free slip boundary conditions on the channel walls.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "FVBCs/inactive='no-slip-u no-slip-v' Outputs/file_base='rc-friction-free-slip' -pc_type lu -pc_factor_shift_type NONZERO"
  issues = '#16765'
[]

[friction-free-slip-action]
  type = 'Exodiff'
  input = 2d-rc-friction-action.i
  exodiff = rc-friction-free-slip-action.e
  method = "!dbg"
  requirement = 'The system shall be able to model porous flow with volumetric friction, using the Darcy and Forchheimer friction models with free slip boundary conditions using the NSFV action.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/momentum_wall_types='slip symmetry' Outputs/file_base='rc-friction-free-slip-action' -pc_type lu -pc_factor_shift_type NONZERO"
  issues = '#19472'
[]

[friction-block]
  type = 'Exodiff'
  input = 2d-rc-friction.i
  exodiff = rc-friction-block.e
  method = "!dbg"
  requirement = 'The system shall be able to model porous flow with block-restricted volumetric friction.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "FVBCs/inactive='no-slip-u no-slip-v' FVKernels/u_friction/block=2  FVKernels/v_friction/block=2 Outputs/file_base='rc-friction-block' -pc_type lu -pc_factor_shift_type NONZERO"
  issues = '#16765'
[]

[friction-block-action]
  type = 'Exodiff'
  input = 2d-rc-friction-action.i
  exodiff = rc-friction-block-action.e
  method = "!dbg"
  requirement = 'The system shall be able to model porous flow with block-restricted volumetric friction using the NSFV action.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  cli_args = "Modules/NavierStokesFV/momentum_wall_types='slip symmetry' Modules/NavierStokesFV/friction_blocks='2' Outputs/file_base='rc-friction-block-action' -pc_type lu -pc_factor_shift_type NONZERO"
  issues = '#19472'
[]

[rz-by-parts]
  type = Exodiff
  input = '2d-rc-rz-by-parts.i'
  exodiff = '2d-rc-rz-by-parts_out.e'
  method = "!dbg"
  requirement = 'The system shall be able to solve the porous flow equations in RZ geometry using an integration by parts of the pressure term.'
  abs_zero = 1e-9
  ad_indexing_type = 'global'
  issues = '#18478'
[]

[check_material]
  type = 'Exodiff'
  input = 2d-rc.i
  exodiff = 2d-rc_out.e
  requirement = 'The system shall be able to compute the speed, eg the norm of the interstitial velocity, in a porous media incompressible flow problem.'
  ad_indexing_type = 'global'
[]

[check_material-action]
  type = 'Exodiff'
  input = 2d-rc-action.i
  exodiff = 2d-rc-action_out.e
  requirement = 'The system shall be able to compute the speed in a porous media incompressible flow problem from within the NSFVAction.'
  ad_indexing_type = 'global'
[]

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

[pressure-corrected]
  type = PythonUnitTest
  input = test.py
  test_case = TestPressureCorrected
  requirement = 'The system shall show second order convergence for all variables when the porosity is interpolated and reconstructed multiple times and when a correction is applied to the pressure interpolation.'
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[pressure-corrected-action]
  type = PythonUnitTest
  input = test.py
  test_case = TestPressureCorrectedAction
  requirement = 'The system shall show second order convergence for all variables when using porosity reconstructions and friction corrections together with NSFV action syntax.'
  method = '!dbg'
  min_parallel = 8
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[pressure-corrected-jac]
  type = 'PetscJacobianTester'
  input = pressure-interpolation-corrected.i
  run_sim = True
  requirement = 'The system shall show a perfect Jacobian when the porosity is interpolated and reconstructed multiple times and when a correction is applied to the pressure interpolation.'
  ad_indexing_type = 'global'
  difference_tol = 1e-6
  cli_args = "Mesh/gen/nx=2 Mesh/gen/ny=2"
  ratio_tol = 2e-8
[]

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

[rc-no-diffusion]
  type = PythonUnitTest
  input = test.py
  test_case = TestRCNoDiffusion
  requirement = 'The system shall be able to solve the incompressible porous flow Euler equations in a 1D channel using a Rhie-Chow interpolation scheme for velocity and upwind interpolation for advected quantities and produce first order convergence for pressure and 1.5 order for velocity.'
  method = '!dbg'
  ad_indexing_type = 'global'
  required_python_packages = 'pandas matplotlib'
[]

[rc-no-diffusion-strong-bc]
  type = PythonUnitTest
  input = test.py
  test_case = TestRCNoDiffusionStrongBC
  requirement = 'The system shall be able to solve the incompressible porous flow Euler equations in a 1D channel using a Rhie-Chow interpolation scheme for velocity and upwind interpolation for advected quantities and produce first order convergence for pressure and 1.5 order for velocity with using flux boundary conditions for both variables.'
  method = '!dbg'
  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'
[]

[transient]
  type = 'Exodiff'
  input = 2d-transient.i
  exodiff = 2d-transient_out.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#16756 #18806'
  requirement = 'The system shall be able to solve transient relaxations within the weakly compressible approximation, with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  design = 'PINSFVEnergyTimeDerivative.md'
  recover = false # see #19126
[]

[transient-action]
  type = 'Exodiff'
  input = 2d-transient-action.i
  exodiff = 2d-transient-action_out.e
  method = "!dbg"
  abs_zero = 1e-9
  issues = '#19472'
  requirement = 'The system shall be able to solve weakly compressible transient problems with the NSFV action syntax.'
  ad_indexing_type = 'global'
  valgrind = HEAVY
  recover = false # see #19126
[]

[transient-jac]
  type = 'PetscJacobianTester'
  input = 2d-transient.i
  run_sim = True
  difference_tol = 1.83
  requirement = 'The system shall be able to solve transient relaxations within the weakly compressible approximation, with fluid energy diffusion, advection and convection with the solid phase in a 2D channel, modeling both fluid and solid temperature and show a perfect Jacobian.'
  ad_indexing_type = 'global'
  cli_args = 'Mesh/gen/nx=2 Mesh/gen/ny=2 Executioner/num_steps=1'
[]

[using_velocity]
  type = 'Exodiff'
  input = dirichlet_bcs_velocity.i
  exodiff = dirichlet_bcs_velocity_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  detail = 'using a velocity postprocessor'
[]

[using_mdot]
  type = 'Exodiff'
  input = dirichlet_bcs_mdot.i
  exodiff = dirichlet_bcs_mdot_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  detail = 'using a mass flow rate postprocessor'
[]

[direct]
  type = 'Exodiff'
  input = flux_bcs_direct.i
  exodiff = flux_bcs_direct_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 3e-9
  recover = false # see #19126
  detail = 'specifying the flux values directly though the flow rates'
[]

[from_velocity]
  type = 'Exodiff'
  input = flux_bcs_velocity.i
  exodiff = flux_bcs_velocity_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  detail = 'using the inlet velocity to compute the fluxes'
[]

[from_velocity_action]
  type = 'Exodiff'
  input = flux_bcs_velocity-action.i
  exodiff = flux_bcs_velocity-action_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  detail = 'using the inlet velocity to compute the fluxes with the NSFV action syntax'
[]

[from_mdot]
  type = 'Exodiff'
  input = flux_bcs_mdot.i
  exodiff = flux_bcs_mdot_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  detail = 'using the inlet mass flow rate to compute the fluxes'
[]

[from_mdot_action]
  type = 'Exodiff'
  input = flux_bcs_mdot-action.i
  exodiff = flux_bcs_mdot-action_out.e
  method = "!dbg"
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  detail = 'using the inlet mass flow rate to compute the fluxes with NSFVAction syntax'
[]

[missing_info_dirichlet_velocity]
  type = RunException
  input = 'dirichlet_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_u FVBCs/inlet_u_bis/type=WCNSFVInletVelocityBC FVBCs/inlet_u_bis/variable=u FVBCs/inlet_u_bis/boundary=left FVBCs/inlet_u_bis/mdot_pp=0'
  expect_err = 'Mass flow rate, area and density should be provided if velocity is not'
  detail = 'the inlet velocity cannot be determined from the mass flow rate due to insufficient parameters'
[]

[missing_info_dirichlet_temperature]
  type = RunException
  input = 'dirichlet_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVInletTemperatureBC FVBCs/inlet_T_bis/variable=T FVBCs/inlet_T_bis/boundary=left'
  expect_err = 'If not providing the temperature, the energy flow rate should be provided'
  detail = 'the inlet temperature cannot be determined because neither the temperature or the energy flow rate have been provided'
[]

[missing_info_dirichlet_temperature_from_energy_flow]
  type = RunException
  input = 'dirichlet_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVInletTemperatureBC FVBCs/inlet_T_bis/variable=T FVBCs/inlet_T_bis/boundary=left FVBCs/inlet_T_bis/energy_pp=0'
  expect_err = 'If not providing the inlet temperature, the inlet velocity or mass flow should'
  detail = 'the inlet temperature cannot be determined from the energy flow rate due to missing fluid velocity information'
[]

[missing_info_dirichlet_temperature_from_velocity]
  type = RunException
  input = 'dirichlet_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVInletTemperatureBC FVBCs/inlet_T_bis/variable=T FVBCs/inlet_T_bis/boundary=left FVBCs/inlet_T_bis/energy_pp=0 FVBCs/inlet_T_bis/velocity_pp=0'
  expect_err = 'If providing the inlet velocity, the density, the area and the specific heat'
  detail = 'the inlet temperature cannot be determined from the velocity due to insufficient parameters'
[]

[missing_info_dirichlet_temperature_from_mdot]
  type = RunException
  input = 'dirichlet_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVInletTemperatureBC FVBCs/inlet_T_bis/variable=T FVBCs/inlet_T_bis/boundary=left FVBCs/inlet_T_bis/energy_pp=0 FVBCs/inlet_T_bis/mdot_pp=0'
  expect_err = 'If providing the inlet mass flow rate, the inlet specific heat capacity should'
  detail = 'the inlet temperature cannot be determined from the mass flow rate due to insufficient parameters'
[]

[both_velocity_and_mdot]
  type = RunException
  input = 'dirichlet_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inlet_u/mdot_pp=0'
  expect_err = 'If setting the velocity directly, no need for inlet mass flow rate or area'
  detail = 'redundant information is provided for setting the inlet velocity'
[]

[both_temperature_and_mdot]
  type = RunException
  input = 'dirichlet_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inlet_T/mdot_pp=0'
  expect_err = 'If setting the temperature directly, no need for inlet velocity, mass flow or energy'
  detail = 'redundant information is provided for setting the inlet temperature'
[]

[mass_mdot_needs_area]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_mass FVBCs/inlet_mass_bis/type=WCNSFVMassFluxBC FVBCs/inlet_mass_bis/variable=pressure FVBCs/inlet_mass_bis/boundary=left FVBCs/inlet_mass_bis/mdot_pp=0'
  expect_err = 'The inlet area should be provided along with the mass flow rate'
  detail = 'only the mass flow rate is provided to compute the inlet mass flux'
[]

[mass_velocity_needs_density]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_mass FVBCs/inlet_mass_bis/type=WCNSFVMassFluxBC FVBCs/inlet_mass_bis/variable=pressure FVBCs/inlet_mass_bis/boundary=left FVBCs/inlet_mass_bis/velocity_pp=0'
  expect_err = 'Velocity and density should be provided if the mass flow rate is not'
  detail = 'only the velocity is provided to compute the inlet mass flux'
[]

[momentum_mdot_needs_area]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_u FVBCs/inlet_u_bis/type=WCNSFVMassFluxBC FVBCs/inlet_u_bis/variable=pressure FVBCs/inlet_u_bis/boundary=left FVBCs/inlet_u_bis/mdot_pp=0'
  expect_err = 'The inlet area should be provided along with the mass flow rate'
  detail = 'only the mass flow rate is provided to compute the inlet momentum flux'
[]

[momentum_velocity_needs_density]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_u FVBCs/inlet_u_bis/type=WCNSFVMassFluxBC FVBCs/inlet_u_bis/variable=pressure FVBCs/inlet_u_bis/boundary=left FVBCs/inlet_u_bis/velocity_pp=0'
  expect_err = 'Velocity and density should be provided if the mass flow rate is not'
  detail = 'only the velocity is provided to compute the inlet momentum flux'
[]

[energy_need_temperature]
  type = RunException
  input = 'flux_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVEnergyFluxBC FVBCs/inlet_T_bis/variable=T_fluid FVBCs/inlet_T_bis/boundary=left'
  expect_err = 'If not providing the energy flow rate, the inlet temperature should be provided'
  detail = 'the inlet energy cannot be determined because neither the temperature or the energy flow rate have been provided'
[]

[energy_temperature_needs_velocity_or_mdot]
  type = RunException
  input = 'flux_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVEnergyFluxBC FVBCs/inlet_T_bis/variable=T_fluid FVBCs/inlet_T_bis/boundary=left FVBCs/inlet_T_bis/temperature_pp=0'
  expect_err = 'If not providing the inlet energy flow rate, the inlet velocity or mass flow'
  detail = 'the inlet energy flux cannot be determined because the temperature was provided but neither the inlet mass flow rate or velocity were provided'
[]

[energy_temperature_velocity_needs_density]
  type = RunException
  input = 'flux_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVEnergyFluxBC FVBCs/inlet_T_bis/variable=T_fluid FVBCs/inlet_T_bis/boundary=left FVBCs/inlet_T_bis/temperature_pp=0 FVBCs/inlet_T_bis/velocity_pp=0'
  expect_err = 'If providing the inlet velocity, the density, the area and the specific heat'
  detail = 'the inlet energy flux cannot be determined because the temperature and velocity were provided but the density or specific heat capacity were not'
[]

[energy_temperature_mdot_needs_cp]
  type = RunException
  input = 'flux_bcs_mdot.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_T FVBCs/inlet_T_bis/type=WCNSFVEnergyFluxBC FVBCs/inlet_T_bis/variable=T_fluid FVBCs/inlet_T_bis/boundary=left FVBCs/inlet_T_bis/temperature_pp=0 FVBCs/inlet_T_bis/mdot_pp=0'
  expect_err = 'If providing the inlet mass flow rate, the inlet specific heat capacity and flow'
  detail = 'the inlet energy flux cannot be determined because the temperature and mass flow rate were provided but the specific heat capacity or flow area were not'
[]

[scalar_need_flux_or_value]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_scalar FVBCs/inlet_scalar_bis/type=WCNSFVScalarFluxBC FVBCs/inlet_scalar_bis/variable=scalar FVBCs/inlet_scalar_bis/boundary=left'
  expect_err = 'If not providing the scalar flow rate, the inlet scalar concentration should be provided'
  detail = 'neither the scalar quantity flux or boundary values are specified when attempting to compute the scalar quantity flux'
[]

[scalar_mdot_needs_area]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_scalar FVBCs/inlet_scalar_bis/type=WCNSFVScalarFluxBC FVBCs/inlet_scalar_bis/variable=scalar FVBCs/inlet_scalar_bis/boundary=left FVBCs/inlet_scalar_bis/scalar_value_pp=0 FVBCs/inlet_scalar_bis/mdot_pp=0'
  expect_err = 'If providing the inlet mass flow rate, the inlet density and flow area should be provided as well'
  detail = 'only the mass flow rate is provided to compute the inlet scalar quantity flux'
[]

[scalar_need_some_velocity_info]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inactive=inlet_scalar FVBCs/inlet_scalar_bis/type=WCNSFVScalarFluxBC FVBCs/inlet_scalar_bis/variable=scalar FVBCs/inlet_scalar_bis/boundary=left FVBCs/inlet_scalar_bis/scalar_value_pp=0'
  expect_err = 'If not providing the scalar flow rate, the inlet velocity or mass flow should be provided'
  detail = 'only the scalar quantity boundary value is provided to compute the inlet scalar quantity flux'
[]

[mass_both_velocity_and_mdot]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inlet_mass/mdot_pp=0'
  expect_err = 'If setting the mass flow rate directly, no need for inlet velocity'
  detail = 'redundant information is provided for setting the inlet mass flux'
[]

[momentum_both_velocity_and_mdot]
  type = RunException
  input = 'flux_bcs_velocity.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inlet_u/mdot_pp=0'
  expect_err = 'If setting the mass flow rate directly, no need for inlet velocity'
  detail = 'redundant information is provided for setting the inlet momentum flux'
[]

[both_temperature_and_mdot]
  type = RunException
  input = 'flux_bcs_direct.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inlet_T/mdot_pp=0'
  expect_err = 'If setting the energy flow rate directly, no need for inlet velocity, mass flow or energy'
  detail = 'redundant information is provided for setting the inlet energy flux'
[]

[both_scalar_and_mdot]
  type = RunException
  input = 'flux_bcs_direct.i'
  ad_indexing_type = 'global'
  cli_args = 'FVBCs/inlet_scalar/mdot_pp=0'
  expect_err = 'If setting the scalar flux directly, no need for inlet velocity, mass flow or scalar '
  detail = 'redundant information is provided for setting the inlet scalar flux'
[]

[transient]
  type = 'Exodiff'
  input = 2d-transient.i
  exodiff = 2d-transient_out.e
  method = "!dbg"
  cli_args = "AuxVariables/inactive='mixing_length'"
  requirement = 'The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation.'
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  valgrind = HEAVY
[]

[transient-action]
  type = 'Exodiff'
  input = 2d-transient-action.i
  exodiff = 2d-transient-action_out.e
  method = "!dbg"
  requirement = 'The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation using the NSFV action syntax.'
  ad_indexing_type = 'global'
  abs_zero = 1e-9
  recover = false # see #19126
  valgrind = HEAVY
[]

[turbulence]
  type = 'Exodiff'
  input = 2d-transient.i
  exodiff = 2d-turbulent.e
  cli_args = "FVKernels/inactive='' AuxKernels/inactive='' Outputs/file_base='2d-turbulent'"
  method = "!dbg"
  requirement = 'The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and a mixing length turbulence model.'
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  abs_zero = 1e-9
  recover = false # see #19126
  valgrind = HEAVY
  min_ad_size = 53 # You can draw out the stencil. For Rhie-Chow we interpolate element 'a' coefficients
  # to a face. Each element 'a' coefficient will be dependent on face computations of rho
  # and the eddy viscosity. For rho that means we need one layer of neighbor dofs from
  # pressure and temperature. For the eddy viscosity we require face gradient of velocity.
  # That means we need one layer of neighbor cell gradients, which in turn means we need
  # *two* layers of neighbor dofs from velocity. When you calculate it all out, you see
  # that we require 52 entries in our sparsity pattern for a face 'a' coefficient (which
  # is used in computation of the Rhie-Chow velocity required for the advection kernels)
  # Unfortunately I need to further bump this by one because I have a bug in my metaphysicl
  # debugging assertions (which I'm fixing)
[]

[turbulence-action]
  type = 'Exodiff'
  input = 2d-transient-action.i
  exodiff = 2d-turbulent-action.e
  cli_args = "Modules/NavierStokesFV/turbulence_handling='mixing-length' Modules/NavierStokesFV/mixing_length_delta='0.5' Modules/NavierStokesFV/mixing_length_walls='top' Modules/NavierStokesFV/mixing_length_aux_execute_on='initial' Modules/NavierStokesFV/mixing_length_two_term_bc_expansion=true Outputs/file_base='2d-turbulent-action'"
  method = "!dbg"
  requirement = 'The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and a mixing length turbulence model using the NSFV action syntax.'
  ad_indexing_type = 'global'
  mesh_mode = REPLICATED
  abs_zero = 1e-9
  recover = false # see #19126
  valgrind = HEAVY
  min_ad_size = 53 # You can draw out the stencil. For Rhie-Chow we interpolate element 'a' coefficients
  # to a face. Each element 'a' coefficient will be dependent on face computations of rho
  # and the eddy viscosity. For rho that means we need one layer of neighbor dofs from
  # pressure and temperature. For the eddy viscosity we require face gradient of velocity.
  # That means we need one layer of neighbor cell gradients, which in turn means we need
  # *two* layers of neighbor dofs from velocity. When you calculate it all out, you see
  # that we require 52 entries in our sparsity pattern for a face 'a' coefficient (which
  # is used in computation of the Rhie-Chow velocity required for the advection kernels)
  # Unfortunately I need to further bump this by one because I have a bug in my metaphysicl
  # debugging assertions (which I'm fixing)
[]

[transient-jac]
  type = 'PetscJacobianTester'
  input = 2d-transient.i
  run_sim = True
  requirement = 'The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and show a perfect Jacobian.'
  ad_indexing_type = 'global'
  difference_tol = 1e-7
  cli_args = 'Mesh/gen/nx=2 Mesh/gen/ny=2 Executioner/num_steps=1'
[]

[turbulence-jac]
  type = 'PetscJacobianTester'
  input = 2d-transient.i
  run_sim = True
  ad_indexing_type = 'global'
  difference_tol = 1e-7
  cli_args = "Mesh/gen/nx=2 Mesh/gen/ny=2 Executioner/num_steps=1 FVKernels/inactive='' AuxKernels/inactive='' Outputs/file_base='2d-turbulent'"
  requirement = 'The system shall be able to solve for a transient 2D channel case with a weakly compressible formulation and a mixing length turbulence model and show a perfect Jacobian.'
  mesh_mode = REPLICATED
  min_ad_size = 53 # You can draw out the stencil. For Rhie-Chow we interpolate element 'a' coefficients
  # to a face. Each element 'a' coefficient will be dependent on face computations of rho
  # and the eddy viscosity. For rho that means we need one layer of neighbor dofs from
  # pressure and temperature. For the eddy viscosity we require face gradient of velocity.
  # That means we need one layer of neighbor cell gradients, which in turn means we need
  # *two* layers of neighbor dofs from velocity. When you calculate it all out, you see
  # that we require 52 entries in our sparsity pattern for a face 'a' coefficient (which
  # is used in computation of the Rhie-Chow velocity required for the advection kernels)
  # Unfortunately I need to further bump this by one because I have a bug in my metaphysicl
  # debugging assertions (which I'm fixing)
[]

[transient]
  type = 'Exodiff'
  input = 2d-transient.i
  exodiff = 2d-transient_out.e
  method = "!dbg"
  requirement = 'The system shall be able to use realistic fluid properties in a weakly compressible flow simulation'
  ad_indexing_type = 'global'
  rel_err = 7e-5
  # The non-linear tolerance is actually fairly tight, but the variable values
  # are small at the no-slip boundary conditions, and this leads to absolute diffs of around 1e-15,
  # triggering the default relative tolerance error criterion
  valgrind = HEAVY
[]

[fluidprops]
  type = 'Exodiff'
  input = functorfluidprops.i
  exodiff = functorfluidprops_out.e
  requirement = 'The system shall be able to output gradeurs, derivatives and adimensional quantities from realistic functor fluid properties'
  ad_indexing_type = 'global'
  recover = false
[]

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

[ns_function_ics]
  type = 'Exodiff'
  input = 'test_function.i'
  exodiff = 'test_function_out.e'
  requirement = 'The system shall be able to set initial conditions for fluid flow variables using functions.'
  design = 'NSFunctionInitialCondition.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'
[]

[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_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' Postprocessors/inactive='inlet_momentum_x inlet_momentum_y'"
  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_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' Postprocessors/inactive='inlet_momentum_x inlet_momentum_y'"
  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'
[]

[rayleigh]
  type = 'CSVDiff'
  input = 'natural_convection.i'
  csvdiff = 'natural_convection_out.csv'
  requirement = 'The system shall be able to compute the Rayleigh number in a natural convection flow simulation'
  ad_indexing_type = 'global'
[]