OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-11-02 20:30:55 -06:00
Code Issues Packages Projects Releases Wiki Activity

update network when group control or gaslift changes

Browse Source
This commit is contained in:
Tor Harald Sandve committed 2025-03-13 08:51:28 +01:00
1 parent 43421a59ec
commit aa64e853bc
1 file changed
+20 -2
Download Patch File Download Diff File Expand all files Collapse all files
opm/simulators/wells/BlackoilWellModel_impl.hpp
+20 -2
View File
@@ -1198,7 +1198,21 @@ namespace Opm {
// after certain number of the iterations, we terminate
const std::size_t max_iteration = param_.network_max_iterations_;
std::size_t network_update_iteration = 0;
const int episodeIdx = simulator_.episodeIndex();
while (do_network_update) {
if (network_update_iteration >= max_iteration ) {
// only output to terminal if we at the last newton iterations where we try to balance the network.
if (this->shouldBalanceNetwork(episodeIdx, iterationIdx + 1)) {
local_deferredLogger.debug("maximum of " + std::to_string(max_iteration) + " network iterations has been used, we stop the update \n "
" and try again after the next newton iteration.");
} else {
if (this->terminal_output_) {
local_deferredLogger.info("maximum of " + std::to_string(max_iteration) + " network iterations has been used and we stop the update. "
"The simulator will continue with unconverged network results.");
}
}
break;
}
if (this->terminal_output_ && (network_update_iteration == iteration_to_relax) ) {
local_deferredLogger.info(" we begin using relaxed tolerance for network update now after " + std::to_string(iteration_to_relax) + " iterations ");
}
@@ -1230,7 +1244,7 @@ namespace Opm {
DeferredLogger& local_deferredLogger)
{
OPM_TIMEFUNCTION();
auto [well_group_control_changed, more_network_update] =
auto [well_group_control_changed, inner_network_not_balanced] =
updateWellControls(mandatory_network_balance,
local_deferredLogger,
relax_network_tolerance);
@@ -1272,7 +1286,11 @@ namespace Opm {
pot,
local_deferredLogger);
}
// we need to re-iterate the network when the well group controls changed or gaslift/alq is changed or
// the inner iterations are did not converge
const int iterationIdx = simulator_.model().newtonMethod().numIterations();
bool more_network_update = this->shouldBalanceNetwork(reportStepIdx, iterationIdx) &&
(inner_network_not_balanced || well_group_control_changed || alq_updated);
return {more_network_update, well_group_control_changed};
}