StiffenedGasTwoPhaseFluidProperties

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

The saturation curve 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): where

  • is the temperature [K],

  • is the critical temperature [K], and

  • [N/m], [-], and [-] 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

References

  1. 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.[BibTeX]
  2. 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.[BibTeX]