StiffenedGasTwoPhaseFluidProperties

This class corresponds to both liquid and vapor phases being approximated using the stiffened gas equation of state (SGEOS).

The saturation curve $var element = document.getElementById("moose-equation-fccea80e-0a6e-4767-b658-0dbb2c2687b4");katex.render("T_\\text{sat}(p)", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is generated by reasoning that at thermodynamic equilibrium, the Gibbs free enthalpy of the phases must be equal at the interface (see Equation (38) of Berry et al. (2010)). A saturation curve is stored at a number of points determined by user-defined parameters, and then linear interpolation is used to get values.

Surface tension is computed using the following form that IAPWS uses for approximating the surface tension of water Water and Steam (2014): $var element = document.getElementById("moose-equation-49bb67ef-fb53-446b-9987-2fc6b87377b6");katex.render(" \\sigma(T) = A \\left(1 - \\frac{T}{T_C}\\right)^B \\left[1 - C \\left(1 - \\frac{T}{T_C}\\right)\\right] \\,,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where

$var element = document.getElementById("moose-equation-0771ac67-b527-409a-8723-d49364f166c4");katex.render("T", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the temperature [K],
$var element = document.getElementById("moose-equation-fda128a1-ded2-4b99-8ee4-350bbeedde53");katex.render("T_C", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the critical temperature [K], and
$var element = document.getElementById("moose-equation-fdfeed89-27a8-4a0e-90d6-69e2fb8fcbec");katex.render("A", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ [N/m], $var element = document.getElementById("moose-equation-7aaa2dda-5105-4055-ba68-cd291e497ba2");katex.render("B", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ [-], and $var element = document.getElementById("moose-equation-0a2bccf7-3d5b-478b-8c19-6f4b8e7ae23b");katex.render("C", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ [-] are a fluid-dependent constants.

Input Parameters

L_fusion0.334Latent heat of fusion [J/kg]
Default:0.334
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Latent heat of fusion [J/kg]
M_liquid0.0180149Liquid molar mass
Default:0.0180149
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid molar mass
M_vapor0.0180149Vapor molar mass
Default:0.0180149
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor molar mass
T_c647.096Critical temperature [K]
Default:647.096
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Critical temperature [K]
T_sat_max647Maximum temperature value in saturation curve [K]
Default:647
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Maximum temperature value in saturation curve [K]
T_sat_min274Minimum temperature value in saturation curve [K]
Default:274
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Minimum temperature value in saturation curve [K]
T_triple273.16Triple-point temperature [K]
Default:273.16
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Triple-point temperature [K]
cv_liquid1816Liquid isochoric specific heat capacity
Default:1816
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid isochoric specific heat capacity
cv_vapor1040Vapor isochoric specific heat capacity
Default:1040
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor isochoric specific heat capacity
e_c2.70298e+06Critical specific internal energy [J/kg]
Default:2.70298e+06
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Critical specific internal energy [J/kg]
emit_on_nannoneWhether to raise a warning, an exception (usually triggering a retry with a smaller time step) or an error (ending the simulation)
Default:none
C++ Type:MooseEnum
Options:none, warning, exception, error
Controllable:No
Description:Whether to raise a warning, an exception (usually triggering a retry with a smaller time step) or an error (ending the simulation)
fp_liquidLiquid single-phase fluid properties user object name
C++ Type:UserObjectName
Controllable:No
Description:Liquid single-phase fluid properties user object name
fp_vaporVapor single-phase fluid properties user object name
C++ Type:UserObjectName
Controllable:No
Description:Vapor single-phase fluid properties user object name
gamma_liquid2.35Liquid heat capacity ratio
Default:2.35
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid heat capacity ratio
gamma_vapor1.43Vapor heat capacity ratio
Default:1.43
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor heat capacity ratio
k_liquid0.5Liquid thermal conductivity
Default:0.5
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid thermal conductivity
k_vapor0.026Vapor thermal conductivity
Default:0.026
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor thermal conductivity
mu_liquid0.0002818Liquid dynamic viscosity
Default:0.0002818
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid dynamic viscosity
mu_vapor1.344e-05Vapor dynamic viscosity
Default:1.344e-05
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor dynamic viscosity
n_sat_samples374Number of samples to take in saturation curve
Default:374
C++ Type:unsigned int
Controllable:No
Description:Number of samples to take in saturation curve
p_c2.209e+07Critical pressure [Pa]
Default:2.209e+07
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Critical pressure [Pa]
p_inf_liquid1e+09Liquid stiffness pressure
Default:1e+09
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid stiffness pressure
p_inf_vapor0Vapor stiffness pressure
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor stiffness pressure
p_sat_guess611Initial guess for saturation pressure Newton solve [Pa]
Default:611
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Initial guess for saturation pressure Newton solve [Pa]
q_liquid-1.167e+06Liquid reference specific internal energy
Default:-1.167e+06
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid reference specific internal energy
q_prime_liquid0Liquid reference specific entropy
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Liquid reference specific entropy
q_prime_vapor-23000Vapor reference specific entropy
Default:-23000
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor reference specific entropy
q_vapor2.03e+06Vapor reference specific internal energy
Default:2.03e+06
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Vapor reference specific internal energy
rho_c322Critical density [kg/m^3]
Default:322
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Critical density [kg/m^3]
sigma_A0.2358'A' constant used in surface tension correlation [N/m]
Default:0.2358
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:'A' constant used in surface tension correlation [N/m]
sigma_B1.256'B' constant used in surface tension correlation
Default:1.256
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:'B' constant used in surface tension correlation
sigma_C0.625'C' constant used in surface tension correlation
Default:0.625
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:'C' constant used in surface tension correlation

Optional Parameters

allow_imperfect_jacobiansFalsetrue to allow unimplemented property derivative terms to be set to zero for the AD API
Default:False
C++ Type:bool
Controllable:No
Description:true to allow unimplemented property derivative terms to be set to zero for the AD API
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
fp_typetwo-phase-fpType of the fluid property object
Default:two-phase-fp
C++ Type:FPType
Controllable:No
Description:Type of the fluid property object

Advanced Parameters

prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

Input Files

(modules/fluid_properties/test/tests/materials/ad_saturation_temperature_material/ad_saturation_temperature_material.i)
(modules/fluid_properties/test/tests/materials/saturation_pressure_material/saturation_pressure_material.i)
(modules/thermal_hydraulics/test/tests/components/flow_channel_1phase/err.wrong_fp.i)
(modules/fluid_properties/test/tests/materials/ad_surface_tension_material/ad_surface_tension_material.i)

References

R. A. Berry, R. Saurel, and O. LeMetayer. The discrete equation method (DEM) for fully compressible, two-phase flows in ducts of spatially varying cross-section. Nuclear Engineering and Design, 240:3797–3818, 2010.

@article{berry2010,
    author = "Berry, R. A. and Saurel, R. and LeMetayer, O.",
    title = "The discrete equation method ({DEM}) for fully compressible, two-phase flows in ducts of spatially varying cross-section",
    journal = "Nuclear Engineering and Design",
    year = "2010",
    volume = "240",
    pages = "3797--3818"
}

The International Association for the Properties of Water and Steam. Revised release on surface tension of ordinary water substance. Technical Report, The International Association for the Properties of Water and Steam, 2014.

@techreport{iapws_surface_tension,
    author = "for the Properties of Water, The International Association and Steam",
    institution = "The International Association for the Properties of Water and Steam",
    title = "Revised Release on Surface Tension of Ordinary Water Substance",
    year = "2014"
}

(modules/fluid_properties/test/tests/materials/ad_saturation_temperature_material/ad_saturation_temperature_material.i)

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 1
  xmin = 0
  xmax = 1
[]

[FluidProperties]
  [./fp_2phase]
    type = StiffenedGasTwoPhaseFluidProperties
  [../]
[]

[Materials]
  [./p_mat]
    type = ADGenericConstantMaterial
    prop_names = 'p_test'
    prop_values = '1e5'
  [../]
  [./T_sat_mat]
    type = ADSaturationTemperatureMaterial
    p = p_test
    T_sat = T_sat_test
    fp_2phase = fp_2phase
  [../]
[]

[Postprocessors]
  [./T_sat_pp]
    type = ADElementIntegralMaterialProperty
    mat_prop = T_sat_test
    execute_on = 'INITIAL'
  [../]
[]

[Problem]
  solve = false
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
  execute_on = 'INITIAL'
[]

(modules/fluid_properties/test/tests/materials/saturation_pressure_material/saturation_pressure_material.i)

# This tests SaturationPressureMaterial, which computes a saturation pressure material
# property from a temperature material property and a TwoPhaseFluidProperties object.

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 1
  xmin = 0
  xmax = 1
[]

[FluidProperties]
  [fp_2phase]
    type = StiffenedGasTwoPhaseFluidProperties
  []
[]

[Materials]
  [T_mat]
    type = ADGenericConstantMaterial
    prop_names = 'T'
    prop_values = '400'
  []
  [p_sat_mat]
    type = ADSaturationPressureMaterial
    T = T
    p_sat = p_sat
    fp_2phase = fp_2phase
  []
[]

[Postprocessors]
  [p_sat_pp]
    type = ADElementAverageMaterialProperty
    mat_prop = p_sat
    execute_on = 'INITIAL'
  []
[]

[Problem]
  solve = false
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
  execute_on = 'INITIAL'
[]

(modules/thermal_hydraulics/test/tests/components/flow_channel_1phase/err.wrong_fp.i)

[GlobalParams]
  closures = simple_closures

  initial_p = 1e5
  initial_T = 300
  initial_vel = 0
[]

[FluidProperties]
  [fp_2phase]
    type = StiffenedGasTwoPhaseFluidProperties
  []
[]

[Closures]
  [simple_closures]
    type = Closures1PhaseSimple
  []
[]

[Components]
  [pipe]
    type = FlowChannel1Phase
    position = '0 0 0'
    orientation = '1 0 0'
    A = 1.
    D_h = 1.12837916709551
    f = 0.01
    length = 1
    n_elems = 100
    fp = fp_2phase   # this is wrong
  []

  [inlet]
    type = InletStagnationPressureTemperature1Phase
    input = 'pipe:in'
    p0 = 1e5
    T0 = 300
  []

  [outlet]
    type = Outlet1Phase
    input = 'pipe:out'
    p = 9.5e4
  []
[]

[Preconditioning]
  [pc]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  scheme = 'bdf2'

  dt = 1e-4
  dtmin = 1.e-7

  solve_type = 'PJFNK'
  nl_rel_tol = 1e-9
  nl_abs_tol = 1e-8
  nl_max_its = 10

  l_tol = 1e-8
  l_max_its = 100

  start_time = 0.0
  num_steps = 1
[]

(modules/fluid_properties/test/tests/materials/ad_surface_tension_material/ad_surface_tension_material.i)

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 1
  xmin = 0
  xmax = 1
[]

[FluidProperties]
  [./fp_2phase]
    type = StiffenedGasTwoPhaseFluidProperties
  [../]
[]

[Materials]
  [./T_mat]
    type = ADGenericConstantMaterial
    prop_names = 'T_test'
    prop_values = '300'
  [../]
  [./sigma_mat]
    type = ADSurfaceTensionMaterial
    T = T_test
    surface_tension = surface_tension_test
    fp_2phase = fp_2phase
  [../]
[]

[Postprocessors]
  [./surface_tension_test_pp]
    type = ADElementIntegralMaterialProperty
    mat_prop = surface_tension_test
    execute_on = 'INITIAL'
  [../]
[]

[Problem]
  solve = false
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
  execute_on = 'INITIAL'
[]

Input Parameters
Input Files
References