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

Enable live oil in initialisation.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2014-02-27 08:31:03 +01:00 committed by Andreas Lauser
parent 074dbbb599
commit 9569ea304f
1 changed files with 28 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
opm/core/simulator/initStateEquil.hpp
View File

@ -257,8 +257,33 @@ namespace Opm
sat_(props.numPhases(),
std::vector<double>(G.number_of_cells))
{
// Get the equilibration records.
const std::vector<EquilRecord> rec = getEquil(deck);
// Create (inverse) region mapping.
const RegionMapping<> eqlmap(equilnum(deck, G));
// Create Rs functions.
rs_func_.reserve(rec.size());
if (deck.hasField("DISGAS")) {
if (deck.hasField("RSVD")) {
// Rs has been specified as a function of depth.
OPM_THROW(std::runtime_error, "Cannot initialise: RSVD field not read by EclipseGridParser class.");
} else {
// Default initialisation: constant Rs below contact, saturated above.
for (size_t i = 0; i < rec.size(); ++i) {
const int cell = *(eqlmap.cells(i + 1).begin());
const double p_contact = rec[i].goc.press;
rs_func_.push_back(std::make_shared<Miscibility::RsSatAtContact>(props, cell, p_contact));
}
}
} else {
for (size_t i = 0; i < rec.size(); ++i) {
rs_func_.push_back(std::make_shared<Miscibility::NoMixing>());
}
}
// Compute phase pressures and saturations.
calcPressSat(eqlmap, rec, props, G, grav);
}
@ -272,6 +297,8 @@ namespace Opm
typedef DensityCalculator<BlackoilPropertiesInterface> RhoCalc;
typedef EquilReg<RhoCalc> EqReg;
std::vector< std::shared_ptr<Miscibility::RsFunction> > rs_func_;
PPress pp_;
PPress sat_;
@ -295,7 +322,7 @@ namespace Opm
const RhoCalc calc(props, repcell);
const EqReg eqreg(rec[r], calc,
std::make_shared<NoMix>(), std::make_shared<NoMix>(),
rs_func_[r], std::make_shared<NoMix>(),
props.phaseUsage());
const PPress press = phasePressures(G, eqreg, cells, grav);