Base material designed to calculate fluid phase porepressure and saturation for the two-phase situation assuming phase porepressures as the nonlinear variables. More...

#include <PorousFlow2PhasePP.h>

Inheritance diagram for PorousFlow2PhasePP:

[legend]

Public Member Functions
	PorousFlow2PhasePP (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

Real	buildQpPPSS ()
	Assemble std::vectors of porepressure and saturation at the nodes and quadpoints, and return the capillary pressure. More...

Protected Attributes
const VariableValue &	_phase0_porepressure
	Nodal or quadpoint value of porepressure of the zero phase (eg, the water phase) More...

const VariableGradient &	_phase0_gradp_qp
	Gradient(phase0_porepressure) at the qps. More...

const unsigned int	_phase0_porepressure_varnum
	Moose variable number of the phase0 porepressure. More...

const unsigned int	_p0var
	PorousFlow variable number of the phase0 porepressure. More...

const VariableValue &	_phase1_porepressure
	Nodal or quadpoint value of porepressure of the one phase (eg, the gas phase) More...

const VariableGradient &	_phase1_gradp_qp
	Gradient(phase1_porepressure) at the qps. More...

const unsigned int	_phase1_porepressure_varnum
	Moose variable number of the phase1 porepressure. More...

const unsigned int	_p1var
	PorousFlow variable number of the phase1 porepressure. More...

const PorousFlowCapillaryPressure &	_pc_uo
	Capillary pressure UserObject. 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

Base material designed to calculate fluid phase porepressure and saturation for the two-phase situation assuming phase porepressures as the nonlinear variables.

Definition at line 24 of file PorousFlow2PhasePP.h.

Constructor & Destructor Documentation

◆ PorousFlow2PhasePP()

PorousFlow2PhasePP::PorousFlow2PhasePP ( const InputParameters & parameters )

Definition at line 34 of file PorousFlow2PhasePP.C.

   : PorousFlowVariableBase(parameters),
  
     _phase0_porepressure(_nodal_material ? coupledDofValues("phase0_porepressure")
                                          : coupledValue("phase0_porepressure")),
     _phase0_gradp_qp(coupledGradient("phase0_porepressure")),
     _phase0_porepressure_varnum(coupled("phase0_porepressure")),
     _p0var(_dictator.isPorousFlowVariable(_phase0_porepressure_varnum)
                ? _dictator.porousFlowVariableNum(_phase0_porepressure_varnum)
                : 0),
  
     _phase1_porepressure(_nodal_material ? coupledDofValues("phase1_porepressure")
                                          : coupledValue("phase1_porepressure")),
     _phase1_gradp_qp(coupledGradient("phase1_porepressure")),
     _phase1_porepressure_varnum(coupled("phase1_porepressure")),
     _p1var(_dictator.isPorousFlowVariable(_phase1_porepressure_varnum)
                ? _dictator.porousFlowVariableNum(_phase1_porepressure_varnum)
                : 0),
     _pc_uo(getUserObject<PorousFlowCapillaryPressure>("capillary_pressure"))
 {
   if (_num_phases != 2)
     mooseError("The Dictator announces that the number of phases is ",
                _dictator.numPhases(),
                " whereas PorousFlow2PhasePP can only be used for 2-phase simulation.  When you "
                "have an efficient government, you have a dictatorship.");
 }

Member Function Documentation

◆ buildQpPPSS()

Real PorousFlow2PhasePP::buildQpPPSS ( )

protected

Assemble std::vectors of porepressure and saturation at the nodes and quadpoints, and return the capillary pressure.

Definition at line 132 of file PorousFlow2PhasePP.C.

 {
   _porepressure[_qp][0] = _phase0_porepressure[_qp];
   _porepressure[_qp][1] = _phase1_porepressure[_qp];
   const Real pc = _phase0_porepressure[_qp] - _phase1_porepressure[_qp]; // this is <= 0
   const Real sat = _pc_uo.saturation(pc);
   _saturation[_qp][0] = sat;
   _saturation[_qp][1] = 1.0 - sat;
   return pc;
 }

Referenced by computeQpProperties(), and initQpStatefulProperties().

◆ computeQpProperties()

void PorousFlow2PhasePP::computeQpProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBase.

Definition at line 69 of file PorousFlow2PhasePP.C.

 {
   // size stuff correctly and prepare the derivative matrices with zeroes
   PorousFlowVariableBase::computeQpProperties();
  
   const Real pc = buildQpPPSS();
   const Real ds = _pc_uo.dSaturation(pc); // dS/d(pc)
  
   if (!_nodal_material)
   {
     (*_gradp_qp)[_qp][0] = _phase0_gradp_qp[_qp];
     (*_gradp_qp)[_qp][1] = _phase1_gradp_qp[_qp];
     (*_grads_qp)[_qp][0] = ds * ((*_gradp_qp)[_qp][0] - (*_gradp_qp)[_qp][1]);
     (*_grads_qp)[_qp][1] = -(*_grads_qp)[_qp][0];
   }
  
   // the derivatives of porepressure with respect to porepressure
   // remain fixed (at unity) throughout the simulation
   if (_dictator.isPorousFlowVariable(_phase0_porepressure_varnum))
   {
     _dporepressure_dvar[_qp][0][_p0var] = 1.0;
     if (!_nodal_material)
       (*_dgradp_qp_dgradv)[_qp][0][_p0var] = 1.0;
   }
   if (_dictator.isPorousFlowVariable(_phase1_porepressure_varnum))
   {
     _dporepressure_dvar[_qp][1][_p1var] = 1.0;
     if (!_nodal_material)
       (*_dgradp_qp_dgradv)[_qp][1][_p1var] = 1.0;
   }
  
   if (_dictator.isPorousFlowVariable(_phase0_porepressure_varnum))
   {
     _dsaturation_dvar[_qp][0][_p0var] = ds;
     _dsaturation_dvar[_qp][1][_p0var] = -ds;
   }
   if (_dictator.isPorousFlowVariable(_phase1_porepressure_varnum))
   {
     _dsaturation_dvar[_qp][0][_p1var] = -ds;
     _dsaturation_dvar[_qp][1][_p1var] = ds;
   }
  
   if (!_nodal_material)
   {
     const Real d2s_qp = _pc_uo.d2Saturation(pc); // d^2(S_qp)/d(pc_qp)^2
     if (_dictator.isPorousFlowVariable(_phase0_porepressure_varnum))
     {
       (*_dgrads_qp_dgradv)[_qp][0][_p0var] = ds;
       (*_dgrads_qp_dv)[_qp][0][_p0var] = d2s_qp * (_phase0_gradp_qp[_qp] - _phase1_gradp_qp[_qp]);
       (*_dgrads_qp_dgradv)[_qp][1][_p0var] = -ds;
       (*_dgrads_qp_dv)[_qp][1][_p0var] = -d2s_qp * (_phase0_gradp_qp[_qp] - _phase1_gradp_qp[_qp]);
     }
     if (_dictator.isPorousFlowVariable(_phase1_porepressure_varnum))
     {
       (*_dgrads_qp_dgradv)[_qp][0][_p1var] = -ds;
       (*_dgrads_qp_dv)[_qp][0][_p1var] = -d2s_qp * (_phase0_gradp_qp[_qp] - _phase1_gradp_qp[_qp]);
       (*_dgrads_qp_dgradv)[_qp][1][_p1var] = ds;
       (*_dgrads_qp_dv)[_qp][1][_p1var] = d2s_qp * (_phase0_gradp_qp[_qp] - _phase1_gradp_qp[_qp]);
     }
   }
 }

◆ initQpStatefulProperties()

void PorousFlow2PhasePP::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBase.

Definition at line 62 of file PorousFlow2PhasePP.C.

 {
   PorousFlowVariableBase::initQpStatefulProperties();
   buildQpPPSS();
 }

Member Data Documentation

◆ _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.

◆ _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(), 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 computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlow1PhaseP::computeQpProperties(), PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), and 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.

◆ _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(), PorousFlowFluidState::PorousFlowFluidState(), PorousFlowFluidStateSingleComponent::PorousFlowFluidStateSingleComponent(), PorousFlowFluidStateSingleComponent::setMaterialVectorSize(), and PorousFlowFluidState::setMaterialVectorSize().

◆ _p0var

const unsigned int PorousFlow2PhasePP::_p0var

protected

PorousFlow variable number of the phase0 porepressure.

Definition at line 46 of file PorousFlow2PhasePP.h.

Referenced by computeQpProperties().

◆ _p1var

const unsigned int PorousFlow2PhasePP::_p1var

protected

PorousFlow variable number of the phase1 porepressure.

Definition at line 54 of file PorousFlow2PhasePP.h.

Referenced by computeQpProperties().

◆ _pc_uo

const PorousFlowCapillaryPressure& PorousFlow2PhasePP::_pc_uo

protected

Capillary pressure UserObject.

Definition at line 56 of file PorousFlow2PhasePP.h.

Referenced by buildQpPPSS(), and computeQpProperties().

◆ _phase0_gradp_qp

const VariableGradient& PorousFlow2PhasePP::_phase0_gradp_qp

protected

Gradient(phase0_porepressure) at the qps.

Definition at line 42 of file PorousFlow2PhasePP.h.

Referenced by computeQpProperties().

◆ _phase0_porepressure

const VariableValue& PorousFlow2PhasePP::_phase0_porepressure

protected

Nodal or quadpoint value of porepressure of the zero phase (eg, the water phase)

Definition at line 40 of file PorousFlow2PhasePP.h.

Referenced by buildQpPPSS().

◆ _phase0_porepressure_varnum

const unsigned int PorousFlow2PhasePP::_phase0_porepressure_varnum

protected

Moose variable number of the phase0 porepressure.

Definition at line 44 of file PorousFlow2PhasePP.h.

Referenced by computeQpProperties().

◆ _phase1_gradp_qp

const VariableGradient& PorousFlow2PhasePP::_phase1_gradp_qp

protected

Gradient(phase1_porepressure) at the qps.

Definition at line 50 of file PorousFlow2PhasePP.h.

Referenced by computeQpProperties().

◆ _phase1_porepressure

const VariableValue& PorousFlow2PhasePP::_phase1_porepressure

protected

Nodal or quadpoint value of porepressure of the one phase (eg, the gas phase)

Definition at line 48 of file PorousFlow2PhasePP.h.

Referenced by buildQpPPSS().

◆ _phase1_porepressure_varnum

const unsigned int PorousFlow2PhasePP::_phase1_porepressure_varnum

protected

Moose variable number of the phase1 porepressure.

Definition at line 52 of file PorousFlow2PhasePP.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(), buildQpPPSS(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), PorousFlowFluidStateSingleComponent::initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ _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(), buildQpPPSS(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), PorousFlowFluidStateSingleComponent::initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlow2PhasePP()

Member Function Documentation

◆ buildQpPPSS()

◆ computeQpProperties()

◆ initQpStatefulProperties()

Member Data Documentation

◆ _dgradp_qp_dgradv

◆ _dgradp_qp_dv

◆ _dgrads_qp_dgradv

◆ _dgrads_qp_dv

◆ _dporepressure_dvar

◆ _dsaturation_dvar

◆ _gradp_qp

◆ _grads_qp

◆ _num_components

◆ _num_pf_vars

◆ _num_phases

◆ _p0var

◆ _p1var

◆ _pc_uo

◆ _phase0_gradp_qp

◆ _phase0_porepressure

◆ _phase0_porepressure_varnum

◆ _phase1_gradp_qp

◆ _phase1_porepressure

◆ _phase1_porepressure_varnum

◆ _porepressure

◆ _saturation