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

using the specified pressure for numerical aquifer

Browse Source 
NOT totally sure how to use this pressure value, do we need to consider
the capillary pressure?
This commit is contained in:
Kai Bao 2020-11-10 14:52:50 +01:00
parent 80c0073a69
commit ef67d9a9e7
1 changed files with 9 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
ebos/equil/initstateequil.hh
View File

@ -1940,7 +1940,7 @@ private:
saturations = psat.deriveSaturations(pos, eqreg, ptable);
if (aquifer.hasCell(global_index)) {
saturations = {0.0, 0.0, 1.0};
const auto& aqu_cell = aquifer.getCell(global_index);
const auto &aqu_cell = aquifer.getCell(global_index);
std::ostringstream ss;
ss << "FOR AQUIFER CELL AT { " << aqu_cell.I + 1 << " " << aqu_cell.J + 1 << " "
<< aqu_cell.J + 1 << " } OF NUMERICAL AQUIFER " << aqu_cell.aquifer_id << " , "
@ -1948,6 +1948,14 @@ private:
OpmLog::info(ss.str());
}
pressures = psat.correctedPhasePressures();
if (aquifer.hasCell(global_index)) {
const auto &aqu_cell = aquifer.getCell(global_index);
// TODO: NOT totally sure what we should do here to empoly the pressure specified by AQUNUM
if (aqu_cell.init_pressure > 0.) {
const double pres = aqu_cell.init_pressure;
pressures = {pres, pres, pres};
}
}
const auto temp = this->temperature_[cell];