27#ifndef OPM_CONSTANT_COMPRESSIBILITY_WATER_PVT_HPP
28#define OPM_CONSTANT_COMPRESSIBILITY_WATER_PVT_HPP
45template <
class Scalar>
61 waterReferenceFormationVolumeFactor_.resize(
numRegions);
66 for (
unsigned regionIdx = 0; regionIdx <
numRegions; ++regionIdx) {
80 { waterReferenceDensity_[regionIdx] = rhoRefWater; }
86 { waterReferencePressure_[regionIdx] = p; }
91 void setViscosity(
unsigned regionIdx, Scalar muw, Scalar waterViscosibility = 0.0)
93 waterViscosity_[regionIdx] = muw;
94 waterViscosibility_[regionIdx] = waterViscosibility;
101 { waterCompressibility_[regionIdx] = waterCompressibility; }
107 { waterReferenceFormationVolumeFactor_[regionIdx] = BwRef; }
113 { waterViscosibility_[regionIdx] = muComp; }
125 {
return waterReferenceDensity_.size(); }
130 template <
class Evaluation>
135 const Evaluation&)
const
137 throw std::runtime_error(
"Requested the enthalpy of water but the thermal option is not enabled");
144 template <
class Evaluation>
146 const Evaluation& temperature,
147 const Evaluation& pressure,
148 const Evaluation& saltconcentration)
const
150 Scalar BwMuwRef = waterViscosity_[regionIdx]*waterReferenceFormationVolumeFactor_[regionIdx];
153 Scalar pRef = waterReferencePressure_[regionIdx];
154 const Evaluation& Y =
155 (waterCompressibility_[regionIdx] - waterViscosibility_[regionIdx])
157 return BwMuwRef*bw/(1 + Y*(1 + Y/2));
162 template <
class Evaluation>
164 const Evaluation& temperature,
165 const Evaluation& pressure,
166 const Evaluation& Rsw,
167 const Evaluation& saltconcentration)
const
169 Scalar BwMuwRef = waterViscosity_[regionIdx]*waterReferenceFormationVolumeFactor_[regionIdx];
172 Scalar pRef = waterReferencePressure_[regionIdx];
173 const Evaluation& Y =
174 (waterCompressibility_[regionIdx] - waterViscosibility_[regionIdx])
176 return BwMuwRef*bw/(1 + Y*(1 + Y/2));
182 template <
class Evaluation>
184 const Evaluation& temperature,
185 const Evaluation& pressure,
186 const Evaluation& saltconcentration)
const
188 Evaluation Rsw = 0.0;
195 template <
class Evaluation>
198 const Evaluation& pressure,
200 const Evaluation& )
const
202 Scalar pRef = waterReferencePressure_[regionIdx];
203 const Evaluation& X = waterCompressibility_[regionIdx]*(pressure - pRef);
205 Scalar BwRef = waterReferenceFormationVolumeFactor_[regionIdx];
208 return (1.0 + X*(1.0 + X/2.0))/BwRef;
211 template <
class Evaluation>
212 void inverseBAndMu(Evaluation& bw, Evaluation& muW,
unsigned regionIdx,
214 const Evaluation& pressure,
216 const Evaluation& )
const{
217 inverseBAndMu(bw, muW, regionIdx,pressure);
220 template <
class Evaluation>
221 void inverseBAndMu(Evaluation& bw, Evaluation& muW,
unsigned regionIdx,
222 const Evaluation& pressure)
const
224 Scalar pRef = waterReferencePressure_[regionIdx];
225 const Evaluation& X = waterCompressibility_[regionIdx]*(pressure - pRef);
227 Scalar BwRef = waterReferenceFormationVolumeFactor_[regionIdx];
230 bw = (1.0 + X*(1.0 + X/2.0))/BwRef;
232 Scalar BwMuwRef = waterViscosity_[regionIdx]*BwRef;
234 const Evaluation& Y =
235 (waterCompressibility_[regionIdx] - waterViscosibility_[regionIdx])
237 muW = BwMuwRef*bw/(1 + Y*(1 + Y/2));
246 template <
class Evaluation>
250 const Evaluation& )
const
253 template <
class Evaluation>
254 Evaluation diffusionCoefficient(
const Evaluation& ,
258 throw std::runtime_error(
"Not implemented: The PVT model does not provide a diffusionCoefficient()");
264 template <
class Evaluation>
268 const Evaluation& )
const
271 Scalar waterReferenceDensity(
unsigned regionIdx)
const
272 {
return waterReferenceDensity_[regionIdx]; }
274 const std::vector<Scalar>& waterReferencePressure()
const
275 {
return waterReferencePressure_; }
277 const std::vector<Scalar>& waterReferenceFormationVolumeFactor()
const
278 {
return waterReferenceFormationVolumeFactor_; }
280 const std::vector<Scalar>& waterCompressibility()
const
281 {
return waterCompressibility_; }
283 const std::vector<Scalar>& waterViscosity()
const
284 {
return waterViscosity_; }
286 const std::vector<Scalar>& waterViscosibility()
const
287 {
return waterViscosibility_; }
290 std::vector<Scalar> waterReferenceDensity_;
291 std::vector<Scalar> waterReferencePressure_;
292 std::vector<Scalar> waterReferenceFormationVolumeFactor_;
293 std::vector<Scalar> waterCompressibility_;
294 std::vector<Scalar> waterViscosity_;
295 std::vector<Scalar> waterViscosibility_;
This class represents the Pressure-Volume-Temperature relations of the gas phase without vaporized oi...
Definition: ConstantCompressibilityWaterPvt.hpp:47
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BwRef)
Set the water reference formation volume factor [-].
Definition: ConstantCompressibilityWaterPvt.hpp:106
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the water "viscosibility" [1/ (Pa s)].
Definition: ConstantCompressibilityWaterPvt.hpp:112
Evaluation saturationPressure(unsigned, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the saturation pressure of the water phase [Pa] depending on its mass fraction of the gas com...
Definition: ConstantCompressibilityWaterPvt.hpp:247
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure, const Evaluation &, const Evaluation &) const
Returns the formation volume factor [-] of the fluid phase.
Definition: ConstantCompressibilityWaterPvt.hpp:196
void setCompressibility(unsigned regionIdx, Scalar waterCompressibility)
Set the compressibility of the water phase.
Definition: ConstantCompressibilityWaterPvt.hpp:100
void setReferenceDensities(unsigned regionIdx, Scalar, Scalar, Scalar rhoRefWater)
Set the water reference density [kg / m^3].
Definition: ConstantCompressibilityWaterPvt.hpp:76
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: ConstantCompressibilityWaterPvt.hpp:124
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the water reference pressure [Pa].
Definition: ConstantCompressibilityWaterPvt.hpp:85
void initEnd()
Finish initializing the water phase PVT properties.
Definition: ConstantCompressibilityWaterPvt.hpp:118
Evaluation viscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &Rsw, const Evaluation &saltconcentration) const
Returns the dynamic viscosity [Pa s] of the fluid phase given a set of parameters.
Definition: ConstantCompressibilityWaterPvt.hpp:163
void setViscosity(unsigned regionIdx, Scalar muw, Scalar waterViscosibility=0.0)
Set the viscosity and "viscosibility" of the water phase.
Definition: ConstantCompressibilityWaterPvt.hpp:91
Evaluation saturatedInverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &saltconcentration) const
Returns the formation volume factor [-] of the fluid phase.
Definition: ConstantCompressibilityWaterPvt.hpp:183
Evaluation saturatedViscosity(unsigned regionIdx, const Evaluation &temperature, const Evaluation &pressure, const Evaluation &saltconcentration) const
Returns the dynamic viscosity [Pa s] of the fluid phase given a set of parameters.
Definition: ConstantCompressibilityWaterPvt.hpp:145
Evaluation internalEnergy(unsigned, const Evaluation &, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the specific enthalpy [J/kg] of water given a set of parameters.
Definition: ConstantCompressibilityWaterPvt.hpp:131
Evaluation saturatedGasDissolutionFactor(unsigned, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the gas dissolution factor [m^3/m^3] of the water phase.
Definition: ConstantCompressibilityWaterPvt.hpp:265
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