StiffenedGasTwoPhaseFluidProperties

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

The saturation curve $T_\text{sat}(p)$ 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): $\sigma(T) = A \left(1 - \frac{T}{T_C}\right)^B \left[1 - C \left(1 - \frac{T}{T_C}\right)\right] \,,$ where

$T$ is the temperature [K],
$T_C$ is the critical temperature [K], and
$A$ [N/m], $B$ [-], and $C$ [-] are a fluid-dependent constants.

Input Parameters

L_fusion0.334Latent heat of fusion [J/kg]
Default:0.334
C++ Type:double
Controllable:No
Description:Latent heat of fusion [J/kg]
M_liquid0.0180149Liquid molar mass
Default:0.0180149
C++ Type:double
Controllable:No
Description:Liquid molar mass
M_vapor0.0180149Vapor molar mass
Default:0.0180149
C++ Type:double
Controllable:No
Description:Vapor molar mass
T_c647.096Critical temperature [K]
Default:647.096
C++ Type:double
Controllable:No
Description:Critical temperature [K]
T_sat_max647Maximum temperature value in saturation curve [K]
Default:647
C++ Type:double
Controllable:No
Description:Maximum temperature value in saturation curve [K]
T_sat_min274Minimum temperature value in saturation curve [K]
Default:274
C++ Type:double
Controllable:No
Description:Minimum temperature value in saturation curve [K]
T_triple273.16Triple-point temperature [K]
Default:273.16
C++ Type:double
Controllable:No
Description:Triple-point temperature [K]
cv_liquid1816Liquid isochoric specific heat capacity
Default:1816
C++ Type:double
Controllable:No
Description:Liquid isochoric specific heat capacity
cv_vapor1040Vapor isochoric specific heat capacity
Default:1040
C++ Type:double
Controllable:No
Description:Vapor isochoric specific heat capacity
e_c2.70298e+06Critical specific internal energy [J/kg]
Default:2.70298e+06
C++ Type:double
Controllable:No
Description:Critical specific internal energy [J/kg]
emit_on_nannoneRaise mooseWarning or mooseError?
Default:none
C++ Type:MooseEnum
Options:none, warning, error
Controllable:No
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, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS
Controllable:No
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, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS.
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
Controllable:No
Description:Liquid heat capacity ratio
gamma_vapor1.43Vapor heat capacity ratio
Default:1.43
C++ Type:double
Controllable:No
Description:Vapor heat capacity ratio
k_liquid0.5Liquid thermal conductivity
Default:0.5
C++ Type:double
Controllable:No
Description:Liquid thermal conductivity
k_vapor0.026Vapor thermal conductivity
Default:0.026
C++ Type:double
Controllable:No
Description:Vapor thermal conductivity
mu_liquid0.0002818Liquid dynamic viscosity
Default:0.0002818
C++ Type:double
Controllable:No
Description:Liquid dynamic viscosity
mu_vapor1.344e-05Vapor dynamic viscosity
Default:1.344e-05
C++ Type:double
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
Controllable:No
Description:Critical pressure [Pa]
p_inf_liquid1e+09Liquid stiffness pressure
Default:1e+09
C++ Type:double
Controllable:No
Description:Liquid stiffness pressure
p_inf_vapor0Vapor stiffness pressure
Default:0
C++ Type:double
Controllable:No
Description:Vapor stiffness pressure
p_sat_guess611Initial guess for saturation pressure Newton solve [Pa]
Default:611
C++ Type:double
Controllable:No
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
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.
q_liquid-1.167e+06Liquid reference specific internal energy
Default:-1.167e+06
C++ Type:double
Controllable:No
Description:Liquid reference specific internal energy
q_prime_liquid0Liquid reference specific entropy
Default:0
C++ Type:double
Controllable:No
Description:Liquid reference specific entropy
q_prime_vapor-23000Vapor reference specific entropy
Default:-23000
C++ Type:double
Controllable:No
Description:Vapor reference specific entropy
q_vapor2.03e+06Vapor reference specific internal energy
Default:2.03e+06
C++ Type:double
Controllable:No
Description:Vapor reference specific internal energy
rho_c322Critical density [kg/m^3]
Default:322
C++ Type:double
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
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
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
Controllable:No
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
Controllable:No
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).
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.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
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
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
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
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
Controllable:No
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

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"
}Copy

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"
}Copy

Input Parameters
References