OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-25 10:40:21 -06:00
Code Issues Packages Projects Releases Wiki Activity

Added default concentration if guess_old_solution is false.

Browse Source
This commit is contained in:
Xavier Raynaud 2012-05-22 17:21:12 +02:00
parent f8461eeec7
commit 6a12fe0208
1 changed files with 2 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/polymer/SinglePointUpwindTwoPhasePolymer.hpp
View File

@ -479,9 +479,11 @@ namespace Opm {
} else {
std::fill(x.begin(), x.end(), 0.0);
const std::vector<double>& s = state.saturation();
const std::vector<double>& c = state.concentration();
for (int cell = 0, ncell = g.number_of_cells; cell < ncell; ++cell) {
// Impose s=0.5 at next time level as an NR initial value.
x[2*cell + 0] = 0.5 - s[2*cell + 0];
x[2*cell + 1] = 1e-5 - c[cell];
}
}
}