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

ECL problem: use spline based two-phase material laws

Browse Source 
... instead of piece-wise linear ones. This improves convergence rates
slightly.
This commit is contained in:
Andreas Lauser 2014-08-06 13:52:53 +02:00
parent 06501bb851
commit 344d613280
1 changed files with 16 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
tests/problems/eclproblem.hh
View File

@ -31,6 +31,7 @@
#include <ewoms/wells/eclwellmanager.hh> #include <ewoms/wells/eclwellmanager.hh>
#include <opm/material/fluidmatrixinteractions/PiecewiseLinearTwoPhaseMaterial.hpp> #include <opm/material/fluidmatrixinteractions/PiecewiseLinearTwoPhaseMaterial.hpp>
#include <opm/material/fluidmatrixinteractions/SplineTwoPhaseMaterial.hpp>
#include <opm/material/fluidmatrixinteractions/EclDefaultMaterial.hpp> #include <opm/material/fluidmatrixinteractions/EclDefaultMaterial.hpp>
#include <opm/material/fluidmatrixinteractions/MaterialTraits.hpp> #include <opm/material/fluidmatrixinteractions/MaterialTraits.hpp>
#include <opm/material/fluidstates/CompositionalFluidState.hpp> #include <opm/material/fluidstates/CompositionalFluidState.hpp>
@ -98,8 +99,11 @@ private:
/*nonWettingPhaseIdx=*/FluidSystem::oilPhaseIdx, /*nonWettingPhaseIdx=*/FluidSystem::oilPhaseIdx,
/*gasPhaseIdx=*/FluidSystem::gasPhaseIdx> Traits; /*gasPhaseIdx=*/FluidSystem::gasPhaseIdx> Traits;
typedef typename Opm::PiecewiseLinearTwoPhaseMaterial<OilWaterTraits> OilWaterLaw; //typedef typename Opm::PiecewiseLinearTwoPhaseMaterial<OilWaterTraits> OilWaterLaw;
typedef typename Opm::PiecewiseLinearTwoPhaseMaterial<GasOilTraits> GasOilLaw; //typedef typename Opm::PiecewiseLinearTwoPhaseMaterial<GasOilTraits> GasOilLaw;
typedef typename Opm::SplineTwoPhaseMaterial<OilWaterTraits> OilWaterLaw;
typedef typename Opm::SplineTwoPhaseMaterial<GasOilTraits> GasOilLaw;
public: public:
typedef Opm::EclDefaultMaterial<Traits, GasOilLaw, OilWaterLaw> type; typedef Opm::EclDefaultMaterial<Traits, GasOilLaw, OilWaterLaw> type;
@ -604,11 +608,12 @@ private:
// must be identical // must be identical
assert(Opm::SwofTable::numTables(swofKeyword) == Opm::SgofTable::numTables(sgofKeyword)); assert(Opm::SwofTable::numTables(swofKeyword) == Opm::SgofTable::numTables(sgofKeyword));
size_t numSatfuncTables = Opm::SwofTable::numTables(swofKeyword);
materialParams_.resize(numSatfuncTables);
typedef typename MaterialLawParams::GasOilParams GasOilParams; typedef typename MaterialLawParams::GasOilParams GasOilParams;
typedef typename MaterialLawParams::OilWaterParams OilWaterParams; typedef typename MaterialLawParams::OilWaterParams OilWaterParams;
size_t numSatfuncTables = Opm::SwofTable::numTables(swofKeyword);
materialParams_.resize(numSatfuncTables);
for (size_t tableIdx = 0; tableIdx < numSatfuncTables; ++ tableIdx) { for (size_t tableIdx = 0; tableIdx < numSatfuncTables; ++ tableIdx) {
// set the parameters of the material law for a given table // set the parameters of the material law for a given table
OilWaterParams owParams; OilWaterParams owParams;
@ -617,20 +622,19 @@ private:
Opm::SwofTable swofTable(swofKeyword, tableIdx); Opm::SwofTable swofTable(swofKeyword, tableIdx);
Opm::SgofTable sgofTable(sgofKeyword, tableIdx); Opm::SgofTable sgofTable(sgofKeyword, tableIdx);
owParams.setSwSamples(swofTable.getSwColumn()); const auto &SwColumn = swofTable.getSwColumn();
owParams.setKrwSamples(swofTable.getKrwColumn()); owParams.setKrwSamples(SwColumn, swofTable.getKrwColumn());
owParams.setKrnSamples(swofTable.getKrowColumn()); owParams.setKrnSamples(SwColumn, swofTable.getKrowColumn());
owParams.setPcnwSamples(swofTable.getPcowColumn()); owParams.setPcnwSamples(SwColumn, swofTable.getPcowColumn());
// convert the saturations from gas to oil saturations // convert the saturations from gas to oil saturations
auto SoSamples = sgofTable.getSgColumn(); auto SoSamples = sgofTable.getSgColumn();
for (size_t sampleIdx = 0; sampleIdx < SoSamples.size(); ++ sampleIdx) { for (size_t sampleIdx = 0; sampleIdx < SoSamples.size(); ++ sampleIdx) {
SoSamples[sampleIdx] = 1 - SoSamples[sampleIdx]; SoSamples[sampleIdx] = 1 - SoSamples[sampleIdx];
} }
goParams.setSwSamples(SoSamples); goParams.setKrwSamples(SoSamples, sgofTable.getKrogColumn());
goParams.setKrwSamples(sgofTable.getKrogColumn()); goParams.setKrnSamples(SoSamples, sgofTable.getKrgColumn());
goParams.setKrnSamples(sgofTable.getKrgColumn()); goParams.setPcnwSamples(SoSamples, sgofTable.getPcogColumn());
goParams.setPcnwSamples(sgofTable.getPcogColumn());
owParams.finalize(); owParams.finalize();
goParams.finalize(); goParams.finalize();