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 #5557 from atgeirr/avoid-cprw-incompatibilities

Browse Source 
Intercept CPRW if MICP or polymer injectivity is used.
This commit is contained in:
Atgeirr Flø Rasmussen 2024-08-26 09:50:57 +02:00 committed by GitHub
commit 2e15a04db5
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
opm/simulators/linalg/ISTLSolver.hpp
View File

@ -35,6 +35,7 @@
#include <opm/models/utils/propertysystem.hh>
#include <opm/simulators/flow/BlackoilModelParameters.hpp>
#include <opm/simulators/flow/FlowBaseVanguard.hpp>
#include <opm/simulators/flow/FlowProblemProperties.hpp>
#include <opm/simulators/linalg/ExtractParallelGridInformationToISTL.hpp>
#include <opm/simulators/linalg/FlowLinearSolverParameters.hpp>
#include <opm/simulators/linalg/matrixblock.hh>
@ -159,6 +160,10 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
using ElementMapper = GetPropType<TypeTag, Properties::ElementMapper>;
constexpr static std::size_t pressureIndex = GetPropType<TypeTag, Properties::Indices>::pressureSwitchIdx;
enum { enableMICP = getPropValue<TypeTag, Properties::EnableMICP>() };
enum { enablePolymerMolarWeight = getPropValue<TypeTag, Properties::EnablePolymerMW>() };
constexpr static bool isIncompatibleWithCprw = enableMICP || enablePolymerMolarWeight;
#if HAVE_MPI
using CommunicationType = Dune::OwnerOverlapCopyCommunication<int,int>;
#else
@ -211,6 +216,21 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
{
OPM_TIMEBLOCK(IstlSolver);
if (isIncompatibleWithCprw) {
// Some model variants are incompatible with the CPRW linear solver.
if (parameters_[0].linsolver_ == "cprw" || parameters_[0].linsolver_ == "hybrid") {
std::string incompatible_model = "Unknown";
if (enableMICP) {
incompatible_model = "MICP";
} else if (enablePolymerMolarWeight) {
incompatible_model = "Polymer injectivity";
}
OPM_THROW(std::runtime_error,
incompatible_model + " model is incompatible with the CPRW linear solver.\n"
"Choose a different option, for example --linear-solver=ilu0");
}
}
if (parameters_[0].linsolver_ == "hybrid") {
// Experimental hybrid configuration.
// When chosen, will set up two solvers, one with CPRW