www.mooseframework.org
PoroFullSatTimeDerivative.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 POROFULLSATTIMEDERIVATIVE_H
9 #define POROFULLSATTIMEDERIVATIVE_H
10 
11 #include "DerivativeMaterialInterface.h"
12 #include "TimeDerivative.h"
13 
14 // Forward Declarations
16 
17 template <>
19 
25 class PoroFullSatTimeDerivative : public DerivativeMaterialInterface<TimeDerivative>
26 {
27 public:
28  PoroFullSatTimeDerivative(const InputParameters & parameters);
29 
30 protected:
31  virtual Real computeQpResidual();
32 
33  virtual Real computeQpJacobian();
34 
35  virtual Real computeQpOffDiagJacobian(unsigned int jvar);
36 
37 private:
39  const VariableValue & _u_old;
40 
42  const MaterialProperty<Real> & _volstrain;
43 
45  const MaterialProperty<Real> & _volstrain_old;
46 
48  unsigned int _ndisp;
49 
51  std::vector<unsigned int> _disp_var_num;
52 
54  const MaterialProperty<Real> & _alpha;
55 
57  const MaterialProperty<Real> & _one_over_biot_modulus;
58 
60  const MaterialProperty<Real> & _done_over_biot_modulus_dP;
61 
63  const MaterialProperty<Real> & _done_over_biot_modulus_dep;
64 };
65 
66 #endif // POROFULLSATTIMEDERIVATIVE_H
const MaterialProperty< Real > & _done_over_biot_modulus_dep
d(1/M)/d(volumetric strain)
const VariableValue & _u_old
old value of porepressure
std::vector< unsigned int > _disp_var_num
variable number of the displacements variables
const MaterialProperty< Real > & _volstrain
volumetric strain
const MaterialProperty< Real > & _alpha
Biot coefficient.
const MaterialProperty< Real > & _volstrain_old
old value of volumetric strain
virtual Real computeQpOffDiagJacobian(unsigned int jvar)
const MaterialProperty< Real > & _one_over_biot_modulus
1/M, where M is the Biot modulus
Kernel = biot_coefficient*d(volumetric_strain)/dt + (1/biot_modulus)*d(porepressure)/dt this is the t...
PoroFullSatTimeDerivative(const InputParameters &parameters)
InputParameters validParams< PoroFullSatTimeDerivative >()
const MaterialProperty< Real > & _done_over_biot_modulus_dP
d(1/M)/d(porepressure)
unsigned int _ndisp
number of displacement variables