Darcy advective flux for a fully-saturated, single phase, single component fluid. More...

#include <PorousFlowFullySaturatedDarcyBase.h>

Inheritance diagram for PorousFlowFullySaturatedDarcyBase:

Public Member Functions
	PorousFlowFullySaturatedDarcyBase (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

virtual Real	mobility () const
	The mobility of the fluid = density / viscosity. More...

virtual Real	dmobility (unsigned pvar) const
	The derivative of the mobility with respect to the PorousFlow variable pvar. More...

Protected Attributes
const bool	_multiply_by_density
	If true then the mobility contains the fluid density, otherwise it doesn't. More...

const MaterialProperty< RealTensorValue > &	_permeability
	Permeability of porous material. More...

const MaterialProperty< std::vector< RealTensorValue > > &	_dpermeability_dvar
	d(permeabiity)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< std::vector< RealTensorValue > > > &	_dpermeability_dgradvar
	d(permeabiity)/d(grad(PorousFlow variable)) More...

const MaterialProperty< std::vector< Real > > &	_density
	Fluid density for each phase (at the qp) More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_ddensity_dvar
	Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp) More...

const MaterialProperty< std::vector< Real > > &	_viscosity
	Viscosity of the fluid at the qp. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dviscosity_dvar
	Derivative of the fluid viscosity wrt PorousFlow variables. More...

const MaterialProperty< std::vector< Real > > &	_pp
	Quadpoint pore pressure in each phase. More...

const MaterialProperty< std::vector< RealGradient > > &	_grad_p
	Gradient of the pore pressure in each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dgrad_p_dgrad_var
	Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables) More...

const MaterialProperty< std::vector< std::vector< RealGradient > > > &	_dgrad_p_dvar
	Derivative of Grad porepressure in each phase wrt PorousFlow variables. More...

const PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...

const RealVectorValue	_gravity
	Gravity pointing downwards. More...

const bool	_perm_derivs
	Flag to check whether permeabiity derivatives are non-zero. More...

Detailed Description

Darcy advective flux for a fully-saturated, single phase, single component fluid.

No upwinding or relative-permeability is used.

Definition at line 25 of file PorousFlowFullySaturatedDarcyBase.h.

Constructor & Destructor Documentation

◆ PorousFlowFullySaturatedDarcyBase()

PorousFlowFullySaturatedDarcyBase::PorousFlowFullySaturatedDarcyBase ( const InputParameters & parameters )

Definition at line 36 of file PorousFlowFullySaturatedDarcyBase.C.

   : Kernel(parameters),
     _multiply_by_density(getParam<bool>("multiply_by_density")),
     _permeability(getMaterialProperty<RealTensorValue>("PorousFlow_permeability_qp")),
     _dpermeability_dvar(
         getMaterialProperty<std::vector<RealTensorValue>>("dPorousFlow_permeability_qp_dvar")),
     _dpermeability_dgradvar(getMaterialProperty<std::vector<std::vector<RealTensorValue>>>(
         "dPorousFlow_permeability_qp_dgradvar")),
     _density(getMaterialProperty<std::vector<Real>>("PorousFlow_fluid_phase_density_qp")),
     _ddensity_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
         "dPorousFlow_fluid_phase_density_qp_dvar")),
     _viscosity(getMaterialProperty<std::vector<Real>>("PorousFlow_viscosity_qp")),
     _dviscosity_dvar(
         getMaterialProperty<std::vector<std::vector<Real>>>("dPorousFlow_viscosity_qp_dvar")),
     _pp(getMaterialProperty<std::vector<Real>>("PorousFlow_porepressure_qp")),
     _grad_p(getMaterialProperty<std::vector<RealGradient>>("PorousFlow_grad_porepressure_qp")),
     _dgrad_p_dgrad_var(getMaterialProperty<std::vector<std::vector<Real>>>(
         "dPorousFlow_grad_porepressure_qp_dgradvar")),
     _dgrad_p_dvar(getMaterialProperty<std::vector<std::vector<RealGradient>>>(
         "dPorousFlow_grad_porepressure_qp_dvar")),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _gravity(getParam<RealVectorValue>("gravity")),
     _perm_derivs(_dictator.usePermDerivs())
 {
   if (_dictator.numPhases() != 1)
     mooseError("PorousFlowFullySaturatedDarcyBase should not be used for multi-phase scenarios as "
                "it does no upwinding and does not include relative-permeability effects");
 }

Member Function Documentation

◆ computeQpJacobian()

Real PorousFlowFullySaturatedDarcyBase::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 77 of file PorousFlowFullySaturatedDarcyBase.C.

 {
   return computeQpOffDiagJacobian(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowFullySaturatedDarcyBase::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 83 of file PorousFlowFullySaturatedDarcyBase.C.

 {
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
  
   const unsigned ph = 0;
   const unsigned pvar = _dictator.porousFlowVariableNum(jvar);
  
   const Real mob = mobility();
   const Real dmob = dmobility(pvar) * _phi[_j][_qp];
  
   const RealVectorValue flow =
       _permeability[_qp] * (_grad_p[_qp][ph] - _density[_qp][ph] * _gravity);
  
   RealVectorValue dflow =
       _permeability[_qp] * (_grad_phi[_j][_qp] * _dgrad_p_dgrad_var[_qp][ph][pvar] -
                             _phi[_j][_qp] * _ddensity_dvar[_qp][ph][pvar] * _gravity);
   dflow += _permeability[_qp] * (_dgrad_p_dvar[_qp][ph][pvar] * _phi[_j][_qp]);
  
   if (_perm_derivs)
   {
     dflow += _dpermeability_dvar[_qp][pvar] * _phi[_j][_qp] *
              (_grad_p[_qp][ph] - _density[_qp][ph] * _gravity);
     for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
       dflow += _dpermeability_dgradvar[_qp][i][pvar] * _grad_phi[_j][_qp](i) *
                (_grad_p[_qp][ph] - _density[_qp][ph] * _gravity);
   }
  
   return _grad_test[_i][_qp] * (dmob * flow + mob * dflow);
 }

Referenced by computeQpJacobian().

◆ computeQpResidual()

Real PorousFlowFullySaturatedDarcyBase::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 67 of file PorousFlowFullySaturatedDarcyBase.C.

 {
   const unsigned ph = 0;
   const Real mob = mobility();
   const RealVectorValue flow =
       _permeability[_qp] * (_grad_p[_qp][ph] - _density[_qp][ph] * _gravity);
   return _grad_test[_i][_qp] * mob * flow;
 }

◆ dmobility()

Real PorousFlowFullySaturatedDarcyBase::dmobility ( unsigned pvar ) const

protectedvirtual

The derivative of the mobility with respect to the PorousFlow variable pvar.

Parameters

pvar	Take the derivative with respect to this PorousFlow variable

Reimplemented in PorousFlowFullySaturatedDarcyFlow, and PorousFlowFullySaturatedHeatAdvection.

Definition at line 125 of file PorousFlowFullySaturatedDarcyBase.C.

 {
   const unsigned ph = 0;
   Real dmob = -_dviscosity_dvar[_qp][ph][pvar] / std::pow(_viscosity[_qp][ph], 2);
   if (_multiply_by_density)
     dmob = _density[_qp][ph] * dmob + _ddensity_dvar[_qp][ph][pvar] / _viscosity[_qp][ph];
   return dmob;
 }

Referenced by computeQpOffDiagJacobian(), PorousFlowFullySaturatedHeatAdvection::dmobility(), and PorousFlowFullySaturatedDarcyFlow::dmobility().

◆ mobility()

Real PorousFlowFullySaturatedDarcyBase::mobility ( ) const

protectedvirtual

The mobility of the fluid = density / viscosity.

Reimplemented in PorousFlowFullySaturatedDarcyFlow, and PorousFlowFullySaturatedHeatAdvection.

Definition at line 115 of file PorousFlowFullySaturatedDarcyBase.C.

 {
   const unsigned ph = 0;
   Real mob = 1.0 / _viscosity[_qp][ph];
   if (_multiply_by_density)
     mob *= _density[_qp][ph];
   return mob;
 }

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), PorousFlowFullySaturatedHeatAdvection::dmobility(), PorousFlowFullySaturatedDarcyFlow::dmobility(), PorousFlowFullySaturatedHeatAdvection::mobility(), and PorousFlowFullySaturatedDarcyFlow::mobility().

Member Data Documentation

◆ _ddensity_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFullySaturatedDarcyBase::_ddensity_dvar

protected

Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp)

Definition at line 62 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian(), and dmobility().

◆ _density

const MaterialProperty<std::vector<Real> >& PorousFlowFullySaturatedDarcyBase::_density

protected

Fluid density for each phase (at the qp)

Definition at line 59 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), dmobility(), and mobility().

◆ _dgrad_p_dgrad_var

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFullySaturatedDarcyBase::_dgrad_p_dgrad_var

protected

Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables)

Definition at line 77 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _dgrad_p_dvar

const MaterialProperty<std::vector<std::vector<RealGradient> > >& PorousFlowFullySaturatedDarcyBase::_dgrad_p_dvar

protected

Derivative of Grad porepressure in each phase wrt PorousFlow variables.

Definition at line 80 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _dictator

const PorousFlowDictator& PorousFlowFullySaturatedDarcyBase::_dictator

protected

PorousFlowDictator UserObject.

Definition at line 83 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian(), PorousFlowFullySaturatedDarcyBase(), and PorousFlowFullySaturatedDarcyFlow::PorousFlowFullySaturatedDarcyFlow().

◆ _dpermeability_dgradvar

const MaterialProperty<std::vector<std::vector<RealTensorValue> > >& PorousFlowFullySaturatedDarcyBase::_dpermeability_dgradvar

protected

d(permeabiity)/d(grad(PorousFlow variable))

Definition at line 56 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _dpermeability_dvar

const MaterialProperty<std::vector<RealTensorValue> >& PorousFlowFullySaturatedDarcyBase::_dpermeability_dvar

protected

d(permeabiity)/d(PorousFlow variable)

Definition at line 53 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _dviscosity_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowFullySaturatedDarcyBase::_dviscosity_dvar

protected

Derivative of the fluid viscosity wrt PorousFlow variables.

Definition at line 68 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by dmobility().

◆ _grad_p

const MaterialProperty<std::vector<RealGradient> >& PorousFlowFullySaturatedDarcyBase::_grad_p

protected

Gradient of the pore pressure in each phase.

Definition at line 74 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _gravity

const RealVectorValue PorousFlowFullySaturatedDarcyBase::_gravity

protected

Gravity pointing downwards.

Definition at line 86 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _multiply_by_density

const bool PorousFlowFullySaturatedDarcyBase::_multiply_by_density

protected

If true then the mobility contains the fluid density, otherwise it doesn't.

Definition at line 47 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by dmobility(), and mobility().

◆ _perm_derivs

const bool PorousFlowFullySaturatedDarcyBase::_perm_derivs

protected

Flag to check whether permeabiity derivatives are non-zero.

Definition at line 89 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _permeability

const MaterialProperty<RealTensorValue>& PorousFlowFullySaturatedDarcyBase::_permeability

protected

Permeability of porous material.

Definition at line 50 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _pp

const MaterialProperty<std::vector<Real> >& PorousFlowFullySaturatedDarcyBase::_pp

protected

Quadpoint pore pressure in each phase.

Definition at line 71 of file PorousFlowFullySaturatedDarcyBase.h.

◆ _viscosity

const MaterialProperty<std::vector<Real> >& PorousFlowFullySaturatedDarcyBase::_viscosity

protected

Viscosity of the fluid at the qp.

Definition at line 65 of file PorousFlowFullySaturatedDarcyBase.h.

Referenced by dmobility(), and mobility().

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

porous_flow/include/kernels/PorousFlowFullySaturatedDarcyBase.h
porous_flow/src/kernels/PorousFlowFullySaturatedDarcyBase.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowFullySaturatedDarcyBase()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ dmobility()

◆ mobility()

Member Data Documentation

◆ _ddensity_dvar

◆ _density

◆ _dgrad_p_dgrad_var

◆ _dgrad_p_dvar

◆ _dictator

◆ _dpermeability_dgradvar

◆ _dpermeability_dvar

◆ _dviscosity_dvar

◆ _grad_p

◆ _gravity

◆ _multiply_by_density

◆ _perm_derivs

◆ _permeability

◆ _pp

◆ _viscosity