My Project
SolventPvt.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_SOLVENT_PVT_HPP
28#define OPM_SOLVENT_PVT_HPP
29
31
32#include <vector>
33
34namespace Opm {
35
36#if HAVE_ECL_INPUT
37class EclipseState;
38class Schedule;
39#endif
40
45template <class Scalar>
47{
48 using SamplingPoints = std::vector<std::pair<Scalar, Scalar>>;
49
50public:
52
53#if HAVE_ECL_INPUT
59 void initFromState(const EclipseState& eclState, const Schedule&);
60#endif
61
62 void setNumRegions(size_t numRegions);
63
67 void setReferenceDensity(unsigned regionIdx, Scalar rhoRefSolvent)
68 { solventReferenceDensity_[regionIdx] = rhoRefSolvent; }
69
75 void setSolventViscosity(unsigned regionIdx, const TabulatedOneDFunction& mug)
76 { solventMu_[regionIdx] = mug; }
77
83 void setSolventFormationVolumeFactor(unsigned regionIdx,
84 const SamplingPoints& samplePoints);
85
89 void initEnd();
90
94 unsigned numRegions() const
95 { return solventReferenceDensity_.size(); }
96
100 template <class Evaluation>
101 Evaluation viscosity(unsigned regionIdx,
102 const Evaluation&,
103 const Evaluation& pressure) const
104 {
105 const Evaluation& invBg = inverseSolventB_[regionIdx].eval(pressure, /*extrapolate=*/true);
106 const Evaluation& invMugBg = inverseSolventBMu_[regionIdx].eval(pressure, /*extrapolate=*/true);
107
108 return invBg/invMugBg;
109 }
110
111 template <class Evaluation>
112 Evaluation diffusionCoefficient(const Evaluation& /*temperature*/,
113 const Evaluation& /*pressure*/,
114 unsigned /*compIdx*/) const
115 {
116 throw std::runtime_error("Not implemented: The PVT model does not provide a diffusionCoefficient()");
117 }
118
122 Scalar referenceDensity(unsigned regionIdx) const
123 { return solventReferenceDensity_[regionIdx]; }
124
128 template <class Evaluation>
129 Evaluation inverseFormationVolumeFactor(unsigned regionIdx,
130 const Evaluation&,
131 const Evaluation& pressure) const
132 { return inverseSolventB_[regionIdx].eval(pressure, /*extrapolate=*/true); }
133
134 const std::vector<Scalar>& solventReferenceDensity() const
135 { return solventReferenceDensity_; }
136
137 const std::vector<TabulatedOneDFunction>& inverseSolventB() const
138 { return inverseSolventB_; }
139
140 const std::vector<TabulatedOneDFunction>& solventMu() const
141 { return solventMu_; }
142
143 const std::vector<TabulatedOneDFunction>& inverseSolventBMu() const
144 { return inverseSolventBMu_; }
145
146private:
147 std::vector<Scalar> solventReferenceDensity_;
148 std::vector<TabulatedOneDFunction> inverseSolventB_;
149 std::vector<TabulatedOneDFunction> solventMu_;
150 std::vector<TabulatedOneDFunction> inverseSolventBMu_;
151};
152
153} // namespace Opm
154
155#endif
Definition: EclipseState.hpp:55
Definition: Schedule.hpp:130
This class represents the Pressure-Volume-Temperature relations of the "second" gas phase in the of E...
Definition: SolventPvt.hpp:47
Evaluation viscosity(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the dynamic viscosity [Pa s] of the fluid phase given a set of parameters.
Definition: SolventPvt.hpp:101
void setReferenceDensity(unsigned regionIdx, Scalar rhoRefSolvent)
Initialize the reference density of the solvent gas for a given PVT region.
Definition: SolventPvt.hpp:67
void setSolventViscosity(unsigned regionIdx, const TabulatedOneDFunction &mug)
Initialize the viscosity of the solvent gas phase.
Definition: SolventPvt.hpp:75
Scalar referenceDensity(unsigned regionIdx) const
Return the reference density the solvent phase for a given PVT region.
Definition: SolventPvt.hpp:122
void initEnd()
Finish initializing the oil phase PVT properties.
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: SolventPvt.hpp:94
void setSolventFormationVolumeFactor(unsigned regionIdx, const SamplingPoints &samplePoints)
Initialize the function for the formation volume factor of solvent gas.
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of the fluid phase.
Definition: SolventPvt.hpp:129
Implements a linearly interpolated scalar function that depends on one variable.
Definition: Tabulated1DFunction.hpp:51
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition: Exceptions.hpp:30