Specialized class for water and a non-condensable gas (NCG) Includes dissolution of gas in liquid water phase using Henry's law. More...

#include <PorousFlowWaterNCG.h>

Inheritance diagram for PorousFlowWaterNCG:

[legend]

Public Member Functions
	PorousFlowWaterNCG (const InputParameters &parameters)

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

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

void	equilibriumMassFractions (const DualReal &pressure, const DualReal &temperature, DualReal &Xncg, DualReal &Yh2o) const
	Mass fractions of NCG in liquid phase and H2O in gas phase at thermodynamic equilibrium. More...

void	massFractions (const DualReal &pressure, const DualReal &temperature, const DualReal &Z, FluidStatePhaseEnum &phase_state, std::vector< FluidStateProperties > &fsp) const
	Mass fractions of NCG and H2O in both phases, as well as derivatives wrt PorousFlow variables. More...

void	gasProperties (const DualReal &pressure, const DualReal &temperature, std::vector< FluidStateProperties > &fsp) const
	Gas properties - density, viscosity and enthalpy. More...

void	liquidProperties (const DualReal &pressure, const DualReal &temperature, std::vector< FluidStateProperties > &fsp) const
	Liquid properties - density, viscosity and enthalpy Note: The pressure here is the liquid pressure. More...

DualReal	liquidDensity (const DualReal &pressure, const DualReal &temperature) const
	Density of the liquid phase Note: The pressure here is the gas pressure. More...

DualReal	gasDensity (const DualReal &pressure, const DualReal &temperature, std::vector< FluidStateProperties > &fsp) const
	Density of the gas phase. More...

DualReal	saturation (const DualReal &pressure, const DualReal &temperature, const DualReal &Z, std::vector< FluidStateProperties > &fsp) const
	Gas saturation in the two-phase region. More...

void	twoPhaseProperties (const DualReal &pressure, const DualReal &temperature, const DualReal &Z, unsigned int qp, std::vector< FluidStateProperties > &fsp) const
	Gas and liquid properties in the two-phase region. More...

DualReal	enthalpyOfDissolution (const DualReal &temperature) const
	Enthalpy of dissolution of NCG in water calculated using Henry's constant From Himmelblau, Partial molal heats and entropies of solution for gases dissolved in water from the freezing to the near critical point, J. More...

virtual Real	totalMassFraction (Real pressure, Real temperature, Real Xnacl, Real saturation, unsigned int qp) const override
	Total mass fraction of fluid component summed over all phases in the two-phase state for a specified gas saturation. More...

void	phaseState (Real Zi, Real Xi, Real Yi, FluidStatePhaseEnum &phase_state) const
	Determines the phase state gven the total mass fraction and equilibrium mass fractions. More...

unsigned int	getPressureIndex () const

unsigned int	getTemperatureIndex () const

unsigned int	getZIndex () const

unsigned int	getXIndex () const

Real	rachfordRice (Real vf, std::vector< Real > &Zi, std::vector< Real > &Ki) const
	Rachford-Rice equation for vapor fraction. More...

Real	rachfordRiceDeriv (Real vf, std::vector< Real > &Zi, std::vector< Real > &Ki) const
	Derivative of Rachford-Rice equation wrt vapor fraction. More...

Real	vaporMassFraction (Real Z0, Real K0, Real K1) const
	Solves Rachford-Rice equation to provide vapor mass fraction. More...

DualReal	vaporMassFraction (const DualReal &Z0, const DualReal &K0, const DualReal &K1) const

Real	vaporMassFraction (std::vector< Real > &Zi, std::vector< Real > &Ki) 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 Member Functions
DualReal	moleFractionToMassFraction (const DualReal &xmol) const
	Convert mole fraction to mass fraction. More...

void	checkVariables (Real temperature) const
	Check that the temperature is between the triple and critical values. More...

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

const Water97FluidProperties &	_water97_fp
	Fluid properties UserObject for water (used to access Henry's law) More...

const SinglePhaseFluidProperties &	_ncg_fp
	Fluid properties UserObject for the NCG. More...

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

const Real	_Mncg
	Molar mass of non-condensable gas (kg/mol) More...

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

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

const std::vector< Real >	_ncg_henry
	Henry's coefficients for the NCG. More...

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

const unsigned int	_Zidx
	Index of derivative wrt total mass fraction Z. More...

const unsigned int	_Tidx
	Index of derivative wrt temperature. More...

const unsigned int	_Xidx
	Index of derivative wrt salt mass fraction X. More...

const Real	_nr_max_its
	Maximum number of iterations for the Newton-Raphson routine. More...

const Real	_nr_tol
	Tolerance for Newton-Raphson iterations. 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 a non-condensable gas (NCG) Includes dissolution of gas in liquid water phase using Henry's law.

Notation convention Throughout this class, both mole fractions and mass fractions will be used. The following notation will be used: yk: mole fraction of component k in the gas phase xk: mole fraction of component k in the liquid phase Yk: mass fraction of component k in the gas phase Xk: mass fraction of component k in the liquid phase

Definition at line 33 of file PorousFlowWaterNCG.h.

Constructor & Destructor Documentation

◆ PorousFlowWaterNCG()

PorousFlowWaterNCG::PorousFlowWaterNCG ( const InputParameters & parameters )

Definition at line 29 of file PorousFlowWaterNCG.C.

   : PorousFlowFluidStateMultiComponentBase(parameters),
     _water_fp(getUserObject<SinglePhaseFluidProperties>("water_fp")),
     _water97_fp(getUserObject<Water97FluidProperties>("water_fp")),
     _ncg_fp(getUserObject<SinglePhaseFluidProperties>("gas_fp")),
     _Mh2o(_water_fp.molarMass()),
     _Mncg(_ncg_fp.molarMass()),
     _water_triple_temperature(_water_fp.triplePointTemperature()),
     _water_critical_temperature(_water_fp.criticalTemperature()),
     _ncg_henry(_ncg_fp.henryCoefficients())
 {
   // 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 = 2;
   _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; };

◆ checkVariables()

void PorousFlowWaterNCG::checkVariables ( Real temperature ) const

protected

Check that the temperature is between the triple and critical values.

Parameters

temperature fluid temperature (K)

Definition at line 384 of file PorousFlowWaterNCG.C.

 {
   // Check whether the input temperature is within the region of validity of this equation
   // of state (T_triple <= T <= T_critical)
   if (temperature < _water_triple_temperature || temperature > _water_critical_temperature)
     mooseException(name() + ": temperature " + Moose::stringify(temperature) +
                    " is outside range 273.16 K <= T <= 647.096 K");
 }

Referenced by thermophysicalProperties(), and totalMassFraction().

◆ 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 thermophysicalProperties(), PorousFlowFluidStateSingleComponent::thermophysicalProperties(), and PorousFlowBrineCO2::thermophysicalProperties().

◆ enthalpyOfDissolution()

DualReal PorousFlowWaterNCG::enthalpyOfDissolution ( const DualReal & temperature ) const

Enthalpy of dissolution of NCG in water calculated using Henry's constant From Himmelblau, Partial molal heats and entropies of solution for gases dissolved in water from the freezing to the near critical point, J.

Phys. Chem. 63 (1959)

Parameters

temperature fluid temperature (K)

Returns: enthalpy of dissolution (J/kg)

Definition at line 394 of file PorousFlowWaterNCG.C.

 {
   // Henry's constant
   const DualReal Kh = _water97_fp.henryConstant(temperature, _ncg_henry);
  
   DualReal hdis = -_R * temperature * temperature * Kh.derivatives()[_Tidx] / Kh / _Mncg;
  
   // Derivative of enthalpy of dissolution wrt temperature requires the second derivative of
   // Henry's constant wrt temperature. For simplicity, approximate this numerically
   const Real dT = temperature.value() * 1.0e-8;
   const DualReal t2 = temperature + dT;
   const DualReal Kh2 = _water97_fp.henryConstant(t2, _ncg_henry);
  
   const Real dhdis_dT =
       (-_R * t2 * t2 * Kh2.derivatives()[_Tidx] / Kh2 / _Mncg - hdis).value() / dT;
  
   hdis.derivatives() = temperature.derivatives() * dhdis_dT;
  
   return hdis;
 }

Referenced by liquidProperties().

◆ equilibriumMassFractions()

void PorousFlowWaterNCG::equilibriumMassFractions	(	const DualReal &	pressure,
		const DualReal &	temperature,
		DualReal &	Xncg,
		DualReal &	Yh2o
	)		const

Mass fractions of NCG in liquid phase and H2O in gas phase at thermodynamic equilibrium.

Calculated using Henry's law (for NCG component), and Raoult's law (for water).

Parameters

	pressure	phase pressure (Pa)
	temperature	phase temperature (K)
[out]	Xncg	mass fraction of NCG in liquid (kg/kg)
[out]	Yh2o	mass fraction of H2O in gas (kg/kg)

Definition at line 354 of file PorousFlowWaterNCG.C.

 {
   // Equilibrium constants for each component (Henry's law for the NCG
   // component, and Raoult's law for water).
   const DualReal Kh = _water97_fp.henryConstant(temperature, _ncg_henry);
   const DualReal psat = _water_fp.vaporPressure(temperature);
  
   const DualReal Kncg = Kh / pressure;
   const DualReal Kh2o = psat / pressure;
  
   // The mole fractions for the NCG component in the two component
   // case can be expressed in terms of the equilibrium constants only
   const DualReal xncg = (1.0 - Kh2o) / (Kncg - Kh2o);
   const DualReal yncg = Kncg * xncg;
  
   // Convert mole fractions to mass fractions
   Xncg = moleFractionToMassFraction(xncg);
   Yh2o = 1.0 - moleFractionToMassFraction(yncg);
 }

Referenced by massFractions(), and totalMassFraction().

◆ 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 PorousFlowWaterNCG::fluidStateName ( ) const

overridevirtual

Name of FluidState.

Implements PorousFlowFluidStateBase.

Definition at line 66 of file PorousFlowWaterNCG.C.

 {
   return "water-ncg";
 }

◆ 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; };

◆ gasDensity()

DualReal PorousFlowWaterNCG::gasDensity	(	const DualReal &	pressure,
		const DualReal &	temperature,
		std::vector< FluidStateProperties > &	fsp
	)		const

Density of the gas phase.

Parameters

pressure	pressure (Pa)
temperature	temperature (K)

Returns: gas density (kg/m^3)

Definition at line 280 of file PorousFlowWaterNCG.C.

 {
   auto & liquid = fsp[_aqueous_phase_number];
   auto & gas = fsp[_gas_phase_number];
  
   DualReal psat = _water_fp.vaporPressure(temperature);
  
   const DualReal Yncg = gas.mass_fraction[_gas_fluid_component];
   const DualReal Xncg = liquid.mass_fraction[_gas_fluid_component];
  
   DualReal ncg_density = _ncg_fp.rho_from_p_T(Yncg * pressure, temperature);
   DualReal vapor_density = _water_fp.rho_from_p_T((1.0 - Xncg) * psat, temperature);
  
   // Density is just the sum of individual component densities
   return ncg_density + vapor_density;
 }

Referenced by saturation(), and totalMassFraction().

◆ 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; };

◆ gasProperties()

void PorousFlowWaterNCG::gasProperties	(	const DualReal &	pressure,
		const DualReal &	temperature,
		std::vector< FluidStateProperties > &	fsp
	)		const

Gas properties - density, viscosity and enthalpy.

Parameters

	pressure	gas pressure (Pa)
	temperature	temperature (K)
[out]	FluidStateProperties	data structure

Definition at line 203 of file PorousFlowWaterNCG.C.

 {
   FluidStateProperties & liquid = fsp[_aqueous_phase_number];
   FluidStateProperties & gas = fsp[_gas_phase_number];
  
   const DualReal psat = _water_fp.vaporPressure(temperature);
  
   const DualReal Yncg = gas.mass_fraction[_gas_fluid_component];
   const DualReal Xncg = liquid.mass_fraction[_gas_fluid_component];
  
   // NCG density, viscosity and enthalpy calculated using partial pressure
   // Yncg * gas_poreressure (Dalton's law)
   DualReal ncg_density, ncg_viscosity;
   _ncg_fp.rho_mu_from_p_T(Yncg * pressure, temperature, ncg_density, ncg_viscosity);
   DualReal ncg_enthalpy = _ncg_fp.h_from_p_T(Yncg * pressure, temperature);
  
   // Vapor density, viscosity and enthalpy calculated using partial pressure
   // X1 * psat (Raoult's law)
   DualReal vapor_density, vapor_viscosity;
  
   _water_fp.rho_mu_from_p_T((1.0 - Xncg) * psat, temperature, vapor_density, vapor_viscosity);
   DualReal vapor_enthalpy = _water_fp.h_from_p_T((1.0 - Xncg) * psat, temperature);
  
   // Density is just the sum of individual component densities
   gas.density = ncg_density + vapor_density;
  
   // Viscosity of the gas phase is a weighted sum of the individual viscosities
   gas.viscosity = Yncg * ncg_viscosity + (1.0 - Yncg) * vapor_viscosity;
  
   // Enthalpy of the gas phase is a weighted sum of the individual enthalpies
   gas.enthalpy = Yncg * ncg_enthalpy + (1.0 - Yncg) * vapor_enthalpy;
  
   //  Internal energy of the gas phase (e = h - pv)
   mooseAssert(gas.density.value() > 0.0, "Gas density must be greater than zero");
   gas.internal_energy = gas.enthalpy - pressure / gas.density;
 }

Referenced by thermophysicalProperties(), and twoPhaseProperties().

◆ getPressureIndex()

unsigned int PorousFlowFluidStateMultiComponentBase::getPressureIndex ( ) const

inlineinherited

Definition at line 70 of file PorousFlowFluidStateMultiComponentBase.h.

70 { return _pidx; };

◆ getTemperatureIndex()

unsigned int PorousFlowFluidStateMultiComponentBase::getTemperatureIndex ( ) const

inlineinherited

Definition at line 71 of file PorousFlowFluidStateMultiComponentBase.h.

71 { return _Tidx; };

◆ getXIndex()

unsigned int PorousFlowFluidStateMultiComponentBase::getXIndex ( ) const

inlineinherited

Definition at line 73 of file PorousFlowFluidStateMultiComponentBase.h.

73 { return _Xidx; };

◆ getZIndex()

unsigned int PorousFlowFluidStateMultiComponentBase::getZIndex ( ) const

inlineinherited

Definition at line 72 of file PorousFlowFluidStateMultiComponentBase.h.

72 { return _Zidx; };

◆ initialize()

void PorousFlowFluidStateBase::initialize ( )

inlinefinalinherited

Definition at line 60 of file PorousFlowFluidStateBase.h.

60 {};

◆ liquidDensity()

DualReal PorousFlowWaterNCG::liquidDensity	(	const DualReal &	pressure,
		const DualReal &	temperature
	)		const

Density of the liquid phase Note: The pressure here is the gas pressure.

As a result, the liquid pressure can include a dependence on saturation due to the capillary pressure, so this method should be called after the saturation is calculated for the two phase case ie: after calling saturation(). For the single phase liquid case, it is ok to call this method by itself, as gas saturation is initialized to zero.

Parameters

pressure	gas pressure (Pa)
temperature	temperature (K)

Returns: liquid density (kg/m^3)

Definition at line 274 of file PorousFlowWaterNCG.C.

 {
   return _water_fp.rho_from_p_T(pressure, temperature);
 }

Referenced by saturation(), and totalMassFraction().

◆ liquidProperties()

void PorousFlowWaterNCG::liquidProperties	(	const DualReal &	pressure,
		const DualReal &	temperature,
		std::vector< FluidStateProperties > &	fsp
	)		const

Liquid properties - density, viscosity and enthalpy Note: The pressure here is the liquid pressure.

In this class, enthalpy includes a contribution due to the enthalpy of dissolution of the NCG into the liquid phase. As a result, the derivatives can include a dependence on the capillary pressure, so this method should be called after the saturation is calculated for the two phase case ie: after calling saturation(). For the single phase liquid case, it is ok to call this method by itself, as gas saturation is initialized to zero.

Parameters

	pressure	liquid pressure (Pa)
	temperature	temperature (K)
[out]	FluidStateProperties	data structure

Definition at line 243 of file PorousFlowWaterNCG.C.

 {
   FluidStateProperties & liquid = fsp[_aqueous_phase_number];
  
   // Calculate liquid density and viscosity if in the two phase or single phase
   // liquid region, assuming they are not affected by the presence of dissolved
   // NCG. Note: the (small) contribution due to derivative of capillary pressure
   // wrt pressure (using the chain rule) is not implemented.
   DualReal liquid_density, liquid_viscosity;
   _water_fp.rho_mu_from_p_T(pressure, temperature, liquid_density, liquid_viscosity);
  
   liquid.density = liquid_density;
   liquid.viscosity = liquid_viscosity;
  
   // Enthalpy does include a contribution due to the enthalpy of dissolution
   const DualReal hdis = enthalpyOfDissolution(temperature);
  
   const DualReal water_enthalpy = _water_fp.h_from_p_T(pressure, temperature);
   const DualReal ncg_enthalpy = _ncg_fp.h_from_p_T(pressure, temperature);
  
   const DualReal Xncg = liquid.mass_fraction[_gas_fluid_component];
   liquid.enthalpy = (1.0 - Xncg) * water_enthalpy + Xncg * (ncg_enthalpy + hdis);
  
   //  Internal energy of the liquid phase (e = h - pv)
   mooseAssert(liquid.density.value() > 0.0, "Liquid density must be greater than zero");
   liquid.internal_energy = liquid.enthalpy - pressure / liquid.density;
 }

Referenced by thermophysicalProperties(), and twoPhaseProperties().

◆ massFractions()

void PorousFlowWaterNCG::massFractions	(	const DualReal &	pressure,
		const DualReal &	temperature,
		const DualReal &	Z,
		FluidStatePhaseEnum &	phase_state,
		std::vector< FluidStateProperties > &	fsp
	)		const

Mass fractions of NCG and H2O in both phases, as well as derivatives wrt PorousFlow variables.

Values depend on the phase state (liquid, gas or two phase)

Parameters

	pressure	phase pressure (Pa)
	temperature	phase temperature (K)
	Z	total mass fraction of NCG component
[out]	PhaseStateEnum	current phase state
[out]	FluidStateMassFractions	data structure

Definition at line 139 of file PorousFlowWaterNCG.C.

 {
   FluidStateProperties & liquid = fsp[_aqueous_phase_number];
   FluidStateProperties & gas = fsp[_gas_phase_number];
  
   // Equilibrium mass fraction of NCG in liquid and H2O in gas phases
   DualReal Xncg, Yh2o;
   equilibriumMassFractions(pressure, temperature, Xncg, Yh2o);
  
   DualReal Yncg = 1.0 - Yh2o;
  
   // Determine which phases are present based on the value of Z
   phaseState(Z.value(), Xncg.value(), Yncg.value(), phase_state);
  
   // The equilibrium mass fractions calculated above are only correct in the two phase
   // state. If only liquid or gas phases are present, the mass fractions are given by
   // the total mass fraction Z.
   DualReal Xh2o = 0.0;
  
   switch (phase_state)
   {
     case FluidStatePhaseEnum::LIQUID:
     {
       Xncg = Z;
       Yncg = 0.0;
       Xh2o = 1.0 - Z;
       Yh2o = 0.0;
       Moose::derivInsert(Xncg.derivatives(), _pidx, 0.0);
       Moose::derivInsert(Xncg.derivatives(), _Tidx, 0.0);
       Moose::derivInsert(Xncg.derivatives(), _Zidx, 1.0);
       break;
     }
  
     case FluidStatePhaseEnum::GAS:
     {
       Xncg = 0.0;
       Yncg = Z;
       Yh2o = 1.0 - Z;
       Moose::derivInsert(Yncg.derivatives(), _pidx, 0.0);
       Moose::derivInsert(Yncg.derivatives(), _Tidx, 0.0);
       Moose::derivInsert(Yncg.derivatives(), _Zidx, 1.0);
       break;
     }
  
     case FluidStatePhaseEnum::TWOPHASE:
     {
       // Keep equilibrium mass fractions
       Xh2o = 1.0 - Xncg;
       break;
     }
   }
  
   // Save the mass fractions in the FluidStateMassFractions object
   liquid.mass_fraction[_aqueous_fluid_component] = Xh2o;
   liquid.mass_fraction[_gas_fluid_component] = Xncg;
   gas.mass_fraction[_aqueous_fluid_component] = Yh2o;
   gas.mass_fraction[_gas_fluid_component] = Yncg;
 }

Referenced by thermophysicalProperties().

◆ moleFractionToMassFraction()

DualReal PorousFlowWaterNCG::moleFractionToMassFraction ( const DualReal & xmol ) const

protected

Convert mole fraction to mass fraction.

Parameters

xmol	mole fraction

Returns: mass fraction

Definition at line 378 of file PorousFlowWaterNCG.C.

 {
   return xmol * _Mncg / (xmol * _Mncg + (1.0 - xmol) * _Mh2o);
 }

Referenced by equilibriumMassFractions().

◆ 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().

◆ phaseState()

void PorousFlowFluidStateMultiComponentBase::phaseState	(	Real	Zi,
		Real	Xi,
		Real	Yi,
		FluidStatePhaseEnum &	phase_state
	)		const

inherited

Determines the phase state gven the total mass fraction and equilibrium mass fractions.

Parameters

	Zi	total mass fraction
	Xi	equilibrium mass fraction in liquid
	Yi	equilibrium mass fraction in gas
[out]	phase_state	the phase state (gas, liquid, two phase)

Definition at line 28 of file PorousFlowFluidStateMultiComponentBase.C.

 {
   if (Zi <= Xi)
   {
     // In this case, there is not enough component i to form a gas phase,
     // so only a liquid phase is present
     phase_state = FluidStatePhaseEnum::LIQUID;
   }
   else if (Zi > Xi && Zi < Yi)
   {
     // Two phases are present
     phase_state = FluidStatePhaseEnum::TWOPHASE;
   }
   else // (Zi >= Yi)
   {
     // In this case, there is not enough water to form a liquid
     // phase, so only a gas phase is present
     phase_state = FluidStatePhaseEnum::GAS;
   }
 }

Referenced by massFractions(), and PorousFlowBrineCO2::massFractions().

◆ rachfordRice()

Real PorousFlowFluidStateFlash::rachfordRice	(	Real	vf,
		std::vector< Real > &	Zi,
		std::vector< Real > &	Ki
	)		const

inherited

Rachford-Rice equation for vapor fraction.

Can be solved analytically for two components in two phases, but must be solved iteratively using a root finding algorithm for more components. This equation has the nice property that it is monotonic in the interval [0,1], so that only a small number of iterations are typically required to find the root.

The Rachford-Rice equation can also be used to check whether the phase state is two phase, single phase gas, or single phase liquid. Evaluate f(v), the Rachford-Rice equation evaluated at the vapor mass fraction.

If f(0) < 0, then the mixture is below the bubble point, and only a single phase liquid can exist

If f(1) > 0, then the mixture is above the dew point, and only a single phase gas exists.

If f(0) >= 0 and f(1) <= 0, the mixture is between the bubble and dew points, and both gas and liquid phases exist.

Parameters

vf	vapor fraction
Zi	mass fractions
Ki	equilibrium constants

Returns: f(x)

Definition at line 27 of file PorousFlowFluidStateFlash.C.

 {
   const std::size_t num_z = Zi.size();
   // Check that the sizes of the mass fractions and equilibrium constant vectors are correct
   if (Ki.size() != num_z + 1)
     mooseError("The number of mass fractions or equilibrium components passed to rachfordRice is "
                "not correct");
  
   Real f = 0.0;
   Real Z_total = 0.0;
  
   for (std::size_t i = 0; i < num_z; ++i)
   {
     f += Zi[i] * (Ki[i] - 1.0) / (1.0 + x * (Ki[i] - 1.0));
     Z_total += Zi[i];
   }
  
   // Add the last component (with total mass fraction = 1 - z_total)
   f += (1.0 - Z_total) * (Ki[num_z] - 1.0) / (1.0 + x * (Ki[num_z] - 1.0));
  
   return f;
 }

Referenced by PorousFlowFluidStateFlash::vaporMassFraction().

◆ rachfordRiceDeriv()

Real PorousFlowFluidStateFlash::rachfordRiceDeriv	(	Real	vf,
		std::vector< Real > &	Zi,
		std::vector< Real > &	Ki
	)		const

inherited

Derivative of Rachford-Rice equation wrt vapor fraction.

Has the nice property that it is strictly negative in the interval [0,1]

Parameters

vf	vapor fraction
Zi	mass fractions
Ki	equilibrium constants

Returns: f'(x)

Definition at line 53 of file PorousFlowFluidStateFlash.C.

 {
   const std::size_t num_Z = Zi.size();
   // Check that the sizes of the mass fractions and equilibrium constant vectors are correct
   if (Ki.size() != num_Z + 1)
     mooseError("The number of mass fractions or equilibrium components passed to rachfordRice is "
                "not correct");
  
   Real df = 0.0;
   Real Z_total = 0.0;
  
   for (std::size_t i = 0; i < num_Z; ++i)
   {
     df -= Zi[i] * (Ki[i] - 1.0) * (Ki[i] - 1.0) / (1.0 + x * (Ki[i] - 1.0)) /
           (1.0 + x * (Ki[i] - 1.0));
     Z_total += Zi[i];
   }
  
   // Add the last component (with total mass fraction = 1 - z_total)
   df -= (1.0 - Z_total) * (Ki[num_Z] - 1.0) * (Ki[num_Z] - 1.0) / (1.0 + x * (Ki[num_Z] - 1.0)) /
         (1.0 + x * (Ki[num_Z] - 1.0));
  
   return df;
 }

Referenced by PorousFlowFluidStateFlash::vaporMassFraction().

◆ 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; };

◆ saturation()

DualReal PorousFlowWaterNCG::saturation	(	const DualReal &	pressure,
		const DualReal &	temperature,
		const DualReal &	Z,
		std::vector< FluidStateProperties > &	fsp
	)		const

Gas saturation in the two-phase region.

Parameters

	pressure	gas pressure (Pa)
	temperature	phase temperature (K)
	Z	total mass fraction of NCG component
[out]	FluidStateProperties	data structure

Returns: gas saturation (-)

Definition at line 300 of file PorousFlowWaterNCG.C.

 {
   auto & gas = fsp[_gas_phase_number];
   auto & liquid = fsp[_aqueous_fluid_component];
  
   // Approximate liquid density as saturation isn't known yet, by using the gas
   // pressure rather than the liquid pressure. This does result in a small error
   // in the calculated saturation, but this is below the error associated with
   // the correlations. A more accurate saturation could be found iteraviely,
   // at the cost of increased computational expense
  
   // The gas and liquid densities
   const DualReal gas_density = gasDensity(pressure, temperature, fsp);
   const DualReal liquid_density = liquidDensity(pressure, temperature);
  
   // Set mass equilibrium constants used in the calculation of vapor mass fraction
   const DualReal Xncg = liquid.mass_fraction[_gas_fluid_component];
   const DualReal Yncg = gas.mass_fraction[_gas_fluid_component];
  
   const DualReal K0 = Yncg / Xncg;
   const DualReal K1 = (1.0 - Yncg) / (1.0 - Xncg);
   const DualReal vapor_mass_fraction = vaporMassFraction(Z, K0, K1);
  
   // The gas saturation in the two phase case
   const DualReal saturation = vapor_mass_fraction * liquid_density /
                               (gas_density + vapor_mass_fraction * (liquid_density - gas_density));
  
   return saturation;
 }

Referenced by totalMassFraction(), and twoPhaseProperties().

◆ thermophysicalProperties()

void PorousFlowWaterNCG::thermophysicalProperties	(	Real	pressure,
		Real	temperature,
		Real	Xnacl,
		Real	Z,
		unsigned int	qp,
		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)
	temperature	fluid temperature (K)
	Xnacl	mass fraction of NaCl
	Z	total mass fraction of fluid component
	qp	quadpoint index
[out]	fsp	the FluidStateProperties struct containing all properties

Implements PorousFlowFluidStateMultiComponentBase.

Definition at line 72 of file PorousFlowWaterNCG.C.

 {
   FluidStateProperties & liquid = fsp[_aqueous_phase_number];
   FluidStateProperties & gas = fsp[_gas_phase_number];
  
   // Check whether the input temperature is within the region of validity
   checkVariables(temperature);
  
   // AD versions of primary variables
   DualReal p = pressure;
   Moose::derivInsert(p.derivatives(), _pidx, 1.0);
   DualReal T = temperature;
   Moose::derivInsert(T.derivatives(), _Tidx, 1.0);
   DualReal Zncg = Z;
   Moose::derivInsert(Zncg.derivatives(), _Zidx, 1.0);
  
   // Clear all of the FluidStateProperties data
   clearFluidStateProperties(fsp);
  
   FluidStatePhaseEnum phase_state;
   massFractions(p, T, Zncg, phase_state, fsp);
  
   switch (phase_state)
   {
     case FluidStatePhaseEnum::GAS:
     {
       // Set the gas saturations
       gas.saturation = 1.0;
  
       // Calculate gas properties
       gasProperties(p, T, fsp);
  
       break;
     }
  
     case FluidStatePhaseEnum::LIQUID:
     {
       // Calculate the liquid properties
       const DualReal liquid_pressure = p - _pc.capillaryPressure(1.0, qp);
       liquidProperties(liquid_pressure, T, fsp);
  
       break;
     }
  
     case FluidStatePhaseEnum::TWOPHASE:
     {
       // Calculate the gas and liquid properties in the two phase region
       twoPhaseProperties(p, T, Zncg, qp, fsp);
  
       break;
     }
   }
  
   // Liquid saturations can now be set
   liquid.saturation = 1.0 - gas.saturation;
  
   // Save pressures to FluidStateProperties object
   gas.pressure = p;
   liquid.pressure = p - _pc.capillaryPressure(liquid.saturation, qp);
 }

◆ totalMassFraction()

Real PorousFlowWaterNCG::totalMassFraction	(	Real	pressure,
		Real	temperature,
		Real	Xnacl,
		Real	saturation,
		unsigned int	qp
	)		const

overridevirtual

Total mass fraction of fluid component summed over all phases in the two-phase state for a specified gas saturation.

Parameters

pressure	gas pressure (Pa)
temperature	temperature (K)
Xnacl	NaCl mass fraction (kg/kg)
saturation	gas saturation (-)
qp	quadpoint index

Returns: total mass fraction Z (-)

Implements PorousFlowFluidStateMultiComponentBase.

Definition at line 416 of file PorousFlowWaterNCG.C.

 {
   // Check whether the input temperature is within the region of validity
   checkVariables(temperature);
  
   // As we do not require derivatives, we can simply ignore their initialisation
   const DualReal p = pressure;
   const DualReal T = temperature;
  
   // FluidStateProperties data structure
   std::vector<FluidStateProperties> fsp(_num_phases, FluidStateProperties(_num_components));
   auto & liquid = fsp[_aqueous_phase_number];
   auto & gas = fsp[_gas_phase_number];
  
   // Calculate equilibrium mass fractions in the two-phase state
   DualReal Xncg, Yh2o;
   equilibriumMassFractions(p, T, Xncg, Yh2o);
  
   // Save the mass fractions in the FluidStateMassFractions object to calculate gas density
   const DualReal Yncg = 1.0 - Yh2o;
   liquid.mass_fraction[_aqueous_fluid_component] = 1.0 - Xncg;
   liquid.mass_fraction[_gas_fluid_component] = Xncg;
   gas.mass_fraction[_aqueous_fluid_component] = Yh2o;
   gas.mass_fraction[_gas_fluid_component] = Yncg;
  
   // Gas density
   const Real gas_density = gasDensity(p, T, fsp).value();
  
   // Liquid density
   const DualReal liquid_pressure = p - _pc.capillaryPressure(1.0 - saturation, qp);
   const Real liquid_density = liquidDensity(liquid_pressure, T).value();
  
   // The total mass fraction of ncg (Z) can now be calculated
   const Real Z = (saturation * gas_density * Yncg.value() +
                   (1.0 - saturation) * liquid_density * Xncg.value()) /
                  (saturation * gas_density + (1.0 - saturation) * liquid_density);
  
   return Z;
 }

◆ twoPhaseProperties()

void PorousFlowWaterNCG::twoPhaseProperties	(	const DualReal &	pressure,
		const DualReal &	temperature,
		const DualReal &	Z,
		unsigned int	qp,
		std::vector< FluidStateProperties > &	fsp
	)		const

Gas and liquid properties in the two-phase region.

Parameters

	pressure	gas pressure (Pa)
	temperature	phase temperature (K)
	Z	total mass fraction of NCG component
	qp	quadpoint for capillary pressure
[out]	FluidStateProperties	data structure

Definition at line 334 of file PorousFlowWaterNCG.C.

 {
   auto & gas = fsp[_gas_phase_number];
  
   // Calculate all of the gas phase properties, as these don't depend on saturation
   gasProperties(pressure, temperature, fsp);
  
   // The gas saturation in the two phase case
   gas.saturation = saturation(pressure, temperature, Z, fsp);
  
   // The liquid pressure and properties can now be calculated
   const DualReal liquid_pressure = pressure - _pc.capillaryPressure(1.0 - gas.saturation, qp);
   liquidProperties(liquid_pressure, temperature, fsp);
 }

Referenced by thermophysicalProperties().

◆ vaporMassFraction() [1/3]

DualReal PorousFlowFluidStateFlash::vaporMassFraction	(	const DualReal &	Z0,
		const DualReal &	K0,
		const DualReal &	K1
	)		const

inherited

Definition at line 87 of file PorousFlowFluidStateFlash.C.

 {
   return (Z0 * (K1 - K0) - (K1 - 1.0)) / ((K0 - 1.0) * (K1 - 1.0));
 }

◆ vaporMassFraction() [2/3]

Real PorousFlowFluidStateFlash::vaporMassFraction	(	Real	Z0,
		Real	K0,
		Real	K1
	)		const

inherited

Solves Rachford-Rice equation to provide vapor mass fraction.

For two components, the analytical solution is used, while for cases with more than two components, a Newton-Raphson iterative solution is calculated.

Parameters

Zi	total mass fraction(s)
Ki	equilibrium constant(s)

Returns: vapor mass fraction

Definition at line 81 of file PorousFlowFluidStateFlash.C.

 {
   return (Z0 * (K1 - K0) - (K1 - 1.0)) / ((K0 - 1.0) * (K1 - 1.0));
 }

Referenced by saturation(), PorousFlowBrineCO2::saturation(), and PorousFlowFluidStateFlash::vaporMassFraction().

◆ vaporMassFraction() [3/3]

Real PorousFlowFluidStateFlash::vaporMassFraction	(	std::vector< Real > &	Zi,
		std::vector< Real > &	Ki
	)		const

inherited

Definition at line 95 of file PorousFlowFluidStateFlash.C.

 {
   // Check that the sizes of the mass fractions and equilibrium constant vectors are correct
   if (Ki.size() != Zi.size() + 1)
     mooseError("The number of mass fractions or equilibrium components passed to rachfordRice is "
                "not correct");
   Real v;
  
   // If there are only two components, an analytical solution is possible
   if (Ki.size() == 2)
     v = vaporMassFraction(Zi[0], Ki[0], Ki[1]);
   else
   {
     // More than two components - solve the Rachford-Rice equation using
     // Newton-Raphson method
     // Initial guess for vapor mass fraction
     Real v0 = 0.5;
     unsigned int iter = 0;
  
     while (std::abs(rachfordRice(v0, Zi, Ki)) > _nr_tol)
     {
       v0 = v0 - rachfordRice(v0, Zi, Ki) / rachfordRiceDeriv(v0, Zi, Ki);
       iter++;
  
       if (iter > _nr_max_its)
         break;
     }
     v = v0;
   }
   return v;
 }

Member Data Documentation

◆ _aqueous_fluid_component

const unsigned int PorousFlowFluidStateBase::_aqueous_fluid_component

protectedinherited

Fluid component number of the aqueous component.

Definition at line 127 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::aqueousComponentIndex(), massFractions(), PorousFlowBrineCO2::massFractions(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG(), PorousFlowWaterVapor::PorousFlowWaterVapor(), saturation(), PorousFlowBrineCO2::saturation(), totalMassFraction(), and PorousFlowBrineCO2::totalMassFraction().

◆ _aqueous_phase_number

const unsigned int PorousFlowFluidStateBase::_aqueous_phase_number

protectedinherited

Phase number of the aqueous phase.

Definition at line 123 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::aqueousPhaseIndex(), gasDensity(), gasProperties(), liquidProperties(), PorousFlowBrineCO2::liquidProperties(), massFractions(), PorousFlowBrineCO2::massFractions(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG(), PorousFlowWaterVapor::PorousFlowWaterVapor(), PorousFlowWaterVapor::thermophysicalProperties(), thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), totalMassFraction(), and PorousFlowBrineCO2::totalMassFraction().

◆ _empty_fsp

FluidStateProperties PorousFlowFluidStateBase::_empty_fsp

protectedinherited

Empty FluidStateProperties object.

Definition at line 139 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::clearFluidStateProperties(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG(), and PorousFlowWaterVapor::PorousFlowWaterVapor().

◆ _gas_fluid_component

unsigned int PorousFlowFluidStateBase::_gas_fluid_component

protectedinherited

Fluid component number of the gas phase.

Definition at line 129 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::gasComponentIndex(), gasDensity(), gasProperties(), liquidProperties(), PorousFlowBrineCO2::liquidProperties(), massFractions(), PorousFlowBrineCO2::massFractions(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG(), PorousFlowWaterVapor::PorousFlowWaterVapor(), saturation(), PorousFlowBrineCO2::saturation(), totalMassFraction(), and PorousFlowBrineCO2::totalMassFraction().

◆ _gas_phase_number

unsigned int PorousFlowFluidStateBase::_gas_phase_number

protectedinherited

◆ _Mh2o

const Real PorousFlowWaterNCG::_Mh2o

protected

Molar mass of water (kg/mol)

Definition at line 192 of file PorousFlowWaterNCG.h.

Referenced by moleFractionToMassFraction().

◆ _Mncg

const Real PorousFlowWaterNCG::_Mncg

protected

Molar mass of non-condensable gas (kg/mol)

Definition at line 194 of file PorousFlowWaterNCG.h.

Referenced by enthalpyOfDissolution(), and moleFractionToMassFraction().

◆ _ncg_fp

const SinglePhaseFluidProperties& PorousFlowWaterNCG::_ncg_fp

protected

Fluid properties UserObject for the NCG.

Definition at line 190 of file PorousFlowWaterNCG.h.

Referenced by gasDensity(), gasProperties(), and liquidProperties().

◆ _ncg_henry

const std::vector<Real> PorousFlowWaterNCG::_ncg_henry

protected

Henry's coefficients for the NCG.

Definition at line 200 of file PorousFlowWaterNCG.h.

Referenced by enthalpyOfDissolution(), and equilibriumMassFractions().

◆ _nr_max_its

const Real PorousFlowFluidStateFlash::_nr_max_its

protectedinherited

Maximum number of iterations for the Newton-Raphson routine.

Definition at line 80 of file PorousFlowFluidStateFlash.h.

Referenced by PorousFlowFluidStateFlash::vaporMassFraction().

◆ _nr_tol

const Real PorousFlowFluidStateFlash::_nr_tol

protectedinherited

Tolerance for Newton-Raphson iterations.

Definition at line 82 of file PorousFlowFluidStateFlash.h.

Referenced by PorousFlowFluidStateFlash::vaporMassFraction().

◆ _num_components

unsigned int PorousFlowFluidStateBase::_num_components

protectedinherited

Number of components.

Definition at line 121 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::numComponents(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG(), PorousFlowWaterVapor::PorousFlowWaterVapor(), totalMassFraction(), and PorousFlowBrineCO2::totalMassFraction().

◆ _num_phases

unsigned int PorousFlowFluidStateBase::_num_phases

protectedinherited

Number of phases.

Definition at line 119 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::numPhases(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG(), PorousFlowWaterVapor::PorousFlowWaterVapor(), totalMassFraction(), and PorousFlowBrineCO2::totalMassFraction().

◆ _pc

const PorousFlowCapillaryPressure& PorousFlowFluidStateBase::_pc

protectedinherited

Capillary pressure UserObject.

Definition at line 137 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowWaterVapor::thermophysicalProperties(), thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), totalMassFraction(), PorousFlowBrineCO2::totalMassFraction(), twoPhaseProperties(), and PorousFlowBrineCO2::twoPhaseProperties().

◆ _pidx

const unsigned int PorousFlowFluidStateMultiComponentBase::_pidx

protectedinherited

Index of derivative wrt pressure.

Definition at line 73 of file PorousFlowFluidStateMultiComponentBase.h.

Referenced by PorousFlowBrineCO2::equilibriumMoleFractions(), PorousFlowFluidStateMultiComponentBase::getPressureIndex(), massFractions(), PorousFlowBrineCO2::massFractions(), thermophysicalProperties(), and PorousFlowBrineCO2::thermophysicalProperties().

◆ _R

const Real PorousFlowFluidStateBase::_R

protectedinherited

Universal gas constant (J/mol/K)

Definition at line 133 of file PorousFlowFluidStateBase.h.

Referenced by enthalpyOfDissolution(), and PorousFlowBrineCO2::enthalpyOfDissolutionGas().

◆ _salt_component

const unsigned int PorousFlowFluidStateBase::_salt_component

protectedinherited

Salt component index.

Definition at line 131 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::saltComponentIndex().

◆ _T_c2k

const Real PorousFlowFluidStateBase::_T_c2k

protectedinherited

Conversion from C to K.

Definition at line 135 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowBrineCO2::equilibriumConstantCO2(), PorousFlowBrineCO2::equilibriumConstantH2O(), PorousFlowBrineCO2::funcABHighTemp(), PorousFlowBrineCO2::henryConstant(), and PorousFlowBrineCO2::partialDensityCO2().

◆ _Tidx

const unsigned int PorousFlowFluidStateMultiComponentBase::_Tidx

protectedinherited

Index of derivative wrt temperature.

Definition at line 81 of file PorousFlowFluidStateMultiComponentBase.h.

Referenced by enthalpyOfDissolution(), PorousFlowBrineCO2::enthalpyOfDissolutionGas(), PorousFlowBrineCO2::equilibriumMoleFractions(), PorousFlowFluidStateMultiComponentBase::getTemperatureIndex(), massFractions(), PorousFlowBrineCO2::massFractions(), thermophysicalProperties(), and PorousFlowBrineCO2::thermophysicalProperties().

◆ _water97_fp

const Water97FluidProperties& PorousFlowWaterNCG::_water97_fp

protected

Fluid properties UserObject for water (used to access Henry's law)

Definition at line 188 of file PorousFlowWaterNCG.h.

Referenced by enthalpyOfDissolution(), and equilibriumMassFractions().

◆ _water_critical_temperature

const Real PorousFlowWaterNCG::_water_critical_temperature

protected

Critical temperature of water (K)

Definition at line 198 of file PorousFlowWaterNCG.h.

Referenced by checkVariables().

◆ _water_fp

const SinglePhaseFluidProperties& PorousFlowWaterNCG::_water_fp

protected

Fluid properties UserObject for water.

Definition at line 186 of file PorousFlowWaterNCG.h.

Referenced by equilibriumMassFractions(), gasDensity(), gasProperties(), liquidDensity(), liquidProperties(), and PorousFlowWaterNCG().

◆ _water_triple_temperature

const Real PorousFlowWaterNCG::_water_triple_temperature

protected

Triple point temperature of water (K)

Definition at line 196 of file PorousFlowWaterNCG.h.

◆ _Xidx

const unsigned int PorousFlowFluidStateMultiComponentBase::_Xidx

protectedinherited

Index of derivative wrt salt mass fraction X.

Definition at line 83 of file PorousFlowFluidStateMultiComponentBase.h.

Referenced by PorousFlowBrineCO2::equilibriumMoleFractions(), PorousFlowFluidStateMultiComponentBase::getXIndex(), PorousFlowBrineCO2::massFractions(), and PorousFlowBrineCO2::thermophysicalProperties().

◆ _Zidx

const unsigned int PorousFlowFluidStateMultiComponentBase::_Zidx

protectedinherited

Index of derivative wrt total mass fraction Z.

Definition at line 79 of file PorousFlowFluidStateMultiComponentBase.h.

Referenced by PorousFlowFluidStateMultiComponentBase::getZIndex(), massFractions(), PorousFlowBrineCO2::massFractions(), thermophysicalProperties(), and PorousFlowBrineCO2::thermophysicalProperties().

The documentation for this class was generated from the following files:

porous_flow/include/userobjects/PorousFlowWaterNCG.h
porous_flow/src/userobjects/PorousFlowWaterNCG.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowWaterNCG()

Member Function Documentation

◆ aqueousComponentIndex()

◆ aqueousPhaseIndex()

◆ checkVariables()

◆ clearFluidStateProperties()

◆ enthalpyOfDissolution()

◆ equilibriumMassFractions()

◆ execute()

◆ finalize()

◆ fluidStateName()

◆ gasComponentIndex()

◆ gasDensity()

◆ gasPhaseIndex()

◆ gasProperties()

◆ getPressureIndex()

◆ getTemperatureIndex()

◆ getXIndex()

◆ getZIndex()

◆ initialize()

◆ liquidDensity()

◆ liquidProperties()

◆ massFractions()

◆ moleFractionToMassFraction()

◆ numComponents()

◆ numPhases()

◆ phaseState()

◆ rachfordRice()

◆ rachfordRiceDeriv()

◆ saltComponentIndex()

◆ saturation()

◆ thermophysicalProperties()

◆ totalMassFraction()

◆ twoPhaseProperties()

◆ vaporMassFraction() [1/3]

◆ vaporMassFraction() [2/3]

◆ vaporMassFraction() [3/3]

Member Data Documentation

◆ _aqueous_fluid_component

◆ _aqueous_phase_number

◆ _empty_fsp

◆ _gas_fluid_component

◆ _gas_phase_number

◆ _Mh2o

◆ _Mncg

◆ _ncg_fp

◆ _ncg_henry

◆ _nr_max_its

◆ _nr_tol

◆ _num_components

◆ _num_phases

◆ _pc

◆ _pidx

◆ _R

◆ _salt_component

◆ _T_c2k

◆ _Tidx

◆ _water97_fp

◆ _water_critical_temperature

◆ _water_fp

◆ _water_triple_temperature

◆ _Xidx

◆ _Zidx