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

Don't set lower bound on network update

Browse Source
This commit is contained in:
Stein Krogstad 2024-01-26 14:58:06 +01:00
parent 2626fbb84b
commit a9ab62e6b8
1 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/simulators/wells/BlackoilWellModelGeneric.cpp
View File

@ -1172,15 +1172,11 @@ updateNetworkPressures(const int reportStepIdx)
// due to the fact our nodal pressure calculation is somewhat explicit // due to the fact our nodal pressure calculation is somewhat explicit
// TODO: the following parameters are subject to adjustment for optimization purpose // TODO: the following parameters are subject to adjustment for optimization purpose
constexpr double upper_update_bound = 5.0 * unit::barsa; constexpr double upper_update_bound = 5.0 * unit::barsa;
constexpr double lower_update_bound = 0.05 * unit::barsa;
// relative dampening factor based on update value // relative dampening factor based on update value
constexpr double damping_factor = 0.1; constexpr double damping_factor = 0.1;
const double allowed_change = std::max(std::min(damping_factor * std::abs(change), upper_update_bound), const double damped_change = std::min(damping_factor * std::abs(change), upper_update_bound);
lower_update_bound); const double sign = change > 0 ? 1. : -1.;
if (std::abs(change) > allowed_change) { node_pressures_[name] = pressure + sign * damped_change;
const double sign = change > 0 ? 1. : -1.;
node_pressures_[name] = pressure + sign * allowed_change;
}
} }
} else { } else {
for (const auto& [name, pressure]: node_pressures_) { for (const auto& [name, pressure]: node_pressures_) {