My Project
Loading...
Searching...
No Matches
ConstantCompressibilityOilPvt.hpp
Go to the documentation of this file.
1// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2// vi: set et ts=4 sw=4 sts=4:
3/*
4 This file is part of the Open Porous Media project (OPM).
5
6 OPM is free software: you can redistribute it and/or modify
7 it under the terms of the GNU General Public License as published by
8 the Free Software Foundation, either version 2 of the License, or
9 (at your option) any later version.
10
11 OPM is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 GNU General Public License for more details.
15
16 You should have received a copy of the GNU General Public License
17 along with OPM. If not, see <http://www.gnu.org/licenses/>.
18
19 Consult the COPYING file in the top-level source directory of this
20 module for the precise wording of the license and the list of
21 copyright holders.
22*/
27#ifndef OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
28#define OPM_CONSTANT_COMPRESSIBILITY_OIL_PVT_HPP
29
30#include <stdexcept>
31#include <vector>
32
33namespace Opm {
34
35#if HAVE_ECL_INPUT
36class EclipseState;
37class Schedule;
38#endif
39
44template <class Scalar>
46{
47public:
48#if HAVE_ECL_INPUT
56 void initFromState(const EclipseState& eclState, const Schedule&);
57#endif
58
59 void setNumRegions(size_t numRegions)
60 {
61 oilReferenceDensity_.resize(numRegions);
62 oilReferencePressure_.resize(numRegions);
63 oilReferenceFormationVolumeFactor_.resize(numRegions);
64 oilCompressibility_.resize(numRegions);
65 oilViscosity_.resize(numRegions);
66 oilViscosibility_.resize(numRegions);
67
68 for (unsigned regionIdx = 0; regionIdx < numRegions; ++regionIdx) {
70 setReferencePressure(regionIdx, 1.03125);
71 }
72 }
73
74 void setVapPars(const Scalar, const Scalar)
75 {
76 }
77
81 void setReferenceDensities(unsigned regionIdx,
82 Scalar rhoRefOil,
83 Scalar /*rhoRefGas*/,
84 Scalar /*rhoRefWater*/)
85 { oilReferenceDensity_[regionIdx] = rhoRefOil; }
86
90 void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility = 0.0)
91 {
92 oilViscosity_[regionIdx] = muo;
93 oilViscosibility_[regionIdx] = oilViscosibility;
94 }
95
99 void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
100 { oilCompressibility_[regionIdx] = oilCompressibility; }
101
105 void setReferencePressure(unsigned regionIdx, Scalar p)
106 { oilReferencePressure_[regionIdx] = p; }
107
111 void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
112 { oilReferenceFormationVolumeFactor_[regionIdx] = BoRef; }
113
117 void setViscosibility(unsigned regionIdx, Scalar muComp)
118 { oilViscosibility_[regionIdx] = muComp; }
119
123 void initEnd()
124 { }
125
129 unsigned numRegions() const
130 { return oilViscosity_.size(); }
131
135 template <class Evaluation>
136 Evaluation internalEnergy(unsigned,
137 const Evaluation&,
138 const Evaluation&,
139 const Evaluation&) const
140 {
141 throw std::runtime_error("Requested the enthalpy of oil but the thermal option is not enabled");
142 }
143
148 template <class Evaluation>
149 Evaluation viscosity(unsigned regionIdx,
150 const Evaluation& temperature,
151 const Evaluation& pressure,
152 const Evaluation& /*Rs*/) const
153 { return saturatedViscosity(regionIdx, temperature, pressure); }
154
158 template <class Evaluation>
159 Evaluation saturatedViscosity(unsigned regionIdx,
160 const Evaluation& temperature,
161 const Evaluation& pressure) const
162 {
163 Scalar BoMuoRef = oilViscosity_[regionIdx]*oilReferenceFormationVolumeFactor_[regionIdx];
164 const Evaluation& bo = saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure);
165
166 Scalar pRef = oilReferencePressure_[regionIdx];
167 const Evaluation& Y =
168 (oilCompressibility_[regionIdx] - oilViscosibility_[regionIdx])
169 * (pressure - pRef);
170 return BoMuoRef*bo/(1.0 + Y*(1.0 + Y/2.0));
171 }
172
176 template <class Evaluation>
177 Evaluation inverseFormationVolumeFactor(unsigned regionIdx,
178 const Evaluation& temperature,
179 const Evaluation& pressure,
180 const Evaluation& /*Rs*/) const
181 { return saturatedInverseFormationVolumeFactor(regionIdx, temperature, pressure); }
182
189 template <class Evaluation>
190 Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx,
191 const Evaluation& /*temperature*/,
192 const Evaluation& pressure) const
193 {
194 // cf. ECLiPSE 2011 technical description, p. 116
195 Scalar pRef = oilReferencePressure_[regionIdx];
196 const Evaluation& X = oilCompressibility_[regionIdx]*(pressure - pRef);
197
198 Scalar BoRef = oilReferenceFormationVolumeFactor_[regionIdx];
199 return (1 + X*(1 + X/2))/BoRef;
200 }
201
205 template <class Evaluation>
206 Evaluation saturatedGasDissolutionFactor(unsigned /*regionIdx*/,
207 const Evaluation& /*temperature*/,
208 const Evaluation& /*pressure*/) const
209 { return 0.0; /* this is dead oil! */ }
210
214 template <class Evaluation>
215 Evaluation saturatedGasDissolutionFactor(unsigned /*regionIdx*/,
216 const Evaluation& /*temperature*/,
217 const Evaluation& /*pressure*/,
218 const Evaluation& /*oilSaturation*/,
219 const Evaluation& /*maxOilSaturation*/) const
220 { return 0.0; /* this is dead oil! */ }
221
228 template <class Evaluation>
229 Evaluation saturationPressure(unsigned /*regionIdx*/,
230 const Evaluation& /*temperature*/,
231 const Evaluation& /*Rs*/) const
232 { return 0.0; /* this is dead oil, so there isn't any meaningful saturation pressure! */ }
233
234 template <class Evaluation>
235 Evaluation diffusionCoefficient(const Evaluation& /*temperature*/,
236 const Evaluation& /*pressure*/,
237 unsigned /*compIdx*/) const
238 {
239 throw std::runtime_error("Not implemented: The PVT model does not provide a diffusionCoefficient()");
240 }
241
242 Scalar oilReferenceDensity(unsigned regionIdx) const
243 { return oilReferenceDensity_[regionIdx]; }
244
245 const std::vector<Scalar>& oilReferenceFormationVolumeFactor() const
246 { return oilReferenceFormationVolumeFactor_; }
247
248 const std::vector<Scalar>& oilCompressibility() const
249 { return oilCompressibility_; }
250
251 const std::vector<Scalar>& oilViscosity() const
252 { return oilViscosity_; }
253
254 const std::vector<Scalar>& oilViscosibility() const
255 { return oilViscosibility_; }
256
257private:
258 std::vector<Scalar> oilReferenceDensity_;
259 std::vector<Scalar> oilReferencePressure_;
260 std::vector<Scalar> oilReferenceFormationVolumeFactor_;
261 std::vector<Scalar> oilCompressibility_;
262 std::vector<Scalar> oilViscosity_;
263 std::vector<Scalar> oilViscosibility_;
264};
265
266} // namespace Opm
267
268#endif
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:136
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:229
void setViscosity(unsigned regionIdx, Scalar muo, Scalar oilViscosibility=0.0)
Set the viscosity and "viscosibility" of the oil phase.
Definition ConstantCompressibilityOilPvt.hpp:90
void initEnd()
Finish initializing the oil phase PVT properties.
Definition ConstantCompressibilityOilPvt.hpp:123
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:215
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the oil reference pressure [Pa].
Definition ConstantCompressibilityOilPvt.hpp:105
void setCompressibility(unsigned regionIdx, Scalar oilCompressibility)
Set the compressibility of the oil phase.
Definition ConstantCompressibilityOilPvt.hpp:99
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:149
void setReferenceDensities(unsigned regionIdx, Scalar rhoRefOil, Scalar, Scalar)
Initialize the reference densities of all fluids for a given PVT region.
Definition ConstantCompressibilityOilPvt.hpp:81
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the oil "viscosibility" [1/ (Pa s)].
Definition ConstantCompressibilityOilPvt.hpp:117
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:177
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of gas saturated oil.
Definition ConstantCompressibilityOilPvt.hpp:190
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BoRef)
Set the oil reference formation volume factor [-].
Definition ConstantCompressibilityOilPvt.hpp:111
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the oil phase.
Definition ConstantCompressibilityOilPvt.hpp:206
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:159
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition ConstantCompressibilityOilPvt.hpp:129
Definition EclipseState.hpp:63
Definition Schedule.hpp:88
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30