27#ifndef OPM_VAN_GENUCHTEN_HPP
28#define OPM_VAN_GENUCHTEN_HPP
54template <
class TraitsT,
class ParamsT = VanGenuchtenParams<TraitsT> >
65 typedef typename Traits::Scalar
Scalar;
70 "The van Genuchten capillary pressure law only "
71 "applies to the case of two fluid phases");
113 template <
class Container,
class Flu
idState>
116 typedef typename std::remove_reference<
decltype(values[0])>::type Evaluation;
118 values[Traits::wettingPhaseIdx] = 0.0;
119 values[Traits::nonWettingPhaseIdx] = pcnw<FluidState, Evaluation>(params, fs);
126 template <
class Container,
class Flu
idState>
129 typedef typename std::remove_reference<
decltype(values[0])>::type Evaluation;
131 values[Traits::wettingPhaseIdx] = Sw<FluidState, Evaluation>(params, fs);
132 values[Traits::nonWettingPhaseIdx] = 1 - values[Traits::wettingPhaseIdx];
145 template <
class Container,
class Flu
idState>
148 typedef typename std::remove_reference<
decltype(values[0])>::type Evaluation;
150 values[Traits::wettingPhaseIdx] = krw<FluidState, Evaluation>(params, fs);
151 values[Traits::nonWettingPhaseIdx] = krn<FluidState, Evaluation>(params, fs);
168 template <
class Flu
idState,
class Evaluation =
typename Flu
idState::Scalar>
169 static Evaluation
pcnw(
const Params& params,
const FluidState& fs)
171 const Evaluation&
Sw =
172 decay<Evaluation>(fs.saturation(Traits::wettingPhaseIdx));
174 assert(0 <=
Sw &&
Sw <= 1);
193 template <
class Evaluation>
196 return pow(pow(
Sw, -1.0/params.vgM()) - 1, 1.0/params.vgN())/params.vgAlpha();
211 template <
class Flu
idState,
class Evaluation =
typename Flu
idState::Scalar>
212 static Evaluation
Sw(
const Params& params,
const FluidState& fs)
215 decay<Evaluation>(fs.pressure(Traits::nonWettingPhaseIdx))
216 - decay<Evaluation>(fs.pressure(Traits::wettingPhaseIdx));
217 return twoPhaseSatSw(params, pC);
220 template <
class Evaluation>
221 static Evaluation twoPhaseSatSw(
const Params& params,
const Evaluation& pC)
225 return pow(pow(params.vgAlpha()*pC, params.vgN()) + 1, -params.vgM());
232 template <
class Flu
idState,
class Evaluation =
typename Flu
idState::Scalar>
233 static Evaluation
Sn(
const Params& params,
const FluidState& fs)
234 {
return 1 - Sw<FluidState, Evaluation>(params, fs); }
236 template <
class Evaluation>
237 static Evaluation twoPhaseSatSn(
const Params& params,
const Evaluation& pC)
238 {
return 1 - twoPhaseSatSw(params, pC); }
250 template <
class Flu
idState,
class Evaluation =
typename Flu
idState::Scalar>
251 static Evaluation
krw(
const Params& params,
const FluidState& fs)
253 const Evaluation&
Sw =
254 decay<Evaluation>(fs.saturation(Traits::wettingPhaseIdx));
256 return twoPhaseSatKrw(params,
Sw);
259 template <
class Evaluation>
260 static Evaluation twoPhaseSatKrw(
const Params& params,
const Evaluation&
Sw)
262 assert(0.0 <=
Sw &&
Sw <= 1.0);
264 Evaluation r = 1.0 - pow(1.0 - pow(
Sw, 1/params.vgM()), params.vgM());
277 template <
class Flu
idState,
class Evaluation =
typename Flu
idState::Scalar>
278 static Evaluation
krn(
const Params& params,
const FluidState& fs)
280 const Evaluation&
Sw =
281 1.0 - decay<Evaluation>(fs.saturation(Traits::nonWettingPhaseIdx));
283 return twoPhaseSatKrn(params,
Sw);
286 template <
class Evaluation>
287 static Evaluation twoPhaseSatKrn(
const Params& params, Evaluation
Sw)
289 assert(0 <=
Sw &&
Sw <= 1);
293 pow(1 - pow(
Sw, 1/params.vgM()), 2*params.vgM());
Specification of the material parameters for the van Genuchten constitutive relations.
Implementation of the van Genuchten capillary pressure - saturation relation.
Definition VanGenuchten.hpp:56
static const int numPhases
The number of fluid phases.
Definition VanGenuchten.hpp:68
static Evaluation Sw(const Params ¶ms, const FluidState &fs)
The saturation-capillary pressure curve according to van Genuchten.
Definition VanGenuchten.hpp:212
ParamsT Params
The type of the parameter objects for this law.
Definition VanGenuchten.hpp:62
static const bool isSaturationDependent
Specify whether the quantities defined by this material law are saturation dependent.
Definition VanGenuchten.hpp:83
static Evaluation twoPhaseSatPcnw(const Params ¶ms, const Evaluation &Sw)
The saturation-capillary pressure curve according to van Genuchten using a material law specific API.
Definition VanGenuchten.hpp:194
static Evaluation pcnw(const Params ¶ms, const FluidState &fs)
The capillary pressure-saturation curve according to van Genuchten.
Definition VanGenuchten.hpp:169
static void saturations(Container &values, const Params ¶ms, const FluidState &fs)
Calculate the saturations of the phases starting from their pressure differences.
Definition VanGenuchten.hpp:127
static const bool implementsTwoPhaseSatApi
Specify whether this material law implements the two-phase convenience API which only depends on the ...
Definition VanGenuchten.hpp:79
static Evaluation Sn(const Params ¶ms, const FluidState &fs)
Calculate the non-wetting phase saturations depending on the phase pressures.
Definition VanGenuchten.hpp:233
static const bool isPressureDependent
Specify whether the quantities defined by this material law are dependent on the absolute pressure.
Definition VanGenuchten.hpp:87
Traits::Scalar Scalar
The type of the scalar values for this law.
Definition VanGenuchten.hpp:65
static void capillaryPressures(Container &values, const Params ¶ms, const FluidState &fs)
The capillary pressure-saturation curves according to van Genuchten.
Definition VanGenuchten.hpp:114
static const bool isCompositionDependent
Specify whether the quantities defined by this material law are dependent on the phase composition.
Definition VanGenuchten.hpp:95
static Evaluation krn(const Params ¶ms, const FluidState &fs)
The relative permeability for the non-wetting phase of the medium according to van Genuchten.
Definition VanGenuchten.hpp:278
static void relativePermeabilities(Container &values, const Params ¶ms, const FluidState &fs)
The relative permeability-saturation curves according to van Genuchten.
Definition VanGenuchten.hpp:146
static Evaluation krw(const Params ¶ms, const FluidState &fs)
The relative permeability for the wetting phase of the medium according to van Genuchten's curve with...
Definition VanGenuchten.hpp:251
static const bool isTemperatureDependent
Specify whether the quantities defined by this material law are temperature dependent.
Definition VanGenuchten.hpp:91
TraitsT Traits
The traits class for this material law.
Definition VanGenuchten.hpp:59
static const bool implementsTwoPhaseApi
Specify whether this material law implements the two-phase convenience API.
Definition VanGenuchten.hpp:75
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30