Fluid state class for brine and CO2. More...

#include <PorousFlowFluidStateBrineCO2.h>

Inheritance diagram for PorousFlowFluidStateBrineCO2:

Public Member Functions
	PorousFlowFluidStateBrineCO2 (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

void	setMaterialVectorSize () const
	Size material property vectors and initialise with zeros. More...

virtual void	thermophysicalProperties ()
	Calculates all required thermophysical properties and derivatives for each phase and fluid component. More...

Protected Attributes
const VariableValue &	_gas_porepressure
	Porepressure. More...

const VariableGradient &	_gas_gradp_qp
	Gradient of porepressure (only defined at the qps) More...

const unsigned int	_gas_porepressure_varnum
	Moose variable number of the gas porepressure. More...

const unsigned int	_pvar
	PorousFlow variable number of the gas porepressure. More...

std::vector< const VariableValue * >	_Z
	Total mass fraction(s) of the gas component(s) summed over all phases. More...

std::vector< const VariableGradient * >	_gradZ_qp
	Gradient(s) of total mass fraction(s) of the gas component(s) (only defined at the qps) More...

std::vector< unsigned int >	_Z_varnum
	Moose variable number of Z. More...

std::vector< unsigned int >	_Zvar
	PorousFlow variable number of Z. More...

const unsigned int	_num_Z_vars
	Number of coupled total mass fractions. Should be _num_phases - 1. More...

const VariableValue &	_Xnacl
	Salt mass fraction (kg/kg) More...

const VariableGradient &	_grad_Xnacl_qp
	Gradient of salt mass fraction (only defined at the qps) More...

const unsigned int	_Xnacl_varnum
	Salt mass fraction variable number. More...

const unsigned int	_Xvar
	Salt mass fraction PorousFlow variable number. More...

const PorousFlowFluidStateMultiComponentBase &	_fs
	FluidState UserObject. More...

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

const 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...

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

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

const MaterialProperty< Real > &	_temperature
	Temperature. More...

const MaterialProperty< RealGradient > &	_gradT_qp
	Gradient of temperature (only defined at the qps) More...

const MaterialProperty< std::vector< Real > > &	_dtemperature_dvar
	Derivative of temperature wrt PorousFlow variables. More...

const unsigned int	_temperature_varnum
	Moose variable number of the temperature. More...

const unsigned int	_Tvar
	PorousFlow variable number of the temperature. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_mass_frac
	Mass fraction matrix. More...

MaterialProperty< std::vector< std::vector< RealGradient > > > *	_grad_mass_frac_qp
	Gradient of the mass fraction matrix (only defined at the qps) More...

MaterialProperty< std::vector< std::vector< std::vector< Real > > > > &	_dmass_frac_dvar
	Derivative of the mass fraction matrix with respect to the Porous Flow variables. More...

MaterialProperty< std::vector< Real > > &	_fluid_density
	Fluid density of each phase. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dfluid_density_dvar
	Derivative of the fluid density for each phase wrt PorousFlow variables. More...

MaterialProperty< std::vector< Real > > &	_fluid_viscosity
	Viscosity of each phase. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dfluid_viscosity_dvar
	Derivative of the fluid viscosity for each phase wrt PorousFlow variables. More...

MaterialProperty< std::vector< Real > > &	_fluid_enthalpy
	Enthalpy of each phase. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dfluid_enthalpy_dvar
	Derivative of the fluid enthalpy for each phase wrt PorousFlow variables. More...

MaterialProperty< std::vector< Real > > &	_fluid_internal_energy
	Internal energy of each phase. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dfluid_internal_energy_dvar
	Derivative of the fluid internal energy for each phase wrt PorousFlow variables. More...

const Real	_T_c2k
	Conversion from degrees Celsius to degrees Kelvin. More...

bool	_is_initqp
	Flag to indicate whether to calculate stateful properties. More...

std::vector< FluidStateProperties >	_fsp
	FluidStateProperties data structure. More...

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

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

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

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

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

const unsigned int	_num_phases
	Number of phases. More...

const unsigned int	_num_components
	Number of components. More...

const unsigned int	_num_pf_vars
	Number of PorousFlow variables. More...

MaterialProperty< std::vector< Real > > &	_porepressure
	Computed nodal or quadpoint values of porepressure of the phases. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dporepressure_dvar
	d(porepressure)/d(PorousFlow variable) More...

MaterialProperty< std::vector< RealGradient > > *const	_gradp_qp
	Grad(p) at the quadpoints. More...

MaterialProperty< std::vector< std::vector< Real > > > *const	_dgradp_qp_dgradv
	d(grad porepressure)/d(grad PorousFlow variable) at the quadpoints More...

MaterialProperty< std::vector< std::vector< RealGradient > > > *const	_dgradp_qp_dv
	d(grad porepressure)/d(PorousFlow variable) at the quadpoints More...

MaterialProperty< std::vector< Real > > &	_saturation
	Computed nodal or qp saturation of the phases. More...

MaterialProperty< std::vector< std::vector< Real > > > &	_dsaturation_dvar
	d(saturation)/d(PorousFlow variable) More...

MaterialProperty< std::vector< RealGradient > > *const	_grads_qp
	Grad(s) at the quadpoints. More...

MaterialProperty< std::vector< std::vector< Real > > > *const	_dgrads_qp_dgradv
	d(grad saturation)/d(grad PorousFlow variable) at the quadpoints More...

MaterialProperty< std::vector< std::vector< RealGradient > > > *const	_dgrads_qp_dv
	d(grad saturation)/d(PorousFlow variable) at the quadpoints More...

Detailed Description

Fluid state class for brine and CO2.

Includes mutual solubility of CO2 and brine using model of Spycher, Pruess and Ennis-King, CO2-H2O mixtures in the geological sequestration of CO2. I. Assessment and calculation of mutual solubilities from 12 to 100C and up to 600 bar, Geochimica et Cosmochimica Acta, 67, 3015-3031 (2003), and Spycher and Pruess, CO2-H2O mixtures in the geological sequestration of CO2. II. Partitioning in chloride brine at 12-100C and up to 600 bar, Geochimica et Cosmochimica Acta, 69, 3309-3320 (2005)

Definition at line 29 of file PorousFlowFluidStateBrineCO2.h.

Constructor & Destructor Documentation

◆ PorousFlowFluidStateBrineCO2()

PorousFlowFluidStateBrineCO2::PorousFlowFluidStateBrineCO2 ( const InputParameters & parameters )

Definition at line 25 of file PorousFlowFluidStateBrineCO2.C.

   : PorousFlowFluidState(parameters)
 {
   mooseDeprecated(name(),
                   ": the PorousFlowFluidStateBrineCO2 material is deprecated. Just use "
                   "PorousFlowFluidState instead");
 }

Member Function Documentation

◆ computeQpProperties()

void PorousFlowFluidState::computeQpProperties ( )

overrideprotectedvirtualinherited

Reimplemented from PorousFlowVariableBase.

Definition at line 195 of file PorousFlowFluidState.C.

 {
   _is_initqp = false;
   // Prepare the derivative vectors
   PorousFlowVariableBase::computeQpProperties();
  
   // Set the size of all other vectors
   setMaterialVectorSize();
  
   // Calculate all required thermophysical properties
   thermophysicalProperties();
  
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     _saturation[_qp][ph] = _fsp[ph].saturation.value();
     _porepressure[_qp][ph] = _fsp[ph].pressure.value();
     _fluid_density[_qp][ph] = _fsp[ph].density.value();
     _fluid_viscosity[_qp][ph] = _fsp[ph].viscosity.value();
     _fluid_enthalpy[_qp][ph] = _fsp[ph].enthalpy.value();
     _fluid_internal_energy[_qp][ph] = _fsp[ph].internal_energy.value();
  
     for (unsigned int comp = 0; comp < _num_components; ++comp)
       _mass_frac[_qp][ph][comp] = _fsp[ph].mass_fraction[comp].value();
   }
  
   // Derivative of properties wrt variables (calculated in fluid state class)
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     // If porepressure is a PorousFlow variable (it usually is), add derivatives wrt porepressure
     if (_dictator.isPorousFlowVariable(_gas_porepressure_varnum))
     {
       _dporepressure_dvar[_qp][ph][_pvar] = _fsp[ph].pressure.derivatives()[_pidx];
       _dsaturation_dvar[_qp][ph][_pvar] = _fsp[ph].saturation.derivatives()[_pidx];
       _dfluid_density_dvar[_qp][ph][_pvar] = _fsp[ph].density.derivatives()[_pidx];
       _dfluid_viscosity_dvar[_qp][ph][_pvar] = _fsp[ph].viscosity.derivatives()[_pidx];
       _dfluid_enthalpy_dvar[_qp][ph][_pvar] = _fsp[ph].enthalpy.derivatives()[_pidx];
       _dfluid_internal_energy_dvar[_qp][ph][_pvar] = _fsp[ph].internal_energy.derivatives()[_pidx];
  
       for (unsigned int comp = 0; comp < _num_components; ++comp)
         _dmass_frac_dvar[_qp][ph][comp][_pvar] = _fsp[ph].mass_fraction[comp].derivatives()[_pidx];
     }
  
     // If Z is a PorousFlow variable (it usually is), add derivatives wrt Z
     if (_dictator.isPorousFlowVariable(_Z_varnum[0]))
     {
       _dporepressure_dvar[_qp][ph][_Zvar[0]] = _fsp[ph].pressure.derivatives()[_Zidx];
       _dsaturation_dvar[_qp][ph][_Zvar[0]] = _fsp[ph].saturation.derivatives()[_Zidx];
       _dfluid_density_dvar[_qp][ph][_Zvar[0]] = _fsp[ph].density.derivatives()[_Zidx];
       _dfluid_viscosity_dvar[_qp][ph][_Zvar[0]] = _fsp[ph].viscosity.derivatives()[_Zidx];
       _dfluid_enthalpy_dvar[_qp][ph][_Zvar[0]] = _fsp[ph].enthalpy.derivatives()[_Zidx];
       _dfluid_internal_energy_dvar[_qp][ph][_Zvar[0]] =
           _fsp[ph].internal_energy.derivatives()[_Zidx];
  
       for (unsigned int comp = 0; comp < _num_components; ++comp)
         _dmass_frac_dvar[_qp][ph][comp][_Zvar[0]] =
             _fsp[ph].mass_fraction[comp].derivatives()[_Zidx];
     }
  
     // If temperature is a PorousFlow variable (nonisothermal case), add derivatives wrt temperature
     if (_dictator.isPorousFlowVariable(_temperature_varnum))
     {
       _dporepressure_dvar[_qp][ph][_Tvar] = _fsp[ph].pressure.derivatives()[_Tidx];
       _dsaturation_dvar[_qp][ph][_Tvar] = _fsp[ph].saturation.derivatives()[_Tidx];
       _dfluid_density_dvar[_qp][ph][_Tvar] = _fsp[ph].density.derivatives()[_Tidx];
       _dfluid_viscosity_dvar[_qp][ph][_Tvar] = _fsp[ph].viscosity.derivatives()[_Tidx];
       _dfluid_enthalpy_dvar[_qp][ph][_Tvar] = _fsp[ph].enthalpy.derivatives()[_Tidx];
       _dfluid_internal_energy_dvar[_qp][ph][_Tvar] = _fsp[ph].internal_energy.derivatives()[_Tidx];
  
       for (unsigned int comp = 0; comp < _num_components; ++comp)
         _dmass_frac_dvar[_qp][ph][comp][_Tvar] = _fsp[ph].mass_fraction[comp].derivatives()[_Tidx];
     }
  
     // If Xnacl is a PorousFlow variable, add derivatives wrt Xnacl
     if (_dictator.isPorousFlowVariable(_Xnacl_varnum))
     {
       _dporepressure_dvar[_qp][ph][_Xvar] = _fsp[ph].pressure.derivatives()[_Xidx];
       _dsaturation_dvar[_qp][ph][_Xvar] = _fsp[ph].saturation.derivatives()[_Xidx];
       _dfluid_density_dvar[_qp][ph][_Xvar] += _fsp[ph].density.derivatives()[_Xidx];
       _dfluid_viscosity_dvar[_qp][ph][_Xvar] += _fsp[ph].viscosity.derivatives()[_Xidx];
       _dfluid_enthalpy_dvar[_qp][ph][_Xvar] = _fsp[ph].enthalpy.derivatives()[_Xidx];
       _dfluid_internal_energy_dvar[_qp][ph][_Xvar] = _fsp[ph].internal_energy.derivatives()[_Xidx];
  
       for (unsigned int comp = 0; comp < _num_components; ++comp)
         _dmass_frac_dvar[_qp][ph][comp][_Xvar] = _fsp[ph].mass_fraction[comp].derivatives()[_Xidx];
     }
   }
  
   // If the material properties are being evaluated at the qps, calculate the gradients as well
   // Note: only nodal properties are evaluated in initQpStatefulProperties(), so no need to check
   // _is_initqp flag for qp properties
   if (!_nodal_material)
   {
     // Derivatives of capillary pressure
     const Real dpc = _pc.dCapillaryPressure(_fsp[_aqueous_phase_number].saturation.value(), _qp);
     const Real d2pc = _pc.d2CapillaryPressure(_fsp[_aqueous_phase_number].saturation.value(), _qp);
  
     // Gradients of saturation and porepressure in all phases
     (*_grads_qp)[_qp][_gas_phase_number] =
         _dsaturation_dvar[_qp][_gas_phase_number][_pvar] * _gas_gradp_qp[_qp] +
         _dsaturation_dvar[_qp][_gas_phase_number][_Zvar[0]] * (*_gradZ_qp[0])[_qp] +
         _dsaturation_dvar[_qp][_gas_phase_number][_Tvar] * _gradT_qp[_qp];
     (*_grads_qp)[_qp][_aqueous_phase_number] = -(*_grads_qp)[_qp][_gas_phase_number];
  
     (*_gradp_qp)[_qp][_gas_phase_number] = _gas_gradp_qp[_qp];
     (*_gradp_qp)[_qp][_aqueous_phase_number] =
         _gas_gradp_qp[_qp] - dpc * (*_grads_qp)[_qp][_aqueous_phase_number];
  
     // Gradients of mass fractions for each component in each phase
     (*_grad_mass_frac_qp)[_qp][_aqueous_phase_number][_aqueous_fluid_component] =
         _fsp[_aqueous_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_pidx] *
             _gas_gradp_qp[_qp] +
         _fsp[_aqueous_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Zidx] *
             (*_gradZ_qp[0])[_qp] +
         _fsp[_aqueous_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Tidx] *
             _gradT_qp[_qp];
     (*_grad_mass_frac_qp)[_qp][_aqueous_phase_number][_gas_fluid_component] =
         -(*_grad_mass_frac_qp)[_qp][_aqueous_phase_number][_aqueous_fluid_component];
  
     (*_grad_mass_frac_qp)[_qp][_gas_phase_number][_aqueous_fluid_component] =
         _fsp[_gas_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_pidx] *
             _gas_gradp_qp[_qp] +
         _fsp[_gas_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Zidx] *
             (*_gradZ_qp[0])[_qp] +
         _fsp[_gas_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Tidx] *
             _gradT_qp[_qp];
     (*_grad_mass_frac_qp)[_qp][_gas_phase_number][_gas_fluid_component] =
         -(*_grad_mass_frac_qp)[_qp][_gas_phase_number][_aqueous_fluid_component];
  
     // Derivatives of gradients wrt variables
     if (_dictator.isPorousFlowVariable(_gas_porepressure_varnum))
     {
       for (unsigned int ph = 0; ph < _num_phases; ++ph)
         (*_dgradp_qp_dgradv)[_qp][ph][_pvar] = 1.0;
  
       (*_dgradp_qp_dgradv)[_qp][_aqueous_phase_number][_pvar] +=
           -dpc * _dsaturation_dvar[_qp][_aqueous_phase_number][_pvar];
  
       (*_dgradp_qp_dv)[_qp][_aqueous_phase_number][_pvar] =
           -d2pc * (*_grads_qp)[_qp][_aqueous_phase_number] *
           _dsaturation_dvar[_qp][_aqueous_phase_number][_pvar];
     }
  
     if (_dictator.isPorousFlowVariable(_Z_varnum[0]))
     {
       (*_dgradp_qp_dgradv)[_qp][_aqueous_phase_number][_Zvar[0]] =
           -dpc * _dsaturation_dvar[_qp][_aqueous_phase_number][_Zvar[0]];
  
       (*_dgradp_qp_dv)[_qp][_aqueous_phase_number][_Zvar[0]] =
           -d2pc * (*_grads_qp)[_qp][_aqueous_phase_number] *
           _dsaturation_dvar[_qp][_aqueous_phase_number][_Zvar[0]];
     }
  
     if (_dictator.isPorousFlowVariable(_temperature_varnum))
     {
       (*_dgradp_qp_dgradv)[_qp][_aqueous_phase_number][_Tvar] =
           -dpc * _dsaturation_dvar[_qp][_aqueous_phase_number][_Tvar];
  
       (*_dgradp_qp_dv)[_qp][_aqueous_phase_number][_Tvar] =
           -d2pc * (*_grads_qp)[_qp][_aqueous_phase_number] *
           _dsaturation_dvar[_qp][_aqueous_phase_number][_Tvar];
     }
  
     // If Xnacl is a PorousFlow variable, add gradients and derivatives wrt Xnacl
     if (_dictator.isPorousFlowVariable(_Xnacl_varnum))
     {
       (*_dgradp_qp_dgradv)[_qp][_aqueous_phase_number][_Xvar] =
           -dpc * _dsaturation_dvar[_qp][_aqueous_phase_number][_Xvar];
  
       (*_grads_qp)[_qp][_aqueous_phase_number] +=
           _dsaturation_dvar[_qp][_aqueous_phase_number][_Xvar] * _grad_Xnacl_qp[_qp];
  
       (*_grads_qp)[_qp][_gas_phase_number] -=
           _dsaturation_dvar[_qp][_aqueous_phase_number][_Xvar] * _grad_Xnacl_qp[_qp];
  
       (*_gradp_qp)[_qp][_aqueous_phase_number] -=
           dpc * _dsaturation_dvar[_qp][_aqueous_phase_number][_Xvar] * _grad_Xnacl_qp[_qp];
  
       (*_dgradp_qp_dv)[_qp][_aqueous_phase_number][_Xvar] =
           -d2pc * (*_grads_qp)[_qp][_aqueous_phase_number] *
           _dsaturation_dvar[_qp][_aqueous_phase_number][_Xvar];
  
       (*_grad_mass_frac_qp)[_qp][_aqueous_phase_number][_salt_component] = _grad_Xnacl_qp[_qp];
       (*_grad_mass_frac_qp)[_qp][_aqueous_phase_number][_aqueous_fluid_component] +=
           _fsp[_aqueous_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Xidx] *
           _grad_Xnacl_qp[_qp];
       (*_grad_mass_frac_qp)[_qp][_aqueous_phase_number][_gas_fluid_component] -=
           _fsp[_aqueous_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Xidx] *
           _grad_Xnacl_qp[_qp];
       (*_grad_mass_frac_qp)[_qp][_gas_phase_number][_aqueous_fluid_component] +=
           _fsp[_gas_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Xidx] *
           _grad_Xnacl_qp[_qp];
       (*_grad_mass_frac_qp)[_qp][_gas_phase_number][_gas_fluid_component] -=
           _fsp[_gas_phase_number].mass_fraction[_aqueous_fluid_component].derivatives()[_Xidx] *
           _grad_Xnacl_qp[_qp];
     }
   }
 }

◆ initQpStatefulProperties()

void PorousFlowFluidState::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented from PorousFlowVariableBase.

Definition at line 166 of file PorousFlowFluidState.C.

 {
   _is_initqp = true;
   // Set the size of pressure and saturation vectors
   PorousFlowVariableBase::initQpStatefulProperties();
  
   // Set the size of all other vectors
   setMaterialVectorSize();
  
   thermophysicalProperties();
  
   // Set the initial values of the properties at the nodes.
   // Note: not required for qp materials as no old values at the qps are requested
   if (_nodal_material)
     for (unsigned int ph = 0; ph < _num_phases; ++ph)
     {
       _saturation[_qp][ph] = _fsp[ph].saturation.value();
       _porepressure[_qp][ph] = _fsp[ph].pressure.value();
       _fluid_density[_qp][ph] = _fsp[ph].density.value();
       _fluid_viscosity[_qp][ph] = _fsp[ph].viscosity.value();
       _fluid_enthalpy[_qp][ph] = _fsp[ph].enthalpy.value();
       _fluid_internal_energy[_qp][ph] = _fsp[ph].internal_energy.value();
  
       for (unsigned int comp = 0; comp < _num_components; ++comp)
         _mass_frac[_qp][ph][comp] = _fsp[ph].mass_fraction[comp].value();
     }
 }

◆ setMaterialVectorSize()

void PorousFlowFluidState::setMaterialVectorSize ( ) const

protectedinherited

Size material property vectors and initialise with zeros.

Definition at line 394 of file PorousFlowFluidState.C.

 {
   _fluid_density[_qp].assign(_num_phases, 0.0);
   _fluid_viscosity[_qp].assign(_num_phases, 0.0);
   _fluid_enthalpy[_qp].assign(_num_phases, 0.0);
   _fluid_internal_energy[_qp].assign(_num_phases, 0.0);
   _mass_frac[_qp].resize(_num_phases);
  
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
     _mass_frac[_qp][ph].resize(_num_components);
  
   // Derivatives and gradients are not required in initQpStatefulProperties
   if (!_is_initqp)
   {
     _dfluid_density_dvar[_qp].resize(_num_phases);
     _dfluid_viscosity_dvar[_qp].resize(_num_phases);
     _dfluid_enthalpy_dvar[_qp].resize(_num_phases);
     _dfluid_internal_energy_dvar[_qp].resize(_num_phases);
     _dmass_frac_dvar[_qp].resize(_num_phases);
  
     if (!_nodal_material)
       (*_grad_mass_frac_qp)[_qp].resize(_num_phases);
  
     for (unsigned int ph = 0; ph < _num_phases; ++ph)
     {
       _dfluid_density_dvar[_qp][ph].assign(_num_pf_vars, 0.0);
       _dfluid_viscosity_dvar[_qp][ph].assign(_num_pf_vars, 0.0);
       _dfluid_enthalpy_dvar[_qp][ph].assign(_num_pf_vars, 0.0);
       _dfluid_internal_energy_dvar[_qp][ph].assign(_num_pf_vars, 0.0);
       _dmass_frac_dvar[_qp][ph].resize(_num_components);
  
       for (unsigned int comp = 0; comp < _num_components; ++comp)
         _dmass_frac_dvar[_qp][ph][comp].assign(_num_pf_vars, 0.0);
  
       if (!_nodal_material)
         (*_grad_mass_frac_qp)[_qp][ph].assign(_num_components, RealGradient());
     }
   }
 }

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ thermophysicalProperties()

void PorousFlowFluidState::thermophysicalProperties ( )

protectedvirtualinherited

Calculates all required thermophysical properties and derivatives for each phase and fluid component.

Definition at line 157 of file PorousFlowFluidState.C.

 {
   // The FluidProperty objects use temperature in K
   Real Tk = _temperature[_qp] + _T_c2k;
  
   _fs.thermophysicalProperties(_gas_porepressure[_qp], Tk, _Xnacl[_qp], (*_Z[0])[_qp], _qp, _fsp);
 }

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::initQpStatefulProperties().

Member Data Documentation

◆ _aqueous_fluid_component

const unsigned int PorousFlowFluidState::_aqueous_fluid_component

protectedinherited

Fluid component number of the aqueous component.

Definition at line 101 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _aqueous_phase_number

const unsigned int PorousFlowFluidState::_aqueous_phase_number

protectedinherited

Phase number of the aqueous phase.

Definition at line 97 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _dfluid_density_dvar

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFluidState::_dfluid_density_dvar

protectedinherited

Derivative of the fluid density for each phase wrt PorousFlow variables.

Definition at line 126 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _dfluid_enthalpy_dvar

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFluidState::_dfluid_enthalpy_dvar

protectedinherited

Derivative of the fluid enthalpy for each phase wrt PorousFlow variables.

Definition at line 134 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _dfluid_internal_energy_dvar

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFluidState::_dfluid_internal_energy_dvar

protectedinherited

Derivative of the fluid internal energy for each phase wrt PorousFlow variables.

Definition at line 138 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _dfluid_viscosity_dvar

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFluidState::_dfluid_viscosity_dvar

protectedinherited

Derivative of the fluid viscosity for each phase wrt PorousFlow variables.

Definition at line 130 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _dgradp_qp_dgradv

MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowVariableBase::_dgradp_qp_dgradv

protectedinherited

d(grad porepressure)/d(grad PorousFlow variable) at the quadpoints

Definition at line 53 of file PorousFlowVariableBase.h.

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

◆ _dgradp_qp_dv

MaterialProperty<std::vector<std::vector<RealGradient> > >* const PorousFlowVariableBase::_dgradp_qp_dv

protectedinherited

d(grad porepressure)/d(PorousFlow variable) at the quadpoints

Definition at line 56 of file PorousFlowVariableBase.h.

◆ _dgrads_qp_dgradv

MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowVariableBase::_dgrads_qp_dgradv

protectedinherited

d(grad saturation)/d(grad PorousFlow variable) at the quadpoints

Definition at line 68 of file PorousFlowVariableBase.h.

◆ _dgrads_qp_dv

MaterialProperty<std::vector<std::vector<RealGradient> > >* const PorousFlowVariableBase::_dgrads_qp_dv

protectedinherited

d(grad saturation)/d(PorousFlow variable) at the quadpoints

Definition at line 71 of file PorousFlowVariableBase.h.

◆ _dmass_frac_dvar

MaterialProperty<std::vector<std::vector<std::vector<Real> > > >& PorousFlowFluidState::_dmass_frac_dvar

protectedinherited

Derivative of the mass fraction matrix with respect to the Porous Flow variables.

Definition at line 121 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _dporepressure_dvar

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowVariableBase::_dporepressure_dvar

protectedinherited

d(porepressure)/d(PorousFlow variable)

Definition at line 47 of file PorousFlowVariableBase.h.

Referenced by PorousFlow1PhaseFullySaturated::computeQpProperties(), PorousFlow2PhasePP::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlow1PhaseP::computeQpProperties(), PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), and PorousFlowFluidState::computeQpProperties().

◆ _dsaturation_dvar

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowVariableBase::_dsaturation_dvar

protectedinherited

d(saturation)/d(PorousFlow variable)

Definition at line 62 of file PorousFlowVariableBase.h.

Referenced by PorousFlow2PhasePP::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlow1PhaseP::computeQpProperties(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), and PorousFlowFluidState::computeQpProperties().

◆ _dtemperature_dvar

const MaterialProperty<std::vector<Real> >& PorousFlowFluidState::_dtemperature_dvar

protectedinherited

Derivative of temperature wrt PorousFlow variables.

Definition at line 111 of file PorousFlowFluidState.h.

◆ _fluid_density

MaterialProperty<std::vector<Real> >& PorousFlowFluidState::_fluid_density

protectedinherited

Fluid density of each phase.

Definition at line 124 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _fluid_enthalpy

MaterialProperty<std::vector<Real> >& PorousFlowFluidState::_fluid_enthalpy

protectedinherited

Enthalpy of each phase.

Definition at line 132 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _fluid_internal_energy

MaterialProperty<std::vector<Real> >& PorousFlowFluidState::_fluid_internal_energy

protectedinherited

Internal energy of each phase.

Definition at line 136 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _fluid_viscosity

MaterialProperty<std::vector<Real> >& PorousFlowFluidState::_fluid_viscosity

protectedinherited

Viscosity of each phase.

Definition at line 128 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _fs

const PorousFlowFluidStateMultiComponentBase& PorousFlowFluidState::_fs

protectedinherited

FluidState UserObject.

Definition at line 95 of file PorousFlowFluidState.h.

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

◆ _fsp

std::vector<FluidStateProperties> PorousFlowFluidState::_fsp

protectedinherited

FluidStateProperties data structure.

Definition at line 145 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), PorousFlowFluidState::PorousFlowFluidState(), and PorousFlowFluidState::thermophysicalProperties().

◆ _gas_fluid_component

const unsigned int PorousFlowFluidState::_gas_fluid_component

protectedinherited

Fluid component number of the gas phase.

Definition at line 103 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _gas_gradp_qp

const VariableGradient& PorousFlowFluidState::_gas_gradp_qp

protectedinherited

Gradient of porepressure (only defined at the qps)

Definition at line 71 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _gas_phase_number

const unsigned int PorousFlowFluidState::_gas_phase_number

protectedinherited

Phase number of the gas phase.

Definition at line 99 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _gas_porepressure

const VariableValue& PorousFlowFluidState::_gas_porepressure

protectedinherited

Porepressure.

Definition at line 69 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::thermophysicalProperties().

◆ _gas_porepressure_varnum

const unsigned int PorousFlowFluidState::_gas_porepressure_varnum

protectedinherited

Moose variable number of the gas porepressure.

Definition at line 73 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _grad_mass_frac_qp

MaterialProperty<std::vector<std::vector<RealGradient> > >* PorousFlowFluidState::_grad_mass_frac_qp

protectedinherited

Gradient of the mass fraction matrix (only defined at the qps)

Definition at line 119 of file PorousFlowFluidState.h.

◆ _grad_Xnacl_qp

const VariableGradient& PorousFlowFluidState::_grad_Xnacl_qp

protectedinherited

Gradient of salt mass fraction (only defined at the qps)

Definition at line 89 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _gradp_qp

MaterialProperty<std::vector<RealGradient> >* const PorousFlowVariableBase::_gradp_qp

protectedinherited

Grad(p) at the quadpoints.

Definition at line 50 of file PorousFlowVariableBase.h.

◆ _grads_qp

MaterialProperty<std::vector<RealGradient> >* const PorousFlowVariableBase::_grads_qp

protectedinherited

Grad(s) at the quadpoints.

Definition at line 65 of file PorousFlowVariableBase.h.

◆ _gradT_qp

const MaterialProperty<RealGradient>& PorousFlowFluidState::_gradT_qp

protectedinherited

Gradient of temperature (only defined at the qps)

Definition at line 109 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _gradZ_qp

std::vector<const VariableGradient *> PorousFlowFluidState::_gradZ_qp

protectedinherited

Gradient(s) of total mass fraction(s) of the gas component(s) (only defined at the qps)

Definition at line 79 of file PorousFlowFluidState.h.

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

◆ _is_initqp

bool PorousFlowFluidState::_is_initqp

protectedinherited

Flag to indicate whether to calculate stateful properties.

Definition at line 143 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _mass_frac

MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFluidState::_mass_frac

protectedinherited

Mass fraction matrix.

Definition at line 117 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _num_components

const unsigned int PorousFlowVariableBase::_num_components

protectedinherited

Number of components.

Definition at line 38 of file PorousFlowVariableBase.h.

Referenced by PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), PorousFlowFluidState::PorousFlowFluidState(), PorousFlowFluidStateSingleComponent::PorousFlowFluidStateSingleComponent(), PorousFlowFluidStateSingleComponent::setMaterialVectorSize(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _num_pf_vars

const unsigned int PorousFlowVariableBase::_num_pf_vars

protectedinherited

Number of PorousFlow variables.

Definition at line 41 of file PorousFlowVariableBase.h.

Referenced by PorousFlowVariableBase::computeQpProperties(), PorousFlowFluidStateSingleComponent::setMaterialVectorSize(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _num_phases

const unsigned int PorousFlowVariableBase::_num_phases

protectedinherited

Number of phases.

Definition at line 35 of file PorousFlowVariableBase.h.

Referenced by PorousFlowVariableBase::computeQpProperties(), PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), PorousFlowFluidStateSingleComponent::initQpStatefulProperties(), PorousFlowFluidState::initQpStatefulProperties(), PorousFlow1PhaseFullySaturated::PorousFlow1PhaseFullySaturated(), PorousFlow1PhaseMD_Gaussian::PorousFlow1PhaseMD_Gaussian(), PorousFlow1PhaseP::PorousFlow1PhaseP(), PorousFlow2PhasePP::PorousFlow2PhasePP(), PorousFlowFluidState::PorousFlowFluidState(), PorousFlowFluidStateSingleComponent::PorousFlowFluidStateSingleComponent(), PorousFlowFluidStateSingleComponent::setMaterialVectorSize(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _num_Z_vars

const unsigned int PorousFlowFluidState::_num_Z_vars

protectedinherited

Number of coupled total mass fractions. Should be _num_phases - 1.

Definition at line 85 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::PorousFlowFluidState().

◆ _pc

const PorousFlowCapillaryPressure& PorousFlowFluidState::_pc

protectedinherited

Capillary pressure UserObject.

Definition at line 147 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _pidx

const unsigned int PorousFlowFluidState::_pidx

protectedinherited

Index of derivative wrt pressure.

Definition at line 149 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _porepressure

MaterialProperty<std::vector<Real> >& PorousFlowVariableBase::_porepressure

protectedinherited

Computed nodal or quadpoint values of porepressure of the phases.

Definition at line 44 of file PorousFlowVariableBase.h.

Referenced by PorousFlow1PhaseMD_Gaussian::buildPS(), PorousFlow2PhasePS::buildQpPPSS(), PorousFlow1PhaseFullySaturated::buildQpPPSS(), PorousFlow1PhaseP::buildQpPPSS(), PorousFlow2PhasePP::buildQpPPSS(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), PorousFlowFluidStateSingleComponent::initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ _pvar

const unsigned int PorousFlowFluidState::_pvar

protectedinherited

PorousFlow variable number of the gas porepressure.

Definition at line 75 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _salt_component

const unsigned int PorousFlowFluidState::_salt_component

protectedinherited

Salt component index.

Definition at line 105 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _saturation

MaterialProperty<std::vector<Real> >& PorousFlowVariableBase::_saturation

protectedinherited

Computed nodal or qp saturation of the phases.

Definition at line 59 of file PorousFlowVariableBase.h.

Referenced by PorousFlow1PhaseMD_Gaussian::buildPS(), PorousFlow2PhasePS::buildQpPPSS(), PorousFlow1PhaseFullySaturated::buildQpPPSS(), PorousFlow1PhaseP::buildQpPPSS(), PorousFlow2PhasePP::buildQpPPSS(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), PorousFlowFluidStateSingleComponent::initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ _T_c2k

const Real PorousFlowFluidState::_T_c2k

protectedinherited

Conversion from degrees Celsius to degrees Kelvin.

Definition at line 141 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::thermophysicalProperties().

◆ _temperature

const MaterialProperty<Real>& PorousFlowFluidState::_temperature

protectedinherited

Temperature.

Definition at line 107 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::thermophysicalProperties().

◆ _temperature_varnum

const unsigned int PorousFlowFluidState::_temperature_varnum

protectedinherited

Moose variable number of the temperature.

Definition at line 113 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Tidx

const unsigned int PorousFlowFluidState::_Tidx

protectedinherited

Index of derivative wrt temperature.

Definition at line 151 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Tvar

const unsigned int PorousFlowFluidState::_Tvar

protectedinherited

PorousFlow variable number of the temperature.

Definition at line 115 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Xidx

const unsigned int PorousFlowFluidState::_Xidx

protectedinherited

Index of derivative wrt salt mass fraction X.

Definition at line 155 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Xnacl

const VariableValue& PorousFlowFluidState::_Xnacl

protectedinherited

Salt mass fraction (kg/kg)

Definition at line 87 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::thermophysicalProperties().

◆ _Xnacl_varnum

const unsigned int PorousFlowFluidState::_Xnacl_varnum

protectedinherited

Salt mass fraction variable number.

Definition at line 91 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Xvar

const unsigned int PorousFlowFluidState::_Xvar

protectedinherited

Salt mass fraction PorousFlow variable number.

Definition at line 93 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Z

std::vector<const VariableValue *> PorousFlowFluidState::_Z

protectedinherited

Total mass fraction(s) of the gas component(s) summed over all phases.

Definition at line 77 of file PorousFlowFluidState.h.

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

◆ _Z_varnum

std::vector<unsigned int> PorousFlowFluidState::_Z_varnum

protectedinherited

Moose variable number of Z.

Definition at line 81 of file PorousFlowFluidState.h.

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

◆ _Zidx

const unsigned int PorousFlowFluidState::_Zidx

protectedinherited

Index of derivative wrt total mass fraction Z.

Definition at line 153 of file PorousFlowFluidState.h.

Referenced by PorousFlowFluidState::computeQpProperties().

◆ _Zvar

std::vector<unsigned int> PorousFlowFluidState::_Zvar

protectedinherited

PorousFlow variable number of Z.

Definition at line 83 of file PorousFlowFluidState.h.

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

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

porous_flow/include/materials/PorousFlowFluidStateBrineCO2.h
porous_flow/src/materials/PorousFlowFluidStateBrineCO2.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowFluidStateBrineCO2()

Member Function Documentation

◆ computeQpProperties()

◆ initQpStatefulProperties()

◆ setMaterialVectorSize()

◆ thermophysicalProperties()

Member Data Documentation

◆ _aqueous_fluid_component

◆ _aqueous_phase_number

◆ _dfluid_density_dvar

◆ _dfluid_enthalpy_dvar

◆ _dfluid_internal_energy_dvar

◆ _dfluid_viscosity_dvar

◆ _dgradp_qp_dgradv

◆ _dgradp_qp_dv

◆ _dgrads_qp_dgradv

◆ _dgrads_qp_dv

◆ _dmass_frac_dvar

◆ _dporepressure_dvar

◆ _dsaturation_dvar

◆ _dtemperature_dvar

◆ _fluid_density

◆ _fluid_enthalpy

◆ _fluid_internal_energy

◆ _fluid_viscosity

◆ _fs

◆ _fsp

◆ _gas_fluid_component

◆ _gas_gradp_qp

◆ _gas_phase_number

◆ _gas_porepressure

◆ _gas_porepressure_varnum

◆ _grad_mass_frac_qp

◆ _grad_Xnacl_qp

◆ _gradp_qp

◆ _grads_qp

◆ _gradT_qp

◆ _gradZ_qp

◆ _is_initqp

◆ _mass_frac

◆ _num_components

◆ _num_pf_vars

◆ _num_phases

◆ _num_Z_vars

◆ _pc

◆ _pidx

◆ _porepressure

◆ _pvar

◆ _salt_component

◆ _saturation

◆ _T_c2k

◆ _temperature

◆ _temperature_varnum

◆ _Tidx

◆ _Tvar

◆ _Xidx

◆ _Xnacl

◆ _Xnacl_varnum

◆ _Xvar

◆ _Z

◆ _Z_varnum

◆ _Zidx

◆ _Zvar