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

Merge pull request #4127 from GitPaean/equilibration_for_zero_volume_cells

Browse Source 
for the cells having zero thickness we use the cell center depth to do the equilibration
This commit is contained in:
Bård Skaflestad
2022-09-23 10:58:48 +02:00
committed by GitHub
parent 7b5ef36cd1 3bed5f41c1
commit 503965ce01
1 changed files with 12 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
ebos/equil/initstateequil.cc
View File

@@ -1864,8 +1864,18 @@ equilibrateHorizontal(const CellRange& cells,
totfrac += frac;
}
saturations /= totfrac;
pressures /= totfrac;
if (totfrac > 0.) {
saturations /= totfrac;
pressures /= totfrac;
} else {
// Fall back to centre point method for zero-thickness cells.
const auto pos = CellPos {
cell, cellCenterDepth_[cell]
};
saturations = psat.deriveSaturations(pos, eqreg, ptable);
pressures = psat.correctedPhasePressures();
}
const auto temp = this->temperature_[cell];
const auto cz = cellCenterDepth_[cell];