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

Apply dpmaxrel to the update of bhp

Browse Source 
And also change the default value for dpmaxrel.
This commit is contained in:
Kai Bao
2014-05-16 13:46:00 +02:00
parent b001c580b9
commit 55b0164d4a
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@@ -1224,7 +1224,7 @@ namespace {
assert(varstart == dx.size());
// Pressure update.
const double dpmaxrel = 0.8;
const double dpmaxrel = 0.2;
const V p_old = Eigen::Map<const V>(&state.pressure()[0], nc, 1);
const V absdpmax = dpmaxrel*p_old.abs();
const V dp_limited = sign(dp) * dp.abs().min(absdpmax);
@@ -1417,7 +1417,8 @@ namespace {
// Bhp update.
const V bhp_old = Eigen::Map<const V>(&well_state.bhp()[0], nw, 1);
const V bhp = bhp_old - dbhp;
const V dbhp_limited = sign(dbhp) * dbhp.abs().min(bhp_old.abs()*dpmaxrel);
const V bhp = bhp_old - dbhp_limited;
std::copy(&bhp[0], &bhp[0] + bhp.size(), well_state.bhp().begin());
}