27#ifndef OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
28#define OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
44template <
class Scalar>
63 oilReferenceFormationVolumeFactor_.resize(
numRegions);
68 for (
unsigned regionIdx = 0; regionIdx <
numRegions; ++regionIdx) {
81 { oilReferenceDensity_[regionIdx] = rhoRefOil; }
86 void setViscosity(
unsigned regionIdx, Scalar muo, Scalar oilViscosibility = 0.0)
88 oilViscosity_[regionIdx] = muo;
89 oilViscosibility_[regionIdx] = oilViscosibility;
96 { oilCompressibility_[regionIdx] = oilCompressibility; }
102 { oilReferencePressure_[regionIdx] = p; }
108 { oilReferenceFormationVolumeFactor_[regionIdx] = BoRef; }
114 { oilViscosibility_[regionIdx] = muComp; }
126 {
return oilViscosity_.size(); }
131 template <
class Evaluation>
135 const Evaluation&)
const
137 throw std::runtime_error(
"Requested the enthalpy of oil but the thermal option is not enabled");
144 template <
class Evaluation>
146 const Evaluation& temperature,
147 const Evaluation& pressure,
148 const Evaluation& )
const
154 template <
class Evaluation>
156 const Evaluation& temperature,
157 const Evaluation& pressure)
const
159 Scalar BoMuoRef = oilViscosity_[regionIdx]*oilReferenceFormationVolumeFactor_[regionIdx];
162 Scalar pRef = oilReferencePressure_[regionIdx];
163 const Evaluation& Y =
164 (oilCompressibility_[regionIdx] - oilViscosibility_[regionIdx])
166 return BoMuoRef*bo/(1.0 + Y*(1.0 + Y/2.0));
172 template <
class Evaluation>
174 const Evaluation& temperature,
175 const Evaluation& pressure,
176 const Evaluation& )
const
185 template <
class Evaluation>
188 const Evaluation& pressure)
const
191 Scalar pRef = oilReferencePressure_[regionIdx];
192 const Evaluation& X = oilCompressibility_[regionIdx]*(pressure - pRef);
194 Scalar BoRef = oilReferenceFormationVolumeFactor_[regionIdx];
195 return (1 + X*(1 + X/2))/BoRef;
201 template <
class Evaluation>
204 const Evaluation& )
const
210 template <
class Evaluation>
215 const Evaluation& )
const
224 template <
class Evaluation>
227 const Evaluation& )
const
230 template <
class Evaluation>
231 Evaluation diffusionCoefficient(
const Evaluation& ,
235 throw std::runtime_error(
"Not implemented: The PVT model does not provide a diffusionCoefficient()");
238 Scalar oilReferenceDensity(
unsigned regionIdx)
const
239 {
return oilReferenceDensity_[regionIdx]; }
241 const std::vector<Scalar>& oilReferenceFormationVolumeFactor()
const
242 {
return oilReferenceFormationVolumeFactor_; }
244 const std::vector<Scalar>& oilCompressibility()
const
245 {
return oilCompressibility_; }
247 const std::vector<Scalar>& oilViscosity()
const
248 {
return oilViscosity_; }
250 const std::vector<Scalar>& oilViscosibility()
const
251 {
return oilViscosibility_; }
254 std::vector<Scalar> oilReferenceDensity_;
255 std::vector<Scalar> oilReferencePressure_;
256 std::vector<Scalar> oilReferenceFormationVolumeFactor_;
257 std::vector<Scalar> oilCompressibility_;
258 std::vector<Scalar> oilViscosity_;
259 std::vector<Scalar> oilViscosibility_;
This class represents the Pressure-Volume-Temperature relations of the oil phase without dissolved ga...
Definition: ConstantCompressibilityOilPvt.hpp:46
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:132
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:225
void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility=0.0)
Set the viscosity and "viscosibility" of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:86
void initEnd()
Finish initializing the oil phase PVT properties.
Definition: ConstantCompressibilityOilPvt.hpp:119
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:211
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the oil reference pressure [Pa].
Definition: ConstantCompressibilityOilPvt.hpp:101
void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
Set the compressibility of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:95
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:145
void setReferenceDensities(unsigned regionIdx, Scalar rhoRefOil, Scalar, Scalar)
Initialize the reference densities of all fluids for a given PVT region.
Definition: ConstantCompressibilityOilPvt.hpp:77
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the oil "viscosibility" [1/ (Pa s)].
Definition: ConstantCompressibilityOilPvt.hpp:113
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:173
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of gas saturated oil.
Definition: ConstantCompressibilityOilPvt.hpp:186
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
Set the oil reference formation volume factor [-].
Definition: ConstantCompressibilityOilPvt.hpp:107
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition: ConstantCompressibilityOilPvt.hpp:202
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:155
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: ConstantCompressibilityOilPvt.hpp:125
Definition: EclipseState.hpp:55
Definition: Schedule.hpp:130
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition: Exceptions.hpp:30