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-a7e2d15c-1b3d-4a09-8b6c-b5243ecd829b");katex.render("T_\\text{sat}(p)", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ 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-6938b051-59be-4467-a27d-71593eb49819");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:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ where

$var element = document.getElementById("moose-equation-e9e54ffe-f9ab-4931-a76e-b413a77905ee");katex.render("T", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is the temperature [K],
$var element = document.getElementById("moose-equation-0df69b45-6b66-4cfc-a312-327c655990ae");katex.render("T_C", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is the critical temperature [K], and
$var element = document.getElementById("moose-equation-57be3717-4bc2-4da6-bf28-79fa8a05079c");katex.render("A", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ [N/m], $var element = document.getElementById("moose-equation-6b218451-2bb9-4fc6-8346-484e500c756a");katex.render("B", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ [-], and $var element = document.getElementById("moose-equation-8f942818-d6b4-4a9a-a9f0-b6a8e4f93a5a");katex.render("C", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ [-] are a fluid-dependent constants.

Input Parameters

L_fusion0.334Latent heat of fusion [J/kg]
Default:0.334
C++ Type:double
Options:
Description:Latent heat of fusion [J/kg]
M_liquid0.0180149Liquid molar mass
Default:0.0180149
C++ Type:double
Options:
Description:Liquid molar mass
M_vapor0.0180149Vapor molar mass
Default:0.0180149
C++ Type:double
Options:
Description:Vapor molar mass
T_c647.096Critical temperature [K]
Default:647.096
C++ Type:double
Options:
Description:Critical temperature [K]
T_sat_max647Maximum temperature value in saturation curve [K]
Default:647
C++ Type:double
Options:
Description:Maximum temperature value in saturation curve [K]
T_sat_min274Minimum temperature value in saturation curve [K]
Default:274
C++ Type:double
Options:
Description:Minimum temperature value in saturation curve [K]
T_triple273.16Triple-point temperature [K]
Default:273.16
C++ Type:double
Options:
Description:Triple-point temperature [K]
allow_imperfect_jacobiansFalsetrue to allow unimplemented property derivative terms to be set to zero for the AD API
Default:False
C++ Type:bool
Options:
Description:true to allow unimplemented property derivative terms to be set to zero for the AD API
cv_liquid1816Liquid isochoric specific heat capacity
Default:1816
C++ Type:double
Options:
Description:Liquid isochoric specific heat capacity
cv_vapor1040Vapor isochoric specific heat capacity
Default:1040
C++ Type:double
Options:
Description:Vapor isochoric specific heat capacity
e_c2.70298e+06Critical specific internal energy [J/kg]
Default:2.70298e+06
C++ Type:double
Options:
Description:Critical specific internal energy [J/kg]
emit_on_nannoneRaise mooseWarning or mooseError?
Default:none
C++ Type:MooseEnum
Options:none, warning, error
Description:Raise mooseWarning or mooseError?
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
fp_liquidLiquid single-phase fluid properties user object name
C++ Type:UserObjectName
Options:
Description:Liquid single-phase fluid properties user object name
fp_typetwo-phase-fpType of the fluid property object
Default:two-phase-fp
C++ Type:FPType
Options:
Description:Type of the fluid property object
fp_vaporVapor single-phase fluid properties user object name
C++ Type:UserObjectName
Options:
Description:Vapor single-phase fluid properties user object name
gamma_liquid2.35Liquid heat capacity ratio
Default:2.35
C++ Type:double
Options:
Description:Liquid heat capacity ratio
gamma_vapor1.43Vapor heat capacity ratio
Default:1.43
C++ Type:double
Options:
Description:Vapor heat capacity ratio
k_liquid0.5Liquid thermal conductivity
Default:0.5
C++ Type:double
Options:
Description:Liquid thermal conductivity
k_vapor0.026Vapor thermal conductivity
Default:0.026
C++ Type:double
Options:
Description:Vapor thermal conductivity
mu_liquid0.0002818Liquid dynamic viscosity
Default:0.0002818
C++ Type:double
Options:
Description:Liquid dynamic viscosity
mu_vapor1.344e-05Vapor dynamic viscosity
Default:1.344e-05
C++ Type:double
Options:
Description:Vapor dynamic viscosity
n_sat_samples374Number of samples to take in saturation curve
Default:374
C++ Type:unsigned int
Options:
Description:Number of samples to take in saturation curve
p_c2.209e+07Critical pressure [Pa]
Default:2.209e+07
C++ Type:double
Options:
Description:Critical pressure [Pa]
p_inf_liquid1e+09Liquid stiffness pressure
Default:1e+09
C++ Type:double
Options:
Description:Liquid stiffness pressure
p_inf_vapor0Vapor stiffness pressure
Default:0
C++ Type:double
Options:
Description:Vapor stiffness pressure
p_sat_guess611Initial guess for saturation pressure Newton solve [Pa]
Default:611
C++ Type:double
Options:
Description:Initial guess for saturation pressure Newton solve [Pa]
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
Options:
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.
q_liquid-1.167e+06Liquid reference specific internal energy
Default:-1.167e+06
C++ Type:double
Options:
Description:Liquid reference specific internal energy
q_prime_liquid0Liquid reference specific entropy
Default:0
C++ Type:double
Options:
Description:Liquid reference specific entropy
q_prime_vapor-23000Vapor reference specific entropy
Default:-23000
C++ Type:double
Options:
Description:Vapor reference specific entropy
q_vapor2.03e+06Vapor reference specific internal energy
Default:2.03e+06
C++ Type:double
Options:
Description:Vapor reference specific internal energy
rho_c322Critical density [kg/m^3]
Default:322
C++ Type:double
Options:
Description:Critical density [kg/m^3]
sigma_A0.2358'A' constant used in surface tension correlation [N/m]
Default:0.2358
C++ Type:double
Options:
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
Options:
Description:'B' constant used in surface tension correlation
sigma_C0.625'C' constant used in surface tension correlation
Default:0.625
C++ Type:double
Options:
Description:'C' constant used in surface tension correlation

Optional Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Options:
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Options:
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Options:
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Options:
Description:Forces the UserObject to be executed in PREIC during initial setup
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Options:
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced 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/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
[]

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

[Modules]
  [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/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
[]

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

Install MOOSE

New Users

Examples and Tutorials

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools

INL Applications and Remote Access