27 #ifndef OPM_CONSTANT_COMPRESSIBILITY_WATER_PVT_HPP
28 #define OPM_CONSTANT_COMPRESSIBILITY_WATER_PVT_HPP
31 #include <opm/input/eclipse/EclipseState/EclipseState.hpp>
41 template <
class Scalar>
47 const std::vector<Scalar>& waterReferencePressure,
48 const std::vector<Scalar>& waterReferenceFormationVolumeFactor,
49 const std::vector<Scalar>& waterCompressibility,
50 const std::vector<Scalar>& waterViscosity,
51 const std::vector<Scalar>& waterViscosibility)
52 : waterReferenceDensity_(waterReferenceDensity)
53 , waterReferencePressure_(waterReferencePressure)
54 , waterReferenceFormationVolumeFactor_(waterReferenceFormationVolumeFactor)
55 , waterCompressibility_(waterCompressibility)
56 , waterViscosity_(waterViscosity)
57 , waterViscosibility_(waterViscosibility)
64 void initFromState(
const EclipseState& eclState,
const Schedule&)
66 const auto& pvtwTable = eclState.getTableManager().getPvtwTable();
67 const auto& densityTable = eclState.getTableManager().getDensityTable();
69 assert(pvtwTable.size() == densityTable.size());
74 for (
unsigned regionIdx = 0; regionIdx <
numRegions; ++ regionIdx) {
75 waterReferenceDensity_[regionIdx] = densityTable[regionIdx].water;
77 waterReferencePressure_[regionIdx] = pvtwTable[regionIdx].reference_pressure;
78 waterReferenceFormationVolumeFactor_[regionIdx] = pvtwTable[regionIdx].volume_factor;
79 waterCompressibility_[regionIdx] = pvtwTable[regionIdx].compressibility;
80 waterViscosity_[regionIdx] = pvtwTable[regionIdx].viscosity;
81 waterViscosibility_[regionIdx] = pvtwTable[regionIdx].viscosibility;
92 waterReferenceFormationVolumeFactor_.resize(
numRegions);
97 for (
unsigned regionIdx = 0; regionIdx <
numRegions; ++regionIdx) {
111 { waterReferenceDensity_[regionIdx] = rhoRefWater; }
117 { waterReferencePressure_[regionIdx] = p; }
122 void setViscosity(
unsigned regionIdx, Scalar muw, Scalar waterViscosibility = 0.0)
124 waterViscosity_[regionIdx] = muw;
125 waterViscosibility_[regionIdx] = waterViscosibility;
132 { waterCompressibility_[regionIdx] = waterCompressibility; }
138 { waterReferenceFormationVolumeFactor_[regionIdx] = BwRef; }
144 { waterViscosibility_[regionIdx] = muComp; }
156 {
return waterReferenceDensity_.size(); }
161 template <
class Evaluation>
165 const Evaluation&)
const
167 throw std::runtime_error(
"Requested the enthalpy of water but the thermal option is not enabled");
174 template <
class Evaluation>
176 const Evaluation& temperature,
177 const Evaluation& pressure,
178 const Evaluation& saltconcentration)
const
180 Scalar BwMuwRef = waterViscosity_[regionIdx]*waterReferenceFormationVolumeFactor_[regionIdx];
183 Scalar pRef = waterReferencePressure_[regionIdx];
184 const Evaluation& Y =
185 (waterCompressibility_[regionIdx] - waterViscosibility_[regionIdx])
187 return BwMuwRef*bw/(1 + Y*(1 + Y/2));
193 template <
class Evaluation>
196 const Evaluation& pressure,
197 const Evaluation& )
const
200 Scalar pRef = waterReferencePressure_[regionIdx];
201 const Evaluation& X = waterCompressibility_[regionIdx]*(pressure - pRef);
203 Scalar BwRef = waterReferenceFormationVolumeFactor_[regionIdx];
206 return (1.0 + X*(1.0 + X/2.0))/BwRef;
209 const Scalar waterReferenceDensity(
unsigned regionIdx)
const
210 {
return waterReferenceDensity_[regionIdx]; }
212 const std::vector<Scalar>& waterReferencePressure()
const
213 {
return waterReferencePressure_; }
215 const std::vector<Scalar>& waterReferenceFormationVolumeFactor()
const
216 {
return waterReferenceFormationVolumeFactor_; }
218 const std::vector<Scalar>& waterCompressibility()
const
219 {
return waterCompressibility_; }
221 const std::vector<Scalar>& waterViscosity()
const
222 {
return waterViscosity_; }
224 const std::vector<Scalar>& waterViscosibility()
const
225 {
return waterViscosibility_; }
227 bool operator==(
const ConstantCompressibilityWaterPvt<Scalar>& data)
const
229 return this->waterReferenceDensity_ == data.waterReferenceDensity_ &&
230 this->waterReferencePressure() == data.waterReferencePressure() &&
231 this->waterReferenceFormationVolumeFactor() == data.waterReferenceFormationVolumeFactor() &&
232 this->waterCompressibility() == data.waterCompressibility() &&
233 this->waterViscosity() == data.waterViscosity() &&
234 this->waterViscosibility() == data.waterViscosibility();
238 std::vector<Scalar> waterReferenceDensity_;
239 std::vector<Scalar> waterReferencePressure_;
240 std::vector<Scalar> waterReferenceFormationVolumeFactor_;
241 std::vector<Scalar> waterCompressibility_;
242 std::vector<Scalar> waterViscosity_;
243 std::vector<Scalar> waterViscosibility_;
This class represents the Pressure-Volume-Temperature relations of the gas phase without vaporized oi...
Definition: ConstantCompressibilityWaterPvt.hpp:43
void setReferenceFormationVolumeFactor(unsigned regionIdx, Scalar BwRef)
Set the water reference formation volume factor [-].
Definition: ConstantCompressibilityWaterPvt.hpp:137
void setViscosibility(unsigned regionIdx, Scalar muComp)
Set the water "viscosibility" [1/ (Pa s)].
Definition: ConstantCompressibilityWaterPvt.hpp:143
Evaluation inverseFormationVolumeFactor(unsigned regionIdx, const Evaluation &, const Evaluation &pressure, const Evaluation &) const
Returns the formation volume factor [-] of the fluid phase.
Definition: ConstantCompressibilityWaterPvt.hpp:194
void setCompressibility(unsigned regionIdx, Scalar waterCompressibility)
Set the compressibility of the water phase.
Definition: ConstantCompressibilityWaterPvt.hpp:131
void setReferenceDensities(unsigned regionIdx, Scalar, Scalar, Scalar rhoRefWater)
Set the water reference density [kg / m^3].
Definition: ConstantCompressibilityWaterPvt.hpp:107
unsigned numRegions() const
Return the number of PVT regions which are considered by this PVT-object.
Definition: ConstantCompressibilityWaterPvt.hpp:155
Evaluation internalEnergy(unsigned, const Evaluation &, const Evaluation &, const Evaluation &) const
Returns the specific enthalpy [J/kg] of water given a set of parameters.
Definition: ConstantCompressibilityWaterPvt.hpp:162
void setReferencePressure(unsigned regionIdx, Scalar p)
Set the water reference pressure [Pa].
Definition: ConstantCompressibilityWaterPvt.hpp:116
void initEnd()
Finish initializing the water phase PVT properties.
Definition: ConstantCompressibilityWaterPvt.hpp:149
Evaluation viscosity(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:175
void setViscosity(unsigned regionIdx, Scalar muw, Scalar waterViscosibility=0.0)
Set the viscosity and "viscosibility" of the water phase.
Definition: ConstantCompressibilityWaterPvt.hpp:122