Material designed to calculate fluid-phase porepressure and saturation for the single-phase situation, assuming a Gaussian capillary suction function and assuming the independent variable is log(mass density) and assuming the fluid has a constant bulk modulus. More...

#include <PorousFlow1PhaseMD_Gaussian.h>

Inheritance diagram for PorousFlow1PhaseMD_Gaussian:

[legend]

Public Member Functions
	PorousFlow1PhaseMD_Gaussian (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

virtual void	buildPS ()

Protected Attributes
const Real	_al
	Gaussian parameter: saturation = exp(-(al*p)^2) More...

const Real	_al2
	_al2 = al*al More...

const Real	_logdens0
	Fluid density = _dens0*exp(P/_bulk) More...

const Real	_bulk
	Fluid density = _dens0*exp(P/_bulk) More...

const Real	_recip_bulk
	1/_bulk/_al More...

const Real	_recip_bulk2
	(1/_bulk)^2 More...

const VariableValue &	_md_var
	Nodal or quadpoint value of mass-density of the fluid phase. More...

const VariableGradient &	_gradmd_qp_var
	Gradient(_mass-density at quadpoints) More...

const unsigned int	_md_varnum
	Moose variable number of the mass-density. More...

const unsigned int	_pvar
	PorousFlow variable number of the mass-density. 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

Material designed to calculate fluid-phase porepressure and saturation for the single-phase situation, assuming a Gaussian capillary suction function and assuming the independent variable is log(mass density) and assuming the fluid has a constant bulk modulus.

Definition at line 25 of file PorousFlow1PhaseMD_Gaussian.h.

Constructor & Destructor Documentation

◆ PorousFlow1PhaseMD_Gaussian()

PorousFlow1PhaseMD_Gaussian::PorousFlow1PhaseMD_Gaussian ( const InputParameters & parameters )

Definition at line 38 of file PorousFlow1PhaseMD_Gaussian.C.

   : PorousFlowVariableBase(parameters),
  
     _al(getParam<Real>("al")),
     _al2(std::pow(_al, 2)),
     _logdens0(std::log(getParam<Real>("density_P0"))),
     _bulk(getParam<Real>("bulk_modulus")),
     _recip_bulk(1.0 / _al / _bulk),
     _recip_bulk2(std::pow(_recip_bulk, 2)),
  
     _md_var(_nodal_material ? coupledDofValues("mass_density") : coupledValue("mass_density")),
     _gradmd_qp_var(coupledGradient("mass_density")),
     _md_varnum(coupled("mass_density")),
     _pvar(_dictator.isPorousFlowVariable(_md_varnum) ? _dictator.porousFlowVariableNum(_md_varnum)
                                                      : 0)
 {
   if (_num_phases != 1)
     mooseError("The Dictator proclaims that the number of phases is ",
                _dictator.numPhases(),
                " whereas PorousFlow1PhaseMD_Gaussian can only be used for 1-phase simulations.  Be "
                "aware that the Dictator has noted your mistake.");
 }

Member Function Documentation

◆ buildPS()

void PorousFlow1PhaseMD_Gaussian::buildPS ( )

protectedvirtual

Definition at line 139 of file PorousFlow1PhaseMD_Gaussian.C.

 {
   if (_md_var[_qp] >= _logdens0)
   {
     // full saturation
     _porepressure[_qp][0] = (_md_var[_qp] - _logdens0) * _bulk;
     _saturation[_qp][0] = 1.0;
   }
   else
   {
     // v = logdens0 + p/bulk - (al p)^2
     // 0 = (v-logdens0) - p/bulk + (al p)^2
     // 2 al p = (1/al/bulk) +/- sqrt((1/al/bulk)^2 - 4(v-logdens0))  (the "minus" sign is chosen)
     // s = exp(-(al*p)^2)
     _porepressure[_qp][0] =
         (_recip_bulk - std::sqrt(_recip_bulk2 + 4.0 * (_logdens0 - _md_var[_qp]))) / (2.0 * _al);
     _saturation[_qp][0] = std::exp(-std::pow(_al * _porepressure[_qp][0], 2.0));
   }
 }

Referenced by computeQpProperties(), and initQpStatefulProperties().

◆ computeQpProperties()

void PorousFlow1PhaseMD_Gaussian::computeQpProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBase.

Definition at line 69 of file PorousFlow1PhaseMD_Gaussian.C.

 {
   // size stuff correctly and prepare the derivative matrices with zeroes
   PorousFlowVariableBase::computeQpProperties();
  
   buildPS();
  
   if (!_nodal_material)
   {
     if (_md_var[_qp] >= _logdens0)
     {
       (*_gradp_qp)[_qp][0] = _gradmd_qp_var[_qp] * _bulk;
       (*_grads_qp)[_qp][0] = 0.0;
     }
     else
     {
       (*_gradp_qp)[_qp][0] =
           _gradmd_qp_var[_qp] / (_recip_bulk - 2.0 * _al * _porepressure[_qp][0]) / _al;
       (*_grads_qp)[_qp][0] =
           -2.0 * _al2 * _porepressure[_qp][0] * _saturation[_qp][0] * (*_gradp_qp)[_qp][0];
     }
   }
  
   if (_dictator.notPorousFlowVariable(_md_varnum))
     return;
  
   if (_md_var[_qp] >= _logdens0)
   {
     // fully saturated at the node or quadpoint
     _dporepressure_dvar[_qp][0][_pvar] = _bulk;
     _dsaturation_dvar[_qp][0][_pvar] = 0.0;
   }
   else
   {
     const Real pp = _porepressure[_qp][0];
     _dporepressure_dvar[_qp][0][_pvar] = 1.0 / (_recip_bulk - 2.0 * _al * pp) / _al;
     const Real sat = _saturation[_qp][0];
     _dsaturation_dvar[_qp][0][_pvar] = -2.0 * _al2 * pp * sat * _dporepressure_dvar[_qp][0][_pvar];
   }
  
   if (!_nodal_material)
   {
     if (_md_var[_qp] >= _logdens0)
     {
       // fully saturated at the quadpoint
       (*_dgradp_qp_dgradv)[_qp][0][_pvar] = _bulk;
       (*_dgradp_qp_dv)[_qp][0][_pvar] = 0.0;
       (*_dgrads_qp_dgradv)[_qp][0][_pvar] = 0.0;
       (*_dgrads_qp_dv)[_qp][0][_pvar] = 0.0;
     }
     else
     {
       const Real pp = _porepressure[_qp][0];
       (*_dgradp_qp_dgradv)[_qp][0][_pvar] = 1.0 / (_recip_bulk - 2.0 * _al * pp) / _al;
       (*_dgradp_qp_dv)[_qp][0][_pvar] = _gradmd_qp_var[_qp] * 2.0 * _al *
                                         _dporepressure_dvar[_qp][0][_pvar] /
                                         std::pow(_recip_bulk - 2.0 * _al * pp, 2.0) / _al;
       const Real sat = _saturation[_qp][0];
       (*_dgrads_qp_dgradv)[_qp][0][_pvar] =
           -2.0 * _al2 * pp * sat * (*_dgradp_qp_dgradv)[_qp][0][_pvar];
       (*_dgrads_qp_dv)[_qp][0][_pvar] =
           -2.0 * _al2 * _dporepressure_dvar[_qp][0][_pvar] * sat * (*_gradp_qp)[_qp][0];
       (*_dgrads_qp_dv)[_qp][0][_pvar] +=
           -2.0 * _al2 * pp * _dsaturation_dvar[_qp][0][_pvar] * (*_gradp_qp)[_qp][0];
       (*_dgrads_qp_dv)[_qp][0][_pvar] += -2.0 * _al2 * pp * sat * (*_dgradp_qp_dv)[_qp][0][_pvar];
     }
   }
 }

◆ initQpStatefulProperties()

void PorousFlow1PhaseMD_Gaussian::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBase.

Definition at line 62 of file PorousFlow1PhaseMD_Gaussian.C.

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

Member Data Documentation

◆ _al

const Real PorousFlow1PhaseMD_Gaussian::_al

protected

Gaussian parameter: saturation = exp(-(al*p)^2)

Definition at line 35 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by buildPS(), and computeQpProperties().

◆ _al2

const Real PorousFlow1PhaseMD_Gaussian::_al2

protected

_al2 = al*al

Definition at line 38 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by computeQpProperties().

◆ _bulk

const Real PorousFlow1PhaseMD_Gaussian::_bulk

protected

Fluid density = _dens0*exp(P/_bulk)

Definition at line 44 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by buildPS(), and computeQpProperties().

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

◆ _gradmd_qp_var

const VariableGradient& PorousFlow1PhaseMD_Gaussian::_gradmd_qp_var

protected

Gradient(_mass-density at quadpoints)

Definition at line 56 of file PorousFlow1PhaseMD_Gaussian.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.

◆ _logdens0

const Real PorousFlow1PhaseMD_Gaussian::_logdens0

protected

Fluid density = _dens0*exp(P/_bulk)

Definition at line 41 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by buildPS(), and computeQpProperties().

◆ _md_var

const VariableValue& PorousFlow1PhaseMD_Gaussian::_md_var

protected

Nodal or quadpoint value of mass-density of the fluid phase.

Definition at line 53 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by buildPS(), and computeQpProperties().

◆ _md_varnum

const unsigned int PorousFlow1PhaseMD_Gaussian::_md_varnum

protected

Moose variable number of the mass-density.

Definition at line 59 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by computeQpProperties().

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

◆ _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 buildPS(), PorousFlow2PhasePS::buildQpPPSS(), PorousFlow1PhaseFullySaturated::buildQpPPSS(), PorousFlow1PhaseP::buildQpPPSS(), PorousFlow2PhasePP::buildQpPPSS(), computeQpProperties(), PorousFlowVariableBase::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), PorousFlowFluidState::computeQpProperties(), PorousFlowVariableBase::initQpStatefulProperties(), PorousFlowFluidStateSingleComponent::initQpStatefulProperties(), and PorousFlowFluidState::initQpStatefulProperties().

◆ _pvar

const unsigned int PorousFlow1PhaseMD_Gaussian::_pvar

protected

PorousFlow variable number of the mass-density.

Definition at line 62 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by computeQpProperties().

◆ _recip_bulk

const Real PorousFlow1PhaseMD_Gaussian::_recip_bulk

protected

1/_bulk/_al

Definition at line 47 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by buildPS(), and computeQpProperties().

◆ _recip_bulk2

const Real PorousFlow1PhaseMD_Gaussian::_recip_bulk2

protected

(1/_bulk)^2

Definition at line 50 of file PorousFlow1PhaseMD_Gaussian.h.

Referenced by buildPS().

◆ _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 buildPS(), PorousFlow2PhasePS::buildQpPPSS(), PorousFlow1PhaseFullySaturated::buildQpPPSS(), PorousFlow1PhaseP::buildQpPPSS(), PorousFlow2PhasePP::buildQpPPSS(), 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/PorousFlow1PhaseMD_Gaussian.h
porous_flow/src/materials/PorousFlow1PhaseMD_Gaussian.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlow1PhaseMD_Gaussian()

Member Function Documentation

◆ buildPS()

◆ computeQpProperties()

◆ initQpStatefulProperties()

Member Data Documentation

◆ _al

◆ _al2

◆ _bulk

◆ _dgradp_qp_dgradv

◆ _dgradp_qp_dv

◆ _dgrads_qp_dgradv

◆ _dgrads_qp_dv

◆ _dporepressure_dvar

◆ _dsaturation_dvar

◆ _gradmd_qp_var

◆ _gradp_qp

◆ _grads_qp

◆ _logdens0

◆ _md_var

◆ _md_varnum

◆ _num_components

◆ _num_pf_vars

◆ _num_phases

◆ _porepressure

◆ _pvar

◆ _recip_bulk

◆ _recip_bulk2

◆ _saturation