27#ifndef OPM_WET_HUMID_GAS_PVT_HPP
28#define OPM_WET_HUMID_GAS_PVT_HPP
31#include <opm/common/OpmLog/OpmLog.hpp>
51template <
class Scalar>
54 using SamplingPoints = std::vector<std::pair<Scalar, Scalar>>;
69 void extendPvtgwTable_(
unsigned regionIdx,
74 void extendPvtgTable_(
unsigned regionIdx,
84 void setVapPars(
const Scalar par1,
const Scalar)
103 const SamplingPoints& samplePoints)
104 { saturatedWaterVaporizationFactorTable_[regionIdx].setContainerOfTuples(samplePoints); }
112 const SamplingPoints& samplePoints)
113 { saturatedOilVaporizationFactorTable_[regionIdx].setContainerOfTuples(samplePoints); }
124 {
return gasReferenceDensity_.size(); }
129 template <
class Evaluation>
134 const Evaluation&)
const
136 throw std::runtime_error(
"Requested the enthalpy of gas but the thermal "
137 "option is not enabled");
140 Scalar hVap(
unsigned)
const
142 throw std::runtime_error(
"Requested the hvap of oil but the thermal "
143 "option is not enabled");
149 template <
class Evaluation>
152 const Evaluation& pressure,
153 const Evaluation& Rv,
154 const Evaluation& Rvw)
const
156 const Evaluation& temperature = 1E30;
158 if (Rv >= (1.0 - 1e-10) * saturatedOilVaporizationFactor(regionIdx, temperature, pressure)) {
159 const Evaluation& invBg = inverseGasBRvSat_[regionIdx].eval(pressure, Rvw,
true);
160 const Evaluation& invMugBg = inverseGasBMuRvSat_[regionIdx].eval(pressure, Rvw,
true);
161 return invBg / invMugBg;
165 const Evaluation& invBg = inverseGasBRvwSat_[regionIdx].eval(pressure, Rv,
true);
166 const Evaluation& invMugBg = inverseGasBMuRvwSat_[regionIdx].eval(pressure, Rv,
true);
167 return invBg / invMugBg;
174 template <
class Evaluation>
177 const Evaluation& pressure)
const
179 const Evaluation& invBg = inverseSaturatedGasB_[regionIdx].eval(pressure,
true);
180 const Evaluation& invMugBg = inverseSaturatedGasBMu_[regionIdx].eval(pressure,
true);
182 return invBg / invMugBg;
195 template <
class Evaluation>
198 const Evaluation& pressure,
199 const Evaluation& Rv,
200 const Evaluation& Rvw)
const
202 const Evaluation& temperature = 1E30;
204 if (Rv >= (1.0 - 1e-10) * saturatedOilVaporizationFactor(regionIdx, temperature, pressure)) {
205 return inverseGasBRvSat_[regionIdx].eval(pressure, Rvw,
true);
209 return inverseGasBRvwSat_[regionIdx].eval(pressure, Rv,
true);
216 template <
class Evaluation>
219 const Evaluation& pressure)
const
220 {
return inverseSaturatedGasB_[regionIdx].eval(pressure,
true); }
225 template <
class Evaluation>
228 const Evaluation& pressure)
const
230 return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure,
true);
236 template <
class Evaluation>
239 const Evaluation& pressure,
240 const Evaluation& saltConcentration)
const
242 if (enableRwgSalt_) {
243 return saturatedWaterVaporizationSaltFactorTable_[regionIdx].eval(pressure, saltConcentration,
true);
246 return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure,
true);
250 template <
class Evaluation>
251 Evaluation saturatedOilVaporizationFactor(
unsigned regionIdx,
253 const Evaluation& pressure)
const
255 return saturatedOilVaporizationFactorTable_[regionIdx].eval(pressure,
true);
265 template <
class Evaluation>
268 const Evaluation& pressure,
269 const Evaluation& oilSaturation,
270 Evaluation maxOilSaturation)
const
273 saturatedOilVaporizationFactorTable_[regionIdx].eval(pressure,
true);
277 maxOilSaturation = min(maxOilSaturation, Scalar(1.0));
278 if (vapPar1_ > 0.0 && maxOilSaturation > 0.01 && oilSaturation < maxOilSaturation) {
279 constexpr const Scalar eps = 0.001;
280 const Evaluation& So = max(oilSaturation, eps);
281 tmp *= max(1e-3, pow(So / maxOilSaturation, vapPar1_));
295 template <
class Evaluation>
298 const Evaluation& Rw)
const
302 const auto& RwTable = saturatedWaterVaporizationFactorTable_[regionIdx];
303 constexpr const Scalar eps = std::numeric_limits<typename Toolbox::Scalar>::epsilon() * 1e6;
306 Evaluation pSat = saturationPressure_[regionIdx].eval(Rw,
true);
311 bool onProbation =
false;
312 for (
unsigned i = 0; i < 20; ++i) {
313 const Evaluation& f = RwTable.eval(pSat,
true) - Rw;
314 const Evaluation& fPrime = RwTable.evalDerivative(pSat,
true);
318 if (std::abs(scalarValue(fPrime)) < 1.0e-30) {
322 const Evaluation& delta = f / fPrime;
337 if (std::abs(scalarValue(delta)) < std::abs(scalarValue(pSat))*eps) {
342 const std::string msg =
343 "Finding saturation pressure did not converge: "
344 "pSat = " + std::to_string(getValue(pSat)) +
345 ", Rw = " + std::to_string(getValue(Rw));
346 OpmLog::debug(
"Wet gas saturation pressure", msg);
350 template <
class Evaluation>
351 Evaluation diffusionCoefficient(
const Evaluation& ,
355 throw std::runtime_error(
"Not implemented: The PVT model does not provide "
356 "a diffusionCoefficient()");
359 Scalar gasReferenceDensity(
unsigned regionIdx)
const
360 {
return gasReferenceDensity_[regionIdx]; }
362 Scalar oilReferenceDensity(
unsigned regionIdx)
const
363 {
return oilReferenceDensity_[regionIdx]; }
365 Scalar waterReferenceDensity(
unsigned regionIdx)
const
366 {
return waterReferenceDensity_[regionIdx]; }
368 const std::vector<TabulatedTwoDFunction>& inverseGasB()
const
369 {
return inverseGasBRvSat_; }
371 const std::vector<TabulatedOneDFunction>& inverseSaturatedGasB()
const
372 {
return inverseSaturatedGasB_; }
374 const std::vector<TabulatedTwoDFunction>& gasMu()
const
375 {
return gasMuRvSat_; }
377 const std::vector<TabulatedTwoDFunction>& inverseGasBMu()
const
378 {
return inverseGasBMuRvSat_; }
380 const std::vector<TabulatedOneDFunction>& inverseSaturatedGasBMu()
const
381 {
return inverseSaturatedGasBMu_; }
383 const std::vector<TabulatedOneDFunction>& saturatedWaterVaporizationFactorTable()
const
384 {
return saturatedWaterVaporizationFactorTable_; }
386 const std::vector<TabulatedTwoDFunction>& saturatedWaterVaporizationSaltFactorTable()
const
387 {
return saturatedWaterVaporizationSaltFactorTable_; }
389 const std::vector<TabulatedOneDFunction>& saturatedOilVaporizationFactorTable()
const
390 {
return saturatedOilVaporizationFactorTable_; }
393 {
return saturationPressure_; }
395 Scalar vapPar1()
const
399 void updateSaturationPressure_(
unsigned regionIdx);
401 std::vector<Scalar> gasReferenceDensity_{};
402 std::vector<Scalar> oilReferenceDensity_{};
403 std::vector<Scalar> waterReferenceDensity_{};
404 std::vector<TabulatedTwoDFunction> inverseGasBRvwSat_{};
405 std::vector<TabulatedTwoDFunction> inverseGasBRvSat_{};
406 std::vector<TabulatedOneDFunction> inverseSaturatedGasB_{};
407 std::vector<TabulatedTwoDFunction> gasMuRvwSat_{};
408 std::vector<TabulatedTwoDFunction> gasMuRvSat_{};
409 std::vector<TabulatedTwoDFunction> inverseGasBMuRvwSat_{};
410 std::vector<TabulatedTwoDFunction> inverseGasBMuRvSat_{};
411 std::vector<TabulatedOneDFunction> inverseSaturatedGasBMu_{};
412 std::vector<TabulatedOneDFunction> saturatedWaterVaporizationFactorTable_{};
413 std::vector<TabulatedTwoDFunction> saturatedWaterVaporizationSaltFactorTable_{};
414 std::vector<TabulatedOneDFunction> saturatedOilVaporizationFactorTable_{};
415 std::vector<TabulatedOneDFunction> saturationPressure_{};
417 bool enableRwgSalt_ =
false;
418 Scalar vapPar1_ = 0.0;
Provides the OPM specific exception classes.
Implements a linearly interpolated scalar function that depends on one variable.
Definition EclipseState.hpp:62
Definition Exceptions.hpp:40
Definition Schedule.hpp:101
Definition SimpleTable.hpp:35
Implements a linearly interpolated scalar function that depends on one variable.
Definition Tabulated1DFunction.hpp:51
This class represents the Pressure-Volume-Temperature relations of the gas phase with vaporized oil a...
Definition WetHumidGasPvt.hpp:53
Evaluation saturationPressure(unsigned regionIdx, const Evaluation &, const Evaluation &Rw) const
Returns the saturation pressure of the gas phase [Pa] depending on its mass fraction of the water com...
Definition WetHumidGasPvt.hpp:296
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition WetHumidGasPvt.hpp:123
void setSaturatedGasOilVaporizationFactor(unsigned regionIdx, const SamplingPoints &samplePoints)
Initialize the function for the oil vaporization factor .
Definition WetHumidGasPvt.hpp:111
Evaluation viscosity(unsigned regionIdx, const Evaluation &, const Evaluation &pressure, const Evaluation &Rv, const Evaluation &Rvw) const
Returns the dynamic viscosity [Pa s] of the fluid phase given a set of parameters.
Definition WetHumidGasPvt.hpp:150
Evaluation saturatedWaterVaporizationFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure, const Evaluation &saltConcentration) const
Returns the water vaporization factor [m^3/m^3] of the water phase.
Definition WetHumidGasPvt.hpp:237
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of water saturated gas at a given pressure.
Definition WetHumidGasPvt.hpp:217
Evaluation saturatedWaterVaporizationFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the water vaporization factor [m^3/m^3] of the water phase.
Definition WetHumidGasPvt.hpp:226
Evaluation saturatedViscosity(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the dynamic viscosity [Pa s] of oil saturated gas at a given pressure.
Definition WetHumidGasPvt.hpp:175
Evaluation saturatedOilVaporizationFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure, const Evaluation &oilSaturation, Evaluation maxOilSaturation) const
Returns the oil vaporization factor [m^3/m^3] of the gas phase.
Definition WetHumidGasPvt.hpp:266
void setReferenceDensities(unsigned regionIdx, Scalar rhoRefOil, Scalar rhoRefGas, Scalar rhoRefWater)
Initialize the reference densities of all fluids for a given PVT region.
Definition WetHumidGasPvt.cpp:422
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure, const Evaluation &Rv, const Evaluation &Rvw) const
Returns the formation volume factor [-] of the fluid phase.
Definition WetHumidGasPvt.hpp:196
Evaluation internalEnergy(unsigned, const Evaluation &, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the specific enthalpy [J/kg] of gas given a set of parameters.
Definition WetHumidGasPvt.hpp:130
void initEnd()
Finish initializing the gas phase PVT properties.
Definition WetHumidGasPvt.cpp:433
void setSaturatedGasWaterVaporizationFactor(unsigned regionIdx, const SamplingPoints &samplePoints)
Initialize the function for the water vaporization factor .
Definition WetHumidGasPvt.hpp:102
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30