OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-25 10:40:21 -06:00
Code Issues Packages Projects Releases Wiki Activity

Change to maximum norm for the well residuals

Browse Source
This commit is contained in:
Tor Harald Sandve 2015-05-06 08:53:37 +02:00
parent b4369cade8
commit 43b8f9851e
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@ -1973,11 +1973,13 @@ namespace detail {
converged_MB = converged_MB && (mass_balance_residual[idx] < tol_mb); converged_MB = converged_MB && (mass_balance_residual[idx] < tol_mb);
converged_CNV = converged_CNV && (CNV[idx] < tol_cnv); converged_CNV = converged_CNV && (CNV[idx] < tol_cnv);
double maxNormWell = 0.0;
for ( int w=0; w<nw; ++w ) for ( int w=0; w<nw; ++w )
{ {
well_flux_residual[idx] += std::abs(residual_.well_flux_eq.value()[nw*idx + w]); maxNormWell = std::max(maxNormWell, std::abs(residual_.well_flux_eq.value()[nw*idx + w]));
} }
well_flux_residual[idx] *= (B_avg[idx] * dt/nw); well_flux_residual[idx] = B_avg[idx] * dt * maxNormWell;
converged_Well = converged_Well && (well_flux_residual[idx] < tol_wells); converged_Well = converged_Well && (well_flux_residual[idx] < tol_wells);
} }