OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-02 05:49:09 -06:00
Code Issues Packages Projects Releases Wiki Activity

set water saturations for numerical aquifer cells to be 1.

Browse Source
This commit is contained in:
Kai Bao 2020-11-09 14:04:01 +01:00
parent 1d5b60a600
commit 80c0073a69
2 changed files with 10 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
ebos/eclproblem.hh
View File

@ -2735,13 +2735,6 @@ private:
auto& elemFluidState = initialFluidStates_[elemIdx];
elemFluidState.assign(equilInitializer.initialFluidState(elemIdx));
}
const auto& aquifer_cells = simulator.vanguard().eclState().aquifer().numericalAquifers().aquiferCells();
const auto& vanguard = simulator.vanguard();
size_t numCartDof = vanguard.cartesianSize();
std::vector<int> cartesianToCompressedElemIdx(numCartDof, -1);
for (unsigned elemIdx = 0; elemIdx < numElems; ++elemIdx)
cartesianToCompressedElemIdx[vanguard.cartesianIndex(elemIdx)] = elemIdx;
}
void readEclRestartSolution_()

10
ebos/equil/initstateequil.hh
View File

@ -62,6 +62,7 @@
#include <type_traits>
#include <utility>
#include <vector>
#include <string>
namespace Opm {
@ -1937,6 +1938,15 @@ 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);
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 << " , "
<< "WATER SATURATION IS SET TO BE UNITY";
OpmLog::info(ss.str());
}
pressures = psat.correctedPhasePressures();
const auto temp = this->temperature_[cell];