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

Fixed problem induced by s=0 in gravitation transport solver.

Browse Source
This commit is contained in:
Xavier Raynaud 2012-05-15 16:09:07 +02:00
parent 7d54d1f31b
commit c0368fa61d
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
opm/polymer/GravityColumnSolverPolymer_impl.hpp
View File

@ -180,11 +180,11 @@ namespace Opm
// not with arbitrary problem models.
int col_size = column_cells.size();
if (col_size == 1) {
sol_vec[2*column_cells[0] + 0] = 0.0;
sol_vec[2*column_cells[0] + 1] = 0.0;
return;
}
// if (col_size == 1) {
// sol_vec[2*column_cells[0] + 0] = 0.0;
// sol_vec[2*column_cells[0] + 1] = 0.0;
// return;
// }
StateWithZeroFlux state(s, c, cmax); // This holds s by reference.
@ -243,7 +243,11 @@ namespace Opm
hm[bmc(2*ci + 0, 2*ci + 0)] += dF[0];
hm[bmc(2*ci + 0, 2*ci + 1)] += dF[1];
hm[bmc(2*ci + 1, 2*ci + 0)] += dF[2];
hm[bmc(2*ci + 1, 2*ci + 1)] += dF[3];
if (std::abs(dF[3]) < 1e-12) {
hm[bmc(2*ci + 1, 2*ci + 1)] += 1e-12;
} else {
hm[bmc(2*ci + 1, 2*ci + 1)] += dF[3];
}
rhs[2*ci + 0] += F[0];
rhs[2*ci + 1] += F[1];