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

black oil fluid system: keep the saturation pressure spline consistent

Browse Source 
It was not updated if the set PvtoTable was used which then lead to
segfaults.
This commit is contained in:
Andreas Lauser
2014-07-25 14:55:49 +02:00
parent a3c8632358
commit 9f3010eac3
1 changed files with 22 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

38
opm/material/fluidsystems/BlackOilFluidSystem.hpp
View File

@@ -128,6 +128,7 @@ public:
saturatedTable->getPressureColumn(),
saturatedTable->getGasSolubilityColumn(),
/*type=*/Spline::Monotonic);
updateSaturationPressureSpline_();
Scalar rhooRef = surfaceDensity(oilPhaseIdx);
Scalar rhogRef = surfaceDensity(gasPhaseIdx);
@@ -227,22 +228,7 @@ public:
gasDissolutionFactorSpline_.setContainerOfTuples(samplePoints, /*type=*/Spline::Monotonic);
assert(gasDissolutionFactorSpline_.monotonic());
// create the spline representing saturation pressure
// depending of the mass fraction in gas
int n = gasDissolutionFactorSpline_.numSamples()*5;
int delta =
(gasDissolutionFactorSpline_.xMax() - gasDissolutionFactorSpline_.xMin())/(n + 1);
SplineSamplingPoints pSatSamplePoints;
Scalar X_oG = 0;
for (int i=0; i <= n; ++ i) {
Scalar pSat = gasDissolutionFactorSpline_.xMin() + i*delta;
X_oG = saturatedOilGasMassFraction(pSat);
std::pair<Scalar, Scalar> val(X_oG, pSat);
pSatSamplePoints.push_back(val);
}
saturationPressureSpline_.setContainerOfTuples(pSatSamplePoints, /*type=*/Spline::Monotonic);
updateSaturationPressureSpline_();
}
/*!
@@ -803,6 +789,26 @@ public:
}
private:
static void updateSaturationPressureSpline_()
{
// create the spline representing saturation pressure
// depending of the mass fraction in gas
int n = gasDissolutionFactorSpline_.numSamples()*5;
int delta =
(gasDissolutionFactorSpline_.xMax() - gasDissolutionFactorSpline_.xMin())/(n + 1);
SplineSamplingPoints pSatSamplePoints;
Scalar X_oG = 0;
for (int i=0; i <= n; ++ i) {
Scalar pSat = gasDissolutionFactorSpline_.xMin() + i*delta;
X_oG = saturatedOilGasMassFraction(pSat);
std::pair<Scalar, Scalar> val(X_oG, pSat);
pSatSamplePoints.push_back(val);
}
saturationPressureSpline_.setContainerOfTuples(pSatSamplePoints, /*type=*/Spline::Monotonic);
}
static Scalar gasDensity_(Scalar pressure)
{
// gas formation volume factor at reservoir pressure