OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-16 14:51:56 -06:00
Code Issues Packages Projects Releases Wiki Activity

Use cell depths based on averaging cell corners

Browse Source
This commit is contained in:
Tor Harald Sandve 2015-11-24 10:17:12 +01:00 committed by Andreas Lauser
parent ca902168ed
commit d77e33b5dc
1 changed files with 6 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/core/simulator/initStateEquil_impl.hpp
View File

@ -313,7 +313,6 @@ namespace Opm
const CellRange& cells,
std::vector<double>& p )
{
const int nd = UgGridHelpers::dimensions(G);
enum { up = 0, down = 1 };
@ -324,7 +323,7 @@ namespace Opm
{
assert (c < p.size());
const double z = UgGridHelpers::cellCentroidCoordinate(G, *ci, nd-1);
const double z = UgGridHelpers::cellCenterDepth(G, *ci);
p[c] = (z < split) ? f[up](z) : f[down](z);
}
}
@ -703,10 +702,8 @@ namespace Opm
double sw = 0.0;
if (water) {
if (isConstPc(props,waterpos,cell)){
const int nd = UgGridHelpers::dimensions(G);
const double cellDepth = UgGridHelpers::cellCentroidCoordinate(G,
cell,
nd-1);
const double cellDepth = UgGridHelpers::cellCenterDepth(G,
cell);
sw = satFromDepth(props,cellDepth,reg.zwoc(),waterpos,cell,false);
phase_saturations[waterpos][local_index] = sw;
}
@ -725,10 +722,8 @@ namespace Opm
double sg = 0.0;
if (gas) {
if (isConstPc(props,gaspos,cell)){
const int nd = UgGridHelpers::dimensions(G);
const double cellDepth = UgGridHelpers::cellCentroidCoordinate(G,
cell,
nd-1);
const double cellDepth = UgGridHelpers::cellCenterDepth(G,
cell);
sg = satFromDepth(props,cellDepth,reg.zgoc(),gaspos,cell,true);
phase_saturations[gaspos][local_index] = sg;
}
@ -829,7 +824,7 @@ namespace Opm
std::vector<double> rs(cells.size());
int count = 0;
for (auto it = cells.begin(); it != cells.end(); ++it, ++count) {
const double depth = UgGridHelpers::cellCentroidCoordinate(grid, *it, 2);
const double depth = UgGridHelpers::cellCenterDepth(grid, *it);
rs[count] = rs_func(depth, oil_pressure[count], temperature[count], gas_saturation[count]);
}
return rs;