36 #ifndef OPENRS_IMPLICITCAPILLARITY_HEADER
37 #define OPENRS_IMPLICITCAPILLARITY_HEADER
40 #include <opm/porsol/euler/EulerUpstreamResidual.hpp>
41 #include <opm/core/utility/parameters/ParameterGroup.hpp>
42 #include <opm/core/utility/SparseVector.hpp>
43 #include <opm/porsol/mimetic/IncompFlowSolverHybrid.hpp>
44 #include <opm/porsol/mimetic/MimeticIPEvaluator.hpp>
50 template <
class GridInterface,
class ReservoirProperties,
class BoundaryConditions,
51 template <
class,
class>
class InnerProd = MimeticIPEvaluator>
67 const ReservoirProperties& resprop,
68 const BoundaryConditions& boundary);
72 void init(
const Opm::ParameterGroup& param);
76 void init(
const Opm::ParameterGroup& param,
77 const GridInterface& grid,
78 const ReservoirProperties& resprop,
79 const BoundaryConditions& boundary);
83 void initObj(
const GridInterface& grid,
84 const ReservoirProperties& resprop,
85 const BoundaryConditions& boundary);
94 template <
class PressureSolution>
97 const typename GridInterface::Vector& gravity,
98 const PressureSolution& pressure_sol,
99 const Opm::SparseVector<double>& injection_rates)
const;
102 typedef typename GridInterface::CellIterator CIt;
103 typedef typename CIt::FaceIterator FIt;
104 typedef typename FIt::Vector Vector;
108 void checkAndPossiblyClampSat(std::vector<double>& s)
const;
112 BoundaryConditions> residual_;
114 bool method_viscous_;
115 bool method_gravity_;
118 double residual_tolerance_;
119 int linsolver_verbosity_;
121 double update_relaxation_;
128 #include "ImplicitCapillarity_impl.hpp"
132 #endif // OPENRS_IMPLICITCAPILLARITY_HEADER
Class for doing simple transport by explicit Euler upstream method for general grid.
Definition: ImplicitCapillarity.hpp:52
void init(const Opm::ParameterGroup ¶m)
Definition: ImplicitCapillarity_impl.hpp:89
Class for doing simple transport by explicit Euler upstream method for general grid.
Definition: EulerUpstreamResidual.hpp:59
Solve mixed formulation of incompressible flow modelled by Darcy's law ] The solver is based on a ...
Definition: IncompFlowSolverHybrid.hpp:369
void transportSolve(std::vector< double > &saturation, const double time, const typename GridInterface::Vector &gravity, const PressureSolution &pressure_sol, const Opm::SparseVector< double > &injection_rates) const
Solve transport equation, evolving.
void initObj(const GridInterface &grid, const ReservoirProperties &resprop, const BoundaryConditions &boundary)
Definition: ImplicitCapillarity_impl.hpp:111
ImplicitCapillarity()
Definition: ImplicitCapillarity_impl.hpp:58