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

[co2]
  type = CSVDiff
  input = 'co2.i'
  csvdiff = 'co2_out.csv'
  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'
  detail = 'spline-based table-lookup air fluid properties,'
  capabilities = 'airapp'
[]

[CO2_pT]
  type = CSVDiff
  input = 'co2.i'
  csvdiff = 'co2_out.csv'
  detail = 'vapor carbon dioxide fluid properties,'
  capabilities = 'carbondioxideapp'
[]

[helium_pT]
  type = CSVDiff
  input = 'helium.i'
  csvdiff = 'helium_out.csv'
  detail = 'spline-based table-lookup Helium fluid properties,'
  capabilities = 'heliumapp'
[]

[nitrogen_pT]
  type = CSVDiff
  input = 'nitrogen.i'
  csvdiff = 'nitrogen_out.csv'
  detail = 'spline-based table-lookup Nitrogen fluid properties,'
  capabilities = 'nitrogenapp'
[]

[potassium_pT]
  type = CSVDiff
  input = 'potassium.i'
  csvdiff = 'potassium_out.csv'
  detail = 'liquid Potassium fluid properties,'
  capabilities = 'potassiumapp'
[]

[sodium_pT]
  type = CSVDiff
  input = 'sodium.i'
  csvdiff = 'sodium_out.csv'
  detail = 'liquid Sodium fluid properties.'
  capabilities = 'sodiumapp'
[]

[error_nonideal_gases]
  type = RunException
  input = 'ideal_gas_mixture.i'
  expect_err = "Each entry in 'component_fluid_properties' must be an ideal gas that implements IdealGasMixtureInterface"
  requirement = 'The system shall report an error if non-ideal gases are used with the ideal gas mixture model.'
[]

[./1ph.p_T]
  type = 'RunApp'
  input = '1ph.p_T.i'
  allow_warnings = true
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
Specific entropy:                      1.4838960040e\+04  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_partial_pressure_p_T]
  type = RunApp
  input = '2ph_ncg_partial_pressure_p_T.i'
  allow_test_objects = true
  capabilities = 'method=opt'
  expect_out = 'Mass fraction 0:\s+0.9930237384'
  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
  capabilities = 'method=opt'
  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]+,"s":14838.[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
  capabilities = 'method=opt'
  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
Specific entropy:                           5966.595696  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
  capabilities = 'method=opt'
  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]+,"s":5966.[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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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
  capabilities = 'method=opt'
  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'
[../]

[test]
  type = CSVDiff
  input = 'hydrostatic_pressure_function.i'
  csvdiff = 'hydrostatic_pressure_function_out_p_vpp_0000.csv'
  requirement = 'The system shall provide a function to compute the hydrostatic pressure distribution.'
[]

[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'
  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'
  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
  capabilities = 'method=dbg'

  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
  capabilities = 'method=dbg'

  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
  capabilities = 'method=dbg'

  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 = 'RunException'
  input = 'nan_interface.i'
  allow_test_objects = True
  capabilities = 'method=opt'
  expect_err = 'Aborting as solve did not converge'

  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'
  # A NaN is produced by the fluid property in the NaNInterfaceTestKernel, preventing the convergence
  # Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 0
  expect_exit_code = 1
  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'
  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
  capabilities = '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.
  capabilities = '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'
  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/input_fp=co2
                FluidProperties/tabulated/fluid_property_output_file='fluid_properties.csv'
                Outputs/file_base=tabulated_pT_out"
  rel_err = 1e-4
  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
  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
  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/input_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
  allow_warnings = true
  prereq = tabulated_bicubic
  detail = 'using a linear specific volume grid and a log-spaced grid in specific internal energy'
[]

[pT]
  type = CSVDiff
  input = 'tabulated.i'
  csvdiff = 'fluid_properties_10_10_out.csv'
  cli_args = "FluidProperties/tabulated/type=TabulatedBicubicFluidProperties
                  FluidProperties/tabulated/num_p=10
                  FluidProperties/tabulated/num_T=10
                  FluidProperties/tabulated/skip_header_tabulation=true
                  FluidProperties/tabulated/fluid_property_output_file='fluid_properties_10_10_out.csv'"
  detail = 'for the (pressure, temperature) interpolation grids.'
  prereq = tabulated_bicubic

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

[ve]
  type = CSVDiff
  input = 'tabulated_v_e.i'
  csvdiff = 'fluid_properties_ve_10_10_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/skip_header_tabulation=true
                  FluidProperties/tabulated/fluid_property_ve_file='gold/fluid_properties_out_ve.csv'
                  FluidProperties/tabulated/create_ve_interpolations=true
                  FluidProperties/tabulated/fluid_property_ve_output_file='fluid_properties_ve_10_10_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
  allow_warnings = true
  prereq = 'file_writing/ve'
[]

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

[missing_energy_bound]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/fluid_property_file='fluid_properties.csv'
                  FluidProperties/tabulated/e_min=100"
  expect_err = 'Either both or none of the min and max values of the specific internal energy should be specified'
  detail = 'if the user only specified one bound for a range of specific internal energy,'
[]

[missing_specific_volume_bound]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/v_min=1e3
                  FluidProperties/tabulated/fluid_property_file='fluid_properties.csv'"
  expect_err = 'Either both or none of the min and max values of the specific volume should be specified'
  detail = 'if the user only specified one bound for a range of specific volume.'
[]

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

[warn_pressure_out_of_bounds]
  type = RunApp
  input = 'tabulated.i'
  cli_args = "AuxVariables/pressure/initial_condition=1e4 FluidProperties/tabulated/out_of_bounds_behavior=warn_invalid Problem/solve=true Kernels/diff/type=MatDiffusion Kernels/diff/diffusivity=viscosity"
  expect_out = 'Pressure out of bounds'
  detail = 'count the number of occurrences if the desired pressure for a fluid property evaluation is outside the user-specified bounds, when prescribed to do so in such conditions,'
[]

[warn_temperature_out_of_bounds]
  type = RunApp
  input = 'tabulated.i'
  cli_args = "AuxVariables/temperature/initial_condition=250 FluidProperties/tabulated/out_of_bounds_behavior=warn_invalid Problem/solve=true Kernels/diff/type=MatDiffusion Kernels/diff/diffusivity=viscosity"
  expect_out = 'Temperature out of bounds'
  detail = 'count the number of occurrences 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'
[]

[num_grids_ve]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/fluid_property_ve_file='gold/fluid_properties_out_ve.csv'
                  FluidProperties/tabulated/num_v=200"
  expect_err = 'User-specified grid sizes in specific volume and internal energy are ignored'
  allow_warnings = false
  detail = 'the user specifies a grid size when reading the interpolation grid from file,'
[]

[min_max_ve]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/fluid_property_ve_file='gold/fluid_properties_out_ve.csv'
                  FluidProperties/tabulated/v_min=20 FluidProperties/tabulated/v_max=20"
  expect_err = 'User-specified bounds in specific volume and internal energy are ignored'
  allow_warnings = false
  detail = 'the user specifies a grid extent when reading the interpolation grid from file,'
[]

[log_ve]
  type = RunException
  input = 'tabulated_v_e.i'
  cli_args = "FluidProperties/tabulated/fluid_property_ve_file='gold/fluid_properties_out_ve.csv'
                  FluidProperties/tabulated/use_log_grid_v=true"
  expect_err = 'User specfied logarithmic grids in specific volume and energy are ignored'
  allow_warnings = false
  detail = 'the user specifies that the grid is logarithmic when reading the interpolation grid from file.'
[]

[exact]
  type = CSVDiff
  input = 'exact.i'
  csvdiff = 'exact_out.csv'
  abs_zero = 1e-9
  ignore_columns = cp_avg
  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
  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
  # 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
  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'
  allow_test_objects = true
  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
  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'
  requirement = 'The system shall be able to compute the fluid properties of water and steam following the IAPWS-IF97 formulations.'
[]

[ad_derivatives_after_value_call]
  type = GoogleTest
  unit_test = BrineFluidPropertiesTest.adDerivativesAfterValueCall
  requirement = 'The system shall include the pressure and temperature derivatives of the water density and enthalpy in the brine density and enthalpy computed with forward-mode automatic differentiation, regardless of which brine property methods were called before.'
[]

[derivatives_near_critical]
  type = GoogleTest
  unit_test = CO2FluidPropertiesTest.derivativesNearCritical
  requirement = 'The system shall compute derivatives of the CO2 density, internal energy and enthalpy with respect to pressure and temperature, and an isobaric heat capacity equal to the temperature derivative of the enthalpy, near the critical point where the nonanalytic terms of the Span and Wagner equation of state are significant.'
[]

[sound_speed_varying_gamma]
  type = GoogleTest
  unit_test = CaloricallyImperfectGasTest.soundSpeedVaryingGamma
  requirement = 'The system shall compute the derivatives of the speed of sound of a calorically imperfect gas with respect to specific volume and internal energy when the specific heat ratio varies with temperature.'
[]

[lookup_table_bounds]
  type = GoogleTest
  unit_test = CaloricallyImperfectGasTest.lookupTableBounds
  requirement = 'The system shall evaluate the temperature and entropy of a calorically imperfect gas at the bounds of its tabulated temperature range, and return values clamped to those bounds outside of it when out-of-bounds errors are disabled.'
[]

[derivatives_at_edge_states]
  type = GoogleTest
  unit_test = CaloricallyImperfectGasTest.derivativesAtEdgeStates
  requirement = 'The system shall compute finite derivatives of the pressure from enthalpy and entropy of a calorically imperfect gas at zero entropy, and of its isobaric heat capacity and viscosity at the maximum tabulated temperature.'
[]

[derivatives_near_critical]
  type = GoogleTest
  unit_test = HydrogenFluidPropertiesTest.derivativesNearCritical
  requirement = 'The system shall compute derivatives of the hydrogen density, internal energy and enthalpy with respect to pressure and temperature, and an isobaric heat capacity equal to the temperature derivative of the enthalpy, near the critical point.'
[]

[gibbs_free_energy]
  type = GoogleTest
  unit_test = IdealGasFluidPropertiesTest.gibbsFreeEnergy
  requirement = 'The system shall compute the Gibbs free energy of an ideal gas consistently with its enthalpy and entropy, including a nonzero reference internal energy.'
[]

[sound_speed_nonphysical_state]
  type = GoogleTest
  unit_test = IdealGasFluidPropertiesTest.soundSpeedNonphysicalState
  requirement = 'The system shall return a finite speed of sound and finite derivatives of it for an ideal gas state with negative temperature, consistently between the methods with and without derivatives.'
[]

[ad_derivatives]
  type = GoogleTest
  unit_test = IdealGasMixtureFluidPropertiesTest.adDerivatives
  requirement = 'The system shall compute the derivatives of the temperature and pressure from the specific volume and internal energy, of the internal energy from the pressure and density, and of the sound speed from the pressure and temperature of an ideal gas mixture with respect to their arguments and the mass fractions with automatic differentiation, in agreement with finite differences.'
[]

[deprecated_specific_heat_ratio]
  type = GoogleTest
  unit_test = IdealGasMixtureFluidPropertiesTest.deprecatedSpecificHeatRatio
  requirement = 'The system shall continue to compute the specific heat ratio of an ideal gas mixture without pressure and temperature arguments when all components have constant specific heats, and shall report an error when a component has temperature-dependent specific heats.'
[]

[nasa7_ad_derivatives]
  type = GoogleTest
  unit_test = IdealGasMixtureFluidPropertiesTest.nasa7AdDerivatives
  requirement = 'The system shall compute the derivatives of the temperature and pressure from the specific volume and internal energy, of the internal energy from the pressure and density, and of the sound speed from the pressure and temperature of a mixture of ideal gases from NASA 7-coefficient polynomials with respect to their arguments and the mass fractions with automatic differentiation, in agreement with finite differences.'
[]

[nasa7_cantera]
  type = GoogleTest
  unit_test = IdealGasMixtureFluidPropertiesTest.nasa7Cantera
  requirement = 'The system shall compute the density, specific internal energy, isobaric and isochoric specific heats, and sound speed of a mixture of ideal gases from NASA 7-coefficient polynomials, reproducing Cantera values, and the temperature and pressure from its specific volume and internal energy.'
[]

[nasa7_inversion]
  type = GoogleTest
  unit_test = IdealGasMixtureFluidPropertiesTest.nasa7Inversion
  requirement = 'The system shall compute the temperature of mixtures of ideal gases from NASA 7-coefficient polynomials, alone or with an ideal gas of constant specific heats, from their specific volume and internal energy, from half the highest lower end to 1.5 times the lowest upper end of the temperature ranges of the data of their components, and where the internal energy is close to zero.'
[]

[cp_volume_dependent_enthalpy]
  type = GoogleTest
  unit_test = IdealRealGasMixtureFluidPropertiesTest.cpVolumeDependentEnthalpy
  requirement = 'The system shall compute the isobaric heat capacity of an ideal-real gas mixture whose primary component enthalpy depends on specific volume.'
[]

[transport_volume_dependent_secondary]
  type = GoogleTest
  unit_test = IdealRealGasMixtureFluidPropertiesTest.transportVolumeDependentSecondary
  requirement = 'The system shall compute the same viscosity and thermal conductivity of an ideal-real gas mixture with and without derivatives when a secondary component internal energy depends on specific volume.'
[]

[three_components]
  type = GoogleTest
  unit_test = IdealRealGasMixtureFluidPropertiesTest.threeComponents
  requirement = 'The system shall compute derivatives of the pressure, internal energy and entropy of an ideal-real gas mixture of three components with respect to temperature, specific volume and the secondary mass fractions.'
[]

[derivatives_near_critical]
  type = GoogleTest
  unit_test = MethaneFluidPropertiesTest.derivativesNearCritical
  requirement = 'The system shall compute derivatives of the methane density, internal energy and enthalpy with respect to pressure and temperature, and an isobaric heat capacity equal to the temperature derivative of the enthalpy, near the critical point where the Gaussian terms of the Setzmann and Wagner equation of state are significant.'
[]

[cantera]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.cantera
  requirement = 'The system shall compute the density, specific internal energy, specific enthalpy, specific entropy, isobaric and isochoric specific heats, and sound speed of an ideal gas from NASA 7-coefficient polynomials, reproducing Cantera values below and above the temperature where the two polynomials meet and at a pressure other than the reference pressure.'
[]

[consistency]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.consistency
  requirement = 'The system shall compute the same state of an ideal gas from NASA 7-coefficient polynomials with the (pressure, temperature), (temperature, specific volume), and (pressure, density) variable sets, and satisfy the definition of enthalpy, the Mayer relation, the ideal-gas sound speed, and the relations between the temperature and pressure derivatives of enthalpy, internal energy, and entropy and the specific heats and specific volume.'
[]

[derivatives]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.derivatives
  requirement = 'The system shall compute the derivatives of the properties of an ideal gas from NASA 7-coefficient polynomials with respect to their arguments, with real and automatic differentiation arguments, in agreement with finite differences.'
[]

[out_of_range]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.outOfRange
  requirement = 'The system shall extrapolate NASA 7-coefficient polynomials outside the temperature range of their data and count each such evaluation, including the evaluation of a temperature found by inverting the internal energy, as a solution warning, and shall flag a non-positive or not-a-number temperature as an invalid solution.'
[]

[inverses]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.inverses
  requirement = 'The system shall compute the temperature of an ideal gas from NASA 7-coefficient polynomials from its specific volume and internal energy and from its pressure and enthalpy by Newton iteration, together with the properties that depend on these variable sets, consistently with the properties from pressure and temperature, and within the mismatch of the two polynomials at the temperature where they meet.'
[]

[energy_near_zero]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.energyNearZero
  requirement = 'The system shall compute the temperature of an ideal gas from NASA 7-coefficient polynomials whose specific internal energy, which includes the enthalpy of formation, is close to zero.'
[]

[inverse_range]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.inverseRange
  requirement = 'The system shall compute the temperature of each built-in NASA 7-coefficient ideal gas from its specific internal energy and from its specific enthalpy from half the lower end to 1.5 times the upper end of the temperature range of its data.'
[]

[inverse_derivatives]
  type = GoogleTest
  unit_test = NASA7IdealGasFluidPropertiesTest.inverseDerivatives
  requirement = 'The system shall compute the derivatives of the properties of an ideal gas from NASA 7-coefficient polynomials from the (specific volume, internal energy), (pressure, enthalpy), (enthalpy, entropy), and (specific volume, enthalpy) variable sets, with real and automatic differentiation arguments, in agreement with finite differences.'
[]

[below_midpoint]
  type = GoogleTest
  unit_test = NASA7PolynomialsTest.canteraBelowMidpoint
  detail = 'below the temperature where the low- and high-temperature polynomials meet,'
[]

[at_midpoint]
  type = GoogleTest
  unit_test = NASA7PolynomialsTest.canteraAtMidpoint
  detail = 'at that temperature, where the low-temperature polynomial applies, and'
[]

[above_midpoint]
  type = GoogleTest
  unit_test = NASA7PolynomialsTest.canteraAboveMidpoint
  detail = 'above that temperature.'
[]

[molar_mass]
  type = GoogleTest
  unit_test = NASA7PolynomialsTest.molarMass
  requirement = 'The system shall provide the molar masses of the built-in NASA 7-coefficient species, matching the values Cantera computes for GRI-Mech 3.0.'
[]

[consistent_derivatives]
  type = GoogleTest
  unit_test = NASA7PolynomialsTest.consistentDerivatives
  requirement = 'The system shall compute NASA 7-coefficient enthalpies and entropies whose temperature derivatives, obtained by automatic differentiation, equal the heat capacity and the heat capacity divided by temperature, respectively.'
[]

[cv]
  type = GoogleTest
  unit_test = NaClFluidPropertiesTest.cv
  requirement = 'The system shall compute the isochoric heat capacity of halite from its isobaric heat capacity, thermal expansion coefficient and isothermal compressibility.'
[]

[energy_derivatives_low_expansion_modulus_product]
  type = GoogleTest
  unit_test = SimpleFluidPropertiesTest.energyDerivativesLowExpansionModulusProduct
  requirement = 'The system shall compute derivatives of the internal energy of a simple fluid with respect to pressure and density when the product of the thermal expansion coefficient and the bulk modulus is small.'
[]

[default_latent_heat_of_fusion]
  type = GoogleTest
  unit_test = StiffenedGasTwoPhaseFluidPropertiesTest.defaultLatentHeatOfFusion
  requirement = 'The system shall use the latent heat of fusion of water as the default latent heat of fusion of the two-phase stiffened gas fluid properties.'
[]

[gibbs_without_temperature_column]
  type = GoogleTest
  unit_test = TabulatedBicubicFluidPropertiesTest.gibbsWithoutTemperatureColumn
  requirement = 'The system shall compute the Gibbs free energy from a tabulation in specific volume and internal energy that does not interpolate temperature.'
[]

[check_nans_midpoint]
  type = GoogleTest
  unit_test = TabulatedBicubicFluidPropertiesTest.checkNaNsMidpoint
  requirement = 'The system shall replace NaN pressures and temperatures from the inversion of tabulated properties in specific volume and internal energy with a bound of the tabulated range at every specific volume, including the midpoint of the specific volume range.'
[]

[heat_capacity_difference]
  type = GoogleTest
  unit_test = TemperaturePressureFunctionFluidPropertiesTest.heatCapacityDifference
  requirement = 'The system shall compute isobaric and isochoric heat capacities from temperature and pressure functions that satisfy cp - cv = T alpha^2 / (rho kappa_T), with alpha the thermal expansion coefficient and kappa_T the isothermal compressibility.'
[]

[derivatives_at_zero_arguments]
  type = GoogleTest
  unit_test = TemperaturePressureFunctionFluidPropertiesTest.derivativesAtZeroArguments
  requirement = 'The system shall compute finite derivatives of the heat capacities and internal energy defined from temperature and pressure functions at zero pressure and zero internal energy.'
[]