opm-simulators/opm/simulators/wells/SimFIBODetails.hpp

84 lines
2.5 KiB
C++
Raw Normal View History

/*
Copyright 2013 SINTEF ICT, Applied Mathematics.
Copyright 2014-2016 IRIS AS
Copyright 2015 Andreas Lauser
This file is part of the Open Porous Media project (OPM).
OPM is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OPM is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with OPM. If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef OPM_SIM_FIBO_DETAILS_HPP
#define OPM_SIM_FIBO_DETAILS_HPP
#include <utility>
#include <algorithm>
#include <locale>
2019-08-27 17:30:12 -05:00
#include <opm/parser/eclipse/EclipseState/Schedule/Events.hpp>
#include <opm/parser/eclipse/EclipseState/Schedule/Well/Well.hpp>
2019-08-27 17:30:12 -05:00
namespace Opm
{
namespace SimFIBODetails {
typedef std::unordered_map<std::string, Well > WellMap;
inline WellMap
mapWells(const std::vector< Well >& wells)
{
WellMap wmap;
for (const auto& w : wells)
{
wmap.insert(std::make_pair(w.name(), w));
}
return wmap;
}
inline void
historyRates(const PhaseUsage& pu,
const Well::ProductionControls& p,
std::vector<double>& rates)
{
assert (! p.prediction_mode);
assert (rates.size() ==
std::vector<double>::size_type(pu.num_phases));
if (pu.phase_used[ BlackoilPhases::Aqua ]) {
const std::vector<double>::size_type
i = pu.phase_pos[ BlackoilPhases::Aqua ];
rates[i] = p.water_rate;
}
if (pu.phase_used[ BlackoilPhases::Liquid ]) {
const std::vector<double>::size_type
i = pu.phase_pos[ BlackoilPhases::Liquid ];
rates[i] = p.oil_rate;
}
if (pu.phase_used[ BlackoilPhases::Vapour ]) {
const std::vector<double>::size_type
i = pu.phase_pos[ BlackoilPhases::Vapour ];
rates[i] = p.gas_rate;
}
}
} // namespace SimFIBODetails
} // namespace Opm
#endif // OPM_SIM_FIBO_DETAILS_HPP