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 #5929 from steink/pass-on-current-well-state-to-potential-calculations

Browse Source 
Pass on current well state to (implicit) potential calculations
This commit is contained in:
Kai Bao
2025-02-03 11:09:31 +01:00
committed by GitHub
parent 2459bb74b2 d6fbbf64f8
commit 236d09c7db
4 changed files with 8 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
opm/simulators/wells/MultisegmentWell.hpp
View File

@@ -250,6 +250,7 @@ namespace Opm {
DeferredLogger& deferred_logger) const;
bool computeWellPotentialsImplicit(const Simulator& simulator,
const WellState<Scalar>& well_state,
std::vector<Scalar>& well_potentials,
DeferredLogger& deferred_logger) const;

5
opm/simulators/wells/MultisegmentWell_impl.hpp
View File

@@ -296,7 +296,7 @@ namespace Opm
debug_cost_counter_ = 0;
bool converged_implicit = false;
if (this->param_.local_well_solver_control_switching_) {
converged_implicit = computeWellPotentialsImplicit(simulator, well_potentials, deferred_logger);
converged_implicit = computeWellPotentialsImplicit(simulator, well_state, well_potentials, deferred_logger);
if (!converged_implicit) {
deferred_logger.debug("Implicit potential calculations failed for well "
+ this->name() + ", reverting to original aproach.");
@@ -516,6 +516,7 @@ namespace Opm
bool
MultisegmentWell<TypeTag>::
computeWellPotentialsImplicit(const Simulator& simulator,
const WellState<Scalar>& well_state,
std::vector<Scalar>& well_potentials,
DeferredLogger& deferred_logger) const
{
@@ -526,7 +527,7 @@ namespace Opm
well_copy.debug_cost_counter_ = 0;
// store a copy of the well state, we don't want to update the real well state
WellState<Scalar> well_state_copy = simulator.problem().wellModel().wellState();
WellState<Scalar> well_state_copy = well_state;
const auto& group_state = simulator.problem().wellModel().groupState();
auto& ws = well_state_copy.well(this->index_of_well_);

1
opm/simulators/wells/StandardWell.hpp
View File

@@ -325,6 +325,7 @@ namespace Opm
const WellState<Scalar>& well_state) const;
bool computeWellPotentialsImplicit(const Simulator& ebos_simulator,
const WellState<Scalar>& well_state,
std::vector<Scalar>& well_potentials,
DeferredLogger& deferred_logger) const;

5
opm/simulators/wells/StandardWell_impl.hpp
View File

@@ -1622,6 +1622,7 @@ namespace Opm
bool
StandardWell<TypeTag>::
computeWellPotentialsImplicit(const Simulator& simulator,
const WellState<Scalar>& well_state,
std::vector<Scalar>& well_potentials,
DeferredLogger& deferred_logger) const
{
@@ -1631,7 +1632,7 @@ namespace Opm
StandardWell<TypeTag> well_copy(*this);
// store a copy of the well state, we don't want to update the real well state
WellState<Scalar> well_state_copy = simulator.problem().wellModel().wellState();
WellState<Scalar> well_state_copy = well_state;
const auto& group_state = simulator.problem().wellModel().groupState();
auto& ws = well_state_copy.well(this->index_of_well_);
@@ -1762,7 +1763,7 @@ namespace Opm
// group controlled wells with defaulted guiderates will have zero targets as
// the potentials are used to compute the well fractions.
if (this->param_.local_well_solver_control_switching_ && !(this->changed_to_open_this_step_ && this->wellUnderZeroRateTarget(simulator, well_state, deferred_logger))) {
converged_implicit = computeWellPotentialsImplicit(simulator, well_potentials, deferred_logger);
converged_implicit = computeWellPotentialsImplicit(simulator, well_state, well_potentials, deferred_logger);
}
if (!converged_implicit) {
// does the well have a THP related constraint?