www.mooseframework.org
PorousFlowPermeabilityExponential.h
Go to the documentation of this file.
1 /****************************************************************/
2 /* MOOSE - Multiphysics Object Oriented Simulation Environment */
3 /* */
4 /* All contents are licensed under LGPL V2.1 */
5 /* See LICENSE for full restrictions */
6 /****************************************************************/
7 
8 #ifndef POROUSFLOWPERMEABILITYEXPONENTIAL_H
9 #define POROUSFLOWPERMEABILITYEXPONENTIAL_H
10 
12 
13 // Forward Declarations
15 
16 template <>
18 
33 {
34 public:
35  PorousFlowPermeabilityExponential(const InputParameters & parameters);
36 
37 protected:
38  void computeQpProperties() override;
39 
41  const Real _A;
42 
44  const Real _B;
45 
47  const RealTensorValue _k_anisotropy;
48 
50  const MaterialProperty<Real> & _porosity_qp;
51 
53  const MaterialProperty<std::vector<Real>> & _dporosity_qp_dvar;
54 
56  const MaterialProperty<std::vector<RealGradient>> & _dporosity_qp_dgradvar;
57 
60 
62  Real _AA;
63 
65  Real _BB;
66 };
67 
68 #endif // POROUSFLOWPERMEABILITYEXPONENTIAL_H
Real _AA
Empirical constant AA in k = k_ijk * BB * exp(AA * phi)
enum PorousFlowPermeabilityExponential::PoropermFunction _poroperm_function
Material designed to provide the permeability tensor which is calculated from porosity using the equa...
PorousFlowPermeabilityExponential(const InputParameters &parameters)
Base class Material designed to provide the permeability tensor.
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)
InputParameters validParams< PorousFlowPermeabilityExponential >()
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.