Loading...
Searching...
No Matches
EclMultiplexerMaterial.hpp
Go to the documentation of this file.
// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*
This file is part of the Open Porous Media project (OPM).
OPM is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 2 of the License, or
(at your option) any later version.
OPM is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with OPM. If not, see <http://www.gnu.org/licenses/>.
Consult the COPYING file in the top-level source directory of this
module for the precise wording of the license and the list of
copyright holders.
*/
#ifndef OPM_ECL_MULTIPLEXER_MATERIAL_HPP
#define OPM_ECL_MULTIPLEXER_MATERIAL_HPP
#include "EclMultiplexerMaterialParams.hpp"
#include "EclDefaultMaterial.hpp"
#include "EclStone1Material.hpp"
#include "EclStone2Material.hpp"
#include "EclTwoPhaseMaterial.hpp"
#include <opm/common/TimingMacros.hpp>
#include <algorithm>
#include <stdexcept>
namespace Opm {
template <class TraitsT,
class GasOilMaterialLawT,
class OilWaterMaterialLawT,
class GasWaterMaterialLawT,
class ParamsT = EclMultiplexerMaterialParams<TraitsT,
GasOilMaterialLawT,
OilWaterMaterialLawT,
GasWaterMaterialLawT> >
class EclMultiplexerMaterial : public TraitsT
{
public:
using GasOilMaterialLaw = GasOilMaterialLawT;
using OilWaterMaterialLaw = OilWaterMaterialLawT;
using GasWaterMaterialLaw = GasWaterMaterialLawT;
using Stone1Material = EclStone1Material<TraitsT, GasOilMaterialLaw, OilWaterMaterialLaw>;
using Stone2Material = EclStone2Material<TraitsT, GasOilMaterialLaw, OilWaterMaterialLaw>;
using DefaultMaterial = EclDefaultMaterial<TraitsT, GasOilMaterialLaw, OilWaterMaterialLaw>;
using TwoPhaseMaterial = EclTwoPhaseMaterial<TraitsT, GasOilMaterialLaw, OilWaterMaterialLaw, GasWaterMaterialLaw>;
// some safety checks
static_assert(TraitsT::numPhases == 3,
"The number of phases considered by this capillary pressure "
"law is always three!");
static_assert(GasOilMaterialLaw::numPhases == 2,
"The number of phases considered by the gas-oil capillary "
"pressure law must be two!");
static_assert(OilWaterMaterialLaw::numPhases == 2,
"The number of phases considered by the oil-water capillary "
"pressure law must be two!");
static_assert(GasWaterMaterialLaw::numPhases == 2,
"The number of phases considered by the gas-water capillary "
"pressure law must be two!");
static_assert(std::is_same<typename GasOilMaterialLaw::Scalar,
typename OilWaterMaterialLaw::Scalar>::value,
"The two two-phase capillary pressure laws must use the same "
"type of floating point values.");
using Traits = TraitsT;
using Params = ParamsT;
using Scalar = typename Traits::Scalar;
static constexpr int numPhases = 3;
static constexpr int waterPhaseIdx = Traits::wettingPhaseIdx;
static constexpr int oilPhaseIdx = Traits::nonWettingPhaseIdx;
static constexpr int gasPhaseIdx = Traits::gasPhaseIdx;
static constexpr bool implementsTwoPhaseApi = false;
static constexpr bool implementsTwoPhaseSatApi = false;
static constexpr bool isSaturationDependent = true;
static constexpr bool isPressureDependent = false;
static constexpr bool isTemperatureDependent = false;
static constexpr bool isCompositionDependent = false;
template <class ContainerT, class FluidState>
static void capillaryPressures(ContainerT& values,
const Params& params,
const FluidState& fluidState)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
Stone1Material::capillaryPressures(values,
params.template getRealParams<EclMultiplexerApproach::Stone1>(),
fluidState);
break;
case EclMultiplexerApproach::Stone2:
Stone2Material::capillaryPressures(values,
params.template getRealParams<EclMultiplexerApproach::Stone2>(),
fluidState);
break;
case EclMultiplexerApproach::Default:
DefaultMaterial::capillaryPressures(values,
params.template getRealParams<EclMultiplexerApproach::Default>(),
fluidState);
break;
case EclMultiplexerApproach::TwoPhase:
TwoPhaseMaterial::capillaryPressures(values,
params.template getRealParams<EclMultiplexerApproach::TwoPhase>(),
fluidState);
break;
case EclMultiplexerApproach::OnePhase:
values[0] = 0.0;
break;
}
}
/*
* Hysteresis parameters for oil-water
* @see EclHysteresisTwoPhaseLawParams::soMax(...)
* @see EclHysteresisTwoPhaseLawParams::swMax(...)
* @see EclHysteresisTwoPhaseLawParams::swMin(...)
* \param params Parameters
*/
static void oilWaterHysteresisParams(Scalar& soMax,
Scalar& swMax,
Scalar& swMin,
const Params& params)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
Stone1Material::oilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::Stone1>());
break;
case EclMultiplexerApproach::Stone2:
Stone2Material::oilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::Stone2>());
break;
case EclMultiplexerApproach::Default:
DefaultMaterial::oilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::Default>());
break;
case EclMultiplexerApproach::TwoPhase:
TwoPhaseMaterial::oilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::TwoPhase>());
break;
case EclMultiplexerApproach::OnePhase:
// Do nothing.
break;
}
}
/*
* Hysteresis parameters for oil-water
* @see EclHysteresisTwoPhaseLawParams::soMax(...)
* @see EclHysteresisTwoPhaseLawParams::swMax(...)
* @see EclHysteresisTwoPhaseLawParams::swMin(...)
* \param params Parameters
*/
static void setOilWaterHysteresisParams(const Scalar& soMax,
const Scalar& swMax,
const Scalar& swMin,
Params& params)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
Stone1Material::setOilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::Stone1>());
break;
case EclMultiplexerApproach::Stone2:
Stone2Material::setOilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::Stone2>());
break;
case EclMultiplexerApproach::Default:
DefaultMaterial::setOilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::Default>());
break;
case EclMultiplexerApproach::TwoPhase:
TwoPhaseMaterial::setOilWaterHysteresisParams(soMax, swMax, swMin,
params.template getRealParams<EclMultiplexerApproach::TwoPhase>());
break;
case EclMultiplexerApproach::OnePhase:
// Do nothing.
break;
}
}
/*
* Hysteresis parameters for gas-oil
* @see EclHysteresisTwoPhaseLawParams::sgmax(...)
* @see EclHysteresisTwoPhaseLawParams::shmax(...)
* @see EclHysteresisTwoPhaseLawParams::somin(...)
* \param params Parameters
*/
static void gasOilHysteresisParams(Scalar& sgmax,
Scalar& shmax,
Scalar& somin,
const Params& params)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
Stone1Material::gasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::Stone1>());
break;
case EclMultiplexerApproach::Stone2:
Stone2Material::gasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::Stone2>());
break;
case EclMultiplexerApproach::Default:
DefaultMaterial::gasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::Default>());
break;
case EclMultiplexerApproach::TwoPhase:
TwoPhaseMaterial::gasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::TwoPhase>());
break;
case EclMultiplexerApproach::OnePhase:
// Do nothing.
break;
}
}
static Scalar trappedGasSaturation(const Params& params, bool maximumTrapping)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::trappedGasSaturation(params.template getRealParams<EclMultiplexerApproach::Stone1>(), maximumTrapping);
case EclMultiplexerApproach::Stone2:
return Stone2Material::trappedGasSaturation(params.template getRealParams<EclMultiplexerApproach::Stone2>(), maximumTrapping);
case EclMultiplexerApproach::Default:
return DefaultMaterial::trappedGasSaturation(params.template getRealParams<EclMultiplexerApproach::Default>(), maximumTrapping);
case EclMultiplexerApproach::TwoPhase:
return TwoPhaseMaterial::trappedGasSaturation(params.template getRealParams<EclMultiplexerApproach::TwoPhase>(), maximumTrapping);
case EclMultiplexerApproach::OnePhase:
return 0.0;
}
return 0.0;
}
static Scalar strandedGasSaturation(const Params& params, Scalar Sg, Scalar Kg)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::strandedGasSaturation(params.template getRealParams<EclMultiplexerApproach::Stone1>(), Sg, Kg);
case EclMultiplexerApproach::Stone2:
return Stone2Material::strandedGasSaturation(params.template getRealParams<EclMultiplexerApproach::Stone2>(), Sg, Kg);
case EclMultiplexerApproach::Default:
return DefaultMaterial::strandedGasSaturation(params.template getRealParams<EclMultiplexerApproach::Default>(), Sg, Kg);
case EclMultiplexerApproach::TwoPhase:
return TwoPhaseMaterial::strandedGasSaturation(params.template getRealParams<EclMultiplexerApproach::TwoPhase>(), Sg, Kg);
case EclMultiplexerApproach::OnePhase:
return 0.0;
}
return 0.0;
}
static Scalar trappedOilSaturation(const Params& params, bool maximumTrapping)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::trappedOilSaturation(params.template getRealParams<EclMultiplexerApproach::Stone1>(), maximumTrapping);
case EclMultiplexerApproach::Stone2:
return Stone2Material::trappedOilSaturation(params.template getRealParams<EclMultiplexerApproach::Stone2>(), maximumTrapping);
case EclMultiplexerApproach::Default:
return DefaultMaterial::trappedOilSaturation(params.template getRealParams<EclMultiplexerApproach::Default>(), maximumTrapping);
case EclMultiplexerApproach::TwoPhase:
return TwoPhaseMaterial::trappedOilSaturation(params.template getRealParams<EclMultiplexerApproach::TwoPhase>(), maximumTrapping);
case EclMultiplexerApproach::OnePhase:
return 0.0;
}
return 0.0;
}
static Scalar trappedWaterSaturation(const Params& params)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::trappedWaterSaturation(params.template getRealParams<EclMultiplexerApproach::Stone1>());
case EclMultiplexerApproach::Stone2:
return Stone2Material::trappedWaterSaturation(params.template getRealParams<EclMultiplexerApproach::Stone2>());
case EclMultiplexerApproach::Default:
return DefaultMaterial::trappedWaterSaturation(params.template getRealParams<EclMultiplexerApproach::Default>());
case EclMultiplexerApproach::TwoPhase:
return TwoPhaseMaterial::trappedWaterSaturation(params.template getRealParams<EclMultiplexerApproach::TwoPhase>());
case EclMultiplexerApproach::OnePhase:
return 0.0;
}
return 0.0;
}
/*
* Hysteresis parameters for gas-oil
* @see EclHysteresisTwoPhaseLawParams::sgmax(...)
* @see EclHysteresisTwoPhaseLawParams::shmax(...)
* @see EclHysteresisTwoPhaseLawParams::somin(...)
* \param params Parameters
*/
static void setGasOilHysteresisParams(const Scalar& sgmax,
const Scalar& shmax,
const Scalar& somin,
Params& params)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
Stone1Material::setGasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::Stone1>());
break;
case EclMultiplexerApproach::Stone2:
Stone2Material::setGasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::Stone2>());
break;
case EclMultiplexerApproach::Default:
DefaultMaterial::setGasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::Default>());
break;
case EclMultiplexerApproach::TwoPhase:
TwoPhaseMaterial::setGasOilHysteresisParams(sgmax, shmax, somin,
params.template getRealParams<EclMultiplexerApproach::TwoPhase>());
break;
case EclMultiplexerApproach::OnePhase:
// Do nothing.
break;
}
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation pcgn(const Params& /* params */,
const FluidState& /* fs */)
{
throw std::logic_error("Not implemented: pcgn()");
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation pcnw(const Params& /* params */,
const FluidState& /* fs */)
{
throw std::logic_error("Not implemented: pcnw()");
}
template <class ContainerT, class FluidState>
static void saturations(ContainerT& /* values */,
const Params& /* params */,
const FluidState& /* fs */)
{
throw std::logic_error("Not implemented: saturations()");
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation Sg(const Params& /* params */,
const FluidState& /* fluidState */)
{
throw std::logic_error("Not implemented: Sg()");
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation Sn(const Params& /* params */,
const FluidState& /* fluidState */)
{
throw std::logic_error("Not implemented: Sn()");
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation Sw(const Params& /* params */,
const FluidState& /* fluidState */)
{
throw std::logic_error("Not implemented: Sw()");
}
template <class ContainerT, class FluidState>
static void relativePermeabilities(ContainerT& values,
const Params& params,
const FluidState& fluidState)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
Stone1Material::relativePermeabilities(values,
params.template getRealParams<EclMultiplexerApproach::Stone1>(),
fluidState);
break;
case EclMultiplexerApproach::Stone2:
Stone2Material::relativePermeabilities(values,
params.template getRealParams<EclMultiplexerApproach::Stone2>(),
fluidState);
break;
case EclMultiplexerApproach::Default:
DefaultMaterial::relativePermeabilities(values,
params.template getRealParams<EclMultiplexerApproach::Default>(),
fluidState);
break;
case EclMultiplexerApproach::TwoPhase:
TwoPhaseMaterial::relativePermeabilities(values,
params.template getRealParams<EclMultiplexerApproach::TwoPhase>(),
fluidState);
break;
case EclMultiplexerApproach::OnePhase:
values[0] = 1.0;
break;
default:
throw std::logic_error("Not implemented: relativePermeabilities() option for unknown EclMultiplexerApproach (="
+ std::to_string(static_cast<int>(params.approach())) + ")");
}
}
template <class Evaluation, class FluidState>
static Evaluation relpermOilInOilGasSystem(const Params& params,
const FluidState& fluidState)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::template relpermOilInOilGasSystem<Evaluation>
(params.template getRealParams<EclMultiplexerApproach::Stone1>(),
fluidState);
case EclMultiplexerApproach::Stone2:
return Stone2Material::template relpermOilInOilGasSystem<Evaluation>
(params.template getRealParams<EclMultiplexerApproach::Stone2>(),
fluidState);
case EclMultiplexerApproach::Default:
return DefaultMaterial::template relpermOilInOilGasSystem<Evaluation>
(params.template getRealParams<EclMultiplexerApproach::Default>(),
fluidState);
default:
throw std::logic_error {
"relpermOilInOilGasSystem() is specific to three phases"
};
}
}
template <class Evaluation, class FluidState>
static Evaluation relpermOilInOilWaterSystem(const Params& params,
const FluidState& fluidState)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::template relpermOilInOilWaterSystem<Evaluation>
(params.template getRealParams<EclMultiplexerApproach::Stone1>(),
fluidState);
case EclMultiplexerApproach::Stone2:
return Stone2Material::template relpermOilInOilWaterSystem<Evaluation>
(params.template getRealParams<EclMultiplexerApproach::Stone2>(),
fluidState);
case EclMultiplexerApproach::Default:
return DefaultMaterial::template relpermOilInOilWaterSystem<Evaluation>
(params.template getRealParams<EclMultiplexerApproach::Default>(),
fluidState);
default:
throw std::logic_error {
"relpermOilInOilWaterSystem() is specific to three phases"
};
}
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation krg(const Params& /* params */,
const FluidState& /* fluidState */)
{
throw std::logic_error("Not implemented: krg()");
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation krw(const Params& /* params */,
const FluidState& /* fluidState */)
{
throw std::logic_error("Not implemented: krw()");
}
template <class FluidState, class Evaluation = typename FluidState::Scalar>
static Evaluation krn(const Params& /* params */,
const FluidState& /* fluidState */)
{
throw std::logic_error("Not implemented: krn()");
}
template <class FluidState>
static bool updateHysteresis(Params& params, const FluidState& fluidState)
{
OPM_TIMEFUNCTION_LOCAL();
switch (params.approach()) {
case EclMultiplexerApproach::Stone1:
return Stone1Material::updateHysteresis(params.template getRealParams<EclMultiplexerApproach::Stone1>(),
fluidState);
break;
case EclMultiplexerApproach::Stone2:
return Stone2Material::updateHysteresis(params.template getRealParams<EclMultiplexerApproach::Stone2>(),
fluidState);
break;
case EclMultiplexerApproach::Default:
return DefaultMaterial::updateHysteresis(params.template getRealParams<EclMultiplexerApproach::Default>(),
fluidState);
break;
case EclMultiplexerApproach::TwoPhase:
return TwoPhaseMaterial::updateHysteresis(params.template getRealParams<EclMultiplexerApproach::TwoPhase>(),
fluidState);
break;
case EclMultiplexerApproach::OnePhase:
return false;
break;
}
return false;
}
};
} // namespace Opm
#endif
EclDefaultMaterial.hpp
Implements the default three phase capillary pressure law used by the ECLipse simulator.
EclMultiplexerMaterialParams.hpp
Multiplexer implementation for the parameters required by the multiplexed three-phase material law.
EclStone1Material.hpp
Implements the second phase capillary pressure/relperm law suggested by Stone as used by the ECLipse ...
EclStone2Material.hpp
Implements the second phase capillary pressure/relperm law suggested by Stone as used by the ECLipse ...
EclTwoPhaseMaterial.hpp
Implements a multiplexer class that provides ECL saturation functions for twophase simulations.
Opm::EclDefaultMaterial
Implements the default three phase capillary pressure law used by the ECLipse simulator.
Definition EclDefaultMaterial.hpp:62
Opm::EclDefaultMaterial::capillaryPressures
static void capillaryPressures(ContainerT &values, const Params &params, const FluidState &state)
Implements the default three phase capillary pressure law used by the ECLipse simulator.
Definition EclDefaultMaterial.hpp:137
Opm::EclDefaultMaterial::updateHysteresis
static bool updateHysteresis(Params &params, const FluidState &fluidState)
Update the hysteresis parameters after a time step.
Definition EclDefaultMaterial.hpp:461
Opm::EclDefaultMaterial::relativePermeabilities
static void relativePermeabilities(ContainerT &values, const Params &params, const FluidState &fluidState)
The relative permeability of all phases.
Definition EclDefaultMaterial.hpp:341
Opm::EclMultiplexerMaterial
Implements a multiplexer class that provides all three phase capillary pressure laws used by the ECLi...
Definition EclMultiplexerMaterial.hpp:58
Opm::EclMultiplexerMaterial::capillaryPressures
static void capillaryPressures(ContainerT &values, const Params &params, const FluidState &fluidState)
Implements the multiplexer three phase capillary pressure law used by the ECLipse simulator.
Definition EclMultiplexerMaterial.hpp:135
Opm::EclMultiplexerMaterial::relpermOilInOilWaterSystem
static Evaluation relpermOilInOilWaterSystem(const Params &params, const FluidState &fluidState)
The relative permeability of oil in oil/water system.
Definition EclMultiplexerMaterial.hpp:566
Opm::EclMultiplexerMaterial::krw
static Evaluation krw(const Params &, const FluidState &)
The relative permeability of the wetting phase.
Definition EclMultiplexerMaterial.hpp:607
Opm::EclMultiplexerMaterial::Sg
static Evaluation Sg(const Params &, const FluidState &)
The saturation of the gas phase.
Definition EclMultiplexerMaterial.hpp:449
Opm::EclMultiplexerMaterial::relpermOilInOilGasSystem
static Evaluation relpermOilInOilGasSystem(const Params &params, const FluidState &fluidState)
The relative permeability of oil in oil/gas system.
Definition EclMultiplexerMaterial.hpp:535
Opm::EclMultiplexerMaterial::pcgn
static Evaluation pcgn(const Params &, const FluidState &)
Capillary pressure between the gas and the non-wetting liquid (i.e., oil) phase.
Definition EclMultiplexerMaterial.hpp:412
Opm::EclMultiplexerMaterial::pcnw
static Evaluation pcnw(const Params &, const FluidState &)
Capillary pressure between the non-wetting liquid (i.e., oil) and the wetting liquid (i....
Definition EclMultiplexerMaterial.hpp:428
Opm::EclMultiplexerMaterial::isPressureDependent
static constexpr bool isPressureDependent
Specify whether the quantities defined by this material law are dependent on the absolute pressure.
Definition EclMultiplexerMaterial.hpp:110
Opm::EclMultiplexerMaterial::implementsTwoPhaseSatApi
static constexpr bool implementsTwoPhaseSatApi
Specify whether this material law implements the two-phase convenience API which only depends on the ...
Definition EclMultiplexerMaterial.hpp:102
Opm::EclMultiplexerMaterial::isSaturationDependent
static constexpr bool isSaturationDependent
Specify whether the quantities defined by this material law are saturation dependent.
Definition EclMultiplexerMaterial.hpp:106
Opm::EclMultiplexerMaterial::isTemperatureDependent
static constexpr bool isTemperatureDependent
Specify whether the quantities defined by this material law are temperature dependent.
Definition EclMultiplexerMaterial.hpp:114
Opm::EclMultiplexerMaterial::implementsTwoPhaseApi
static constexpr bool implementsTwoPhaseApi
Specify whether this material law implements the two-phase convenience API.
Definition EclMultiplexerMaterial.hpp:98
Opm::EclMultiplexerMaterial::relativePermeabilities
static void relativePermeabilities(ContainerT &values, const Params &params, const FluidState &fluidState)
The relative permeability of all phases.
Definition EclMultiplexerMaterial.hpp:491
Opm::EclMultiplexerMaterial::Sn
static Evaluation Sn(const Params &, const FluidState &)
The saturation of the non-wetting (i.e., oil) phase.
Definition EclMultiplexerMaterial.hpp:459
Opm::EclMultiplexerMaterial::krg
static Evaluation krg(const Params &, const FluidState &)
The relative permeability of the gas phase.
Definition EclMultiplexerMaterial.hpp:597
Opm::EclMultiplexerMaterial::isCompositionDependent
static constexpr bool isCompositionDependent
Specify whether the quantities defined by this material law are dependent on the phase composition.
Definition EclMultiplexerMaterial.hpp:118
Opm::EclMultiplexerMaterial::saturations
static void saturations(ContainerT &, const Params &, const FluidState &)
The inverse of the capillary pressure.
Definition EclMultiplexerMaterial.hpp:438
Opm::EclMultiplexerMaterial::krn
static Evaluation krn(const Params &, const FluidState &)
The relative permeability of the non-wetting (i.e., oil) phase.
Definition EclMultiplexerMaterial.hpp:617
Opm::EclMultiplexerMaterial::Sw
static Evaluation Sw(const Params &, const FluidState &)
The saturation of the wetting (i.e., water) phase.
Definition EclMultiplexerMaterial.hpp:469
Opm::EclMultiplexerMaterial::updateHysteresis
static bool updateHysteresis(Params &params, const FluidState &fluidState)
Update the hysteresis parameters after a time step.
Definition EclMultiplexerMaterial.hpp:632
Opm::EclStone1Material
Implements the second phase capillary pressure/relperm law suggested by Stone as used by the ECLipse ...
Definition EclStone1Material.hpp:60
Opm::EclStone1Material::relativePermeabilities
static void relativePermeabilities(ContainerT &values, const Params &params, const FluidState &fluidState)
The relative permeability of all phases.
Definition EclStone1Material.hpp:339
Opm::EclStone1Material::capillaryPressures
static void capillaryPressures(ContainerT &values, const Params &params, const FluidState &state)
Implements the default three phase capillary pressure law used by the ECLipse simulator.
Definition EclStone1Material.hpp:135
Opm::EclStone1Material::updateHysteresis
static bool updateHysteresis(Params &params, const FluidState &fluidState)
Update the hysteresis parameters after a time step.
Definition EclStone1Material.hpp:446
Opm::EclStone2Material
Implements the second phase capillary pressure/relperm law suggested by Stone as used by the ECLipse ...
Definition EclStone2Material.hpp:61
Opm::EclStone2Material::updateHysteresis
static bool updateHysteresis(Params &params, const FluidState &fluidState)
Update the hysteresis parameters after a time step.
Definition EclStone2Material.hpp:431
Opm::EclStone2Material::capillaryPressures
static void capillaryPressures(ContainerT &values, const Params &params, const FluidState &state)
Implements the default three phase capillary pressure law used by the ECLipse simulator.
Definition EclStone2Material.hpp:136
Opm::EclStone2Material::relativePermeabilities
static void relativePermeabilities(ContainerT &values, const Params &params, const FluidState &fluidState)
The relative permeability of all phases.
Definition EclStone2Material.hpp:341
Opm::EclTwoPhaseMaterial
Implements a multiplexer class that provides ECL saturation functions for twophase simulations.
Definition EclTwoPhaseMaterial.hpp:57
Opm::EclTwoPhaseMaterial::relativePermeabilities
static void relativePermeabilities(ContainerT &values, const Params &params, const FluidState &fluidState)
The relative permeability of all phases.
Definition EclTwoPhaseMaterial.hpp:375
Opm::EclTwoPhaseMaterial::updateHysteresis
static bool updateHysteresis(Params &params, const FluidState &fluidState)
Update the hysteresis parameters after a time step.
Definition EclTwoPhaseMaterial.hpp:452
Opm::EclTwoPhaseMaterial::capillaryPressures
static void capillaryPressures(ContainerT &values, const Params &params, const FluidState &fluidState)
Implements the multiplexer three phase capillary pressure law used by the ECLipse simulator.
Definition EclTwoPhaseMaterial.hpp:145
Opm
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30