OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

THPRES: change the sign of the gravity correction

Browse Source 
we're correcting the pressure at the cell center depths to get the
pressure at the face depth, not the other way around. This is
confusing...

thanks to [at]totto82 for discovering this.
This commit is contained in:
Andreas Lauser 2015-11-05 12:14:35 +01:00
parent 6a9e11a4e1
commit f4fcb616be
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/core/utility/thresholdPressures.hpp
View File

@ -278,8 +278,8 @@ namespace Opm
double p1 = phasePressure[phaseIdx][c1];
double p2 = phasePressure[phaseIdx][c2];
p1 += rhoAvg*gravity*(z1 - zAvg);
p2 += rhoAvg*gravity*(z2 - zAvg);
p1 += rhoAvg*gravity*(zAvg - z1);
p2 += rhoAvg*gravity*(zAvg - z2);
if ((p1 > p2 && s1 > 0.0) || (p2 > p1 && s2 > 0.0))
maxDp[barrierId] = std::max(maxDp[barrierId], std::abs(p1 - p2));