Loading...
Searching...
No Matches
VFPProperties.hpp
/*
Copyright 2015 SINTEF ICT, Applied Mathematics.
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_AUTODIFF_VFPPROPERTIES_HPP_
#define OPM_AUTODIFF_VFPPROPERTIES_HPP_
#include <opm/simulators/wells/VFPInjProperties.hpp>
#include <opm/simulators/wells/VFPProdProperties.hpp>
#include <opm/simulators/wells/WellState.hpp>
#include <opm/simulators/wells/VFPHelpers.hpp>
#include <cstddef>
namespace Opm {
class VFPInjTable;
class VFPProdTable;
template<typename Scalar, typename IndexTraits>
class VFPProperties {
public:
VFPProperties(const std::vector<std::reference_wrapper<const VFPInjTable>>& inj_tables,
const std::vector<std::reference_wrapper<const VFPProdTable>>& prod_tables,
const WellState<Scalar, IndexTraits>& well_state)
: well_state_(well_state)
{
for (const auto& vfpinj : inj_tables)
this->m_inj.addTable( vfpinj );
for (const auto& vfpprod : prod_tables)
this->m_prod.addTable( vfpprod );
};
const VFPInjProperties<Scalar>* getInj() const
{
return &m_inj;
}
const VFPProdProperties<Scalar>* getProd() const
{
return &m_prod;
}
Scalar getExplicitWFR(const int table_id, const std::size_t well_index) const
{
const auto& pu = well_state_.phaseUsageInfo();
const bool has_water = pu.phaseIsActive(IndexTraits::waterPhaseIdx);
const bool has_oil = pu.phaseIsActive(IndexTraits::oilPhaseIdx);
const bool has_gas = pu.phaseIsActive(IndexTraits::gasPhaseIdx);
const auto& rates = well_state_.well(well_index).prev_surface_rates;
const auto& aqua = has_water ? rates[pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx)] : 0.0;
const auto& liquid = has_oil ? rates[pu.canonicalToActivePhaseIdx(IndexTraits::oilPhaseIdx)] : 0.0;
const auto& vapour = has_gas ? rates[pu.canonicalToActivePhaseIdx(IndexTraits::gasPhaseIdx)] : 0.0;
const VFPProdTable& table = this->m_prod.getTable(table_id);
return detail::getWFR(table, aqua, liquid, vapour);
}
Scalar getExplicitGFR(const int table_id, const std::size_t well_index) const
{
const auto& pu = well_state_.phaseUsageInfo();
const auto& rates = well_state_.well(well_index).prev_surface_rates;
const bool has_water = pu.phaseIsActive(IndexTraits::waterPhaseIdx);
const bool has_oil = pu.phaseIsActive(IndexTraits::oilPhaseIdx);
const bool has_gas = pu.phaseIsActive(IndexTraits::gasPhaseIdx);
const auto& aqua = has_water ? rates[pu.canonicalToActivePhaseIdx(IndexTraits::waterPhaseIdx)] : 0.0;
const auto& liquid = has_oil ? rates[pu.canonicalToActivePhaseIdx(IndexTraits::oilPhaseIdx)] : 0.0;
const auto& vapour = has_gas ? rates[pu.canonicalToActivePhaseIdx(IndexTraits::gasPhaseIdx)] : 0.0;
const VFPProdTable& table = this->m_prod.getTable(table_id);
return detail::getGFR(table, aqua, liquid, vapour);
}
private:
VFPInjProperties<Scalar> m_inj;
VFPProdProperties<Scalar> m_prod;
const WellState<Scalar, IndexTraits>& well_state_;
};
} // namespace Opm
#endif /* OPM_AUTODIFF_VFPPROPERTIES_HPP_ */
Opm::VFPInjProperties
Definition VFPInjProperties.hpp:34
Opm::VFPProdProperties
Class which linearly interpolates BHP as a function of rate, tubing head pressure,...
Definition VFPProdProperties.hpp:38
Opm::VFPProperties::getInj
const VFPInjProperties< Scalar > * getInj() const
Returns the VFP properties for injection wells.
Definition VFPProperties.hpp:65
Opm::VFPProperties::VFPProperties
VFPProperties(const std::vector< std::reference_wrapper< const VFPInjTable > > &inj_tables, const std::vector< std::reference_wrapper< const VFPProdTable > > &prod_tables, const WellState< Scalar, IndexTraits > &well_state)
Constructor Takes no ownership of data.
Definition VFPProperties.hpp:50
Opm::VFPProperties::getProd
const VFPProdProperties< Scalar > * getProd() const
Returns the VFP properties for production wells.
Definition VFPProperties.hpp:73
Opm::WellState
The state of a set of wells, tailored for use by the fully implicit blackoil simulator.
Definition WellState.hpp:66
Opm
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition blackoilbioeffectsmodules.hh:43