PorousFlowPermeabilityKozenyCarman Class Reference

Material designed to provide the permeability tensor which is calculated from porosity using a form of the Kozeny-Carman equation (e.g. More...

#include <PorousFlowPermeabilityKozenyCarman.h>

Inheritance diagram for PorousFlowPermeabilityKozenyCarman:
[legend]

Public Member Functions

PorousFlowPermeabilityKozenyCarman (const InputParameters &parameters)

Protected Types

enum  PoropermFunction { PoropermFunction::kozeny_carman_fd2, PoropermFunction::kozeny_carman_phi0 }
Name of porosity-permeability relationship. More...

Protected Member Functions

void computeQpProperties () override

Protected Attributes

const Real _k0
Reference scalar permeability in A = k0 * (1 - phi0)^m / phi0^n. More...

const Real _phi0
Reference porosity in A = k0 * (1 - phi0)^m / phi0^n. More...

const Real _f
Multiplying factor in A = f * d^2. More...

const Real _d
Grain diameter in A = f * d^2. More...

const Real _m
Exponent in k = k_ijk * A * phi^n / (1 - phi)^m. More...

const Real _n
Exponent in k = k_ijk * A * phi^n / (1 - phi)^m. 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

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

enum PorousFlowPermeabilityKozenyCarman::PoropermFunction _poroperm_function

Real _A
Multiplying factor in k = k_ijk * A * phi^n / (1 - phi)^m. More...

MaterialProperty< RealTensorValue > & _permeability_qp

MaterialProperty< std::vector< RealTensorValue > > & _dpermeability_qp_dvar

MaterialProperty< std::vector< std::vector< RealTensorValue > > > & _dpermeability_qp_dgradvar

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 a form of the Kozeny-Carman equation (e.g.

Oelkers 1996: Reviews in Mineralogy v. 34, p. 131-192): k = k_ijk * A * phi^n / (1 - phi)^m where k_ijk is a tensor providing the anisotropy, phi is porosity, n and m are positive scalar constants and A is given in one of the following forms: A = k0 * (1 - phi0)^m / phi0^n where k0 and phi0 are a reference permeability and porosity, or A = f * d^2 where f is a scalar constant and d is grain diameter.

Definition at line 34 of file PorousFlowPermeabilityKozenyCarman.h.

◆ PoropermFunction

 strongprotected

Name of porosity-permeability relationship.

Enumerator
kozeny_carman_fd2
kozeny_carman_phi0

Definition at line 73 of file PorousFlowPermeabilityKozenyCarman.h.

73 { kozeny_carman_fd2, kozeny_carman_phi0 } _poroperm_function;
enum PorousFlowPermeabilityKozenyCarman::PoropermFunction _poroperm_function

◆ PorousFlowPermeabilityKozenyCarman()

 PorousFlowPermeabilityKozenyCarman::PorousFlowPermeabilityKozenyCarman ( const InputParameters & parameters )

Definition at line 54 of file PorousFlowPermeabilityKozenyCarman.C.

56  : PorousFlowPermeabilityBase(parameters),
57  _k0(parameters.isParamValid("k0") ? getParam<Real>("k0") : -1),
58  _phi0(parameters.isParamValid("phi0") ? getParam<Real>("phi0") : -1),
59  _f(parameters.isParamValid("f") ? getParam<Real>("f") : -1),
60  _d(parameters.isParamValid("d") ? getParam<Real>("d") : -1),
61  _m(getParam<Real>("m")),
62  _n(getParam<Real>("n")),
63  _k_anisotropy(parameters.isParamValid("k_anisotropy")
64  ? getParam<RealTensorValue>("k_anisotropy")
65  : RealTensorValue(1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0)),
66  _porosity_qp(getMaterialProperty<Real>("PorousFlow_porosity_qp")),
67  _dporosity_qp_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_qp_dvar")),
70  _poroperm_function(getParam<MooseEnum>("poroperm_function").getEnum<PoropermFunction>())
71 {
72  switch (_poroperm_function)
73  {
75  if (!(parameters.isParamValid("f") && parameters.isParamValid("d")))
76  mooseError("You must specify f and d in order to use kozeny_carman_fd2 in "
77  "PorousFlowPermeabilityKozenyCarman");
78  _A = _f * _d * _d;
79  break;
80
82  if (!(parameters.isParamValid("k0") && parameters.isParamValid("phi0")))
83  mooseError("You must specify k0 and phi0 in order to use kozeny_carman_phi0 in "
84  "PorousFlowPermeabilityKozenyCarman");
85  _A = _k0 * std::pow(1.0 - _phi0, _m) / std::pow(_phi0, _n);
86  break;
87  }
88 }
const Real _k0
Reference scalar permeability in A = k0 * (1 - phi0)^m / phi0^n.
Real _A
Multiplying factor in k = k_ijk * A * phi^n / (1 - phi)^m.
const Real _phi0
Reference porosity in A = k0 * (1 - phi0)^m / phi0^n.
const MaterialProperty< Real > & _porosity_qp
const RealTensorValue _k_anisotropy
Tensor multiplier k_ijk in k = k_ijk * A * phi^n / (1 - phi)^m.
const Real _f
Multiplying factor in A = f * d^2.
PorousFlowPermeabilityBase(const InputParameters &parameters)
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
const Real _d
Grain diameter in A = f * d^2.
PoropermFunction
Name of porosity-permeability relationship.
const Real _n
Exponent in k = k_ijk * A * phi^n / (1 - phi)^m.
enum PorousFlowPermeabilityKozenyCarman::PoropermFunction _poroperm_function
const MaterialProperty< std::vector< Real > > & _dporosity_qp_dvar
const Real _m
Exponent in k = k_ijk * A * phi^n / (1 - phi)^m.

◆ computeQpProperties()

 void PorousFlowPermeabilityKozenyCarman::computeQpProperties ( )
overrideprotected

Definition at line 91 of file PorousFlowPermeabilityKozenyCarman.C.

92 {
93  _permeability_qp[_qp] =
94  _k_anisotropy * _A * std::pow(_porosity_qp[_qp], _n) / std::pow(1.0 - _porosity_qp[_qp], _m);
95
96  _dpermeability_qp_dvar[_qp].resize(_num_var, RealTensorValue());
97  for (unsigned int v = 0; v < _num_var; ++v)
99  (_n / _porosity_qp[_qp] + _m / (1.0 - _porosity_qp[_qp]));
100
102  for (unsigned i = 0; i < LIBMESH_DIM; ++i)
103  {
105  for (unsigned int v = 0; v < _num_var; ++v)
108  (_n / _porosity_qp[_qp] + _m / (1.0 - _porosity_qp[_qp]));
109  }
110 }
Real _A
Multiplying factor in k = k_ijk * A * phi^n / (1 - phi)^m.
MaterialProperty< std::vector< std::vector< RealTensorValue > > > & _dpermeability_qp_dgradvar
const MaterialProperty< Real > & _porosity_qp
const RealTensorValue _k_anisotropy
Tensor multiplier k_ijk in k = k_ijk * A * phi^n / (1 - phi)^m.
const unsigned int _num_var
Number of PorousFlow variables.
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
MaterialProperty< std::vector< RealTensorValue > > & _dpermeability_qp_dvar
const Real _n
Exponent in k = k_ijk * A * phi^n / (1 - phi)^m.
MaterialProperty< RealTensorValue > & _permeability_qp
const MaterialProperty< std::vector< Real > > & _dporosity_qp_dvar
const Real _m
Exponent in k = k_ijk * A * phi^n / (1 - phi)^m.

◆ _A

 Real PorousFlowPermeabilityKozenyCarman::_A
protected

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

Definition at line 76 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by computeQpProperties(), and PorousFlowPermeabilityKozenyCarman().

◆ _d

 const Real PorousFlowPermeabilityKozenyCarman::_d
protected

Grain diameter in A = f * d^2.

Definition at line 52 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by PorousFlowPermeabilityKozenyCarman().

protectedinherited

Definition at line 36 of file PorousFlowPermeabilityBase.h.

◆ _dpermeability_qp_dvar

 MaterialProperty >& PorousFlowPermeabilityBase::_dpermeability_qp_dvar
protectedinherited

Definition at line 33 of file PorousFlowPermeabilityBase.h.

protected

Definition at line 70 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by computeQpProperties().

◆ _dporosity_qp_dvar

 const MaterialProperty >& PorousFlowPermeabilityKozenyCarman::_dporosity_qp_dvar
protected

Definition at line 67 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by computeQpProperties().

◆ _f

 const Real PorousFlowPermeabilityKozenyCarman::_f
protected

Multiplying factor in A = f * d^2.

Definition at line 49 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by PorousFlowPermeabilityKozenyCarman().

◆ _k0

 const Real PorousFlowPermeabilityKozenyCarman::_k0
protected

Reference scalar permeability in A = k0 * (1 - phi0)^m / phi0^n.

Definition at line 43 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by PorousFlowPermeabilityKozenyCarman().

◆ _k_anisotropy

 const RealTensorValue PorousFlowPermeabilityKozenyCarman::_k_anisotropy
protected

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

Definition at line 61 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by computeQpProperties().

◆ _m

 const Real PorousFlowPermeabilityKozenyCarman::_m
protected

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

Definition at line 55 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by computeQpProperties(), and PorousFlowPermeabilityKozenyCarman().

◆ _n

 const Real PorousFlowPermeabilityKozenyCarman::_n
protected

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

Definition at line 58 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by computeQpProperties(), and PorousFlowPermeabilityKozenyCarman().

◆ _num_components

 const unsigned int PorousFlowMaterialVectorBase::_num_components
protectedinherited

◆ _permeability_qp

 MaterialProperty& PorousFlowPermeabilityBase::_permeability_qp
protectedinherited

Definition at line 30 of file PorousFlowPermeabilityBase.h.

◆ _phi0

 const Real PorousFlowPermeabilityKozenyCarman::_phi0
protected

Reference porosity in A = k0 * (1 - phi0)^m / phi0^n.

Definition at line 46 of file PorousFlowPermeabilityKozenyCarman.h.

Referenced by PorousFlowPermeabilityKozenyCarman().

◆ _poroperm_function

 enum PorousFlowPermeabilityKozenyCarman::PoropermFunction PorousFlowPermeabilityKozenyCarman::_poroperm_function
protected

◆ _porosity_qp

 const MaterialProperty& PorousFlowPermeabilityKozenyCarman::_porosity_qp
protected