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

Compute phase pressures in subset of cells

Browse Source 
This commit adds support for assigning the initial phase pressure
distribution to a subset of the total grid cells.  This is needed in
order to fully support equilibration regions.  The existing region
support (template parameter 'Region' in function 'phasePressures()')
was only used/needed to define PVT property (specifically, the fluid
phase density) calculator pertaining to a particular equilibration
region.
This commit is contained in:
Bård Skaflestad
2014-01-17 17:43:27 +01:00
parent 4ea3e7ed53
commit e84e6ee4bb
3 changed files with 87 additions and 28 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
tests/test_equil.cpp
View File

@@ -75,8 +75,11 @@ BOOST_AUTO_TEST_CASE (PhasePressure)
Opm::equil::miscibility::NoMixing(),
props.phaseUsage());
std::vector<int> cells(G->number_of_cells);
std::iota(cells.begin(), cells.end(), 0);
const double grav = 10;
const PPress ppress = Opm::equil::phasePressures(*G, region, grav);
const PPress ppress = Opm::equil::phasePressures(*G, region, cells, grav);
const int first = 0, last = G->number_of_cells - 1;
const double reltol = 1.0e-8;