Fluid state class using a persistent set of primary variables for the mutliphase, single component case. More...

#include <PorousFlowFluidStateSingleComponent.h>

Inheritance diagram for PorousFlowFluidStateSingleComponent:

Public Member Functions
	PorousFlowFluidStateSingleComponent (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 &	_liquid_porepressure
	Porepressure. More...

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

const unsigned int	_liquid_porepressure_varnum
	Moose variable number of the porepressure. More...

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

const VariableValue &	_enthalpy
	Enthalpy. More...

const VariableGradient &	_gradh_qp
	Gradient of enthalpy (only defined at the qps) More...

const unsigned int	_enthalpy_varnum
	Moose variable number of the enthalpy. More...

const unsigned int	_hvar
	PorousFlow variable number of the enthalpy. More...

const PorousFlowFluidStateSingleComponentBase &	_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...

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

MaterialProperty< RealGradient > *	_grad_temperature_qp
	Gradient of temperature (only defined at the qps) More...

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

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

MaterialProperty< std::vector< RealGradient > > *const	_dgrad_temperature_dv
	d(grad temperature)/d(PorousFlow variable) at the quadpoints 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...

FluidStatePhaseEnum	_phase_state
	FluidStatePhaseEnum. More...

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

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

const unsigned int	_hidx
	Index of derivative wrt enthalpy. 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 using a persistent set of primary variables for the mutliphase, single component case.

Primary variables are: liquid pressure and enthalpy.

The PorousFlow kernels expect saturation and mass fractions (as well as pressure and temperature), so these must be calculated from the primary variables once the state of the system is determined.

Definition at line 31 of file PorousFlowFluidStateSingleComponent.h.

Constructor & Destructor Documentation

◆ PorousFlowFluidStateSingleComponent()

PorousFlowFluidStateSingleComponent::PorousFlowFluidStateSingleComponent ( const InputParameters & parameters )

Definition at line 35 of file PorousFlowFluidStateSingleComponent.C.

   : PorousFlowVariableBase(parameters),
  
     _liquid_porepressure(_nodal_material ? coupledDofValues("porepressure")
                                          : coupledValue("porepressure")),
     _liquid_gradp_qp(coupledGradient("porepressure")),
     _liquid_porepressure_varnum(coupled("porepressure")),
     _pvar(_dictator.isPorousFlowVariable(_liquid_porepressure_varnum)
               ? _dictator.porousFlowVariableNum(_liquid_porepressure_varnum)
               : 0),
  
     _enthalpy(_nodal_material ? coupledDofValues("enthalpy") : coupledValue("enthalpy")),
     _gradh_qp(coupledGradient("enthalpy")),
     _enthalpy_varnum(coupled("enthalpy")),
     _hvar(_dictator.isPorousFlowVariable(_enthalpy_varnum)
               ? _dictator.porousFlowVariableNum(_enthalpy_varnum)
               : 0),
  
     _fs(getUserObject<PorousFlowFluidStateSingleComponentBase>("fluid_state")),
     _aqueous_phase_number(_fs.aqueousPhaseIndex()),
     _gas_phase_number(_fs.gasPhaseIndex()),
  
     _temperature(_nodal_material ? declareProperty<Real>("PorousFlow_temperature_nodal")
                                  : declareProperty<Real>("PorousFlow_temperature_qp")),
     _grad_temperature_qp(_nodal_material
                              ? nullptr
                              : &declareProperty<RealGradient>("PorousFlow_grad_temperature_qp")),
     _dtemperature_dvar(
         _nodal_material ? declareProperty<std::vector<Real>>("dPorousFlow_temperature_nodal_dvar")
                         : declareProperty<std::vector<Real>>("dPorousFlow_temperature_qp_dvar")),
     _dgrad_temperature_dgradv(_nodal_material ? nullptr
                                               : &declareProperty<std::vector<Real>>(
                                                     "dPorousFlow_grad_temperature_qp_dgradvar")),
     _dgrad_temperature_dv(_nodal_material ? nullptr
                                           : &declareProperty<std::vector<RealGradient>>(
                                                 "dPorousFlow_grad_temperature_qp_dvar")),
     _mass_frac(_nodal_material
                    ? declareProperty<std::vector<std::vector<Real>>>("PorousFlow_mass_frac_nodal")
                    : declareProperty<std::vector<std::vector<Real>>>("PorousFlow_mass_frac_qp")),
     _grad_mass_frac_qp(_nodal_material ? nullptr
                                        : &declareProperty<std::vector<std::vector<RealGradient>>>(
                                              "PorousFlow_grad_mass_frac_qp")),
     _dmass_frac_dvar(_nodal_material ? declareProperty<std::vector<std::vector<std::vector<Real>>>>(
                                            "dPorousFlow_mass_frac_nodal_dvar")
                                      : declareProperty<std::vector<std::vector<std::vector<Real>>>>(
                                            "dPorousFlow_mass_frac_qp_dvar")),
  
     _fluid_density(_nodal_material
                        ? declareProperty<std::vector<Real>>("PorousFlow_fluid_phase_density_nodal")
                        : declareProperty<std::vector<Real>>("PorousFlow_fluid_phase_density_qp")),
     _dfluid_density_dvar(_nodal_material ? declareProperty<std::vector<std::vector<Real>>>(
                                                "dPorousFlow_fluid_phase_density_nodal_dvar")
                                          : declareProperty<std::vector<std::vector<Real>>>(
                                                "dPorousFlow_fluid_phase_density_qp_dvar")),
     _fluid_viscosity(_nodal_material
                          ? declareProperty<std::vector<Real>>("PorousFlow_viscosity_nodal")
                          : declareProperty<std::vector<Real>>("PorousFlow_viscosity_qp")),
     _dfluid_viscosity_dvar(
         _nodal_material
             ? declareProperty<std::vector<std::vector<Real>>>("dPorousFlow_viscosity_nodal_dvar")
             : declareProperty<std::vector<std::vector<Real>>>("dPorousFlow_viscosity_qp_dvar")),
  
     _fluid_enthalpy(
         _nodal_material
             ? declareProperty<std::vector<Real>>("PorousFlow_fluid_phase_enthalpy_nodal")
             : declareProperty<std::vector<Real>>("PorousFlow_fluid_phase_enthalpy_qp")),
     _dfluid_enthalpy_dvar(_nodal_material ? declareProperty<std::vector<std::vector<Real>>>(
                                                 "dPorousFlow_fluid_phase_enthalpy_nodal_dvar")
                                           : declareProperty<std::vector<std::vector<Real>>>(
                                                 "dPorousFlow_fluid_phase_enthalpy_qp_dvar")),
  
     _fluid_internal_energy(
         _nodal_material
             ? declareProperty<std::vector<Real>>("PorousFlow_fluid_phase_internal_energy_nodal")
             : declareProperty<std::vector<Real>>("PorousFlow_fluid_phase_internal_energy_qp")),
     _dfluid_internal_energy_dvar(_nodal_material
                                      ? declareProperty<std::vector<std::vector<Real>>>(
                                            "dPorousFlow_fluid_phase_internal_energy_nodal_dvar")
                                      : declareProperty<std::vector<std::vector<Real>>>(
                                            "dPorousFlow_fluid_phase_internal_qp_dvar")),
  
     _T_c2k(getParam<MooseEnum>("temperature_unit") == 0 ? 0.0 : 273.15),
     _is_initqp(false),
     _pc(getUserObject<PorousFlowCapillaryPressure>("capillary_pressure")),
     _pidx(_fs.getPressureIndex()),
     _hidx(_fs.getEnthalpyIndex())
 {
   // Check that the number of phases in the fluidstate class is also provided in the Dictator
   if (_fs.numPhases() != _num_phases)
     mooseError(name(),
                ": only ",
                _fs.numPhases(),
                " phases are allowed. Please check the number of phases entered in the dictator is "
                "correct");
  
   // Set the size of the FluidStateProperties vector
   _fsp.resize(_num_phases, FluidStateProperties(_num_components));
 }

Member Function Documentation

◆ computeQpProperties()

void PorousFlowFluidStateSingleComponent::computeQpProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBase.

Definition at line 176 of file PorousFlowFluidStateSingleComponent.C.

 {
   _is_initqp = false;
   // Prepare the derivative vectors
   PorousFlowVariableBase::computeQpProperties();
  
   // Set the size of all other vectors
   setMaterialVectorSize();
  
   // Calculate all required thermophysical properties
   thermophysicalProperties();
  
   // Temperature doesn't depend on fluid phase
   _temperature[_qp] = _fsp[_aqueous_phase_number].temperature.value() - _T_c2k;
   _dtemperature_dvar[_qp][_pvar] = _fsp[_aqueous_phase_number].temperature.derivatives()[_pidx];
   _dtemperature_dvar[_qp][_hvar] = _fsp[_aqueous_phase_number].temperature.derivatives()[_hidx];
  
   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();
   }
  
   // Derivative of pressure, saturation and fluid properties wrt variables
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     _dporepressure_dvar[_qp][ph][_pvar] = _fsp[ph].pressure.derivatives()[_pidx];
     _dporepressure_dvar[_qp][ph][_hvar] = _fsp[ph].pressure.derivatives()[_hidx];
  
     _dsaturation_dvar[_qp][ph][_pvar] = _fsp[ph].saturation.derivatives()[_pidx];
     _dsaturation_dvar[_qp][ph][_hvar] = _fsp[ph].saturation.derivatives()[_hidx];
  
     _dfluid_density_dvar[_qp][ph][_pvar] = _fsp[ph].density.derivatives()[_pidx];
     _dfluid_density_dvar[_qp][ph][_hvar] = _fsp[ph].density.derivatives()[_hidx];
  
     _dfluid_viscosity_dvar[_qp][ph][_pvar] = _fsp[ph].viscosity.derivatives()[_pidx];
     _dfluid_viscosity_dvar[_qp][ph][_hvar] = _fsp[ph].viscosity.derivatives()[_hidx];
  
     _dfluid_enthalpy_dvar[_qp][ph][_pvar] = _fsp[ph].enthalpy.derivatives()[_pidx];
     _dfluid_enthalpy_dvar[_qp][ph][_hvar] = _fsp[ph].enthalpy.derivatives()[_hidx];
  
     _dfluid_internal_energy_dvar[_qp][ph][_pvar] = _fsp[ph].internal_energy.derivatives()[_pidx];
     _dfluid_internal_energy_dvar[_qp][ph][_hvar] = _fsp[ph].internal_energy.derivatives()[_hidx];
   }
  
   // 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)
   {
     // Need to compute second derivatives of properties wrt variables for some of
     // the gradient derivatives. Use finite differences for now
     const Real dp = 1.0e-2;
     const Real dh = 1.0e-2;
  
     std::vector<FluidStateProperties> fsp_dp(_num_phases, FluidStateProperties(_num_components));
     _fs.thermophysicalProperties(
         _liquid_porepressure[_qp] + dp, _enthalpy[_qp], _qp, _phase_state, fsp_dp);
  
     std::vector<FluidStateProperties> fsp_dh(_num_phases, FluidStateProperties(_num_components));
     _fs.thermophysicalProperties(
         _liquid_porepressure[_qp], _enthalpy[_qp] + dh, _qp, _phase_state, fsp_dh);
  
     // Gradient of temperature (non-zero in all phases)
     (*_grad_temperature_qp)[_qp] = _dtemperature_dvar[_qp][_pvar] * _liquid_gradp_qp[_qp] +
                                    _dtemperature_dvar[_qp][_hvar] * _gradh_qp[_qp];
     (*_dgrad_temperature_dgradv)[_qp][_pvar] = _dtemperature_dvar[_qp][_pvar];
     (*_dgrad_temperature_dgradv)[_qp][_hvar] = _dtemperature_dvar[_qp][_hvar];
  
     const Real d2T_dp2 = (fsp_dp[_aqueous_phase_number].temperature.derivatives()[_pidx] -
                           _fsp[_aqueous_phase_number].temperature.derivatives()[_pidx]) /
                          dp;
  
     const Real d2T_dh2 = (fsp_dh[_aqueous_phase_number].temperature.derivatives()[_hidx] -
                           _fsp[_aqueous_phase_number].temperature.derivatives()[_hidx]) /
                          dh;
  
     const Real d2T_dph = (fsp_dp[_aqueous_phase_number].temperature.derivatives()[_hidx] -
                           _fsp[_aqueous_phase_number].temperature.derivatives()[_hidx]) /
                              (2.0 * dp) +
                          (fsp_dh[_aqueous_phase_number].temperature.derivatives()[_pidx] -
                           _fsp[_aqueous_phase_number].temperature.derivatives()[_pidx]) /
                              (2.0 * dh);
  
     (*_dgrad_temperature_dv)[_qp][_pvar] =
         d2T_dp2 * _liquid_gradp_qp[_qp] + d2T_dph * _gradh_qp[_qp];
     (*_dgrad_temperature_dv)[_qp][_hvar] =
         d2T_dph * _liquid_gradp_qp[_qp] + d2T_dh2 * _gradh_qp[_qp];
  
     // Gradient of saturation and derivatives
     (*_grads_qp)[_qp][_gas_phase_number] =
         _dsaturation_dvar[_qp][_gas_phase_number][_pvar] * _liquid_gradp_qp[_qp] +
         _dsaturation_dvar[_qp][_gas_phase_number][_hvar] * _gradh_qp[_qp];
     (*_grads_qp)[_qp][_aqueous_phase_number] = -(*_grads_qp)[_qp][_gas_phase_number];
  
     (*_dgrads_qp_dgradv)[_qp][_gas_phase_number][_pvar] =
         _dsaturation_dvar[_qp][_gas_phase_number][_pvar];
     (*_dgrads_qp_dgradv)[_qp][_aqueous_phase_number][_pvar] =
         -(*_dgrads_qp_dgradv)[_qp][_gas_phase_number][_pvar];
  
     (*_dgrads_qp_dgradv)[_qp][_gas_phase_number][_hvar] =
         _dsaturation_dvar[_qp][_gas_phase_number][_hvar];
     (*_dgrads_qp_dgradv)[_qp][_aqueous_phase_number][_hvar] =
         -(*_dgrads_qp_dgradv)[_qp][_gas_phase_number][_hvar];
  
     const Real d2s_dp2 = (fsp_dp[_gas_phase_number].saturation.derivatives()[_pidx] -
                           _fsp[_gas_phase_number].saturation.derivatives()[_pidx]) /
                          dp;
  
     const Real d2s_dh2 = (fsp_dh[_gas_phase_number].saturation.derivatives()[_hidx] -
                           _fsp[_gas_phase_number].saturation.derivatives()[_hidx]) /
                          dh;
  
     const Real d2s_dph = (fsp_dp[_gas_phase_number].saturation.derivatives()[_hidx] -
                           _fsp[_gas_phase_number].saturation.derivatives()[_hidx]) /
                              (2.0 * dp) +
                          (fsp_dh[_gas_phase_number].saturation.derivatives()[_pidx] -
                           _fsp[_gas_phase_number].saturation.derivatives()[_pidx]) /
                              (2.0 * dh);
  
     (*_dgrads_qp_dv)[_qp][_gas_phase_number][_pvar] =
         d2s_dp2 * _liquid_gradp_qp[_qp] + d2s_dph * _gradh_qp[_qp];
     (*_dgrads_qp_dv)[_qp][_aqueous_phase_number][_pvar] =
         -(*_dgrads_qp_dv)[_qp][_gas_phase_number][_pvar];
  
     (*_dgrads_qp_dv)[_qp][_gas_phase_number][_hvar] =
         d2s_dh2 * _gradh_qp[_qp] + d2s_dph * _liquid_gradp_qp[_qp];
     (*_dgrads_qp_dv)[_qp][_aqueous_phase_number][_hvar] =
         -(*_dgrads_qp_dv)[_qp][_gas_phase_number][_hvar];
  
     // Gradient of porepressure and derivatives
     // Note: need first and second derivativea of capillary pressure
     const Real dpc = _pc.dCapillaryPressure(_fsp[_aqueous_phase_number].saturation.value());
     const Real d2pc = _pc.d2CapillaryPressure(_fsp[_aqueous_phase_number].saturation.value());
  
     (*_gradp_qp)[_qp][_aqueous_phase_number] = _liquid_gradp_qp[_qp];
     (*_gradp_qp)[_qp][_gas_phase_number] =
         _liquid_gradp_qp[_qp] + dpc * (*_grads_qp)[_qp][_aqueous_phase_number];
  
     for (unsigned int ph = 0; ph < _num_phases; ++ph)
       (*_dgradp_qp_dgradv)[_qp][ph][_pvar] = 1.0;
  
     (*_dgradp_qp_dgradv)[_qp][_gas_phase_number][_pvar] +=
         dpc * (*_dgrads_qp_dgradv)[_qp][_aqueous_phase_number][_pvar];
     (*_dgradp_qp_dgradv)[_qp][_gas_phase_number][_hvar] =
         dpc * (*_dgrads_qp_dgradv)[_qp][_aqueous_phase_number][_hvar];
  
     (*_dgradp_qp_dv)[_qp][_gas_phase_number][_pvar] =
         d2pc * (*_grads_qp)[_qp][_aqueous_phase_number] *
             _dsaturation_dvar[_qp][_aqueous_phase_number][_pvar] +
         dpc * (*_dgrads_qp_dv)[_qp][_aqueous_phase_number][_pvar];
  
     (*_dgradp_qp_dv)[_qp][_gas_phase_number][_hvar] =
         d2pc * (*_grads_qp)[_qp][_aqueous_phase_number] *
             _dsaturation_dvar[_qp][_aqueous_phase_number][_hvar] +
         dpc * (*_dgrads_qp_dv)[_qp][_aqueous_phase_number][_hvar];
   }
 }

◆ initQpStatefulProperties()

void PorousFlowFluidStateSingleComponent::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBase.

Definition at line 143 of file PorousFlowFluidStateSingleComponent.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)
   {
     thermophysicalProperties();
  
     // Temperature doesn't depend on fluid phase
     _temperature[_qp] = _fsp[_aqueous_phase_number].temperature.value() - _T_c2k;
  
     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();
     }
   }
 }

◆ setMaterialVectorSize()

void PorousFlowFluidStateSingleComponent::setMaterialVectorSize ( ) const

protected

Size material property vectors and initialise with zeros.

Definition at line 341 of file PorousFlowFluidStateSingleComponent.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].assign(_num_components, 1.0);
  
   // 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);
  
     // Temperature doesn't depend of fluid phase
     _dtemperature_dvar[_qp].assign(_num_pf_vars, 0.0);
  
     if (!_nodal_material)
     {
       (*_grad_mass_frac_qp)[_qp].resize(_num_phases);
       (*_grad_temperature_qp)[_qp] = RealGradient();
       (*_dgrad_temperature_dgradv)[_qp].assign(_num_pf_vars, 0.0);
       (*_dgrad_temperature_dv)[_qp].assign(_num_pf_vars, RealGradient());
     }
  
     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 computeQpProperties(), and initQpStatefulProperties().

◆ thermophysicalProperties()

void PorousFlowFluidStateSingleComponent::thermophysicalProperties ( )

protectedvirtual

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

Definition at line 136 of file PorousFlowFluidStateSingleComponent.C.

 {
   _fs.clearFluidStateProperties(_fsp);
   _fs.thermophysicalProperties(_liquid_porepressure[_qp], _enthalpy[_qp], _qp, _phase_state, _fsp);
 }

Referenced by computeQpProperties(), and initQpStatefulProperties().

Member Data Documentation

◆ _aqueous_phase_number

const unsigned int PorousFlowFluidStateSingleComponent::_aqueous_phase_number

protected

Phase number of the aqueous phase.

Definition at line 68 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

◆ _dfluid_density_dvar

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

protected

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

Definition at line 90 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and setMaterialVectorSize().

◆ _dfluid_enthalpy_dvar

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

protected

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

Definition at line 98 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and setMaterialVectorSize().

◆ _dfluid_internal_energy_dvar

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

protected

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

Definition at line 102 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and setMaterialVectorSize().

◆ _dfluid_viscosity_dvar

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

protected

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

Definition at line 94 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and setMaterialVectorSize().

◆ _dgrad_temperature_dgradv

MaterialProperty<std::vector<Real> >* const PorousFlowFluidStateSingleComponent::_dgrad_temperature_dgradv

protected

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

Definition at line 78 of file PorousFlowFluidStateSingleComponent.h.

◆ _dgrad_temperature_dv

MaterialProperty<std::vector<RealGradient> >* const PorousFlowFluidStateSingleComponent::_dgrad_temperature_dv

protected

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

Definition at line 80 of file PorousFlowFluidStateSingleComponent.h.

◆ _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 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> > > >& PorousFlowFluidStateSingleComponent::_dmass_frac_dvar

protected

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

Definition at line 86 of file PorousFlowFluidStateSingleComponent.h.

Referenced by 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(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlow1PhaseP::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlow2PhasePS::computeQpProperties(), 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(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlow1PhaseP::computeQpProperties(), PorousFlow2PhasePS::computeQpProperties(), computeQpProperties(), and PorousFlowFluidState::computeQpProperties().

◆ _dtemperature_dvar

MaterialProperty<std::vector<Real> >& PorousFlowFluidStateSingleComponent::_dtemperature_dvar

protected

Derivative of temperature wrt PorousFlow variables.

Definition at line 76 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and setMaterialVectorSize().

◆ _enthalpy

const VariableValue& PorousFlowFluidStateSingleComponent::_enthalpy

protected

Enthalpy.

Definition at line 58 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and thermophysicalProperties().

◆ _enthalpy_varnum

const unsigned int PorousFlowFluidStateSingleComponent::_enthalpy_varnum

protected

Moose variable number of the enthalpy.

Definition at line 62 of file PorousFlowFluidStateSingleComponent.h.

◆ _fluid_density

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

protected

Fluid density of each phase.

Definition at line 88 of file PorousFlowFluidStateSingleComponent.h.

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

◆ _fluid_enthalpy

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

protected

Enthalpy of each phase.

Definition at line 96 of file PorousFlowFluidStateSingleComponent.h.

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

◆ _fluid_internal_energy

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

protected

Internal energy of each phase.

Definition at line 100 of file PorousFlowFluidStateSingleComponent.h.

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

◆ _fluid_viscosity

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

protected

Viscosity of each phase.

Definition at line 92 of file PorousFlowFluidStateSingleComponent.h.

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

◆ _fs

const PorousFlowFluidStateSingleComponentBase& PorousFlowFluidStateSingleComponent::_fs

protected

FluidState UserObject.

Definition at line 66 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), PorousFlowFluidStateSingleComponent(), and thermophysicalProperties().

◆ _fsp

std::vector<FluidStateProperties> PorousFlowFluidStateSingleComponent::_fsp

protected

FluidStateProperties data structure.

Definition at line 108 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), initQpStatefulProperties(), PorousFlowFluidStateSingleComponent(), and thermophysicalProperties().

◆ _gas_phase_number

const unsigned int PorousFlowFluidStateSingleComponent::_gas_phase_number

protected

Phase number of the gas phase.

Definition at line 70 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties().

◆ _grad_mass_frac_qp

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

protected

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

Definition at line 84 of file PorousFlowFluidStateSingleComponent.h.

◆ _grad_temperature_qp

MaterialProperty<RealGradient>* PorousFlowFluidStateSingleComponent::_grad_temperature_qp

protected

Gradient of temperature (only defined at the qps)

Definition at line 74 of file PorousFlowFluidStateSingleComponent.h.

◆ _gradh_qp

const VariableGradient& PorousFlowFluidStateSingleComponent::_gradh_qp

protected

Gradient of enthalpy (only defined at the qps)

Definition at line 60 of file PorousFlowFluidStateSingleComponent.h.

Referenced by 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.

◆ _hidx

const unsigned int PorousFlowFluidStateSingleComponent::_hidx

protected

Index of derivative wrt enthalpy.

Definition at line 116 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties().

◆ _hvar

const unsigned int PorousFlowFluidStateSingleComponent::_hvar

protected

PorousFlow variable number of the enthalpy.

Definition at line 64 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties().

◆ _is_initqp

bool PorousFlowFluidStateSingleComponent::_is_initqp

protected

Flag to indicate whether to calculate stateful properties.

Definition at line 106 of file PorousFlowFluidStateSingleComponent.h.

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

◆ _liquid_gradp_qp

const VariableGradient& PorousFlowFluidStateSingleComponent::_liquid_gradp_qp

protected

Gradient of porepressure (only defined at the qps)

Definition at line 52 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties().

◆ _liquid_porepressure

const VariableValue& PorousFlowFluidStateSingleComponent::_liquid_porepressure

protected

Porepressure.

Definition at line 50 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and thermophysicalProperties().

◆ _liquid_porepressure_varnum

const unsigned int PorousFlowFluidStateSingleComponent::_liquid_porepressure_varnum

protected

Moose variable number of the porepressure.

Definition at line 54 of file PorousFlowFluidStateSingleComponent.h.

◆ _mass_frac

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

protected

Mass fraction matrix.

Definition at line 82 of file PorousFlowFluidStateSingleComponent.h.

Referenced by setMaterialVectorSize().

◆ _num_components

const unsigned int PorousFlowVariableBase::_num_components

protectedinherited

Number of components.

Definition at line 38 of file PorousFlowVariableBase.h.

Referenced by computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowFluidState::initQpStatefulProperties(), PorousFlowFluidState::PorousFlowFluidState(), 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(), 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(), computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), initQpStatefulProperties(), PorousFlowFluidState::initQpStatefulProperties(), PorousFlow1PhaseFullySaturated::PorousFlow1PhaseFullySaturated(), PorousFlow1PhaseMD_Gaussian::PorousFlow1PhaseMD_Gaussian(), PorousFlow1PhaseP::PorousFlow1PhaseP(), PorousFlow2PhasePP::PorousFlow2PhasePP(), PorousFlowFluidState::PorousFlowFluidState(), PorousFlowFluidStateSingleComponent(), setMaterialVectorSize(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _pc

const PorousFlowCapillaryPressure& PorousFlowFluidStateSingleComponent::_pc

protected

Capillary pressure UserObject.

Definition at line 112 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties().

◆ _phase_state

FluidStatePhaseEnum PorousFlowFluidStateSingleComponent::_phase_state

protected

FluidStatePhaseEnum.

Definition at line 110 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and thermophysicalProperties().

◆ _pidx

const unsigned int PorousFlowFluidStateSingleComponent::_pidx

protected

Index of derivative wrt pressure.

Definition at line 114 of file PorousFlowFluidStateSingleComponent.h.

Referenced by 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(), computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ _pvar

const unsigned int PorousFlowFluidStateSingleComponent::_pvar

protected

PorousFlow variable number of the porepressure.

Definition at line 56 of file PorousFlowFluidStateSingleComponent.h.

Referenced by 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(), computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ _T_c2k

const Real PorousFlowFluidStateSingleComponent::_T_c2k

protected

Conversion from degrees Celsius to degrees Kelvin.

Definition at line 104 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

◆ _temperature

MaterialProperty<Real>& PorousFlowFluidStateSingleComponent::_temperature

protected

Temperature.

Definition at line 72 of file PorousFlowFluidStateSingleComponent.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowFluidStateSingleComponent()

Member Function Documentation

◆ computeQpProperties()

◆ initQpStatefulProperties()

◆ setMaterialVectorSize()

◆ thermophysicalProperties()

Member Data Documentation

◆ _aqueous_phase_number

◆ _dfluid_density_dvar

◆ _dfluid_enthalpy_dvar

◆ _dfluid_internal_energy_dvar

◆ _dfluid_viscosity_dvar

◆ _dgrad_temperature_dgradv

◆ _dgrad_temperature_dv

◆ _dgradp_qp_dgradv

◆ _dgradp_qp_dv

◆ _dgrads_qp_dgradv

◆ _dgrads_qp_dv

◆ _dmass_frac_dvar

◆ _dporepressure_dvar

◆ _dsaturation_dvar

◆ _dtemperature_dvar

◆ _enthalpy

◆ _enthalpy_varnum

◆ _fluid_density

◆ _fluid_enthalpy

◆ _fluid_internal_energy

◆ _fluid_viscosity

◆ _fs

◆ _fsp

◆ _gas_phase_number

◆ _grad_mass_frac_qp

◆ _grad_temperature_qp

◆ _gradh_qp

◆ _gradp_qp

◆ _grads_qp

◆ _hidx

◆ _hvar

◆ _is_initqp

◆ _liquid_gradp_qp

◆ _liquid_porepressure

◆ _liquid_porepressure_varnum

◆ _mass_frac

◆ _num_components

◆ _num_pf_vars

◆ _num_phases

◆ _pc

◆ _phase_state

◆ _pidx

◆ _porepressure

◆ _pvar

◆ _saturation

◆ _T_c2k

◆ _temperature