27#ifndef OPM_WET_HUMID_GAS_PVT_HPP
28#define OPM_WET_HUMID_GAS_PVT_HPP
31#include <opm/common/OpmLog/OpmLog.hpp>
49template <
class Scalar>
52 using SamplingPoints = std::vector<std::pair<Scalar, Scalar>>;
67 void extendPvtgwTable_(
unsigned regionIdx,
72 void extendPvtgTable_(
unsigned regionIdx,
82 void setVapPars(
const Scalar par1,
const Scalar)
101 { saturatedWaterVaporizationFactorTable_[regionIdx].setContainerOfTuples(samplePoints); }
109 { saturatedOilVaporizationFactorTable_[regionIdx].setContainerOfTuples(samplePoints); }
121 {
return gasReferenceDensity_.size(); }
126 template <
class Evaluation>
131 const Evaluation&)
const
133 throw std::runtime_error(
"Requested the enthalpy of gas but the thermal option is not enabled");
139 template <
class Evaluation>
142 const Evaluation& pressure,
143 const Evaluation& Rv,
144 const Evaluation& Rvw)
const
146 const Evaluation& temperature = 1E30;
148 if (Rv >= (1.0 - 1e-10)*saturatedOilVaporizationFactor(regionIdx, temperature, pressure)) {
149 const Evaluation& invBg = inverseGasBRvSat_[regionIdx].eval(pressure, Rvw,
true);
150 const Evaluation& invMugBg = inverseGasBMuRvSat_[regionIdx].eval(pressure, Rvw,
true);
151 return invBg/invMugBg;
155 const Evaluation& invBg = inverseGasBRvwSat_[regionIdx].eval(pressure, Rv,
true);
156 const Evaluation& invMugBg = inverseGasBMuRvwSat_[regionIdx].eval(pressure, Rv,
true);
157 return invBg/invMugBg;
164 template <
class Evaluation>
167 const Evaluation& pressure)
const
169 const Evaluation& invBg = inverseSaturatedGasB_[regionIdx].eval(pressure,
true);
170 const Evaluation& invMugBg = inverseSaturatedGasBMu_[regionIdx].eval(pressure,
true);
172 return invBg/invMugBg;
185 template <
class Evaluation>
188 const Evaluation& pressure,
189 const Evaluation& Rv,
190 const Evaluation& Rvw)
const
192 const Evaluation& temperature = 1E30;
194 if (Rv >= (1.0 - 1e-10)*saturatedOilVaporizationFactor(regionIdx, temperature, pressure)) {
195 return inverseGasBRvSat_[regionIdx].eval(pressure, Rvw,
true);
199 return inverseGasBRvwSat_[regionIdx].eval(pressure, Rv,
true);
208 template <
class Evaluation>
211 const Evaluation& pressure)
const
212 {
return inverseSaturatedGasB_[regionIdx].eval(pressure,
true); }
217 template <
class Evaluation>
220 const Evaluation& pressure)
const
222 return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure,
true);
228 template <
class Evaluation>
231 const Evaluation& pressure,
232 const Evaluation& saltConcentration)
const
235 return saturatedWaterVaporizationSaltFactorTable_[regionIdx].eval(pressure, saltConcentration,
true);
237 return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure,
true);
242 template <
class Evaluation>
243 Evaluation saturatedOilVaporizationFactor(
unsigned regionIdx,
245 const Evaluation& pressure)
const
247 return saturatedOilVaporizationFactorTable_[regionIdx].eval(pressure,
true);
257 template <
class Evaluation>
260 const Evaluation& pressure,
261 const Evaluation& oilSaturation,
262 Evaluation maxOilSaturation)
const
265 saturatedOilVaporizationFactorTable_[regionIdx].eval(pressure,
true);
269 maxOilSaturation = min(maxOilSaturation, Scalar(1.0));
270 if (vapPar1_ > 0.0 && maxOilSaturation > 0.01 && oilSaturation < maxOilSaturation) {
271 constexpr const Scalar eps = 0.001;
272 const Evaluation& So = max(oilSaturation, eps);
273 tmp *= max(1e-3, pow(So/maxOilSaturation, vapPar1_));
287 template <
class Evaluation>
290 const Evaluation& Rw)
const
294 const auto& RwTable = saturatedWaterVaporizationFactorTable_[regionIdx];
295 constexpr const Scalar eps = std::numeric_limits<typename Toolbox::Scalar>::epsilon()*1e6;
298 Evaluation pSat = saturationPressure_[regionIdx].eval(Rw,
true);
303 bool onProbation =
false;
304 for (
unsigned i = 0; i < 20; ++i) {
305 const Evaluation& f = RwTable.eval(pSat,
true) - Rw;
306 const Evaluation& fPrime = RwTable.evalDerivative(pSat,
true);
310 if (std::abs(scalarValue(fPrime)) < 1.0e-30) {
314 const Evaluation& delta = f/fPrime;
328 if (std::abs(scalarValue(delta)) < std::abs(scalarValue(pSat))*eps)
332 const std::string msg =
333 "Finding saturation pressure did not converge: "
334 "pSat = " + std::to_string(getValue(pSat)) +
335 ", Rw = " + std::to_string(getValue(Rw));
336 OpmLog::debug(
"Wet gas saturation pressure", msg);
340 template <
class Evaluation>
341 Evaluation diffusionCoefficient(
const Evaluation& ,
345 throw std::runtime_error(
"Not implemented: The PVT model does not provide a diffusionCoefficient()");
348 Scalar gasReferenceDensity(
unsigned regionIdx)
const
349 {
return gasReferenceDensity_[regionIdx]; }
351 Scalar oilReferenceDensity(
unsigned regionIdx)
const
352 {
return oilReferenceDensity_[regionIdx]; }
354 Scalar waterReferenceDensity(
unsigned regionIdx)
const
355 {
return waterReferenceDensity_[regionIdx]; }
357 const std::vector<TabulatedTwoDFunction>& inverseGasB()
const {
358 return inverseGasBRvSat_;
361 const std::vector<TabulatedOneDFunction>& inverseSaturatedGasB()
const {
362 return inverseSaturatedGasB_;
365 const std::vector<TabulatedTwoDFunction>& gasMu()
const {
369 const std::vector<TabulatedTwoDFunction>& inverseGasBMu()
const {
370 return inverseGasBMuRvSat_;
373 const std::vector<TabulatedOneDFunction>& inverseSaturatedGasBMu()
const {
374 return inverseSaturatedGasBMu_;
377 const std::vector<TabulatedOneDFunction>& saturatedWaterVaporizationFactorTable()
const {
378 return saturatedWaterVaporizationFactorTable_;
381 const std::vector<TabulatedTwoDFunction>& saturatedWaterVaporizationSaltFactorTable()
const {
382 return saturatedWaterVaporizationSaltFactorTable_;
385 const std::vector<TabulatedOneDFunction>& saturatedOilVaporizationFactorTable()
const {
386 return saturatedOilVaporizationFactorTable_;
390 return saturationPressure_;
393 Scalar vapPar1()
const {
398 void updateSaturationPressure_(
unsigned regionIdx);
400 std::vector<Scalar> gasReferenceDensity_;
401 std::vector<Scalar> oilReferenceDensity_;
402 std::vector<Scalar> waterReferenceDensity_;
403 std::vector<TabulatedTwoDFunction> inverseGasBRvwSat_;
404 std::vector<TabulatedTwoDFunction> inverseGasBRvSat_;
405 std::vector<TabulatedOneDFunction> inverseSaturatedGasB_;
406 std::vector<TabulatedTwoDFunction> gasMuRvwSat_;
407 std::vector<TabulatedTwoDFunction> gasMuRvSat_;
408 std::vector<TabulatedTwoDFunction> inverseGasBMuRvwSat_;
409 std::vector<TabulatedTwoDFunction> inverseGasBMuRvSat_;
410 std::vector<TabulatedOneDFunction> inverseSaturatedGasBMu_;
411 std::vector<TabulatedOneDFunction> saturatedWaterVaporizationFactorTable_;
412 std::vector<TabulatedTwoDFunction> saturatedWaterVaporizationSaltFactorTable_;
413 std::vector<TabulatedOneDFunction> saturatedOilVaporizationFactorTable_;
414 std::vector<TabulatedOneDFunction> saturationPressure_;
416 bool enableRwgSalt_ =
false;
417 Scalar vapPar1_ = 0.0;
Provides the OPM specific exception classes.
Implements a linearly interpolated scalar function that depends on one variable.
Definition EclipseState.hpp:63
Definition Exceptions.hpp:40
Definition Schedule.hpp:88
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:51
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:288
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition WetHumidGasPvt.hpp:120
void setSaturatedGasOilVaporizationFactor(unsigned regionIdx, const SamplingPoints &samplePoints)
Initialize the function for the oil vaporization factor .
Definition WetHumidGasPvt.hpp:108
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:140
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:229
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:209
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:218
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:165
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:258
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:186
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:127
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:100
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30