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

#include <PorousFlowDarcyVelocityComponentLowerDimensional.h>

Inheritance diagram for PorousFlowDarcyVelocityComponentLowerDimensional:

[legend]

Public Member Functions
	PorousFlowDarcyVelocityComponentLowerDimensional (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeValue () override

Protected Attributes
const VariableValue &	_aperture
	Fracture aperture (width) More...

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 a) (nabla_j P - w_j) where k_ij is the permeability tensor, krel is the relative permeaility, mu is the fluid viscosity, a is the fracture aperture, P is the fluid pressure and w_j is the fluid weight tensor that is projected in the tangent direction of this element This is measured in m^3 .

s^-1 . m^-2

Definition at line 30 of file PorousFlowDarcyVelocityComponentLowerDimensional.h.

Constructor & Destructor Documentation

◆ PorousFlowDarcyVelocityComponentLowerDimensional()

PorousFlowDarcyVelocityComponentLowerDimensional::PorousFlowDarcyVelocityComponentLowerDimensional ( const InputParameters & parameters )

Definition at line 34 of file PorousFlowDarcyVelocityComponentLowerDimensional.C.

   : PorousFlowDarcyVelocityComponent(parameters), _aperture(coupledValue("aperture"))
 {
 }

Member Function Documentation

◆ computeValue()

Real PorousFlowDarcyVelocityComponentLowerDimensional::computeValue ( )

overrideprotectedvirtual

Reimplemented from PorousFlowDarcyVelocityComponent.

Definition at line 41 of file PorousFlowDarcyVelocityComponentLowerDimensional.C.

 {
   const unsigned elem_dim = _current_elem->dim();
   if (elem_dim == _mesh.dimension())
     mooseError("The variable ",
                _var.name(),
                " must must be defined on lower-dimensional elements "
                "only since it employs "
                "PorousFlowDarcyVelocityComponentLowerDimensional\n");
  
   // MOOSE will automatically make grad(P) lie along the element's tangent direction
   // but we need to project gravity in that direction too
  
   // tang_xi is the element's tangent vector in xi direction
   const std::vector<RealGradient> & tang_xi = _assembly.getFE(FEType(), elem_dim)->get_dxyzdxi();
   RealVectorValue tangential_gravity =
       (_gravity * tang_xi[_qp] / tang_xi[_qp].norm_sq()) * tang_xi[_qp];
   if (elem_dim == 2)
   {
     // tang_eta is the element's tangent vector in eta direction
     const std::vector<RealGradient> & tang_eta =
         _assembly.getFE(FEType(), elem_dim)->get_dxyzdeta();
     tangential_gravity += (_gravity * tang_eta[_qp] / tang_eta[_qp].norm_sq()) * tang_eta[_qp];
   }
  
   return -(_permeability[_qp] *
            (_grad_p[_qp][_ph] - _fluid_density_qp[_qp][_ph] * tangential_gravity) *
            _relative_permeability[_qp][_ph] / _fluid_viscosity[_qp][_ph])(_component) /
          _aperture[_qp];
 }