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

Merge pull request #4712 from vkip/network_solver_params

Browse Source 
Parameterized iteration limits in network solver
This commit is contained in:
Kai Bao 2023-06-16 13:43:07 +02:00 committed by GitHub
commit 8bfec82fdd
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 30 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
opm/simulators/flow/BlackoilModelParametersEbos.hpp
View File

@ -168,6 +168,16 @@ template<class TypeTag, class MyTypeTag>
struct UseAverageDensityMsWells { struct UseAverageDensityMsWells {
using type = UndefinedProperty; using type = UndefinedProperty;
}; };
// Network solver parameters
template<class TypeTag, class MyTypeTag>
struct NetworkMaxStrictIterations {
using type = UndefinedProperty;
};
template<class TypeTag, class MyTypeTag>
struct NetworkMaxIterations {
using type = UndefinedProperty;
};
template<class TypeTag> template<class TypeTag>
struct DbhpMaxRel<TypeTag, TTag::FlowModelParameters> { struct DbhpMaxRel<TypeTag, TTag::FlowModelParameters> {
@ -316,6 +326,15 @@ template<class TypeTag>
struct UseAverageDensityMsWells<TypeTag, TTag::FlowModelParameters> { struct UseAverageDensityMsWells<TypeTag, TTag::FlowModelParameters> {
static constexpr bool value = false; static constexpr bool value = false;
}; };
// Network solver parameters
template<class TypeTag>
struct NetworkMaxStrictIterations<TypeTag, TTag::FlowModelParameters> {
static constexpr int value = 100;
};
template<class TypeTag>
struct NetworkMaxIterations<TypeTag, TTag::FlowModelParameters> {
static constexpr int value = 200;
};
@ -436,6 +455,11 @@ namespace Opm
/// Whether to approximate segment densities by averaging over segment and its outlet /// Whether to approximate segment densities by averaging over segment and its outlet
bool use_average_density_ms_wells_; bool use_average_density_ms_wells_;
/// Maximum number of iterations in the network solver before relaxing tolerance
int network_max_strict_iterations_;
/// Maximum number of iterations in the network solver before giving up
int network_max_iterations_;
/// Construct from user parameters or defaults. /// Construct from user parameters or defaults.
@ -474,6 +498,8 @@ namespace Opm
max_number_of_well_switches_ = EWOMS_GET_PARAM(TypeTag, int, MaximumNumberOfWellSwitches); max_number_of_well_switches_ = EWOMS_GET_PARAM(TypeTag, int, MaximumNumberOfWellSwitches);
use_average_density_ms_wells_ = EWOMS_GET_PARAM(TypeTag, bool, UseAverageDensityMsWells); use_average_density_ms_wells_ = EWOMS_GET_PARAM(TypeTag, bool, UseAverageDensityMsWells);
deck_file_name_ = EWOMS_GET_PARAM(TypeTag, std::string, EclDeckFileName); deck_file_name_ = EWOMS_GET_PARAM(TypeTag, std::string, EclDeckFileName);
network_max_strict_iterations_ = EWOMS_GET_PARAM(TypeTag, int, NetworkMaxStrictIterations);
network_max_iterations_ = EWOMS_GET_PARAM(TypeTag, int, NetworkMaxIterations);
} }
static void registerParameters() static void registerParameters()
@ -512,6 +538,8 @@ namespace Opm
EWOMS_REGISTER_PARAM(TypeTag, bool, EnableWellOperabilityCheckIter, "Enable the well operability checking during iterations"); EWOMS_REGISTER_PARAM(TypeTag, bool, EnableWellOperabilityCheckIter, "Enable the well operability checking during iterations");
EWOMS_REGISTER_PARAM(TypeTag, int, MaximumNumberOfWellSwitches, "Maximum number of times a well can switch to the same control"); EWOMS_REGISTER_PARAM(TypeTag, int, MaximumNumberOfWellSwitches, "Maximum number of times a well can switch to the same control");
EWOMS_REGISTER_PARAM(TypeTag, bool, UseAverageDensityMsWells, "Approximate segment densitities by averaging over segment and its outlet"); EWOMS_REGISTER_PARAM(TypeTag, bool, UseAverageDensityMsWells, "Approximate segment densitities by averaging over segment and its outlet");
EWOMS_REGISTER_PARAM(TypeTag, int, NetworkMaxStrictIterations, "Maximum iterations in network solver before relaxing tolerance");
EWOMS_REGISTER_PARAM(TypeTag, int, NetworkMaxIterations, "Maximum number of iterations in the network solver before giving up");
} }
}; };
} // namespace Opm } // namespace Opm

4
opm/simulators/wells/BlackoilWellModel_impl.hpp
View File

@ -972,9 +972,9 @@ namespace Opm {
bool do_network_update = true; bool do_network_update = true;
bool well_group_control_changed = false; bool well_group_control_changed = false;
// after certain number of the iterations, we use relaxed tolerance for the network update // after certain number of the iterations, we use relaxed tolerance for the network update
constexpr size_t iteration_to_relax = 100; const size_t iteration_to_relax = param_.network_max_strict_iterations_;
// after certain number of the iterations, we terminate // after certain number of the iterations, we terminate
constexpr size_t max_iteration = 200; const size_t max_iteration = param_.network_max_iterations_;
std::size_t network_update_iteration = 0; std::size_t network_update_iteration = 0;
while (do_network_update) { while (do_network_update) {
if (network_update_iteration == iteration_to_relax) { if (network_update_iteration == iteration_to_relax) {