27 #ifndef OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
28 #define OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
31 #include <opm/input/eclipse/EclipseState/EclipseState.hpp>
32 #include <opm/input/eclipse/Schedule/Schedule.hpp>
40 template <
class Scalar>
46 const std::vector<Scalar>& oilReferencePressure,
47 const std::vector<Scalar>& oilReferenceFormationVolumeFactor,
48 const std::vector<Scalar>& oilCompressibility,
49 const std::vector<Scalar>& oilViscosity,
50 const std::vector<Scalar>& oilViscosibility)
51 : oilReferenceDensity_(oilReferenceDensity)
52 , oilReferencePressure_(oilReferencePressure)
53 , oilReferenceFormationVolumeFactor_(oilReferenceFormationVolumeFactor)
54 , oilCompressibility_(oilCompressibility)
55 , oilViscosity_(oilViscosity)
56 , oilViscosibility_(oilViscosibility)
68 void initFromState(
const EclipseState& eclState,
const Schedule&)
70 const auto& pvcdoTable = eclState.getTableManager().getPvcdoTable();
71 const auto& densityTable = eclState.getTableManager().getDensityTable();
73 assert(pvcdoTable.size() == densityTable.size());
78 for (
unsigned regionIdx = 0; regionIdx <
numRegions; ++ regionIdx) {
79 Scalar rhoRefO = densityTable[regionIdx].oil;
80 Scalar rhoRefG = densityTable[regionIdx].gas;
81 Scalar rhoRefW = densityTable[regionIdx].water;
85 oilReferencePressure_[regionIdx] = pvcdoTable[regionIdx].reference_pressure;
86 oilReferenceFormationVolumeFactor_[regionIdx] = pvcdoTable[regionIdx].volume_factor;
87 oilCompressibility_[regionIdx] = pvcdoTable[regionIdx].compressibility;
88 oilViscosity_[regionIdx] = pvcdoTable[regionIdx].viscosity;
89 oilViscosibility_[regionIdx] = pvcdoTable[regionIdx].viscosibility;
100 oilReferenceFormationVolumeFactor_.resize(
numRegions);
105 for (
unsigned regionIdx = 0; regionIdx <
numRegions; ++regionIdx) {
118 { oilReferenceDensity_[regionIdx] = rhoRefOil; }
123 void setViscosity(
unsigned regionIdx, Scalar muo, Scalar oilViscosibility = 0.0)
125 oilViscosity_[regionIdx] = muo;
126 oilViscosibility_[regionIdx] = oilViscosibility;
133 { oilCompressibility_[regionIdx] = oilCompressibility; }
139 { oilReferencePressure_[regionIdx] = p; }
145 { oilReferenceFormationVolumeFactor_[regionIdx] = BoRef; }
151 { oilViscosibility_[regionIdx] = muComp; }
163 {
return oilViscosity_.size(); }
168 template <
class Evaluation>
172 const Evaluation&)
const
174 throw std::runtime_error(
"Requested the enthalpy of oil but the thermal option is not enabled");
181 template <
class Evaluation>
183 const Evaluation& temperature,
184 const Evaluation& pressure,
185 const Evaluation& )
const
191 template <
class Evaluation>
193 const Evaluation& temperature,
194 const Evaluation& pressure)
const
196 Scalar BoMuoRef = oilViscosity_[regionIdx]*oilReferenceFormationVolumeFactor_[regionIdx];
199 Scalar pRef = oilReferencePressure_[regionIdx];
200 const Evaluation& Y =
201 (oilCompressibility_[regionIdx] - oilViscosibility_[regionIdx])
203 return BoMuoRef*bo/(1.0 + Y*(1.0 + Y/2.0));
209 template <
class Evaluation>
211 const Evaluation& temperature,
212 const Evaluation& pressure,
213 const Evaluation& )
const
222 template <
class Evaluation>
225 const Evaluation& pressure)
const
228 Scalar pRef = oilReferencePressure_[regionIdx];
229 const Evaluation& X = oilCompressibility_[regionIdx]*(pressure - pRef);
231 Scalar BoRef = oilReferenceFormationVolumeFactor_[regionIdx];
232 return (1 + X*(1 + X/2))/BoRef;
238 template <
class Evaluation>
241 const Evaluation& )
const
247 template <
class Evaluation>
252 const Evaluation& )
const
261 template <
class Evaluation>
264 const Evaluation& )
const
267 template <
class Evaluation>
268 Evaluation diffusionCoefficient(
const Evaluation& ,
272 throw std::runtime_error(
"Not implemented: The PVT model does not provide a diffusionCoefficient()");
275 const Scalar oilReferenceDensity(
unsigned regionIdx)
const
276 {
return oilReferenceDensity_[regionIdx]; }
278 const std::vector<Scalar>& oilReferenceFormationVolumeFactor()
const
279 {
return oilReferenceFormationVolumeFactor_; }
281 const std::vector<Scalar>& oilCompressibility()
const
282 {
return oilCompressibility_; }
284 const std::vector<Scalar>& oilViscosity()
const
285 {
return oilViscosity_; }
287 const std::vector<Scalar>& oilViscosibility()
const
288 {
return oilViscosibility_; }
290 bool operator==(
const ConstantCompressibilityOilPvt<Scalar>& data)
const
292 return this->oilReferenceDensity_ == data.oilReferenceDensity_ &&
293 this->oilReferencePressure_ == data.oilReferencePressure_ &&
294 this->oilReferenceFormationVolumeFactor() == data.oilReferenceFormationVolumeFactor() &&
295 this->oilCompressibility() == data.oilCompressibility() &&
296 this->oilViscosity() == data.oilViscosity() &&
297 this->oilViscosibility() == data.oilViscosibility();
301 std::vector<Scalar> oilReferenceDensity_;
302 std::vector<Scalar> oilReferencePressure_;
303 std::vector<Scalar> oilReferenceFormationVolumeFactor_;
304 std::vector<Scalar> oilCompressibility_;
305 std::vector<Scalar> oilViscosity_;
306 std::vector<Scalar> oilViscosibility_;
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
Definition: ConstantCompressibilityOilPvt.hpp:42
Evaluation internalEnergy(unsigned, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the specific enthalpy [J/kg] of oil given a set of parameters.
Definition: ConstantCompressibilityOilPvt.hpp:169
Evaluation saturationPressure(unsigned, const Evaluation &, const Evaluation &) const
Returns the saturation pressure of the oil phase [Pa] depending on its mass fraction of the gas compo...
Definition: ConstantCompressibilityOilPvt.hpp:262
void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility=0.0)
Set the viscosity and "viscosibility" of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:123
void initEnd()
Finish initializing the oil phase PVT properties.
Definition: ConstantCompressibilityOilPvt.hpp:156
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:248
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the oil reference pressure [Pa].
Definition: ConstantCompressibilityOilPvt.hpp:138
void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
Set the compressibility of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:132
Evaluation viscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &) const
Returns the dynamic viscosity [Pa s] of gas saturated oil given a pressure and a phase composition.
Definition: ConstantCompressibilityOilPvt.hpp:182
void setReferenceDensities(unsigned regionIdx, Scalar rhoRefOil, Scalar, Scalar)
Initialize the reference densities of all fluids for a given PVT region.
Definition: ConstantCompressibilityOilPvt.hpp:114
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the oil "viscosibility" [1/ (Pa s)].
Definition: ConstantCompressibilityOilPvt.hpp:150
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &) const
Returns the formation volume factor [-] of the fluid phase.
Definition: ConstantCompressibilityOilPvt.hpp:210
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of gas saturated oil.
Definition: ConstantCompressibilityOilPvt.hpp:223
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
Set the oil reference formation volume factor [-].
Definition: ConstantCompressibilityOilPvt.hpp:144
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:239
Evaluation saturatedViscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure) const
Returns the dynamic viscosity [Pa s] of gas saturated oil given a pressure.
Definition: ConstantCompressibilityOilPvt.hpp:192
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: ConstantCompressibilityOilPvt.hpp:162