https://mooseframework.inl.gov
Loading...
Searching...
No Matches
PorousFlowCapillaryPressureVG.h
Go to the documentation of this file.
1//* This file is part of the MOOSE framework
2//* https://mooseframework.inl.gov
3//*
4//* All rights reserved, see COPYRIGHT for full restrictions
5//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
6//*
7//* Licensed under LGPL 2.1, please see LICENSE for details
8//* https://www.gnu.org/licenses/lgpl-2.1.html
9
10#pragma once
11
13
21{
22public:
24
26
27 virtual Real capillaryPressureCurve(Real saturation, unsigned qp = 0) const override;
28 virtual Real dCapillaryPressureCurve(Real saturation, unsigned qp = 0) const override;
29 virtual Real d2CapillaryPressureCurve(Real saturation, unsigned qp = 0) const override;
30
31 virtual Real effectiveSaturation(Real pc, unsigned qp = 0) const override;
32 virtual Real dEffectiveSaturation(Real pc, unsigned qp = 0) const override;
33 virtual Real d2EffectiveSaturation(Real pc, unsigned qp = 0) const override;
34
35protected:
37 const Real _m;
39 const Real _alpha;
41 const Real _s_scale;
43 const Real _pc_sscale;
44};
const InputParameters & parameters() const
van Genuchten form of capillary pressure.
virtual Real dCapillaryPressureCurve(Real saturation, unsigned qp=0) const override
Derivative of raw capillary pressure wrt true saturation.
virtual Real d2EffectiveSaturation(Real pc, unsigned qp=0) const override
Second derivative of effective saturation wrt capillary pressure.
virtual Real d2CapillaryPressureCurve(Real saturation, unsigned qp=0) const override
Second derivative of raw capillary pressure wrt true saturation.
const Real _m
van Genuchten exponent m
virtual Real dEffectiveSaturation(Real pc, unsigned qp=0) const override
Derivative of effective saturation wrt capillary pressure.
virtual Real capillaryPressureCurve(Real saturation, unsigned qp=0) const override
Raw capillary pressure curve (does not include logarithmic extension)
virtual Real effectiveSaturation(Real pc, unsigned qp=0) const override
Effective saturation as a function of capillary pressure.
const Real _s_scale
Capillary pressure = f(Seff * s_scale) - pc_sscale, where f is the van Genuchten function....
const Real _alpha
van Genuchten capillary coefficient alpha
const Real _pc_sscale
pc_sscale = f(s_scale), where f is the van Genuchten function
Base class for capillary pressure for multiphase flow in porous media.
Real saturation(Real pc, unsigned qp=0) const
Saturation as a function of capillary pressure.