36#ifndef OPENRS_IMPLICITCAPILLARITY_HEADER
37#define OPENRS_IMPLICITCAPILLARITY_HEADER
40#include <opm/porsol/euler/EulerUpstreamResidual.hpp>
41#include <opm/common/utility/parameters/ParameterGroup.hpp>
42#include <opm/common/utility/numeric/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"
Class for doing simple transport by explicit Euler upstream method for general grid.
Definition EulerUpstreamResidual.hpp:58
Class for doing simple transport by explicit Euler upstream method for general grid.
Definition ImplicitCapillarity.hpp:53
void initObj(const GridInterface &grid, const ReservoirProperties &resprop, const BoundaryConditions &boundary)
Definition ImplicitCapillarity_impl.hpp:111
ImplicitCapillarity()
Definition ImplicitCapillarity_impl.hpp:58
void init(const Opm::ParameterGroup ¶m, const GridInterface &grid, const ReservoirProperties &resprop, const BoundaryConditions &boundary)
void init(const Opm::ParameterGroup ¶m)
Definition ImplicitCapillarity_impl.hpp:89
ImplicitCapillarity(const GridInterface &grid, const ReservoirProperties &resprop, const BoundaryConditions &boundary)
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.
Solve mixed formulation of incompressible flow modelled by Darcy's law.
Definition IncompFlowSolverHybrid.hpp:365
Inverting small matrices.
Definition ImplicitAssembly.hpp:43