My Project
Loading...
Searching...
No Matches
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 <cstddef>
33#include <vector>
34
35namespace Opm {
36
37#if HAVE_ECL_INPUT
38class EclipseState;
39class Schedule;
40#endif
41
46template <class Scalar>
48{
49 using SamplingPoints = std::vector<std::pair<Scalar, Scalar>>;
50
51public:
53
54#if HAVE_ECL_INPUT
60 void initFromState(const EclipseState& eclState, const Schedule&);
61#endif
62
63 void setNumRegions(std::size_t numRegions);
64
65 void setVapPars(const Scalar, const Scalar)
66 {
67 }
68
72 void setReferenceDensity(unsigned regionIdx, Scalar rhoRefSolvent)
73 { solventReferenceDensity_[regionIdx] = rhoRefSolvent; }
74
80 void setSolventViscosity(unsigned regionIdx, const TabulatedOneDFunction& mug)
81 { solventMu_[regionIdx] = mug; }
82
88 void setSolventFormationVolumeFactor(unsigned regionIdx,
89 const SamplingPoints& samplePoints);
90
94 void initEnd();
95
99 unsigned numRegions() const
100 { return solventReferenceDensity_.size(); }
101
105 template <class Evaluation>
106 Evaluation viscosity(unsigned regionIdx,
107 const Evaluation&,
108 const Evaluation& pressure) const
109 {
110 const Evaluation& invBg = inverseSolventB_[regionIdx].eval(pressure, /*extrapolate=*/true);
111 const Evaluation& invMugBg = inverseSolventBMu_[regionIdx].eval(pressure, /*extrapolate=*/true);
112
113 return invBg / invMugBg;
114 }
115
116 template <class Evaluation>
117 Evaluation diffusionCoefficient(const Evaluation& /*temperature*/,
118 const Evaluation& /*pressure*/,
119 unsigned /*compIdx*/) const
120 {
121 throw std::runtime_error("Not implemented: The PVT model does not provide "
122 "a diffusionCoefficient()");
123 }
124
128 Scalar referenceDensity(unsigned regionIdx) const
129 { return solventReferenceDensity_[regionIdx]; }
130
134 template <class Evaluation>
135 Evaluation inverseFormationVolumeFactor(unsigned regionIdx,
136 const Evaluation&,
137 const Evaluation& pressure) const
138 { return inverseSolventB_[regionIdx].eval(pressure, /*extrapolate=*/true); }
139
140 const std::vector<Scalar>& solventReferenceDensity() const
141 { return solventReferenceDensity_; }
142
143 const std::vector<TabulatedOneDFunction>& inverseSolventB() const
144 { return inverseSolventB_; }
145
146 const std::vector<TabulatedOneDFunction>& solventMu() const
147 { return solventMu_; }
148
149 const std::vector<TabulatedOneDFunction>& inverseSolventBMu() const
150 { return inverseSolventBMu_; }
151
152private:
153 std::vector<Scalar> solventReferenceDensity_{};
154 std::vector<TabulatedOneDFunction> inverseSolventB_{};
155 std::vector<TabulatedOneDFunction> solventMu_{};
156 std::vector<TabulatedOneDFunction> inverseSolventBMu_{};
157};
158
159} // namespace Opm
160
161#endif
Implements a linearly interpolated scalar function that depends on one variable.
Definition EclipseState.hpp:63
Definition Schedule.hpp:101
This class represents the Pressure-Volume-Temperature relations of the "second" gas phase in the of E...
Definition SolventPvt.hpp:48
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:106
void setReferenceDensity(unsigned regionIdx, Scalar rhoRefSolvent)
Initialize the reference density of the solvent gas for a given PVT region.
Definition SolventPvt.hpp:72
void setSolventViscosity(unsigned regionIdx, const TabulatedOneDFunction &mug)
Initialize the viscosity of the solvent gas phase.
Definition SolventPvt.hpp:80
Scalar referenceDensity(unsigned regionIdx) const
Return the reference density the solvent phase for a given PVT region.
Definition SolventPvt.hpp:128
void initEnd()
Finish initializing the oil phase PVT properties.
Definition SolventPvt.cpp:105
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition SolventPvt.hpp:99
void setSolventFormationVolumeFactor(unsigned regionIdx, const SamplingPoints &samplePoints)
Initialize the function for the formation volume factor of solvent gas.
Definition SolventPvt.cpp:91
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure) const
Returns the formation volume factor [-] of the fluid phase.
Definition SolventPvt.hpp:135
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