TwoPhaseNCGPartialPressureFluidProperties

This fluid properties class is used for two-phase fluids that have a single, non-condensable gas (NCG) mixed with their vapor phase. The vapor mixture model of IdealRealGasMixtureFluidProperties is used, which corresponds to a partial pressure mixture model for real gases.

In addition to the interfaces of its parent class, this class provides an interface for computing the mass fraction of the NCG $var element = document.getElementById("moose-equation-91f1266e-22f5-46a9-a9f2-1f70caef4989");katex.render("\\xi_\\text{NCG}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ assuming saturation of the main fluid at the temperature $var element = document.getElementById("moose-equation-bbab91ab-c9ef-4acd-9319-4efee568ad1a");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}"}});$ and mixture pressure $var element = document.getElementById("moose-equation-46d723e6-c9f0-4305-9322-5dfdced9155b");katex.render("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}"}});$ . The molar fraction of the NCG $var element = document.getElementById("moose-equation-5c2d3389-8541-4fed-a0b3-5addbd67526b");katex.render("\\psi_\\text{NCG}", 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 computed by combining the two relations:

var element = document.getElementById("moose-equation-84a4e5b2-89cd-43f3-ac46-65ae93e6cdde");katex.render("p = p_{v,\\text{sat}}(T) + p_\\text{NCG} \\,,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

var element = document.getElementById("moose-equation-50e1bbc3-1aee-4dce-9e86-ac63afaff87c");katex.render("p_\\text{NCG} = \\psi_\\text{NCG} p \\,,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

which with the relation between molar fractions and mass fractions, yields

var element = document.getElementById("moose-equation-297be0b0-8e13-4cfe-9f7e-461922788bd9");katex.render("\\psi_\\text{NCG} = \\frac{p - p_{v,\\text{sat}}(T)}{p} \\,,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

var element = document.getElementById("moose-equation-f814820e-bf0f-4129-aa64-d2d39fee3a8b");katex.render("\\xi_\\text{NCG} = \\frac{M_\\text{NCG}}{M} \\,,", 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-ec86c11c-623c-43c3-862f-4e9de4bf28cc");katex.render("M_\\text{NCG}", 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 NCG molar mass, and $var element = document.getElementById("moose-equation-b5dbc284-9292-4aa6-9f7b-056a12729248");katex.render("M", 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 mixture molar mass,

var element = document.getElementById("moose-equation-a5be2777-2cbf-406a-bb5f-e6ef814fd78e");katex.render("M = (1 - \\psi_\\text{NCG}) M_v + \\psi_\\text{NCG} M_\\text{NCG} \\,,", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});

with $var element = document.getElementById("moose-equation-5ecf19ac-3688-4253-a54e-f85266354c3f");katex.render("M_v", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ denoting the main fluid molar mass.

Input Parameters

fp_ncgNCG fluid properties object
C++ Type:UserObjectName
Controllable:No
Description:NCG fluid properties object

Required Parameters

fp_2phaseName of fluid properties user object(s) for two-phase model
C++ Type:UserObjectName
Controllable:No
Description:Name of fluid properties user object(s) for two-phase model
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

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-ncg-fpType of the fluid property object
Default:two-phase-ncg-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/fp_interrogator/2ph_ncg_partial_pressure_p_T.i)
(modules/fluid_properties/test/tests/two_phase_ncg_partial_pressure/test.i)

References

No citations exist within this document.

(modules/fluid_properties/test/tests/fp_interrogator/2ph_ncg_partial_pressure_p_T.i)

[FluidPropertiesInterrogator]
  fp = fp_2phase_ncg_partial_pressure
  p = 1e5
  T = 372.7559289
[]

[FluidProperties]
  [fp_water]
    type = StiffenedGasFluidProperties
    gamma = 2.35
    cv = 1816.0
    q = -1.167e6
    p_inf = 1.0e9
    q_prime = 0
    k = 0.5
    mu = 281.8e-6
    M = 0.01801488
  []
  [fp_steam]
    type = IdealGasFluidProperties
    gamma = 1.43
    molar_mass = 0.01801488
  []
  [fp_2phase]
    type = TestTwoPhaseFluidProperties
    fp_liquid = fp_water
    fp_vapor = fp_steam
  []

  [fp_air]
    type = IdealGasFluidProperties
  []

  [fp_2phase_ncg_partial_pressure]
    type = TwoPhaseNCGPartialPressureFluidProperties
    fp_2phase = fp_2phase
    fp_ncg = fp_air
  []
[]

(modules/fluid_properties/test/tests/two_phase_ncg_partial_pressure/test.i)

p = 40e3
T = 300

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 1
[]

[FluidProperties]
  [fp_water]
    type = StiffenedGasFluidProperties
    gamma = 2.35
    cv = 1816.0
    q = -1.167e6
    p_inf = 1.0e9
    q_prime = 0
    k = 0.5
    mu = 281.8e-6
    M = 0.01801488
  []
  [fp_steam]
    type = IdealGasFluidProperties
    gamma = 1.43
    molar_mass = 0.01801488
  []
  [fp_air]
    type = IdealGasFluidProperties
  []
  [fp_2phase]
    type = TestTwoPhaseFluidProperties
    fp_liquid = fp_water
    fp_vapor = fp_steam
  []
  [fp_2phase_ncg]
    type = TwoPhaseNCGPartialPressureFluidProperties
    fp_2phase = fp_2phase
    fp_ncg = fp_air
  []
[]

[Functions]
  [p_sat_fn]
    type = TwoPhaseNCGPartialPressureFunction
    fluid_properties = fp_2phase_ncg
    property_call = p_sat
    arg1 = ${T}
  []
  [x_sat_ncg_fn]
    type = TwoPhaseNCGPartialPressureFunction
    fluid_properties = fp_2phase_ncg
    property_call = x_sat_ncg_from_p_T
    arg1 = ${p}
    arg2 = ${T}
  []
[]

[Postprocessors]
  [p_sat]
    type = FunctionValuePostprocessor
    function = p_sat_fn
    execute_on = 'INITIAL'
  []
  [x_sat_ncg]
    type = FunctionValuePostprocessor
    function = x_sat_ncg_fn
    execute_on = 'INITIAL'
  []
[]

[Problem]
  solve = false
[]


[Executioner]
  type = Steady
[]

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

Input Parameters
Input Files
References