OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-29 12:33:49 -06:00
Code Issues Packages Projects Releases Wiki Activity

fix for gaswater system with water evaporation

Browse Source
This commit is contained in:
Paul Egberts 2022-08-22 17:48:43 +02:00
parent 6fe6fbf141
commit 45a968fe46
1 changed files with 15 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/models/blackoil/blackoilprimaryvariables.hh
View File

@ -498,14 +498,23 @@ public:
Scalar So = 1.0 - Sw - Sg - solventSaturation_();
Scalar So3 = 1.0 - Sg - solventSaturation_();
//water disappears
if(Sw < -eps && So3 > 0.0 && Sg > 0.0 && FluidSystem::enableVaporizedWater()) {
Scalar po = (*this)[Indices::pressureSwitchIdx];
if(Sw < -eps && FluidSystem::enableVaporizedWater()) {
Scalar pg = 0.0;
Scalar T = asImp_().temperature_();
std::array<Scalar, numPhases> pC = { 0.0 };
const MaterialLawParams& matParams = problem.materialLawParams(globalDofIdx);
computeCapillaryPressures_(pC, So3, Sg + solventSaturation_(), /*Sw=*/ 0.0, matParams);
Scalar pg = po + (pC[gasPhaseIdx] - pC[oilPhaseIdx]);
if constexpr (waterEnabled && gasEnabled && !oilEnabled) {
// twophase water-gas system
pg = (*this)[Indices::pressureSwitchIdx];
}
else if (So3 > 0.0 && Sg > 0.0) {
// threephase case
std::array<Scalar, numPhases> pC = { 0.0 };
const MaterialLawParams& matParams = problem.materialLawParams(globalDofIdx);
computeCapillaryPressures_(pC, So3, Sg + solventSaturation_(), /*Sw=*/ 0.0, matParams);
Scalar po = (*this)[Indices::pressureSwitchIdx];
pg = po + (pC[gasPhaseIdx] - pC[oilPhaseIdx]);
}
Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
T,
pg);