Specialized class for water and vapor mixture using pressure and enthalpy. More...

#include <PorousFlowWaterVapor.h>

Inheritance diagram for PorousFlowWaterVapor:

Public Member Functions
	PorousFlowWaterVapor (const InputParameters &parameters)

virtual std::string	fluidStateName () const override
	Name of FluidState. More...

void	thermophysicalProperties (Real pressure, Real enthalpy, unsigned int qp, FluidStatePhaseEnum &phase_state, std::vector< FluidStateProperties > &fsp) const override
	Determines the complete thermophysical state of the system for a given set of primary variables. More...

unsigned int	getPressureIndex () const

unsigned int	getEnthalpyIndex () const

void	initialize () final

void	execute () final

void	finalize () final

unsigned int	numPhases () const
	The maximum number of phases in this model. More...

unsigned int	numComponents () const
	The maximum number of components in this model. More...

unsigned int	aqueousPhaseIndex () const
	The index of the aqueous phase. More...

unsigned int	gasPhaseIndex () const
	The index of the gas phase. More...

unsigned int	aqueousComponentIndex () const
	The index of the aqueous fluid component. More...

unsigned int	gasComponentIndex () const
	The index of the gas fluid component. More...

unsigned int	saltComponentIndex () const
	The index of the salt component. More...

void	clearFluidStateProperties (std::vector< FluidStateProperties > &fsp) const
	Clears the contents of the FluidStateProperties data structure. More...

Protected Attributes
const SinglePhaseFluidProperties &	_water_fp
	Fluid properties UserObject for water. More...

const Real	_Mh2o
	Molar mass of water (kg/mol) More...

const Real	_p_triple
	Triple point pressure of water (Pa) More...

const Real	_p_critical
	Critical pressure of water (Pa) More...

const Real	_T_triple
	Triple point temperature of water (K) More...

const Real	_T_critical
	Critical temperature of water (K) More...

const unsigned int	_pidx
	Index of derivative wrt pressure. More...

const unsigned int	_hidx
	Index of derivative wrt enthalpy. More...

const Real	dT
	Perturbation applied to saturation temperature to move to gas/liquid phase. More...

unsigned int	_num_phases
	Number of phases. More...

unsigned int	_num_components
	Number of components. More...

const unsigned int	_aqueous_phase_number
	Phase number of the aqueous phase. More...

unsigned int	_gas_phase_number
	Phase number of the gas phase. More...

const unsigned int	_aqueous_fluid_component
	Fluid component number of the aqueous component. More...

unsigned int	_gas_fluid_component
	Fluid component number of the gas phase. More...

const unsigned int	_salt_component
	Salt component index. More...

const Real	_R
	Universal gas constant (J/mol/K) More...

const Real	_T_c2k
	Conversion from C to K. More...

const PorousFlowCapillaryPressure &	_pc
	Capillary pressure UserObject. More...

FluidStateProperties	_empty_fsp
	Empty FluidStateProperties object. More...

Detailed Description

Specialized class for water and vapor mixture using pressure and enthalpy.

Note: It is important to note that this class must be used in a non-isothernal setting (i.e., with both fluid and heat transport) in order for the problem to be well posed.

Definition at line 26 of file PorousFlowWaterVapor.h.

Constructor & Destructor Documentation

◆ PorousFlowWaterVapor()

PorousFlowWaterVapor::PorousFlowWaterVapor ( const InputParameters & parameters )

Definition at line 27 of file PorousFlowWaterVapor.C.

   : PorousFlowFluidStateSingleComponentBase(parameters),
     _water_fp(getUserObject<SinglePhaseFluidProperties>("water_fp")),
     _Mh2o(_water_fp.molarMass()),
     _p_triple(_water_fp.triplePointPressure()),
     _p_critical(_water_fp.criticalPressure()),
     _T_triple(_water_fp.triplePointTemperature()),
     _T_critical(_water_fp.criticalTemperature())
 {
   // Check that the correct FluidProperties UserObjects have been provided
   if (_water_fp.fluidName() != "water")
     paramError("water_fp", "A valid water FluidProperties UserObject must be provided in water_fp");
  
   // Set the number of phases and components, and their indexes
   _num_phases = 2;
   _num_components = 1;
   _gas_phase_number = 1 - _aqueous_phase_number;
   _gas_fluid_component = 1 - _aqueous_fluid_component;
  
   // Check that _aqueous_phase_number is <= total number of phases
   if (_aqueous_phase_number >= _num_phases)
     paramError("liquid_phase_number",
                "This value is larger than the possible number of phases ",
                _num_phases);
  
   // Check that _aqueous_fluid_component is <= total number of fluid components
   if (_aqueous_fluid_component >= _num_components)
     paramError("liquid_fluid_component",
                "This value is larger than the possible number of fluid components",
                _num_components);
  
   _empty_fsp = FluidStateProperties(_num_components);
 }

Member Function Documentation

◆ aqueousComponentIndex()

unsigned int PorousFlowFluidStateBase::aqueousComponentIndex ( ) const

inlineinherited

The index of the aqueous fluid component.

Returns: aqueous fluid component number

Definition at line 92 of file PorousFlowFluidStateBase.h.

92 { return _aqueous_fluid_component; };

◆ aqueousPhaseIndex()

unsigned int PorousFlowFluidStateBase::aqueousPhaseIndex ( ) const

inlineinherited

The index of the aqueous phase.

Returns: aqueous phase number

Definition at line 80 of file PorousFlowFluidStateBase.h.

80 { return _aqueous_phase_number; };

◆ clearFluidStateProperties()

void PorousFlowFluidStateBase::clearFluidStateProperties ( std::vector< FluidStateProperties > & fsp ) const

inherited

Clears the contents of the FluidStateProperties data structure.

Parameters

[out] fsp FluidStateProperties data structure with all data initialized to 0

Definition at line 39 of file PorousFlowFluidStateBase.C.

 {
   std::fill(fsp.begin(), fsp.end(), _empty_fsp);
 }

Referenced by PorousFlowWaterNCG::thermophysicalProperties(), PorousFlowFluidStateSingleComponent::thermophysicalProperties(), and PorousFlowBrineCO2::thermophysicalProperties().

◆ execute()

void PorousFlowFluidStateBase::execute ( )

inlinefinalinherited

Definition at line 61 of file PorousFlowFluidStateBase.h.

61 {};

◆ finalize()

void PorousFlowFluidStateBase::finalize ( )

inlinefinalinherited

Definition at line 62 of file PorousFlowFluidStateBase.h.

62 {};

◆ fluidStateName()

std::string PorousFlowWaterVapor::fluidStateName ( ) const

overridevirtual

Name of FluidState.

Implements PorousFlowFluidStateBase.

Definition at line 62 of file PorousFlowWaterVapor.C.

 {
   return "water-vapor";
 }

◆ gasComponentIndex()

unsigned int PorousFlowFluidStateBase::gasComponentIndex ( ) const

inlineinherited

The index of the gas fluid component.

Returns: gas fluid component number

Definition at line 98 of file PorousFlowFluidStateBase.h.

98 { return _gas_fluid_component; };

◆ gasPhaseIndex()

unsigned int PorousFlowFluidStateBase::gasPhaseIndex ( ) const

inlineinherited

The index of the gas phase.

Returns: gas phase number

Definition at line 86 of file PorousFlowFluidStateBase.h.

86 { return _gas_phase_number; };

◆ getEnthalpyIndex()

unsigned int PorousFlowFluidStateSingleComponentBase::getEnthalpyIndex ( ) const

inlineinherited

Definition at line 45 of file PorousFlowFluidStateSingleComponentBase.h.

45 { return _hidx; };

◆ getPressureIndex()

unsigned int PorousFlowFluidStateSingleComponentBase::getPressureIndex ( ) const

inlineinherited

Definition at line 44 of file PorousFlowFluidStateSingleComponentBase.h.

44 { return _pidx; };

◆ initialize()

void PorousFlowFluidStateBase::initialize ( )

inlinefinalinherited

Definition at line 60 of file PorousFlowFluidStateBase.h.

60 {};

◆ numComponents()

unsigned int PorousFlowFluidStateBase::numComponents ( ) const

inlineinherited

The maximum number of components in this model.

Returns: number of components

Definition at line 74 of file PorousFlowFluidStateBase.h.

74 { return _num_components; };

◆ numPhases()

unsigned int PorousFlowFluidStateBase::numPhases ( ) const

inlineinherited

The maximum number of phases in this model.

Returns: number of phases

Definition at line 68 of file PorousFlowFluidStateBase.h.

68 { return _num_phases; };

Referenced by PorousFlowFluidState::PorousFlowFluidState(), and PorousFlowFluidStateSingleComponent::PorousFlowFluidStateSingleComponent().

◆ saltComponentIndex()

unsigned int PorousFlowFluidStateBase::saltComponentIndex ( ) const

inlineinherited

The index of the salt component.

Returns: salt component number

Definition at line 104 of file PorousFlowFluidStateBase.h.

104 { return _salt_component; };

◆ thermophysicalProperties()

void PorousFlowWaterVapor::thermophysicalProperties	(	Real	pressure,
		Real	enthalpy,
		unsigned int	qp,
		FluidStatePhaseEnum &	phase_state,
		std::vector< FluidStateProperties > &	fsp
	)		const

overridevirtual

Determines the complete thermophysical state of the system for a given set of primary variables.

Parameters

	pressure	gas phase pressure (Pa)
	enthalpy	fluid enthalpy (J/kg)
	qp	quadpoint index
[out]	phase_state	the current phase (gas, liquid or two-phase)
[out]	fsp	the FluidStateProperties struct containing all properties

Implements PorousFlowFluidStateSingleComponentBase.

Definition at line 68 of file PorousFlowWaterVapor.C.

 {
   FluidStateProperties & liquid = fsp[_aqueous_phase_number];
   FluidStateProperties & gas = fsp[_gas_phase_number];
  
   // AD versions of primary variables
   DualReal p = pressure;
   Moose::derivInsert(p.derivatives(), _pidx, 1.0);
   DualReal h = enthalpy;
   Moose::derivInsert(h.derivatives(), _hidx, 1.0);
  
   DualReal Tsat = 0.0;
   DualReal hl = 0.0;
   DualReal hv = 0.0;
  
   // Determine the phase state of the system
   if (p.value() >= _p_triple && p.value() <= _p_critical)
   {
     // Saturation temperature at the given pressure
     Tsat = _water_fp.vaporTemperature(p);
  
     // Enthalpy of saturated liquid and saturated vapor
     hl = _water_fp.h_from_p_T(p, Tsat - dT);
     hv = _water_fp.h_from_p_T(p, Tsat + dT);
  
     if (h.value() < hl)
       phase_state = FluidStatePhaseEnum::LIQUID;
  
     else if (h.value() >= hl && h.value() <= hv)
       phase_state = FluidStatePhaseEnum::TWOPHASE;
  
     else // h > hv
       phase_state = FluidStatePhaseEnum::GAS;
   }
   else // p.value() > _p_critical
   {
     // Check whether the phase point is in the liquid or vapor state
     const DualReal T = _water_fp.T_from_p_h(p, h);
  
     if (T.value() <= _T_critical)
       phase_state = FluidStatePhaseEnum::LIQUID;
     else
     {
       // The supercritical state is treated as a gas
       phase_state = FluidStatePhaseEnum::GAS;
     }
   }
  
   // Calculate the properties for each phase as required
   switch (phase_state)
   {
     case FluidStatePhaseEnum::GAS:
     {
       gas.pressure = p + _pc.capillaryPressure(0.0, qp);
       gas.saturation = 1.0;
  
       const DualReal T = _water_fp.T_from_p_h(gas.pressure, h);
  
       gas.temperature = T;
       liquid.temperature = T;
  
       gas.density = _water_fp.rho_from_p_T(gas.pressure, T);
       gas.viscosity = _water_fp.mu_from_p_T(gas.pressure, T);
       gas.enthalpy = h;
       gas.internal_energy = _water_fp.e_from_p_T(gas.pressure, T);
  
       break;
     }
  
     case FluidStatePhaseEnum::LIQUID:
     {
       const DualReal T = _water_fp.T_from_p_h(p, h);
  
       liquid.pressure = p;
       liquid.temperature = T;
       liquid.density = _water_fp.rho_from_p_T(p, T);
       liquid.viscosity = _water_fp.mu_from_p_T(p, T);
       liquid.enthalpy = h;
       liquid.internal_energy = _water_fp.e_from_p_T(p, T);
       liquid.saturation = 1.0;
  
       break;
     }
  
     case FluidStatePhaseEnum::TWOPHASE:
     {
       // Latent heat of vaporization
       const DualReal hvl = hv - hl;
  
       // Vapor quality
       const DualReal X = (h - hl) / hvl;
  
       // Perturbed saturation temperature to ensure that the correct
       // phase properties are calculated
       const DualReal Tsatl = Tsat - dT;
       const DualReal Tsatv = Tsat + dT;
  
       // Density
       const DualReal rhol = _water_fp.rho_from_p_T(p, Tsatl);
       const DualReal rhov = _water_fp.rho_from_p_T(p, Tsatv);
  
       // Vapor (gas) saturation
       const DualReal satv = X * rhol / (rhov + X * (rhol - rhov));
  
       gas.temperature = Tsat;
       gas.density = rhov;
       gas.viscosity = _water_fp.mu_from_p_T(p, Tsatv);
       gas.enthalpy = hv;
       gas.internal_energy = _water_fp.e_from_p_T(p, Tsatv);
       gas.saturation = satv;
  
       liquid.temperature = Tsat;
       liquid.density = rhol;
       liquid.viscosity = _water_fp.mu_from_p_T(p, Tsatl);
       liquid.enthalpy = hl;
       liquid.internal_energy = _water_fp.e_from_p_T(p, Tsatl);
       liquid.saturation = 1.0 - satv;
  
       liquid.pressure = p;
       gas.pressure = p + _pc.capillaryPressure(liquid.saturation, qp);
  
       break;
     }
   }
 }