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

Change non-convergence monitor to checking if number of non-converged has increased

Browse Source
This commit is contained in:
jakobtorben 2024-09-10 15:08:41 +02:00
parent cda47a6387
commit eef0ba50bf
4 changed files with 19 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
opm/simulators/flow/BlackoilModel.hpp
View File

@ -1064,7 +1064,7 @@ namespace Opm {
report.setReservoirFailed({types[ii], CR::Severity::Normal, compIdx});
}
report.setReservoirConvergenceMetric(types[ii], compIdx, res[ii]);
report.setReservoirConvergenceMetric(types[ii], compIdx, res[ii], tol[ii]);
}
}
@ -1113,22 +1113,23 @@ namespace Opm {
void checkCardPenalty(ConvergenceReport& report, int iteration)
{
if (iteration > 0) {
const auto& prev_report = convergence_reports_.back().report.back();
const auto& prev_metrics = prev_report.reservoirConvergence();
const auto& current_metrics = report.reservoirConvergence();
auto distances = std::vector<double>(current_metrics.size(), 0.0);
int current_above_tolerance = 0;
for (size_t i = 0; i < current_metrics.size(); ++i) {
distances[i] = std::max(std::log(current_metrics[i].value()), 0.0);
if (current_metrics[i].type() == prev_metrics[i].type() && current_metrics[i].phase() == prev_metrics[i].phase()) {
if (current_metrics[i].value() > prev_metrics[i].value()) {
// Metric value has increased, add a penalty card
report.addNonConvergedPenalty();
// Count number of metrics above tolerance
if (current_metrics[i].value() > current_metrics[i].tolerance()) {
current_above_tolerance++;
}
}
// Add penalty if number of metrics above tolerance has increased
if (current_above_tolerance > prev_above_tolerance_) {
report.addNonConvergedPenalty();
}
// use L1 norm of the distances vector
double current_distance = std::accumulate(distances.begin(), distances.end(), 0.0);
@ -1136,6 +1137,7 @@ namespace Opm {
report.addDistanceDecayPenalty();
}
prev_distance_ = current_distance;
prev_above_tolerance_ = current_above_tolerance;
}
if (report.wellFailures().size() > 0) {
@ -1445,7 +1447,7 @@ namespace Opm {
double linear_solve_setup_time_;
ConvergenceReport::PenaltyCard total_penaltyCard_;
double prev_distance_ = std::numeric_limits<double>::infinity();
int prev_above_tolerance_ = 0;
public:
std::vector<bool> wasSwitched_;
};

2
opm/simulators/flow/BlackoilModelNldd.hpp
View File

@ -678,7 +678,7 @@ private:
report.setReservoirFailed({types[ii], CR::Severity::Normal, compIdx});
}
report.setReservoirConvergenceMetric(types[ii], compIdx, res[ii]);
report.setReservoirConvergenceMetric(types[ii], compIdx, res[ii], tol[ii]);
}
}

2
opm/simulators/flow/BlackoilModelParameters.hpp
View File

@ -138,7 +138,7 @@ struct LocalDomainsPartitioningMethod { static constexpr auto value = "zoltan";
struct LocalDomainsOrderingMeasure { static constexpr auto value = "maxpressure"; };
struct ConvergenceMonitoring { static constexpr bool value = false; };
struct ConvergenceMonitoringCutOff { static constexpr int value = 30; };
struct ConvergenceMonitoringCutOff { static constexpr int value = 6; };
template<class Scalar>
struct ConvergenceMonitoringDecayFactor { static constexpr Scalar value = 0.75; };

7
opm/simulators/timestepping/ConvergenceReport.hpp
View File

@ -129,13 +129,14 @@ namespace Opm
// use this for anything else.
ReservoirConvergenceMetric() = default;
ReservoirConvergenceMetric(ReservoirFailure::Type t, int phase, double value)
: type_(t), phase_(phase), value_(value)
ReservoirConvergenceMetric(ReservoirFailure::Type t, int phase, double value, double tolerance)
: type_(t), phase_(phase), value_(value), tolerance_(tolerance)
{}
ReservoirFailure::Type type() const { return type_; }
int phase() const { return phase_; }
double value() const { return value_; }
double tolerance() const { return tolerance_; }
template <typename Serializer>
void serializeOp(Serializer& serializer)
@ -143,6 +144,7 @@ namespace Opm
serializer(this->type_);
serializer(this->phase_);
serializer(this->value_);
serializer(this->tolerance_);
}
private:
@ -151,6 +153,7 @@ namespace Opm
ReservoirFailure::Type type_ { ReservoirFailure::Type::Invalid };
int phase_ { -1 };
double value_ { 0.0 };
double tolerance_ { 0.0 };
};
class WellFailure