mirror of
https://github.com/OPM/opm-simulators.git
synced 2024-11-28 20:13:49 -06:00
72ae444566
The PreconditionerFactory requires 2.6, so this is addressed by simply not using that code at all for Dune < 2.6. No FlexibleSolver etc.
104 lines
4.5 KiB
C++
104 lines
4.5 KiB
C++
/*
|
|
Copyright 2013, 2014, 2015, 2019 SINTEF Digital, Mathematics and Cybernetics.
|
|
Copyright 2014 Dr. Blatt - HPC-Simulation-Software & Services
|
|
Copyright 2015, 2017 IRIS AS
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
#include "config.h"
|
|
#include "flow/flow_tag.hpp"
|
|
|
|
#if DUNE_VERSION_NEWER(DUNE_ISTL, 2, 6)
|
|
#include <opm/simulators/linalg/ISTLSolverEbosFlexible.hpp>
|
|
#else
|
|
#include <opm/simulators/linalg/ISTLSolverEbosCpr.hpp>
|
|
#endif
|
|
|
|
BEGIN_PROPERTIES
|
|
NEW_TYPE_TAG(EclFlowProblemSimple, INHERITS_FROM(EclFlowProblem));
|
|
NEW_PROP_TAG(FluidState);
|
|
//SET_TYPE_PROP(EclBaseProblem, Problem, Ewoms::EclProblem<TypeTag>);
|
|
SET_PROP(EclFlowProblemSimple, FluidState)
|
|
{
|
|
private:
|
|
typedef typename GET_PROP_TYPE(TypeTag, FluidSystem) FluidSystem;
|
|
typedef typename GET_PROP_TYPE(TypeTag, Indices) Indices;
|
|
enum { enableTemperature = GET_PROP_VALUE(TypeTag, EnableTemperature) };
|
|
enum { enableSolvent = GET_PROP_VALUE(TypeTag, EnableSolvent) };
|
|
enum { enableEnergy = GET_PROP_VALUE(TypeTag, EnableEnergy) };
|
|
enum { numPhases = GET_PROP_VALUE(TypeTag, NumPhases) };
|
|
typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
|
|
typedef typename GET_PROP_TYPE(TypeTag, Evaluation) Evaluation;
|
|
static const bool compositionSwitchEnabled = Indices::gasEnabled;
|
|
|
|
public:
|
|
//typedef Opm::BlackOilFluidSystemSimple<Scalar> type;
|
|
typedef Opm::BlackOilFluidState<Evaluation, FluidSystem, enableTemperature, enableEnergy, compositionSwitchEnabled, Indices::numPhases > type;
|
|
};
|
|
|
|
SET_BOOL_PROP(EclFlowProblemSimple, MatrixAddWellContributions, true);
|
|
SET_INT_PROP(EclFlowProblemSimple, LinearSolverVerbosity,0);
|
|
SET_SCALAR_PROP(EclFlowProblemSimple, LinearSolverReduction, 1e-2);
|
|
SET_INT_PROP(EclFlowProblemSimple, LinearSolverMaxIter, 100);
|
|
SET_BOOL_PROP(EclFlowProblemSimple, UseAmg, true);//probably not used
|
|
SET_BOOL_PROP(EclFlowProblemSimple, UseCpr, true);
|
|
SET_INT_PROP(EclFlowProblemSimple, CprMaxEllIter, 1);
|
|
SET_INT_PROP(EclFlowProblemSimple, CprEllSolvetype, 3);
|
|
SET_INT_PROP(EclFlowProblemSimple, CprReuseSetup, 3);
|
|
SET_INT_PROP(EclFlowProblemSimple, CprSolverVerbose, 0);
|
|
SET_STRING_PROP(EclFlowProblemSimple, SystemStrategy, "quasiimpes");
|
|
END_PROPERTIES
|
|
|
|
namespace Ewoms {
|
|
namespace Properties {
|
|
|
|
SET_PROP(EclFlowProblemSimple, FluidSystem)
|
|
{
|
|
private:
|
|
//typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
|
|
typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
|
|
typedef typename GET_PROP_TYPE(TypeTag, Evaluation) Evaluation;
|
|
|
|
public:
|
|
typedef Opm::BlackOilFluidSystem<Scalar> type;
|
|
};
|
|
//NEW_TYPE_TAG(EclFlowProblem, INHERITS_FROM(BlackOilModel, EclBaseProblem));
|
|
SET_TYPE_PROP(EclFlowProblemSimple, IntensiveQuantities, Ewoms::BlackOilIntensiveQuantities<TypeTag>);
|
|
//SET_TYPE_PROP(EclFlowProblemSimple, LinearSolverBackend, Opm::ISTLSolverEbos<TypeTag>);
|
|
//SET_TAG_PROP(EclFlowProblemSimple, LinearSolverSplice, ParallelBiCGStabLinearSolver);
|
|
//SET_TYPE_PROP(EclFlowProblemSimple, LinearSolverBackend, Ewoms::Linear::ParallelBiCGStabSolverBackend<TypeTag>);//not work
|
|
//SET_TYPE_PROP(EclFlowProblemSimple, LinearSolverBackend, Ewoms::Linear::SuperLUBackend<TypeTag>)//not work
|
|
//SET_TAG_PROP(EclFlowProblem, FluidState, Opm::BlackOilFluidState);
|
|
#if DUNE_VERSION_NEWER(DUNE_ISTL, 2, 6)
|
|
SET_TYPE_PROP(EclFlowProblemSimple, LinearSolverBackend, Opm::ISTLSolverEbosFlexible<TypeTag>);
|
|
#else
|
|
SET_TYPE_PROP(EclFlowProblemSimple, LinearSolverBackend, Opm::ISTLSolverEbosCpr<TypeTag>);
|
|
#endif
|
|
SET_BOOL_PROP(EclFlowProblemSimple, EnableStorageCache, true);
|
|
SET_BOOL_PROP(EclFlowProblemSimple, EnableIntensiveQuantityCache, true);
|
|
|
|
//SET_INT_PROP(EclFlowProblemSimple, NumWellAdjoint, 1);
|
|
//SET_BOOL_PROP(EclFlowProblem, EnableStorageCache, true);
|
|
//SET_BOOL_PROP(EclFlowProblem, EnableIntensiveQuantityCache, true);
|
|
}
|
|
}
|
|
|
|
int main(int argc, char** argv)
|
|
{
|
|
typedef TTAG(EclFlowProblemSimple) TypeTag;
|
|
return mainFlow<TypeTag>(argc, argv);
|
|
}
|