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

Fix bug: exactly the wrong potentials were nullified.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2014-08-27 15:10:48 +02:00
parent efbf42120b
commit a48ab0188a
1 changed files with 5 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@@ -1635,11 +1635,12 @@ namespace {
// Identify the set of faces where the potential is under the // Identify the set of faces where the potential is under the
// threshold, that shall have zero flow. Storing the bool // threshold, that shall have zero flow. Storing the bool
// Array as a V (a double Array) with 1 and 0 elements. // Array as a V (a double Array) with 1 and 0 elements, a
const V low_potential = (dp.value().abs() < threshold_pressures_by_interior_face_).template cast<double>(); // 1 where flow is allowed, a 0 where it is not.
const V high_potential = (dp.value().abs() >= threshold_pressures_by_interior_face_).template cast<double>();
// Create a sparse vector that nullifies the low potential elements. // Create a sparse vector that nullifies the low potential elements.
const M nullify_low_potential = spdiag(low_potential); const M keep_high_potential = spdiag(high_potential);
// The threshold modification must have the sign that reduces the // The threshold modification must have the sign that reduces the
// absolute value of the potential. // absolute value of the potential.
@@ -1647,7 +1648,7 @@ namespace {
const V threshold_modification = sign_for_mod * threshold_pressures_by_interior_face_; const V threshold_modification = sign_for_mod * threshold_pressures_by_interior_face_;
// Modify potential and nullify where appropriate. // Modify potential and nullify where appropriate.
dp = nullify_low_potential * (dp + threshold_modification); dp = keep_high_potential * (dp + threshold_modification);
} }