Fluid Properties Requirements Traceability Matrix

[./specific_enthalpy_aux]
  type = 'Exodiff'
  input = 'specific_enthalpy_aux.i'
  exodiff = 'specific_enthalpy_aux_out.e'
  requirement = 'The system shall compute specific enthalpy from pressure and temperature'
  threading = '!pthreads'
  issues = '#19225'
  design = 'SpecificEnthalpyAux.md'
[../]

[./stagnation_pressure_aux]
  type = 'Exodiff'
  input = 'stagnation_pressure_aux.i'
  exodiff = 'stagnation_pressure_aux_out.e'
  requirement = 'The system shall compute stagnation pressure from specific volume, specific internal energy, and velocit'
  threading = '!pthreads'
  issues = '#19225'
  design = 'StagnationPressureAux.md'
[../]

[./stagnation_temperature_aux]
  type = 'Exodiff'
  input = 'stagnation_temperature_aux.i'
  exodiff = 'stagnation_temperature_aux_out.e'
  requirement = 'The system shall compute stagnation temperature from specific volume, specific internal energy, and velocity'
  threading = '!pthreads'
  issues = '#19225'
  design = 'StagnationTemperatureAux.md'
[../]

[./fluid_density_aux]
  type = 'Exodiff'
  input = 'fluid_density_aux.i'
  exodiff = 'fluid_density_aux_out.e'
  requirement = 'The system shall compute fluid density from pressure and temperature.'
  issues = '#17546'
  design = 'FluidDensityAux.md'
[../]

[./brine]
  type = CSVDiff
  input = 'brine.i'
  csvdiff = 'brine_out.csv'
  allow_test_objects = true
  threading = '!pthreads'
  requirement = 'The system shall compute properties of brine'
  design = 'BrineFluidProperties.md'
  issues = '#6972'
[../]

[./brine_tabulated]
  type = CSVDiff
  input = 'brine_tabulated.i'
  csvdiff = 'brine_out.csv'
  allow_test_objects = true
  threading = '!pthreads'
  prereq = brine
  requirement = 'The system shall compute properties of brine using tabulated water properties'
  design = 'BrineFluidProperties.md'
  issues = '#6972 #12812'
[../]

[calorically_imperfect_gas]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  requirement = 'The system shall compute properties for a calorically imperfect but otherwise ideal gas'
  design = 'CaloricallyImperfectGas.md'
  issues = '#20101'
  threading = '!pthreads'
[]

[co2]
  type = CSVDiff
  input = 'co2.i'
  csvdiff = 'co2_out.csv'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute fluid properties for carbon dioxide in a gaseous phase.'
  design = 'CO2FluidProperties.md'
  issues = '#6972'
[]

[air_pT]
  type = CSVDiff
  input = 'air.i'
  csvdiff = 'air_out.csv'
  threading = '!pthreads'
  detail = 'spline-based table-lookup air fluid properties,'
  required_submodule = 'contrib/air'
[]

[CO2_pT]
  type = CSVDiff
  input = 'co2.i'
  csvdiff = 'co2_out.csv'
  threading = '!pthreads'
  detail = 'vapor carbon dioxide fluid properties,'
  required_submodule = 'contrib/carbon_dioxide'
[]

[helium_pT]
  type = CSVDiff
  input = 'helium.i'
  csvdiff = 'helium_out.csv'
  threading = '!pthreads'
  detail = 'spline-based table-lookup Helium fluid properties,'
  required_submodule = 'contrib/helium'
[]

[nitrogen_pT]
  type = CSVDiff
  input = 'nitrogen.i'
  csvdiff = 'nitrogen_out.csv'
  threading = '!pthreads'
  detail = 'spline-based table-lookup Nitrogen fluid properties,'
  required_submodule = 'contrib/nitrogen'
[]

[potassium_pT]
  type = CSVDiff
  input = 'potassium.i'
  csvdiff = 'potassium_out.csv'
  threading = '!pthreads'
  detail = 'liquid Potassium fluid properties,'
  required_submodule = 'contrib/potassium'
[]

[sodium_pT]
  type = CSVDiff
  input = 'sodium.i'
  csvdiff = 'sodium_out.csv'
  threading = '!pthreads'
  detail = 'liquid Sodium fluid properties.'
  required_submodule = 'contrib/sodium'
[]

[./1ph.p_T]
  type = 'RunApp'
  input = '1ph.p_T.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'Single-phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                              1.0000000000e\+05  Pa
Temperature:                                        300  K
Density:                                    1.162655505  kg/m\^3
Specific volume:                                 0.8601  m\^3/kg
Specific internal energy:              2.1502500000e\+05  J/kg
Specific enthalpy:                     3.0103500000e\+05  J/kg
Specific entropy:                           2422.704793  J/\(kg-K\)

Sound speed:                                347.0072045  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:             1003.45  J/\(kg-K\)
Specific heat at constant volume:                716.75  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:        0.003333333333  1/K'

  requirement = "The fluid properties interrogator shall output static-state fluid properties for (p, T) input."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.p_T.json]
  type = 'RunApp'
  input = '1ph.p_T.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"static":{"T":300.0,"beta":0.[0-9]+,"c":347.[0-9]+,"cp":1003.45,"cv":716.[0-9]+,"e":215025.[0-9]+,"h":301035.0,"k":0.[0-9]+,"mu":1.823e-05,"p":100000.0,"rho":1.[0-9]+,"s":2422.[0-9]+,"v":0.86[0-9]+}}'

  requirement = "The fluid properties interrogator shall output static-state fluid properties for (p, T) input in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.rho_e]
  type = 'RunApp'
  input = '1ph.rho_e.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'Single-phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                              8.6010000000e\+04  Pa
Temperature:                                        300  K
Density:                                              1  kg/m\^3
Specific volume:                                      1  m\^3/kg
Specific internal energy:              2.1502500000e\+05  J/kg
Specific enthalpy:                     3.0103500000e\+05  J/kg
Specific entropy:                           2465.912381  J/\(kg-K\)

Sound speed:                                347.0072045  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:             1003.45  J/\(kg-K\)
Specific heat at constant volume:                716.75  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:        0.003333333333  1/K'

  requirement = "The fluid properties interrogator shall output static-state fluid properties for (rho, e) input."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.rho_e.json]
  type = 'RunApp'
  input = '1ph.rho_e.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"static":{"T":299.[0-9]+,"beta":0.[0-9]+,"c":347.[0-9]+,"cp":1003.45,"cv":716.[0-9]+,"e":215025.0,"h":301034.[0-9]+,"k":0.[0-9]+,"mu":1.823e-05,"p":86009.[0-9]+,"rho":1.0,"s":2465.[0-9]+,"v":1.0}}'

  requirement = "The fluid properties interrogator shall output static-state fluid properties for (rho, e) input in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.rho_p]
  type = 'RunApp'
  input = '1ph.rho_p.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'Single-phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                              1.0000000000e\+05  Pa
Temperature:                                348.7966516  K
Density:                                              1  kg/m\^3
Specific volume:                                      1  m\^3/kg
Specific internal energy:              2.5000000000e\+05  J/kg
Specific enthalpy:                     3.5000000000e\+05  J/kg
Specific entropy:                           2573.931349  J/\(kg-K\)

Sound speed:                                374.1657387  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:             1003.45  J/\(kg-K\)
Specific heat at constant volume:                716.75  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:              0.002867  1/K'

  requirement = "The fluid properties interrogator shall output static-state fluid properties for (rho, p) input."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.rho_p.json]
  type = 'RunApp'
  input = '1ph.rho_p.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"static":{"T":348.[0-9]+,"beta":0.[0-9]+,"c":374.[0-9]+,"cp":1003.45,"cv":716.[0-9]+,"e":250000.[0-9]+,"h":350000.[0-9]+,"k":0.[0-9]+,"mu":1.823e-05,"p":100000.0,"rho":1.0,"s":2573.[0-9]+,"v":1.0}}'

  requirement = "The fluid properties interrogator shall output static-state fluid properties for (rho, p) input in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.rho_rhou_rhoE]
  type = 'RunApp'
  input = '1ph.rho_rhou_rhoE.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'Single-phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                                             1  Pa
Temperature:                                        2.8  K
Density:                                            0.5  kg/m\^3
Specific volume:                                      2  m\^3/kg
Specific internal energy:                             5  J/kg
Specific enthalpy:                                    7  J/kg
Specific entropy:                           2.574048543  J/\(kg-K\)

Sound speed:                                1.673320053  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:                 2.5  J/\(kg-K\)
Specific heat at constant volume:           1.785714286  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:          0.3571428571  1/K



Single-phase STAGNATION properties:
--------------------------------------------------------------------------------
Pressure:                                    1.27312569  Pa
Temperature:                                          3  K
Density:                                   0.5941253221  kg/m\^3
Specific volume:                            1.683146573  m\^3/kg
Specific internal energy:                   5.357142857  J/kg
Specific enthalpy:                                  7.5  J/kg
Specific entropy:                           2.574048543  J/\(kg-K\)

Sound speed:                                1.732050808  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:                 2.5  J/\(kg-K\)
Specific heat at constant volume:           1.785714286  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:          0.3333333333  1/K'

  requirement = "The fluid properties interrogator shall output static-state and stagnation-state fluid properties for (rho, rhou, rhoE) input with a single-phase fluid properties object."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./1ph.rho_rhou_rhoE.json]
  type = 'RunApp'
  input = '1ph.rho_rhou_rhoE.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"stagnation":{"T0":2.[0-9]+,"beta0":0.[0-9]+,"c0":1.[0-9]+,"cp0":2.[0-9]+,"cv0":1.[0-9]+,"e0":5.[0-9]+,"h0":7.[0-9]+,"k0":0.[0-9]+,"mu0":1.823e-05,"p0":1.[0-9]+,"rho0":0.[0-9]+,"s0":2.[0-9]+,"v0":1.[0-9]+,"vel":1.0},"static":{"T":2.[0-9]+,"beta":0.[0-9]+,"c":1.[0-9]+,"cp":2.[0-9]+,"cv":1.[0-9]+,"e":5.0,"h":6.[0-9]+,"k":0.[0-9]+,"mu":1.823e-05,"p":0.[0-9]+,"rho":0.5,"s":2.[0-9]+,"v":2.0}}'

  requirement = "The fluid properties interrogator shall output static-state and stagnation-state fluid properties for (rho, rhou, rhoE) input with a single-phase fluid properties object in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph.p_T]
  type = 'RunApp'
  input = '2ph.p_T.i'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'TWO-PHASE properties:
--------------------------------------------------------------------------------
Critical pressure:                                   25  Pa
Latent heat of vaporization:           6.8896500000e\+05  J/kg



LIQUID phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                              1.0000000000e\+05  Pa
Temperature:                                        300  K
Density:                                    1.162655505  kg/m\^3
Specific volume:                                 0.8601  m\^3/kg
Specific internal energy:              2.1502500000e\+05  J/kg
Specific enthalpy:                     3.0103500000e\+05  J/kg
Specific entropy:                           2422.704793  J/\(kg-K\)

Sound speed:                                347.0072045  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:             1003.45  J/\(kg-K\)
Specific heat at constant volume:                716.75  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:        0.003333333333  1/K



VAPOR phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                              1.0000000000e\+05  Pa
Temperature:                                        300  K
Density:                                    1.111111111  kg/m\^3
Specific volume:                                    0.9  m\^3/kg
Specific internal energy:              9.0000000000e\+05  J/kg
Specific enthalpy:                     9.9000000000e\+05  J/kg
Specific entropy:                      1.5368604527e\+04  J/\(kg-K\)

Sound speed:                                314.6426545  m/s
Dynamic viscosity:                     1.7000000000e-05  Pa-s
Specific heat at constant pressure:                3300  J/\(kg-K\)
Specific heat at constant volume:                  3000  J/\(kg-K\)
Thermal conductivity:                              0.05  W/\(m-K\)
Volumetric expansion coefficient:        0.003333333333  1/K'

  requirement = "The fluid properties interrogator shall output two-phase and static-state, single-phase fluid properties for (p, T) input with a two-phase fluid properties object."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph.p_T.json]
  type = 'RunApp'
  input = '2ph.p_T.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"2-phase":{"h_lat":688964.[0-9]+,"p_critical":25.0},"liquid":{"static":{"T":300.0,"beta":0.[0-9]+,"c":347.[0-9]+,"cp":1003.45,"cv":716.[0-9]+,"e":215025.[0-9]+,"h":301035.0,"k":0.[0-9]+,"mu":1.823e-05,"p":100000.0,"rho":1.[0-9]+,"s":2422.[0-9]+,"v":0.[0-9]+}},"vapor":{"static":{"T":300.0,"beta":0.[0-9]+,"c":314.[0-9]+,"cp":3299.[0-9]+,"cv":2999.[0-9]+,"e":899999.[0-9]+,"h":989999.[0-9]+,"k":0.[0-9]+,"mu":1.7e-05,"p":100000.0,"rho":1.[0-9]+,"s":15368.[0-9]+,"v":0.[0-9]+}}}'

  requirement = "The fluid properties interrogator shall output two-phase and static-state, single-phase fluid properties for (p, T) input with a two-phase fluid properties object in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph_ncg_p_T]
  type = 'RunApp'
  input = '2ph_ncg_p_T.i'
  cli_args = 'FluidPropertiesInterrogator/x_ncg=0.1'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'TWO-PHASE properties:
--------------------------------------------------------------------------------
Critical pressure:                                   25  Pa
Latent heat of vaporization:           8.7784021256e\+05  J/kg



LIQUID phase STATIC properties:
--------------------------------------------------------------------------------
Pressure:                              1.0000000000e\+05  Pa
Temperature:                                372.7559289  K
Density:                                    0.993600213  kg/m\^3
Specific volume:                            1.006441008  m\^3/kg
Specific internal energy:              2.5161025201e\+05  J/kg
Specific enthalpy:                     3.5225435281e\+05  J/kg
Specific entropy:                           2486.783172  J/\(kg-K\)

Sound speed:                                375.3688068  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:                 945  J/\(kg-K\)
Specific heat at constant volume:                   675  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)
Volumetric expansion coefficient:        0.002682720575  1/K



Vapor mixture STATIC properties:
--------------------------------------------------------------------------------
Mass fraction 0:                                    0.1  -
Pressure:                              1.0000000000e\+05  Pa
Temperature:                                372.7559289  K
Density:                                   0.8972309616  kg/m\^3
Specific volume:                            1.114540227  m\^3/kg
Specific internal energy:              1.0334658129e\+06  J/kg
Specific enthalpy:                     1.1449198356e\+06  J/kg

Sound speed:                                351.3883482  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:              3071.5  J/\(kg-K\)
Specific heat at constant volume:                2772.5  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)'

  requirement = "The fluid properties interrogator shall output two-phase and static-state, single-phase fluid properties for (p, T) input with a two-phase NCG fluid properties object when the NCG mass fraction is provided."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[2ph_ncg_p_T_sat]
  type = RunApp
  input = '2ph_ncg_p_T.i'
  cli_args = "
      FluidPropertiesInterrogator/fp=fp_2phase_ncg_partial_pressure
      FluidPropertiesInterrogator/p=40e3
      FluidPropertiesInterrogator/T=300"
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'Mass fraction 0:\s+0.8463343673'
  requirement = "The fluid properties interrogator shall output two-phase and static-state, single-phase fluid properties for (p, T) input with a two-phase NCG fluid properties object when the NCG mass fraction is not provided."
  issues = '#28335'
  design = 'FluidPropertiesInterrogator.md'
[]

[./2ph_ncg_p_T.json]
  type = 'RunApp'
  input = '2ph_ncg_p_T.i'
  cli_args = 'FluidPropertiesInterrogator/x_ncg=0.1 FluidPropertiesInterrogator/json=true'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"2-phase":{"h_lat":877840.[0-9]+,"p_critical":25.0},"liquid":{"static":{"T":372.[0-9]+,"beta":0.[0-9]+,"c":375.[0-9]+,"cp":945.[0-9]+,"cv":675.[0-9]+,"e":251610.[0-9]+,"h":352254.[0-9]+,"k":0.[0-9]+,"mu":1.823e-05,"p":100000.0,"rho":0.[0-9]+,"s":2486.[0-9]+,"v":1.[0-9]+}},"vapor-mixture":{"static":{"T":372.[0-9]+,"c":351.[0-9]+,"cp":3071.[0-9]+,"cv":2772.[0-9]+,"e":1033465.[0-9]+,"h":1144919.[0-9]+,"k":0.[0-9]+,"mu":1.[0-9]+e-05,"p":100000.0,"rho":0.[0-9]+,"v":1.[0-9]+}}}'

  requirement = "The fluid properties interrogator shall output two-phase and static-state, single-phase fluid properties for (p, T) input with a two-phase NCG fluid properties object in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./vapor_mixture_rho_e]
  type = 'RunApp'
  input = 'vapor_mixture_rho_e.i'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'Vapor mixture STATIC properties:
--------------------------------------------------------------------------------
Mass fraction 0:                                    0.1  -
Pressure:                              9.0404345635e\+05  Pa
Temperature:                                2547.214296  K
Density:                                    1.187005237  kg/m\^3
Specific volume:                           0.8424562661  m\^3/kg
Specific internal energy:              2.4771659033e\+06  J/kg
Specific enthalpy:                     3.2387829780e\+06  J/kg

Sound speed:                                997.8878004  m/s
Dynamic viscosity:                     1.8230000000e-05  Pa-s
Specific heat at constant pressure:              1271.5  J/\(kg-K\)
Specific heat at constant volume:                 972.5  J/\(kg-K\)
Thermal conductivity:                           0.02568  W/\(m-K\)'

  requirement = "The fluid properties interrogator shall output static-state, single-phase fluid properties for (rho, e) input with a vapor mixture fluid properties object."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./vapor_mixture_rho_e.json]
  type = 'RunApp'
  input = 'vapor_mixture_rho_e.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"static":{"T":2547.[0-9]+,"c":997.[0-9]+,"cp":1271.[0-9]+,"cv":972.[0-9]+,"e":2477165.[0-9]+,"h":3238782.[0-9]+,"k":0.[0-9]+,"mu":1.[0-9]+e-05,"p":904043.[0-9]+,"rho":1.[0-9]+,"v":0.[0-9]+}}'

  requirement = "The fluid properties interrogator shall output static-state, single-phase fluid properties for (rho, e) input with a vapor mixture fluid properties object in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph_p]
  type = 'RunApp'
  input = '2ph.p.i'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'TWO-PHASE properties:
--------------------------------------------------------------------------------
Critical pressure:                                   25  Pa
Saturation temperature:                2.0000000000e\+05  K
Latent heat of vaporization:           4.5931000000e\+08  J/kg'

  requirement = "The fluid properties interrogator shall output two-phase fluid properties for (p) input with a two-phase fluid properties object."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph_p.json]
  type = 'RunApp'
  input = '2ph.p.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"2-phase":{"T_sat":200000.0,"h_lat":459309999.[0-9]+,"p_critical":25.0}}'

  requirement = "The fluid properties interrogator shall output two-phase fluid properties for (p) input with a two-phase fluid properties object in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph_T]
  type = 'RunApp'
  input = '2ph.T.i'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = 'TWO-PHASE properties:
--------------------------------------------------------------------------------
Critical pressure:                                   25  Pa
Saturation pressure:                                900  Pa
Latent heat of vaporization:           6.8896500000e\+05  J/kg
Surface tension:                                   1500  N/m'

  requirement = "The fluid properties interrogator shall output two-phase fluid properties for (T) input with a two-phase fluid properties object."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./2ph_T.json]
  type = 'RunApp'
  input = '2ph.T.i'
  cli_args = 'FluidPropertiesInterrogator/json=true'
  allow_test_objects = true
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_out = '{"2-phase":{"h_lat":688964.[0-9]+,"p_critical":25.0,"p_sat":900.0,"sigma":1500.0}}'

  requirement = "The fluid properties interrogator shall output two-phase fluid properties for (T) input with a two-phase fluid properties object in JSON format."
  issues = '#13741'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./err.wrong_fp_type]
  type = RunException
  input = 'err.wrong_fp_type.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_err = "The type of the parameter 'fp' must be derived from type 'SinglePhaseFluidProperties', 'VaporMixtureFluidProperties', or 'TwoPhaseFluidProperties'"

  requirement = "The fluid properties interrogator shall throw an error if an
      incompatible fluid properties object is supplied."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./err.extraneous_parameter]
  type = RunException
  input = '1ph.rho_e.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_err = "\(rho\,e\) has been specified\; T cannot be specified\."
  cli_args = "FluidPropertiesInterrogator/T=300"

  requirement = "The fluid properties interrogator shall throw an error if an
      extraneous parameter(s) are supplied."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[./err.no_params]
  type = RunException
  input = 'err.no_params.i'
  allow_warnings = true
  method = 'OPT'
  threading = '!pthreads'
  expect_err = "For one-phase fluid properties\, you must provide one of the following
combinations of thermodynamic properties:"

  requirement = "The fluid properties interrogator shall throw an error if an
      no valid input sets were supplied."
  issues = '#12995'
  design = 'FluidPropertiesInterrogator.md'
[../]

[liquid]
  type = 'CSVDiff'
  input = 'saturation_density_function.i'
  csvdiff = 'liquid.csv'
  allow_test_objects = True
  detail = 'for the liquid phase'
[]

[vapor]
  type = 'CSVDiff'
  input = 'saturation_density_function.i'
  csvdiff = 'vapor.csv'
  cli_args = "Functions/rho_sat_fn/use_liquid=false Outputs/file_base=vapor"
  allow_test_objects = True
  detail = 'for the vapor phase'
[]

[./test]
  type = 'CSVDiff'
  input = 'saturation_pressure_function.i'
  csvdiff = 'saturation_pressure_function_out.csv'
  allow_test_objects = True

  requirement = 'The system shall provide a function that computes saturation pressure from a temperature function'
  design = '/SaturationPressureFunction.md'
  issues = '#14755'
[../]

[./test]
  type = 'CSVDiff'
  input = 'saturation_temperature_function.i'
  csvdiff = 'saturation_temperature_function_out.csv'
  allow_test_objects = True

  requirement = 'The system shall provide a function that computes saturation temperature from a pressure function'
  design = '/SaturationTemperatureFunction.md'
  issues = '#14755'
[../]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-10
  requirement = 'The system shall be able to set an initial condition for density given pressure and temperature as variables'
  design = 'RhoFromPressureTemperatureIC.md'
  issues = '#15524'
[]

[./test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-10
  requirement = 'The system shall be able to set an initial condition for density of vapor mixture given pressure and temperature as variables'
  design = 'RhoVaporMixtureFromPressureTemperatureIC.md'
  issues = '#15524'
[../]

[test]
  type = 'CSVDiff'
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-10
  requirement = 'The system shall be able to set an initial condition for specific enthalpy given pressure and temperature as variables'
  design = 'SpecificEnthalpyFromPressureTemperatureIC.md'
  issues = '#15561'
[]

[conservative]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute the fluid properties of an ideal gas with the specific energy and volume variables.'
[]

[primitive]
  type = Exodiff
  input = 'test2.i'
  exodiff = 'test2_out.e'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute the fluid properties of an ideal gas with the pressure and temperature variables.'
[]

[quiet_nan_scalar]
  type = 'RunException'
  input = 'nan_interface.i'
  cli_args = 'FluidProperties/fp/emit_on_nan=warning'
  allow_test_objects = True
  method = 'DBG'
  threading = '!pthreads'

  requirement = 'The system should produce a warning when a scalar NaN is produced and user required that the execution would not terminate'
  expect_err = "fp: A NaN was produced."
  issues = '#12234 #12350'
[]

[quiet_nan_vector]
  type = 'RunException'
  input = 'nan_interface.i'
  cli_args = 'FluidProperties/fp/emit_on_nan=warning Kernels/test_kernel/test_vector_version=true'
  allow_test_objects = True
  method = 'DBG'
  threading = '!pthreads'

  requirement = 'The system should produce a warning when a vector NaN is produced and user required that the execution would not terminate'
  expect_err = "fp: A NaN was produced."
  issues = '#12234 #12350'
[]

[signaling_nan_dbg]
  type = 'RunException'
  input = 'nan_interface.i'
  allow_test_objects = True
  method = 'DBG'
  threading = '!pthreads'

  requirement = 'The system should report an error when a NaN is produced by a computation in DEBUG mode, by default'
  expect_err = "fp: A NaN was produced."
  issues = '#12234 #12350'
[]

[signaling_nan_opt]
  type = 'RunApp'
  input = 'nan_interface.i'
  allow_test_objects = True
  method = 'OPT'
  threading = '!pthreads'

  requirement = 'The system should not report an error when a NaN is produced by a computation in OPT mode, by default'
  issues = '#12234 #12350'

  # Make sure this test does not run in --recover mode because the default
  # value of execute_on for checkpoint is TIMESTEP_END. This test will fail
  # with a nan before it solves a timestep and writes a checkpoint, causing
  # the test to fail in --recover mode.
  recover = False
[]

[signaling_nan_warning]
  type = 'RunApp'
  input = 'nan_interface.i'
  cli_args = 'FluidProperties/fp/emit_on_nan=warning'
  allow_test_objects = True
  allow_warnings = true
  expect_out = 'fp: A NaN was produced.'

  requirement = 'The system should be able to return a warning when a fluid property evaluation encounters a NaN.'
  issues = '#12234 #12350'
  recover = False
[]

[signaling_nan_exception]
  type = 'RunException'
  input = 'nan_interface.i'
  cli_args = 'FluidProperties/fp/emit_on_nan=exception Executioner/type=Transient'
  allow_test_objects = True
  allow_warnings = true
  expect_err = 'Solve failed and timestep already at or below dtmin'

  requirement = 'The system should be able to throw an exception when a fluid property evaluation encounters a NaN.'
  issues = '#12234 #12350'
  recover = False
[]

[./test]
  type = 'CSVDiff'
  input = 'ad_saturation_temperature_material.i'
  csvdiff = 'ad_saturation_temperature_material_out.csv'
  allow_test_objects = True
  recover = false

  requirement = 'The system shall provide an AD material that computes saturation temperature.'
  design = '/ADSaturationTemperatureMaterial.md'
  issues = '#15308'
[../]

[./test]
  type = 'CSVDiff'
  input = 'ad_surface_tension_material.i'
  csvdiff = 'ad_surface_tension_material_out.csv'
  allow_test_objects = True
  recover = false

  requirement = 'The system shall provide an AD material that computes surface tension.'
  design = '/ADSurfaceTensionMaterial.md'
  issues = '#15308'
[../]

[pressure_temperature]
  type = Exodiff
  input = 'test_pt.i'
  exodiff = 'test_pt_out.e'
  detail = 'using a (pressure, temperature) variable set'
[]

[pressure_enthalpy]
  type = Exodiff
  input = 'test_ph.i'
  exodiff = 'test_ph_out.e'
  detail = 'using a (pressure, enthalpy) variable set'
  allow_test_objects = true
[]

[volume_energy]
  type = Exodiff
  input = 'test_ve.i'
  exodiff = 'test_ve_out.e'
  detail = 'using a (specific volume, specific internal energy) variable set'
[]

[ad]
  type = 'CSVDiff'
  input = 'saturation_pressure_material.i'
  csvdiff = 'saturation_pressure_material_out.csv'
  recover = false

  requirement = 'The system shall provide a material that computes saturation pressure using automatic differentiation material properties.'
[]

[nonad]
  type = 'CSVDiff'
  input = 'saturation_pressure_material.i'
  cli_args = "
      Materials/T_mat/type=GenericConstantMaterial
      Materials/p_sat_mat/type=SaturationPressureMaterial
      Postprocessors/p_sat_pp/type=ElementAverageMaterialProperty"
  csvdiff = 'saturation_pressure_material_out.csv'
  prereq = 'ad'
  recover = false

  requirement = 'The system shall provide a material that computes saturation pressure using non-automatic differentiation material properties.'
[]

[methane]
  type = Exodiff
  input = 'methane.i'
  exodiff = 'methane_out.e'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute fluid properties for methane in a gaseous phase.'
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-8
  allow_test_objects = true
  threading = '!pthreads'
  required_objects = 'SalineMoltenSaltFluidProperties'
  required_submodule = 'contrib/saline'
  requirement = "The system shall be able to leverage the Saline submodule to compute composed salt fluid properties"
  design = "SalineMoltenSaltFluidProperties.md"
  issues = "#22724"
[]

[exact]
  type = CSVDiff
  input = 'exact.i'
  csvdiff = 'exact_out.csv'
  abs_zero = 1e-9
  # At the initial time, an assertion is triggered regarding the temperature
  # exceeding the critical temperature. The temperature variable stays within
  # this critical temperature, but temperature computed from the given enthalpy
  # does not.
  method = 'OPT'
  requirement = 'The system shall be able to compute liquid sodium properties and compare exactly to analytical expressions.'
[]

[constant]
  type = CSVDiff
  input = 'constant.i'
  csvdiff = 'constant_out.csv'
  requirement = 'The system shall be able to compute liquid sodium properties given constant thermal conductivity and specific heat values.'
[]

[test]
  type = Exodiff
  input = 'test.i'
  exodiff = 'test_out.e'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute fluid properties with a stiffened gas model using a conservative variable set.'
[]

[tabulated_bicubic]
  type = CSVDiff
  input = 'tabulated.i'
  csvdiff = 'tabulated_pT_out.csv'
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                FluidProperties/tabulated/fp=co2
                FluidProperties/tabulated/fluid_property_output_file='fluid_properties.csv'
                Outputs/file_base=tabulated_pT_out"
  rel_err = 1e-4
  threading = '!pthreads'
  requirement = 'The system shall be able to generate a tabulation of fluid properties based on pressure and temperature, and use a tabulated interplation of these properties.'
  # Steady solve
  recover = false
[]

[regular_grids]
  type = CSVDiff
  input = 'tabulated_v_e.i'
  csvdiff = 'tabulated_v_e_out.csv'
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                  FluidProperties/tabulated/fluid_property_file='fluid_properties.csv'
                  Outputs/file_base=tabulated_v_e_out"
  rel_err = 1e-4
  threading = '!pthreads'
  allow_warnings = true
  prereq = tabulated_bicubic
  detail = 'using a regularly spaced interpolation grid'
[]

[log_v]
  type = CSVDiff
  input = 'tabulated_v_e.i'
  csvdiff = 'tabulated_v_e_logv_out.csv'
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                  FluidProperties/tabulated/fluid_property_file='fluid_properties.csv'
                  Outputs/file_base=tabulated_v_e_logv_out
                  FluidProperties/tabulated/use_log_grid_v=true"
  rel_err = 1e-4
  threading = '!pthreads'
  allow_warnings = true
  prereq = tabulated_bicubic
  detail = 'using a log-spaced specific volume grid and a regularly spaced grid in specific internal energy'
[]

[log_e]
  type = CSVDiff
  input = 'tabulated_v_e.i'
  csvdiff = 'tabulated_v_e_loge_out.csv'
  # raise the minimum temperature to avoid negative internal energy
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                  FluidProperties/tabulated/fp=co2
                  FluidProperties/tabulated/temperature_min=400
                  FluidProperties/tabulated/temperature_max=800
                  FluidProperties/tabulated/pressure_min=2e5
                  FluidProperties/tabulated/use_log_grid_e=true
                  Outputs/file_base=tabulated_v_e_loge_out"
  rel_err = 1e-4
  threading = '!pthreads'
  allow_warnings = true
  prereq = tabulated_bicubic
  detail = 'using a linear specific volume grid and a log-spaced grid in specific internal energy'
[]

[ve]
  type = CSVDiff
  input = 'tabulated_v_e.i'
  csvdiff = 'fluid_properties_out_ve.csv'
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                  FluidProperties/tabulated/num_v=10
                  FluidProperties/tabulated/num_e=10
                  FluidProperties/tabulated/fluid_property_file='fluid_properties.csv'
                  FluidProperties/tabulated/fluid_property_output_file='fluid_properties_out.csv'"
  detail = 'for the (specific volume, internal energy) conversion grids.'
  prereq = tabulated_bicubic

  # Some conversions are going to fail
  allow_warnings = true
[]

[tabulated_v_e_bicubic_direct]
  # Note: the output of this test depends on the gold file from the file_writing/ve
  type = CSVDiff
  requirement = 'The system shall be able to read a tabulation of fluid properties with a (specific volume, specific internal energy) grid and use this tabulation directly.'
  input = 'tabulated_v_e.i'
  csvdiff = 'tabulated_direct_v_e_out.csv'
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                FluidProperties/tabulated/create_pT_interpolations=false
                FluidProperties/tabulated/construct_pT_from_ve=false
                FluidProperties/tabulated/construct_pT_from_vh=false
                FluidProperties/tabulated/create_ve_interpolations=true
                FluidProperties/tabulated/fluid_property_ve_file='gold/fluid_properties_out_ve.csv'
                Outputs/file_base=tabulated_direct_v_e_out"
  rel_err = 1e-4
  threading = '!pthreads'
  allow_warnings = true
[]

[bad_pressure_bounds]
  type = RunException
  input = 'tabulated.i'
  cli_args = 'FluidProperties/tabulated/pressure_min=1e6 FluidProperties/tabulated/pressure_max=1e5'
  expect_err = 'pressure_max must be greater than pressure_min'
  detail = 'if the user-specified pressure boundaries are inverted'
[]

[bad_temperature_bounds]
  type = RunException
  input = 'tabulated.i'
  cli_args = 'FluidProperties/tabulated/temperature_min=1e3 FluidProperties/tabulated/temperature_max=1e2'
  expect_err = 'temperature_max must be greater than temperature_min'
  detail = 'if the user-specified temperature boundaries are inverted'
[]

[bad_pressure_guess]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/p_initial_guess=10 FluidProperties/tabulated/save_file=false FluidProperties/tabulated/fp=co2"
  expect_err = 'Pressure initial guess for \(p,T\), \(v,e\) conversions 10 is outside the range of tabulated pressure'
  detail = 'if the user-specified pressure initial guess for variable set inversions is out of bounds'
[]

[bad_temperature_guess]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/T_initial_guess=1e5 FluidProperties/tabulated/save_file=false FluidProperties/tabulated/fp=co2"
  expect_err = 'Temperature initial guess for \(p,T\), \(v,e\) conversions 100000 is outside the range of tabulated temperature'
  detail = 'if the user-specified temperature initial guess for variable set inversions is out of bounds'
[]

[error_pressure_out_of_bounds]
  type = RunException
  input = 'tabulated.i'
  cli_args = 'AuxVariables/pressure/initial_condition=1e4 FluidProperties/tabulated/error_on_out_of_bounds=true'
  expect_err = 'is outside the range of tabulated pressure'
  detail = 'error if the desired pressure for a fluid property evaluation is outside the user-specified bounds, when prescribed to error in such conditions'
[]

[error_temperature_out_of_bounds]
  type = RunException
  input = 'tabulated.i'
  cli_args = 'AuxVariables/temperature/initial_condition=250 FluidProperties/tabulated/error_on_out_of_bounds=true'
  expect_err = 'is outside the range of tabulated temperature'
  detail = 'error if the desired temperature for a fluid property evaluation is outside the user-specified bounds, when prescribed to error in such conditions'
[]

[invalid_pressure_out_of_bounds]
  type = RunException
  input = 'tabulated.i'
  cli_args = "AuxVariables/pressure/initial_condition=1e4 FluidProperties/tabulated/out_of_bounds_behavior=declare_invalid Problem/solve=true Kernels/diff/type=MatDiffusion Kernels/diff/diffusivity=viscosity"
  expect_err = 'Pressure out of bounds'
  detail = 'declare the solution invalid if the desired pressure for a fluid property evaluation is outside the user-specified bounds, when prescribed to do so in such conditions'
[]

[invalid_temperature_out_of_bounds]
  type = RunException
  input = 'tabulated.i'
  cli_args = "AuxVariables/temperature/initial_condition=250 FluidProperties/tabulated/out_of_bounds_behavior=declare_invalid Problem/solve=true Kernels/diff/type=MatDiffusion Kernels/diff/diffusivity=viscosity"
  expect_err = 'Temperature out of bounds'
  detail = 'declare the solution invalid if the desired temperature for a fluid property evaluation is outside the user-specified bounds, when prescribed to do so in such conditions'
[]

[ignore_pressure_out_of_bounds]
  type = RunApp
  input = 'tabulated.i'
  cli_args = 'AuxVariables/pressure/initial_condition=1e4 FluidProperties/tabulated/out_of_bounds_behavior=ignore '
  detail = 'ignore if the desired pressure for a fluid property evaluation is outside the user-specified bounds, when prescribed to ignore in such conditions'
[]

[ignore_temperature_out_of_bounds]
  type = RunApp
  input = 'tabulated.i'
  cli_args = 'AuxVariables/temperature/initial_condition=250 FluidProperties/tabulated/out_of_bounds_behavior=ignore'
  detail = 'ignore if the desired temperature for a fluid property evaluation is outside the user-specified bounds, when prescribed to ignore in such conditions'
[]

[exact]
  type = CSVDiff
  input = 'exact.i'
  csvdiff = 'exact_out.csv'
  abs_zero = 1e-9
  requirement = 'The system shall be able to compute functionalized fluid properties.'
[]

[example]
  type = CSVDiff
  input = 'example.i'
  csvdiff = 'example_out.csv'
  abs_zero = 1e-9
  requirement = "The system shall be able to compute fluid properties defined through functions of pressure and temperature."
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  rel_err = 1e-8
  allow_test_objects = true
  threading = '!pthreads'
  requirement = 'The system shall be able to mix two independent single-phase fluid properties to form a basic two-phase model.'
[]

[no_error_on_unimplemented]
  type = RunApp
  input = 'test.i'
  cli_args = "
      FluidProperties/fp_2phase/error_on_unimplemented=false
      AuxVariables/T_sat/family=LAGRANGE
      AuxVariables/T_sat/order=FIRST
      AuxKernels/T_sat_aux/type=SaturationTemperatureAux
      AuxKernels/T_sat_aux/variable=T_sat
      AuxKernels/T_sat_aux/p=p
      AuxKernels/T_sat_aux/fp_2phase=fp_2phase
      AuxKernels/T_sat_aux/execute_on=initial
      Outputs/csv=false"
  allow_test_objects = true
  threading = '!pthreads'
  # this test modifies the checkpoint file causing recovery of 'test' to fail
  prereq = test
  requirement = 'The system shall be able to not error on unimplemented two-phase fluid properties routines if the user requested so.'
[]

[error:error_on_unimplemented]
  type = 'RunException'
  input = 'test.i'
  cli_args = "
      AuxVariables/T_sat/family=LAGRANGE
      AuxVariables/T_sat/order=FIRST
      AuxKernels/T_sat_aux/type=SaturationTemperatureAux
      AuxKernels/T_sat_aux/variable=T_sat
      AuxKernels/T_sat_aux/p=p
      AuxKernels/T_sat_aux/fp_2phase=fp_2phase
      AuxKernels/T_sat_aux/execute_on=initial"
  expect_err = "The 2-phase fluid properties class 'TwoPhaseFluidPropertiesIndependent' does not allow calling any 2-phase property interfaces"
  allow_test_objects = true
  threading = '!pthreads'
  requirement = 'The system shall report an error if 2-phase property interfaces are called on a class designed to handle two single-phase fluid properties.'
[]

[test]
  type = CSVDiff
  input = 'test.i'
  csvdiff = 'test_out.csv'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute properties from a two-phase with non-condensable gas user object.'
[]

[error_wrong_number_of_args]
  type = RunException
  input = 'test.i'
  cli_args = "Functions/p_sat_fn/arg2=40000"
  expect_err = "The property call 'p_sat' expects the parameter\(s\) {'arg1'} to be provided, but the provided argument\(s\) were {'arg1', 'arg2'}"
  allow_test_objects = true
  threading = '!pthreads'
  requirement = 'The system shall report an error if an incorrect number of arguments are provided when computing properties from a two-phase with non-condensable gas user object.'
[]

[water]
  type = CSVDiff
  input = 'water.i'
  csvdiff = 'water_out.csv'
  threading = '!pthreads'
  requirement = 'The system shall be able to compute the fluid properties of water and steam following the IAPWS-IF97 formulations.'
[]