www.mooseframework.org
Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | List of all members
PorousFlowPermeabilityExponential Class Reference

Material designed to provide the permeability tensor which is calculated from porosity using the equation: permeability = k_ijk * k, with k = BB * exp(AA * phi) where k_ijk is a tensor providing the anisotropy, phi is porosity, and A and B are empirical constants. More...

#include <PorousFlowPermeabilityExponential.h>

Inheritance diagram for PorousFlowPermeabilityExponential:
[legend]

Public Member Functions

 PorousFlowPermeabilityExponential (const InputParameters &parameters)
 

Protected Types

enum  PoropermFunction { PoropermFunction::log_k, PoropermFunction::ln_k, PoropermFunction::exp_k }
 Name of porosity-permeability relationship. More...
 

Protected Member Functions

void computeQpProperties () override
 

Protected Attributes

const Real _A
 Empirical constant A. More...
 
const Real _B
 Empirical constant B. More...
 
const RealTensorValue _k_anisotropy
 Tensor multiplier k_ijk in k = k_ijk * A * phi^n / (1 - phi)^m. More...
 
const MaterialProperty< Real > & _porosity_qp
 quadpoint porosity More...
 
const MaterialProperty< std::vector< Real > > & _dporosity_qp_dvar
 d(quadpoint porosity)/d(PorousFlow variable) More...
 
const MaterialProperty< std::vector< RealGradient > > & _dporosity_qp_dgradvar
 d(quadpoint porosity)/d(grad(PorousFlow variable)) More...
 
enum PorousFlowPermeabilityExponential::PoropermFunction _poroperm_function
 
Real _AA
 Empirical constant AA in k = k_ijk * BB * exp(AA * phi) More...
 
Real _BB
 Empirical constant BB in k = k_ijk * BB * exp(AA * phi) More...
 
MaterialProperty< RealTensorValue > & _permeability_qp
 quadpoint permeability More...
 
MaterialProperty< std::vector< RealTensorValue > > & _dpermeability_qp_dvar
 d(quadpoint permeability)/d(PorousFlow variable) More...
 
MaterialProperty< std::vector< std::vector< RealTensorValue > > > & _dpermeability_qp_dgradvar
 d(quadpoint permeability)/d(grad(PorousFlow variable)) More...
 
const unsigned int _num_phases
 Number of phases. More...
 
const unsigned int _num_components
 Number of fluid components. More...
 
const unsigned int _num_var
 Number of PorousFlow variables. More...
 

Detailed Description

Material designed to provide the permeability tensor which is calculated from porosity using the equation: permeability = k_ijk * k, with k = BB * exp(AA * phi) where k_ijk is a tensor providing the anisotropy, phi is porosity, and A and B are empirical constants.

The user can provide input for the equation expressed in any of the following 3 forms: log k = A * phi + B ln k = A * phi + B k = B * exp(A * phi) A and B are then converted to AA and BB.

Definition at line 32 of file PorousFlowPermeabilityExponential.h.

Member Enumeration Documentation

Name of porosity-permeability relationship.

Enumerator
log_k 
ln_k 
exp_k 

Definition at line 59 of file PorousFlowPermeabilityExponential.h.

59 { log_k, ln_k, exp_k } _poroperm_function;
enum PorousFlowPermeabilityExponential::PoropermFunction _poroperm_function

Constructor & Destructor Documentation

PorousFlowPermeabilityExponential::PorousFlowPermeabilityExponential ( const InputParameters &  parameters)

Definition at line 37 of file PorousFlowPermeabilityExponential.C.

39  : PorousFlowPermeabilityBase(parameters),
40  _A(getParam<Real>("A")),
41  _B(getParam<Real>("B")),
42  _k_anisotropy(parameters.isParamValid("k_anisotropy")
43  ? getParam<RealTensorValue>("k_anisotropy")
44  : RealTensorValue(1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0)),
45  _porosity_qp(getMaterialProperty<Real>("PorousFlow_porosity_qp")),
46  _dporosity_qp_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_qp_dvar")),
48  getMaterialProperty<std::vector<RealGradient>>("dPorousFlow_porosity_qp_dgradvar")),
49  _poroperm_function(getParam<MooseEnum>("poroperm_function").getEnum<PoropermFunction>())
50 {
51  switch (_poroperm_function)
52  {
54  _AA = _A * std::log(10.0);
55  _BB = std::pow(10.0, _B);
56  break;
57 
59  _AA = _A;
60  _BB = std::exp(_B);
61  break;
62 
64  _AA = _A;
65  _BB = _B;
66  break;
67  }
68 }
Real _AA
Empirical constant AA in k = k_ijk * BB * exp(AA * phi)
enum PorousFlowPermeabilityExponential::PoropermFunction _poroperm_function
const MaterialProperty< std::vector< RealGradient > > & _dporosity_qp_dgradvar
d(quadpoint porosity)/d(grad(PorousFlow variable))
const MaterialProperty< std::vector< Real > > & _dporosity_qp_dvar
d(quadpoint porosity)/d(PorousFlow variable)
PorousFlowPermeabilityBase(const InputParameters &parameters)
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
Real _BB
Empirical constant BB in k = k_ijk * BB * exp(AA * phi)
PoropermFunction
Name of porosity-permeability relationship.
const MaterialProperty< Real > & _porosity_qp
quadpoint porosity
const RealTensorValue _k_anisotropy
Tensor multiplier k_ijk in k = k_ijk * A * phi^n / (1 - phi)^m.

Member Function Documentation

void PorousFlowPermeabilityExponential::computeQpProperties ( )
overrideprotected

Definition at line 71 of file PorousFlowPermeabilityExponential.C.

72 {
73  _permeability_qp[_qp] = _k_anisotropy * _BB * std::exp(_porosity_qp[_qp] * _AA);
74 
75  _dpermeability_qp_dvar[_qp].resize(_num_var, RealTensorValue());
76  for (unsigned int v = 0; v < _num_var; ++v)
78 
79  _dpermeability_qp_dgradvar[_qp].resize(LIBMESH_DIM);
80  for (unsigned i = 0; i < LIBMESH_DIM; ++i)
81  {
82  _dpermeability_qp_dgradvar[_qp][i].resize(_num_var, RealTensorValue());
83  for (unsigned int v = 0; v < _num_var; ++v)
84  _dpermeability_qp_dgradvar[_qp][i][v] =
85  _AA * _permeability_qp[_qp] * _dporosity_qp_dgradvar[_qp][v](i);
86  }
87 }
Real _AA
Empirical constant AA in k = k_ijk * BB * exp(AA * phi)
MaterialProperty< std::vector< std::vector< RealTensorValue > > > & _dpermeability_qp_dgradvar
d(quadpoint permeability)/d(grad(PorousFlow variable))
const MaterialProperty< std::vector< RealGradient > > & _dporosity_qp_dgradvar
d(quadpoint porosity)/d(grad(PorousFlow variable))
const MaterialProperty< std::vector< Real > > & _dporosity_qp_dvar
d(quadpoint porosity)/d(PorousFlow variable)
const unsigned int _num_var
Number of PorousFlow variables.
MaterialProperty< std::vector< RealTensorValue > > & _dpermeability_qp_dvar
d(quadpoint permeability)/d(PorousFlow variable)
Real _BB
Empirical constant BB in k = k_ijk * BB * exp(AA * phi)
MaterialProperty< RealTensorValue > & _permeability_qp
quadpoint permeability
const MaterialProperty< Real > & _porosity_qp
quadpoint porosity
const RealTensorValue _k_anisotropy
Tensor multiplier k_ijk in k = k_ijk * A * phi^n / (1 - phi)^m.

Member Data Documentation

const Real PorousFlowPermeabilityExponential::_A
protected

Empirical constant A.

Definition at line 41 of file PorousFlowPermeabilityExponential.h.

Referenced by PorousFlowPermeabilityExponential().

Real PorousFlowPermeabilityExponential::_AA
protected

Empirical constant AA in k = k_ijk * BB * exp(AA * phi)

Definition at line 62 of file PorousFlowPermeabilityExponential.h.

Referenced by computeQpProperties(), and PorousFlowPermeabilityExponential().

const Real PorousFlowPermeabilityExponential::_B
protected

Empirical constant B.

Definition at line 44 of file PorousFlowPermeabilityExponential.h.

Referenced by PorousFlowPermeabilityExponential().

Real PorousFlowPermeabilityExponential::_BB
protected

Empirical constant BB in k = k_ijk * BB * exp(AA * phi)

Definition at line 65 of file PorousFlowPermeabilityExponential.h.

Referenced by computeQpProperties(), and PorousFlowPermeabilityExponential().

MaterialProperty<std::vector<std::vector<RealTensorValue> > >& PorousFlowPermeabilityBase::_dpermeability_qp_dgradvar
protectedinherited
MaterialProperty<std::vector<RealTensorValue> >& PorousFlowPermeabilityBase::_dpermeability_qp_dvar
protectedinherited
const MaterialProperty<std::vector<RealGradient> >& PorousFlowPermeabilityExponential::_dporosity_qp_dgradvar
protected

d(quadpoint porosity)/d(grad(PorousFlow variable))

Definition at line 56 of file PorousFlowPermeabilityExponential.h.

Referenced by computeQpProperties().

const MaterialProperty<std::vector<Real> >& PorousFlowPermeabilityExponential::_dporosity_qp_dvar
protected

d(quadpoint porosity)/d(PorousFlow variable)

Definition at line 53 of file PorousFlowPermeabilityExponential.h.

Referenced by computeQpProperties().

const RealTensorValue PorousFlowPermeabilityExponential::_k_anisotropy
protected

Tensor multiplier k_ijk in k = k_ijk * A * phi^n / (1 - phi)^m.

Definition at line 47 of file PorousFlowPermeabilityExponential.h.

Referenced by computeQpProperties().

const unsigned int PorousFlowMaterialVectorBase::_num_components
protectedinherited
const unsigned int PorousFlowMaterialVectorBase::_num_phases
protectedinherited
const unsigned int PorousFlowMaterialVectorBase::_num_var
protectedinherited
MaterialProperty<RealTensorValue>& PorousFlowPermeabilityBase::_permeability_qp
protectedinherited
enum PorousFlowPermeabilityExponential::PoropermFunction PorousFlowPermeabilityExponential::_poroperm_function
protected
const MaterialProperty<Real>& PorousFlowPermeabilityExponential::_porosity_qp
protected

quadpoint porosity

Definition at line 50 of file PorousFlowPermeabilityExponential.h.

Referenced by computeQpProperties().


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