27#ifndef OPM_OIL_PVT_MULTIPLEXER_HPP
28#define OPM_OIL_PVT_MULTIPLEXER_HPP
44#define OPM_OIL_PVT_MULTIPLEXER_CALL(codeToCall) \
45 switch (approach_) { \
46 case OilPvtApproach::ConstantCompressibilityOil: { \
47 auto& pvtImpl = getRealPvt<OilPvtApproach::ConstantCompressibilityOil>(); \
51 case OilPvtApproach::DeadOil: { \
52 auto& pvtImpl = getRealPvt<OilPvtApproach::DeadOil>(); \
56 case OilPvtApproach::LiveOil: { \
57 auto& pvtImpl = getRealPvt<OilPvtApproach::LiveOil>(); \
61 case OilPvtApproach::ThermalOil: { \
62 auto& pvtImpl = getRealPvt<OilPvtApproach::ThermalOil>(); \
66 case OilPvtApproach::BrineCo2: { \
67 auto& pvtImpl = getRealPvt<OilPvtApproach::BrineCo2>(); \
71 case OilPvtApproach::BrineH2: { \
72 auto& pvtImpl = getRealPvt<OilPvtApproach::BrineH2>(); \
76 case OilPvtApproach::NoOil: \
77 throw std::logic_error("Not implemented: Oil PVT of this deck!"); \
80enum class OilPvtApproach {
84 ConstantCompressibilityOil,
102template <
class Scalar,
bool enableThermal = true>
108 approach_ = OilPvtApproach::NoOil;
109 realOilPvt_ =
nullptr;
114 , realOilPvt_(realOilPvt)
125 case OilPvtApproach::LiveOil: {
126 delete &getRealPvt<OilPvtApproach::LiveOil>();
129 case OilPvtApproach::DeadOil: {
130 delete &getRealPvt<OilPvtApproach::DeadOil>();
133 case OilPvtApproach::ConstantCompressibilityOil: {
134 delete &getRealPvt<OilPvtApproach::ConstantCompressibilityOil>();
137 case OilPvtApproach::ThermalOil: {
138 delete &getRealPvt<OilPvtApproach::ThermalOil>();
141 case OilPvtApproach::BrineCo2: {
142 delete &getRealPvt<OilPvtApproach::BrineCo2>();
145 case OilPvtApproach::BrineH2: {
146 delete &getRealPvt<OilPvtApproach::BrineH2>();
149 case OilPvtApproach::NoOil:
165 { OPM_OIL_PVT_MULTIPLEXER_CALL(pvtImpl.initEnd()); }
171 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.numRegions());
return 1; }
173 void setVapPars(
const Scalar par1,
const Scalar par2)
175 OPM_OIL_PVT_MULTIPLEXER_CALL(pvtImpl.setVapPars(par1, par2));
182 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.oilReferenceDensity(regionIdx));
return 700.; }
187 template <
class Evaluation>
189 const Evaluation& temperature,
190 const Evaluation& pressure,
191 const Evaluation& Rs)
const
192 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.internalEnergy(regionIdx, temperature, pressure, Rs));
return 0; }
197 template <
class Evaluation>
199 const Evaluation& temperature,
200 const Evaluation& pressure,
201 const Evaluation& Rs)
const
202 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.viscosity(regionIdx, temperature, pressure, Rs));
return 0; }
207 template <
class Evaluation>
209 const Evaluation& temperature,
210 const Evaluation& pressure)
const
211 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.saturatedViscosity(regionIdx, temperature, pressure));
return 0; }
216 template <
class Evaluation>
218 const Evaluation& temperature,
219 const Evaluation& pressure,
220 const Evaluation& Rs)
const
221 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.inverseFormationVolumeFactor(regionIdx, temperature, pressure, Rs));
return 0; }
226 template <
class Evaluation>
228 const Evaluation& temperature,
229 const Evaluation& pressure)
const
230 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure));
return 0; }
235 template <
class Evaluation>
237 const Evaluation& temperature,
238 const Evaluation& pressure)
const
239 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.saturatedGasDissolutionFactor(regionIdx, temperature, pressure));
return 0; }
244 template <
class Evaluation>
246 const Evaluation& temperature,
247 const Evaluation& pressure,
248 const Evaluation& oilSaturation,
249 const Evaluation& maxOilSaturation)
const
250 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.saturatedGasDissolutionFactor(regionIdx, temperature, pressure, oilSaturation, maxOilSaturation));
return 0; }
259 template <
class Evaluation>
261 const Evaluation& temperature,
262 const Evaluation& Rs)
const
263 { OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.saturationPressure(regionIdx, temperature, Rs));
return 0; }
268 template <
class Evaluation>
270 const Evaluation& pressure,
271 unsigned compIdx)
const
273 OPM_OIL_PVT_MULTIPLEXER_CALL(
return pvtImpl.diffusionCoefficient(temperature, pressure, compIdx));
return 0;
276 void setApproach(OilPvtApproach appr)
279 case OilPvtApproach::LiveOil:
283 case OilPvtApproach::DeadOil:
287 case OilPvtApproach::ConstantCompressibilityOil:
291 case OilPvtApproach::ThermalOil:
295 case OilPvtApproach::BrineCo2:
299 case OilPvtApproach::BrineH2:
303 case OilPvtApproach::NoOil:
304 throw std::logic_error(
"Not implemented: Oil PVT of this deck!");
316 {
return approach_; }
319 template <OilPvtApproach approachV>
320 typename std::enable_if<approachV == OilPvtApproach::LiveOil, LiveOilPvt<Scalar> >::type& getRealPvt()
326 template <OilPvtApproach approachV>
327 typename std::enable_if<approachV == OilPvtApproach::LiveOil, const LiveOilPvt<Scalar> >::type& getRealPvt()
const
330 return *
static_cast<LiveOilPvt<Scalar>*
>(realOilPvt_);
333 template <OilPvtApproach approachV>
334 typename std::enable_if<approachV == OilPvtApproach::DeadOil, DeadOilPvt<Scalar> >::type& getRealPvt()
337 return *
static_cast<DeadOilPvt<Scalar>*
>(realOilPvt_);
340 template <OilPvtApproach approachV>
341 typename std::enable_if<approachV == OilPvtApproach::DeadOil, const DeadOilPvt<Scalar> >::type& getRealPvt()
const
344 return *
static_cast<DeadOilPvt<Scalar>*
>(realOilPvt_);
347 template <OilPvtApproach approachV>
348 typename std::enable_if<approachV == OilPvtApproach::ConstantCompressibilityOil, ConstantCompressibilityOilPvt<Scalar> >::type& getRealPvt()
351 return *
static_cast<ConstantCompressibilityOilPvt<Scalar>*
>(realOilPvt_);
354 template <OilPvtApproach approachV>
355 typename std::enable_if<approachV == OilPvtApproach::ConstantCompressibilityOil, const ConstantCompressibilityOilPvt<Scalar> >::type& getRealPvt()
const
358 return *
static_cast<ConstantCompressibilityOilPvt<Scalar>*
>(realOilPvt_);
361 template <OilPvtApproach approachV>
362 typename std::enable_if<approachV == OilPvtApproach::ThermalOil, OilPvtThermal<Scalar> >::type& getRealPvt()
365 return *
static_cast<OilPvtThermal<Scalar>*
>(realOilPvt_);
368 template <OilPvtApproach approachV>
369 typename std::enable_if<approachV == OilPvtApproach::ThermalOil, const OilPvtThermal<Scalar> >::type& getRealPvt()
const
372 return *
static_cast<const OilPvtThermal<Scalar>*
>(realOilPvt_);
375 template <OilPvtApproach approachV>
376 typename std::enable_if<approachV == OilPvtApproach::BrineCo2, BrineCo2Pvt<Scalar> >::type& getRealPvt()
379 return *
static_cast<BrineCo2Pvt<Scalar>*
>(realOilPvt_);
382 template <OilPvtApproach approachV>
383 typename std::enable_if<approachV == OilPvtApproach::BrineCo2, const BrineCo2Pvt<Scalar> >::type& getRealPvt()
const
386 return *
static_cast<const BrineCo2Pvt<Scalar>*
>(realOilPvt_);
389 const void* realOilPvt()
const {
return realOilPvt_; }
391 template <OilPvtApproach approachV>
392 typename std::enable_if<approachV == OilPvtApproach::BrineH2, BrineH2Pvt<Scalar> >::type& getRealPvt()
395 return *
static_cast<BrineH2Pvt<Scalar>*
>(realOilPvt_);
398 template <OilPvtApproach approachV>
399 typename std::enable_if<approachV == OilPvtApproach::BrineH2, const BrineH2Pvt<Scalar> >::type& getRealPvt()
const
402 return *
static_cast<const BrineH2Pvt<Scalar>*
>(realOilPvt_);
405 OilPvtMultiplexer<Scalar,enableThermal>& operator=(
const OilPvtMultiplexer<Scalar,enableThermal>& data)
407 approach_ = data.approach_;
409 case OilPvtApproach::ConstantCompressibilityOil:
410 realOilPvt_ =
new ConstantCompressibilityOilPvt<Scalar>(*
static_cast<const ConstantCompressibilityOilPvt<Scalar>*
>(data.realOilPvt_));
412 case OilPvtApproach::DeadOil:
413 realOilPvt_ =
new DeadOilPvt<Scalar>(*
static_cast<const DeadOilPvt<Scalar>*
>(data.realOilPvt_));
415 case OilPvtApproach::LiveOil:
416 realOilPvt_ =
new LiveOilPvt<Scalar>(*
static_cast<const LiveOilPvt<Scalar>*
>(data.realOilPvt_));
418 case OilPvtApproach::ThermalOil:
419 realOilPvt_ =
new OilPvtThermal<Scalar>(*
static_cast<const OilPvtThermal<Scalar>*
>(data.realOilPvt_));
421 case OilPvtApproach::BrineCo2:
422 realOilPvt_ =
new BrineCo2Pvt<Scalar>(*
static_cast<const BrineCo2Pvt<Scalar>*
>(data.realOilPvt_));
424 case OilPvtApproach::BrineH2:
425 realOilPvt_ =
new BrineH2Pvt<Scalar>(*
static_cast<const BrineH2Pvt<Scalar>*
>(data.realOilPvt_));
435 OilPvtApproach approach_;
This class represents the Pressure-Volume-Temperature relations of the liquid phase for a CO2-Brine s...
This class represents the Pressure-Volume-Temperature relations of the liquid phase for a H2-Brine sy...
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
This class represents the Pressure-Volume-Temperature relations of the oil phas with dissolved gas.
This class implements temperature dependence of the PVT properties of oil.
This class represents the Pressure-Volume-Temperature relations of the liquid phase for a CO2-Brine s...
Definition BrineCo2Pvt.hpp:61
This class represents the Pressure-Volume-Temperature relations of the liquid phase for a H2-Brine sy...
Definition BrineH2Pvt.hpp:53
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
Definition ConstantCompressibilityOilPvt.hpp:46
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
Definition DeadOilPvt.hpp:45
Definition EclipseState.hpp:63
This class represents the Pressure-Volume-Temperature relations of the oil phas with dissolved gas.
Definition LiveOilPvt.hpp:51
This class represents the Pressure-Volume-Temperature relations of the oil phase in the black-oil mod...
Definition OilPvtMultiplexer.hpp:104
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition OilPvtMultiplexer.hpp:170
Evaluation internalEnergy(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &Rs) const
Returns the specific enthalpy [J/kg] oil given a set of parameters.
Definition OilPvtMultiplexer.hpp:188
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &Rs) const
Returns the formation volume factor [-] of the fluid phase.
Definition OilPvtMultiplexer.hpp:217
Evaluation saturatedViscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure) const
Returns the dynamic viscosity [Pa s] of the fluid phase given a set of parameters.
Definition OilPvtMultiplexer.hpp:208
OilPvtApproach approach() const
Returns the concrete approach for calculating the PVT relations.
Definition OilPvtMultiplexer.hpp:315
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure) const
Returns the formation volume factor [-] of the fluid phase.
Definition OilPvtMultiplexer.hpp:227
Evaluation saturatedGasDissolutionFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure) const
Returns the gas dissolution factor [m^3/m^3] of saturated oil.
Definition OilPvtMultiplexer.hpp:236
Evaluation viscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &Rs) const
Returns the dynamic viscosity [Pa s] of the fluid phase given a set of parameters.
Definition OilPvtMultiplexer.hpp:198
Evaluation saturationPressure(unsigned regionIdx, const Evaluation &temperature, const Evaluation &Rs) const
Returns the saturation pressure [Pa] of oil given the mass fraction of the gas component in the oil p...
Definition OilPvtMultiplexer.hpp:260
Evaluation saturatedGasDissolutionFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &oilSaturation, const Evaluation &maxOilSaturation) const
Returns the gas dissolution factor [m^3/m^3] of saturated oil.
Definition OilPvtMultiplexer.hpp:245
const Scalar oilReferenceDensity(unsigned regionIdx) const
Return the reference density which are considered by this PVT-object.
Definition OilPvtMultiplexer.hpp:181
Evaluation diffusionCoefficient(const Evaluation &temperature, const Evaluation &pressure, unsigned compIdx) const
Calculate the binary molecular diffusion coefficient for a component in a fluid phase [mol^2 * s / (k...
Definition OilPvtMultiplexer.hpp:269
This class implements temperature dependence of the PVT properties of oil.
Definition OilPvtThermal.hpp:50
Definition Schedule.hpp:88
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30