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

clean-up relaxed tolerances

Browse Source
This commit is contained in:
Tor Harald Sandve
2024-02-15 09:22:42 +01:00
parent c546d7fc30
commit c73337186e
1 changed files with 26 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -909,25 +909,34 @@ namespace Opm {
auto cnvErrorPvFraction = computeCnvErrorPv(B_avg, dt); auto cnvErrorPvFraction = computeCnvErrorPv(B_avg, dt);
cnvErrorPvFraction /= (pvSum - numAquiferPvSum); cnvErrorPvFraction /= (pvSum - numAquiferPvSum);
// Default value of relaxed_max_pv_fraction_ is 0.03 and min_strict_cnv_iter_ is 0.
// For each iteration, we need to determine whether to use the relaxed CNV tolerance.
// To disable the usage of relaxed CNV tolerance, you can set the relaxed_max_pv_fraction_ to be 0.
bool use_relaxed_cnv = cnvErrorPvFraction < param_.relaxed_max_pv_fraction_ && iteration >= param_.min_strict_cnv_iter_;
bool use_relaxed_mb = param_.min_strict_mb_iter_ < 0? false : iteration >= param_.min_strict_mb_iter_;
if ((iteration + 1) == maxIter) { // For each iteration, we need to determine whether to use the relaxed tolerances.
use_relaxed_mb = true; // To disable the usage of relaxed tolerances, you can set the relaxed tolerances as the strict tolerances.
use_relaxed_cnv = true;
// If min_strict_mb_iter = -1 (default) we use a relaxed tolerance for the mass balance for the last iterations
// For positive values we use the relaxed tolerance after the given number of iterations
const bool relax_final_iteration_mb = (param_.min_strict_mb_iter_ < 0) && (iteration == maxIter);
const bool use_relaxed_mb = relax_final_iteration_mb || (param_.min_strict_mb_iter_ >= 0 && iteration >= param_.min_strict_mb_iter_);
// If min_strict_cnv_iter = -1 we use a relaxed tolerance for the cnv for the last iterations
// For positive values we use the relaxed tolerance after the given number of iterations
// We also use relaxed tolerances for cells with total poro volume less than relaxed_max_pv_fraction_
// Default value of relaxed_max_pv_fraction_ is 0.03
const bool relax_final_iteration_cnv = (param_.min_strict_cnv_iter_ < 0) && (iteration == maxIter);
const bool relax_iter_cnv = param_.min_strict_mb_iter_ >= 0 && iteration >= param_.min_strict_mb_iter_;
const bool relax_pv_fraction_cnv = cnvErrorPvFraction < param_.relaxed_max_pv_fraction_;
const bool use_relaxed_cnv = relax_final_iteration_cnv || relax_pv_fraction_cnv || relax_iter_cnv;
if (relax_final_iteration_mb || relax_final_iteration_cnv) {
if ( terminal_output_ ) { if ( terminal_output_ ) {
std::ostringstream ss; std::string message = "Number of newton iterations reached its maximum try to continue with relaxed tolerances:";
ss << "Number of newton iterations reached its maximum try to continue with relaxed tolerances:"; if (relax_final_iteration_mb)
const std::streamsize oprec = ss.precision(3); message += fmt::format(" MB: {:.1e}", param_.tolerance_mb_relaxed_);
const std::ios::fmtflags oflags = ss.setf(std::ios::scientific); if (relax_final_iteration_cnv)
ss << " CNV TOL: " << param_.tolerance_cnv_relaxed_; message += fmt::format(" CNV: {:.1e}", param_.tolerance_cnv_relaxed_);
ss << " MB TOL: " << param_.tolerance_mb_relaxed_;
ss.precision(oprec); OpmLog::debug(message);
ss.flags(oflags);
OpmLog::debug(ss.str());
} }
} }
const double tol_cnv = use_relaxed_cnv ? param_.tolerance_cnv_relaxed_ : param_.tolerance_cnv_; const double tol_cnv = use_relaxed_cnv ? param_.tolerance_cnv_relaxed_ : param_.tolerance_cnv_;