Computes a component of the Darcy velocity: -k_ij * krel /mu (nabla_j P - w_j) where k_ij is the permeability tensor, krel is the relative permeaility, mu is the fluid viscosity, P is the fluid pressure and w_j is the fluid weight This is measured in m^3 . More...

#include <PorousFlowDarcyVelocityComponent.h>

Inheritance diagram for PorousFlowDarcyVelocityComponent:

[legend]

Public Member Functions
	PorousFlowDarcyVelocityComponent (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeValue ()

Protected Attributes
const MaterialProperty< std::vector< Real > > &	_relative_permeability
	Relative permeability of each phase. More...

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

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

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

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

const PorousFlowDictator &	_dictator
	PorousFlowDicatator UserObject. More...

const unsigned int	_ph
	Index of the fluid phase. More...

unsigned int	_component
	Desired spatial component. More...

const RealVectorValue	_gravity
	Gravitational acceleration. More...

Detailed Description

Computes a component of the Darcy velocity: -k_ij * krel /mu (nabla_j P - w_j) where k_ij is the permeability tensor, krel is the relative permeaility, mu is the fluid viscosity, P is the fluid pressure and w_j is the fluid weight This is measured in m^3 .

s^-1 . m^-2

Definition at line 30 of file PorousFlowDarcyVelocityComponent.h.

Constructor & Destructor Documentation

◆ PorousFlowDarcyVelocityComponent()

PorousFlowDarcyVelocityComponent::PorousFlowDarcyVelocityComponent ( const InputParameters & parameters )

Definition at line 34 of file PorousFlowDarcyVelocityComponent.C.

   : AuxKernel(parameters),
     _relative_permeability(
         getMaterialProperty<std::vector<Real>>("PorousFlow_relative_permeability_qp")),
     _fluid_viscosity(getMaterialProperty<std::vector<Real>>("PorousFlow_viscosity_qp")),
     _permeability(getMaterialProperty<RealTensorValue>("PorousFlow_permeability_qp")),
     _grad_p(getMaterialProperty<std::vector<RealGradient>>("PorousFlow_grad_porepressure_qp")),
     _fluid_density_qp(getMaterialProperty<std::vector<Real>>("PorousFlow_fluid_phase_density_qp")),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _ph(getParam<unsigned int>("fluid_phase")),
     _component(getParam<MooseEnum>("component")),
     _gravity(getParam<RealVectorValue>("gravity"))
 {
   if (_ph >= _dictator.numPhases())
     paramError("fluid_phase",
                "The Dictator proclaims that the maximum phase index in this simulation is ",
                _dictator.numPhases() - 1,
                " whereas you have used ",
                _ph,
                ". Remember that indexing starts at 0. The Dictator is watching you, to "
                "ensure your wellbeing.");
 }

Member Function Documentation

◆ computeValue()

Real PorousFlowDarcyVelocityComponent::computeValue ( )

protectedvirtual

Reimplemented in PorousFlowDarcyVelocityComponentLowerDimensional.

Definition at line 59 of file PorousFlowDarcyVelocityComponent.C.

 {
   return -(_permeability[_qp] * (_grad_p[_qp][_ph] - _fluid_density_qp[_qp][_ph] * _gravity) *
            _relative_permeability[_qp][_ph] / _fluid_viscosity[_qp][_ph])(_component);
 }