/*
Copyright 2017 SINTEF Digital, Mathematics and Cybernetics.
Copyright 2017 Statoil ASA.
Copyright 2016 - 2017 IRIS AS.
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 .
*/
#ifndef OPM_STANDARDWELL_CONNECTIONS_HEADER_INCLUDED
#define OPM_STANDARDWELL_CONNECTIONS_HEADER_INCLUDED
#include
#include
namespace Opm
{
class DeferredLogger;
template class WellInterfaceIndices;
class WellState;
template
class StandardWellConnections
{
public:
StandardWellConnections(const WellInterfaceIndices& well);
void computePropertiesForPressures(const WellState& well_state,
const std::function& getTemperature,
const std::function& getSaltConcentration,
const std::function& pvtRegionIdx,
const std::function& solventInverseFormationVolumeFactor,
const std::function& solventRefDensity,
std::vector& b_perf,
std::vector& rsmax_perf,
std::vector& rvmax_perf,
std::vector& rvwmax_perf,
std::vector& rswmax_perf,
std::vector& surf_dens_perf) const;
//! \brief Compute connection properties (densities, pressure drop, ...)
void computeProperties(const WellState& well_state,
const std::function& invB,
const std::function& mobility,
const std::function& solventInverseFormationVolumeFactor,
const std::function& solventMobility,
const std::vector& b_perf,
const std::vector& rsmax_perf,
const std::vector& rvmax_perf,
const std::vector& rvwmax_perf,
const std::vector& rswmax_perf,
const std::vector& surf_dens_perf,
DeferredLogger& deferred_logger);
//! \brief Returns density for first perforation.
Scalar rho() const
{
return this->perf_densities_.empty() ? 0.0 : perf_densities_[0];
}
//! \brief Returns pressure drop for a given perforation.
Scalar pressure_diff(const unsigned perf) const
{ return perf_pressure_diffs_[perf]; }
private:
void computePressureDelta();
// TODO: not total sure whether it is a good idea to put this function here
// the major reason to put here is to avoid the usage of Wells struct
void computeDensities(const std::vector& perfComponentRates,
const std::vector& b_perf,
const std::vector& rsmax_perf,
const std::vector& rvmax_perf,
const std::vector& rvwmax_perf,
const std::vector& rswmax_perf,
const std::vector& surf_dens_perf,
DeferredLogger& deferred_logger);
const WellInterfaceIndices& well_; //!< Reference to well interface
std::vector perf_densities_; //!< densities of the fluid in each perforation
std::vector perf_pressure_diffs_; //!< // pressure drop between different perforations
};
}
#endif // OPM_STANDARDWELL_CONNECTIONS_HEADER_INCLUDED