OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #707 from plgbrts/gas-water-fix

Browse Source 
fix for gaswater system with water evaporation
This commit is contained in:
Bård Skaflestad 2022-08-23 15:14:23 +02:00 committed by GitHub
commit 758473ffd8
1 changed files with 17 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -498,14 +498,25 @@ 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);