My Project
Loading...
Searching...
No Matches
SteadyStateUpscalerImplicit.hpp
1//===========================================================================
2//
3// File: SteadyStateUpscaler.hpp
4//
5// Created: Fri Aug 28 14:01:19 2009
6//
7// Author(s): Atgeirr F Rasmussen <atgeirr@sintef.no>
8// B�rd Skaflestad <bard.skaflestad@sintef.no>
9//
10// $Date$
11//
12// $Revision$
13//
14//===========================================================================
15
16/*
17 Copyright 2009, 2010 SINTEF ICT, Applied Mathematics.
18 Copyright 2009, 2010 Statoil ASA.
19
20 This file is part of The Open Reservoir Simulator Project (OpenRS).
21
22 OpenRS is free software: you can redistribute it and/or modify
23 it under the terms of the GNU General Public License as published by
24 the Free Software Foundation, either version 3 of the License, or
25 (at your option) any later version.
26
27 OpenRS is distributed in the hope that it will be useful,
28 but WITHOUT ANY WARRANTY; without even the implied warranty of
29 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
30 GNU General Public License for more details.
31
32 You should have received a copy of the GNU General Public License
33 along with OpenRS. If not, see <http://www.gnu.org/licenses/>.
34*/
35
36#ifndef OPM_STEADYSTATEUPSCALERIMPLICIT_HEADER
37#define OPM_STEADYSTATEUPSCALERIMPLICIT_HEADER
38
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>
43
44namespace Opm
45{
50 template <class Traits>
52 {
53 public:
54 // ------- Typedefs and enums -------
55
57 typedef typename Super::permtensor_t permtensor_t;
59 typedef typename UpscalerBase<Traits>::GridType GridType;
60 enum { Dimension = UpscalerBase<Traits>::Dimension };
61
62 // ------- Methods -------
63
66
78 std::pair<permtensor_t, permtensor_t> upscaleSteadyState(const int flow_direction,
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);
83
87 const std::vector<double>& lastSaturationState() const;
88
92 double lastSaturationUpscaled() const;
93
95 void initSatLimits(std::vector<double>& s) const;
96
97 void setToCapillaryLimit(double average_s, std::vector<double>& s) const;
98
99
100 protected:
101 // ------- Typedefs -------
102 typedef typename Traits::template TransportSolver<GridInterface, typename Super::BCs>::Type TransportSolver;
103
104 // ------- Methods -------
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;
112
113 // ------- Data members -------
114 std::vector<double> last_saturation_state_;
115 bool use_gravity_;
116 bool output_vtk_;
117 bool output_ecl_;
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_;
124 int max_it_;
125 double max_stepsize_;
126 double dt_sat_tol_;
127 bool use_maxdiff_;
128 TransportSolver transport_solver_;
129 GridAdapter grid_adapter_;
130 };
131
132} // namespace Opm
133
134#include "SteadyStateUpscalerImplicit_impl.hpp"
135
136
137#endif // OPM_STEADYSTATEUPSCALERIMPLICIT_HEADER
A class for doing steady state upscaling.
Definition SteadyStateUpscalerImplicit.hpp:52
void initImpl(const Opm::ParameterGroup &param) 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