www.mooseframework.org
RichardsDensityConstBulkCut.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 RICHARDSDENSITYCONSTBULKCUT_H
9 #define RICHARDSDENSITYCONSTBULKCUT_H
10 
11 #include "RichardsDensity.h"
12 
14 
15 template <>
17 
25 {
26 public:
27  RichardsDensityConstBulkCut(const InputParameters & parameters);
28 
33  Real density(Real p) const;
34 
39  Real ddensity(Real p) const;
40 
45  Real d2density(Real p) const;
46 
47 protected:
49  Real _dens0;
50 
52  Real _bulk;
53 
55  Real _cut_limit;
56 
59 
61  Real _c3;
62 };
63 
64 #endif // RICHARDSDENSITYCONSTBULKCUT_H
Real density(Real p) const
fluid density as a function of porepressure
InputParameters validParams< RichardsDensityConstBulkCut >()
Real _cut_limit
where the cubic starts
Real ddensity(Real p) const
derivative of fluid density wrt porepressure
Real _dens0
density = _dens0*exp(p/_bulk), modified by cubic
Real _bulk
density = _dens0*exp(p/_bulk), modified by cubic
Real _c3
(cut_limit-zero_point)^3
Base class for fluid density as a function of porepressure The functions density, ddensity and d2dens...
Fluid density assuming constant bulk modulus, for p>cut_limit Then following a cubic for zero_point <...
Real d2density(Real p) const
second derivative of fluid density wrt porepressure
Real _zero_point
where the density is zero
RichardsDensityConstBulkCut(const InputParameters &parameters)