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

fix some bugs in WetGasPvt

Browse Source 
so far, this has never been tested on a deck which features vaporized oil.
This commit is contained in:
Andreas Lauser 2015-08-27 19:07:08 +02:00
parent 5edc12bef2
commit 543c3cb603
1 changed files with 9 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/material/fluidsystems/blackoilpvt/WetGasPvt.hpp
View File

@ -344,8 +344,8 @@ private:
{
const LhsEval& Rv =
XgO/(1 - XgO)
* (BlackOilFluidSystem::referenceDensity(oilPhaseIdx, regionIdx)
/ BlackOilFluidSystem::referenceDensity(gasPhaseIdx, regionIdx));
* (BlackOilFluidSystem::referenceDensity(gasPhaseIdx, regionIdx)
/ BlackOilFluidSystem::referenceDensity(oilPhaseIdx, regionIdx));
const LhsEval& invBg = inverseGasB_[regionIdx].eval(pressure, Rv, /*extrapolate=*/true);
const LhsEval& invMugBg = inverseGasBMu_[regionIdx].eval(pressure, Rv, /*extrapolate=*/true);
@ -366,15 +366,16 @@ private:
Scalar rhogRef = BlackOilFluidSystem::referenceDensity(gasPhaseIdx, regionIdx);
const LhsEval& Bg = formationVolumeFactor_(regionIdx, temperature, pressure, XgO);
LhsEval rhoo = rhooRef/Bg;
LhsEval rhog = rhogRef/Bg;
// the oil formation volume factor just represents the partial density of the gas
// component in the gas phase. to get the total density of the phase, we have to
// add the partial density of the oil component.
const LhsEval& Rv = XgO/(1 - XgO) * (rhooRef/rhogRef);
rhoo += (rhogRef*Rv)/Bg;
const LhsEval& Rv = XgO/(1 - XgO) * (rhogRef/rhooRef);
rhog += (rhogRef*Rv)/Bg;
return rhoo;
return rhog;
}
/*!
@ -388,8 +389,8 @@ private:
{
const LhsEval& Rv =
XgO/(1-XgO)
*BlackOilFluidSystem::referenceDensity(oilPhaseIdx, regionIdx)
/BlackOilFluidSystem::referenceDensity(gasPhaseIdx, regionIdx);
*(BlackOilFluidSystem::referenceDensity(gasPhaseIdx, regionIdx)
/BlackOilFluidSystem::referenceDensity(oilPhaseIdx, regionIdx));
return 1.0 / inverseGasB_[regionIdx].eval(pressure, Rv, /*extrapolate=*/true);
}