36 #ifndef OPM_STEADYSTATEUPSCALERIMPLICIT_HEADER
37 #define OPM_STEADYSTATEUPSCALERIMPLICIT_HEADER
39 #include <opm/upscaling/UpscalerBase.hpp>
40 #include <opm/porsol/euler/EulerUpstream.hpp>
41 #include <opm/porsol/euler/ImplicitCapillarity.hpp>
42 #include <opm/grid/common/GridAdapter.hpp>
51 template <
class Traits>
60 typedef typename UpscalerBase<Traits>::GridType GridType;
80 const std::vector<double>& initial_saturation,
81 const double boundary_saturation,
82 const double pressure_drop,
83 const permtensor_t& upscaled_perm,
bool& success);
98 void setToCapillaryLimit(
double average_s, std::vector<double>& s)
const;
103 typedef typename Traits::template TransportSolver<GridInterface, typename Super::BCs>::Type TransportSolver;
106 template <
class FlowSol>
107 void computeInOutFlows(std::pair<double, double>& water_inout,
108 std::pair<double, double>& oil_inout,
109 const FlowSol& flow_solution,
110 const std::vector<double>& saturations)
const;
112 virtual void initImpl(
const Opm::ParameterGroup& param);
116 std::vector<double> last_saturation_state_;
120 bool print_inoutflows_;
121 int simulation_steps_;
122 double init_stepsize_;
123 double relperm_threshold_;
124 double maximum_mobility_contrast_;
125 double sat_change_year_;
127 double max_stepsize_;
130 TransportSolver transport_solver_;
131 GridAdapter grid_adapter_;
136 #include "SteadyStateUpscalerImplicit_impl.hpp"
A class for doing steady state upscaling.
Definition: SteadyStateUpscalerImplicit.hpp:53
void initSatLimits(std::vector< double > &s) const
Ensure saturations are not outside table.
Definition: SteadyStateUpscalerImplicit_impl.hpp:419
virtual void initImpl(const Opm::ParameterGroup ¶m)
Override from superclass.
Definition: SteadyStateUpscalerImplicit_impl.hpp:78
double lastSaturationUpscaled() const
Computes the upscaled saturation corresponding to the saturation field returned by lastSaturationStat...
Definition: SteadyStateUpscalerImplicit_impl.hpp:403
std::pair< permtensor_t, permtensor_t > upscaleSteadyState(const int flow_direction, const std::vector< double > &initial_saturation, const double boundary_saturation, const double pressure_drop, const permtensor_t &upscaled_perm, bool &success)
Does a steady-state upscaling.
Definition: SteadyStateUpscalerImplicit_impl.hpp:175
SteadyStateUpscalerImplicit()
Default constructor.
Definition: SteadyStateUpscalerImplicit_impl.hpp:56
const std::vector< double > & lastSaturationState() const
Accessor for the steady state saturation field.
Definition: SteadyStateUpscalerImplicit_impl.hpp:394
A base class for upscaling.
Definition: UpscalerBase.hpp:56
ResProp::MutablePermTensor permtensor_t
A type for the upscaled permeability.
Definition: UpscalerBase.hpp:66
Inverting small matrices.
Definition: ImplicitAssembly.hpp:43