DarcyFluxComponent Class Reference

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

#include <DarcyFluxComponent.h>

Inheritance diagram for DarcyFluxComponent:

Public Member Functions
	DarcyFluxComponent (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeValue ()

Protected Attributes
const VariableGradient &	_grad_pp
	gradient of the pressure More...
RealVectorValue	_fluid_weight
	fluid weight (gravity*density) as a vector pointing downwards, eg '0 0 -10000' More...
Real	_fluid_viscosity
	fluid dynamic viscosity More...
Real	_poro_recip
	(1/velocity_scaling) More...
const MaterialProperty< RealTensorValue > &	_permeability
	Material permeability. More...
unsigned int	_component
	Desired component. More...

Detailed Description

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

s^-1 . m^-2

Sometimes the fluid velocity is required, rather than the flux. In this case velocity_scaling may be used, and the result quoted above is multiplied by (1/velocity_scaling)

Definition at line 35 of file DarcyFluxComponent.h.

Constructor & Destructor Documentation

DarcyFluxComponent()

DarcyFluxComponent::DarcyFluxComponent

(

const InputParameters &

parameters

)

Definition at line 40 of file DarcyFluxComponent.C.

  : AuxKernel(parameters),
    _grad_pp(coupledGradient("porepressure")),
    _fluid_weight(getParam<RealVectorValue>("fluid_weight")),
    _fluid_viscosity(getParam<Real>("fluid_viscosity")),
    _poro_recip(1.0 / getParam<Real>("velocity_scaling")),
    _permeability(getMaterialProperty<RealTensorValue>("permeability")),
    _component(getParam<MooseEnum>("component"))
{
}

Member Function Documentation

computeValue()

Real DarcyFluxComponent::computeValue ( )

protectedvirtual

Definition at line 52 of file DarcyFluxComponent.C.

{
  return -_poro_recip *
         (_permeability[_qp] * (_grad_pp[_qp] - _fluid_weight) / _fluid_viscosity)(_component);
}

Member Data Documentation

_component

unsigned int DarcyFluxComponent::_component

protected

Desired component.

Definition at line 59 of file DarcyFluxComponent.h.

Referenced by computeValue().

_fluid_viscosity

Real DarcyFluxComponent::_fluid_viscosity

protected

fluid dynamic viscosity

Definition at line 50 of file DarcyFluxComponent.h.

Referenced by computeValue().

_fluid_weight

RealVectorValue DarcyFluxComponent::_fluid_weight

protected

fluid weight (gravity*density) as a vector pointing downwards, eg '0 0 -10000'

Definition at line 47 of file DarcyFluxComponent.h.

Referenced by computeValue().

_grad_pp

const VariableGradient& DarcyFluxComponent::_grad_pp

protected

gradient of the pressure

Definition at line 44 of file DarcyFluxComponent.h.

Referenced by computeValue().

_permeability

const MaterialProperty<RealTensorValue>& DarcyFluxComponent::_permeability

protected

Material permeability.

Definition at line 56 of file DarcyFluxComponent.h.

Referenced by computeValue().

_poro_recip

Real DarcyFluxComponent::_poro_recip

protected

(1/velocity_scaling)

Definition at line 53 of file DarcyFluxComponent.h.

Referenced by computeValue().

---

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

• richards/include/auxkernels/DarcyFluxComponent.h
• richards/src/auxkernels/DarcyFluxComponent.C