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

Fixed a bug in non-linear Gauss-Seidel column segregation solver.

Browse Source 
With columns containing an odd number of cells, middle cell would
get solved twice, and the second time with incorrect s0.
We still solve twice, but with correct s0.
This commit is contained in:
Atgeirr Flø Rasmussen
2012-05-21 11:11:06 +02:00
parent 5abf755e2d
commit fb5113836c
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
opm/core/transport/reorder/TransportModelTwophase.cpp
View File

@@ -584,8 +584,8 @@ namespace Opm
double old_s[2] = { saturation_[cells[ci]], double old_s[2] = { saturation_[cells[ci]],
saturation_[cells[ci2]] }; saturation_[cells[ci2]] };
saturation_[cells[ci]] = s0_[ci]; saturation_[cells[ci]] = s0_[ci];
saturation_[cells[ci2]] = s0_[ci2];
solveSingleCellGravity(cells, ci, &col_gravflux[0]); solveSingleCellGravity(cells, ci, &col_gravflux[0]);
saturation_[cells[ci2]] = s0_[ci2];
solveSingleCellGravity(cells, ci2, &col_gravflux[0]); solveSingleCellGravity(cells, ci2, &col_gravflux[0]);
max_s_change = std::max(max_s_change, std::max(std::fabs(saturation_[cells[ci]] - old_s[0]), max_s_change = std::max(max_s_change, std::max(std::fabs(saturation_[cells[ci]] - old_s[0]),
std::fabs(saturation_[cells[ci2]] - old_s[1]))); std::fabs(saturation_[cells[ci2]] - old_s[1])));