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

Bugfix: ensure all cells have a valid saturation initially.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-03-29 21:10:14 +02:00
parent 5e78fc1c9f
commit ed5a7802b3
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
opm/core/utility/initState_impl.hpp
View File

@ -161,6 +161,12 @@ namespace Opm
} }
state.init(grid); state.init(grid);
const int num_cells = props.numCells(); const int num_cells = props.numCells();
// By default: initialise water saturation to minimum everywhere.
std::vector<int> all_cells(num_cells);
for (int i = 0; i < num_cells; ++i) {
all_cells[i] = i;
}
state.setWaterSat(all_cells, props, State::MinSat);
const bool convection_testcase = param.getDefault("convection_testcase", false); const bool convection_testcase = param.getDefault("convection_testcase", false);
const bool segregation_testcase = param.getDefault("segregation_testcase", false); const bool segregation_testcase = param.getDefault("segregation_testcase", false);
if (convection_testcase) { if (convection_testcase) {
@ -222,12 +228,7 @@ namespace Opm
const double ref_z = grid.cell_centroids[0 + grid.dimensions - 1]; const double ref_z = grid.cell_centroids[0 + grid.dimensions - 1];
initHydrostaticPressure(grid, dens, ref_z, gravity, ref_z, ref_p, state); initHydrostaticPressure(grid, dens, ref_z, gravity, ref_z, ref_p, state);
} else { } else {
// By default: initialise water saturation to minimum everywhere. // Use default: water saturation is minimum everywhere.
std::vector<int> all_cells(num_cells);
for (int i = 0; i < num_cells; ++i) {
all_cells[i] = i;
}
state.setWaterSat(all_cells, props, State::MinSat);
// Initialise pressure to hydrostatic state. // Initialise pressure to hydrostatic state.
const double ref_p = param.getDefault("ref_pressure", 100)*unit::barsa; const double ref_p = param.getDefault("ref_pressure", 100)*unit::barsa;
const double rho = props.density()[1]; const double rho = props.density()[1];