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>
50 template <
class Traits>
59 typedef typename UpscalerBase<Traits>::GridType GridType;
79 const std::vector<double>& initial_saturation,
80 const double boundary_saturation,
81 const double pressure_drop,
82 const permtensor_t& upscaled_perm,
bool& success);
97 void setToCapillaryLimit(
double average_s, std::vector<double>& s)
const;
102 typedef typename Traits::template TransportSolver<GridInterface, typename Super::BCs>::Type TransportSolver;
105 template <
class FlowSol>
106 void computeInOutFlows(std::pair<double, double>& water_inout,
107 std::pair<double, double>& oil_inout,
108 const FlowSol& flow_solution,
109 const std::vector<double>& saturations)
const;
111 void initImpl(
const Opm::ParameterGroup& param)
override;
114 std::vector<double> last_saturation_state_;
118 bool print_inoutflows_;
119 int simulation_steps_;
120 double init_stepsize_;
121 double relperm_threshold_;
122 double maximum_mobility_contrast_;
123 double sat_change_year_;
125 double max_stepsize_;
128 TransportSolver transport_solver_;
129 GridAdapter grid_adapter_;
134#include "SteadyStateUpscalerImplicit_impl.hpp"
A class for doing steady state upscaling.
Definition SteadyStateUpscalerImplicit.hpp:52
void initImpl(const Opm::ParameterGroup ¶m) override
Override from superclass.
Definition SteadyStateUpscalerImplicit_impl.hpp:78
void initSatLimits(std::vector< double > &s) const
Ensure saturations are not outside table.
Definition SteadyStateUpscalerImplicit_impl.hpp:416
double lastSaturationUpscaled() const
Computes the upscaled saturation corresponding to the saturation field returned by lastSaturationStat...
Definition SteadyStateUpscalerImplicit_impl.hpp:400
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:391
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