www.mooseframework.org
PorousFlowPorosityHM.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 POROUSFLOWPOROSITYHM_H
9 #define POROUSFLOWPOROSITYHM_H
10 
12 
13 // Forward Declarations
15 
16 template <>
17 InputParameters validParams<PorousFlowPorosityHM>();
18 
25 {
26 public:
27  PorousFlowPorosityHM(const InputParameters & parameters);
28 
29 protected:
30  virtual Real atNegInfinityQp() const override;
31  virtual Real atZeroQp() const override;
32  virtual Real decayQp() const override;
33  virtual Real ddecayQp_dvar(unsigned pvar) const override;
34  virtual RealGradient ddecayQp_dgradvar(unsigned pvar) const override;
35 
37  const VariableValue & _phi0;
38 
40  const Real _biot;
41 
43  const Real _solid_bulk;
44 
46  const Real _coeff;
47 
49  const VariableValue & _p_reference;
50 
52  const unsigned int _ndisp;
53 
55  std::vector<unsigned int> _disp_var_num;
56 
58  const MaterialProperty<Real> & _vol_strain_qp;
59 
61  const MaterialProperty<std::vector<RealGradient>> & _dvol_strain_qp_dvar;
62 
64  const MaterialProperty<Real> & _pf;
65 
67  const MaterialProperty<std::vector<Real>> & _dpf_dvar;
68 };
69 
70 #endif // POROUSFLOWPOROSITYHM_H
virtual Real decayQp() const override
Returns "decay" at the quadpoint (porosity = a + (b - a) * exp(decay))
const Real _solid_bulk
drained bulk modulus of the porous skeleton
virtual Real atNegInfinityQp() const override
Returns "a" at the quadpoint (porosity = a + (b - a) * exp(decay))
Base class Material designed to provide the porosity.
std::vector< unsigned int > _disp_var_num
variable number of the displacements variables
const VariableValue & _phi0
porosity at zero strain and zero porepressure
const VariableValue & _p_reference
reference porepressure
Material designed to provide the porosity in hydro-mechanical simulations biot + (phi0 - biot)*exp(-v...
InputParameters validParams< PorousFlowPorosityHM >()
const MaterialProperty< std::vector< Real > > & _dpf_dvar
d(effective porepressure)/(d porflow variable)
const Real _coeff
short-hand number (biot-1)/solid_bulk
const unsigned int _ndisp
number of displacement variables
virtual RealGradient ddecayQp_dgradvar(unsigned pvar) const override
d(decay)/d(grad(porous-flow variable pvar))
const Real _biot
biot coefficient
const MaterialProperty< std::vector< RealGradient > > & _dvol_strain_qp_dvar
d(strain)/(dvar)
const MaterialProperty< Real > & _vol_strain_qp
strain
virtual Real atZeroQp() const override
Returns "b" at the quadpoint (porosity = a + (b - a) * exp(decay))
virtual Real ddecayQp_dvar(unsigned pvar) const override
d(decay)/d(porous-flow variable pvar)
PorousFlowPorosityHM(const InputParameters &parameters)
const MaterialProperty< Real > & _pf
effective porepressure