36 #ifndef OPENRS_EULERUPSTREAM_HEADER 37 #define OPENRS_EULERUPSTREAM_HEADER 39 #include <opm/porsol/euler/EulerUpstreamResidual.hpp> 41 #include <boost/unordered_map.hpp> 43 #include <opm/core/utility/parameters/ParameterGroup.hpp> 44 #include <opm/core/utility/SparseVector.hpp> 51 template <
class Gr
idInterface,
class ReservoirProperties,
class BoundaryConditions>
62 const ReservoirProperties& resprop,
63 const BoundaryConditions& boundary);
67 void init(
const Opm::ParameterGroup& param);
71 void init(
const Opm::ParameterGroup& param,
72 const GridInterface& grid,
73 const ReservoirProperties& resprop,
74 const BoundaryConditions& boundary);
78 void initObj(
const GridInterface& grid,
79 const ReservoirProperties& resprop,
80 const BoundaryConditions& boundary);
97 template <
class PressureSolution>
100 const typename GridInterface::Vector& gravity,
101 const PressureSolution& pressure_sol,
102 const Opm::SparseVector<double>& injection_rates)
const;
105 typedef typename GridInterface::CellIterator CIt;
106 typedef typename CIt::FaceIterator FIt;
107 typedef typename FIt::Vector Vector;
110 template <
class PressureSolution>
111 double computeCflTime(
const std::vector<double>& saturation,
113 const typename GridInterface::Vector& gravity,
114 const PressureSolution& pressure_sol)
const;
116 template <
class PressureSolution>
117 void smallTimeStep(std::vector<double>& saturation,
119 const typename GridInterface::Vector& gravity,
120 const PressureSolution& pressure_sol,
121 const Opm::SparseVector<double>& injection_rates)
const;
123 void checkAndPossiblyClampSat(std::vector<double>& s)
const;
128 BoundaryConditions> residual_computer_;
130 bool method_viscous_;
131 bool method_gravity_;
132 bool method_capillary_;
133 bool use_cfl_viscous_;
134 bool use_cfl_gravity_;
135 bool use_cfl_capillary_;
137 double courant_number_;
138 int minimum_small_steps_;
139 int maximum_small_steps_;
142 std::vector<double> porevol_;
145 mutable std::vector<double> residual_;
150 #include "EulerUpstream_impl.hpp" 152 #endif // OPENRS_EULERUPSTREAM_HEADER void setCourantNumber(double cn)
Set the Courant number.
Definition: EulerUpstream_impl.hpp:144
Class for immiscible dead oil and dry gas.
Definition: applier.hpp:18
Class for doing simple transport by explicit Euler upstream method for general grid.
Definition: EulerUpstream.hpp:52
void initObj(const GridInterface &grid, const ReservoirProperties &resprop, const BoundaryConditions &boundary)
Definition: EulerUpstream_impl.hpp:120
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 display()
Definition: EulerUpstream_impl.hpp:132
EulerUpstream()
Definition: EulerUpstream_impl.hpp:60
Class for doing simple transport by explicit Euler upstream method for general grid.
Definition: EulerUpstreamResidual.hpp:59
void init(const Opm::ParameterGroup ¶m)
Definition: EulerUpstream_impl.hpp:95