StandardWellConnections.hpp

Go to the documentation of this file.

/*
  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 <http://www.gnu.org/licenses/>.
*/
 
 
#ifndef OPM_STANDARDWELL_CONNECTIONS_HEADER_INCLUDED
#define OPM_STANDARDWELL_CONNECTIONS_HEADER_INCLUDED
 
#include <opm/simulators/wells/StandardWellPrimaryVariables.hpp>
 
#include <functional>
#include <variant>
#include <vector>
 
namespace Opm
{
 
class DeferredLogger;
enum class Phase;
template<class FluidSystem, class Indices> class WellInterfaceIndices;
template<class Scalar> class WellState;
 
template<class FluidSystem, class Indices>
class StandardWellConnections
{
public:
    using Scalar = typename FluidSystem::Scalar;
    StandardWellConnections(const WellInterfaceIndices<FluidSystem,Indices>& well);
 
    struct Properties
    {
        std::vector<Scalar> b_perf;
        std::vector<Scalar> rsmax_perf;
        std::vector<Scalar> rvmax_perf;
        std::vector<Scalar> rvwmax_perf;
        std::vector<Scalar> rswmax_perf;
        std::vector<Scalar> surf_dens_perf;
    };
 
    void computePropertiesForPressures(const WellState<Scalar>& well_state,
                                       const std::function<Scalar(int,int)>& getTemperature,
                                       const std::function<Scalar(int)>& getSaltConcentration,
                                       const std::function<int(int)>& pvtRegionIdx,
                                       const std::function<Scalar(int)>& solventInverseFormationVolumeFactor,
                                       const std::function<Scalar(int)>& solventRefDensity,
                                       Properties& props) const;
 
    void computeProperties(const WellState<Scalar>& well_state,
                           const std::function<Scalar(int,int)>& invB,
                           const std::function<Scalar(int,int)>& mobility,
                           const std::function<Scalar(int)>& solventInverseFormationVolumeFactor,
                           const std::function<Scalar(int)>& solventMobility,
                           const Properties& props,
                           DeferredLogger& deferred_logger);
 
    Scalar rho() const
    {
        return this->rho(0);
    }
 
    Scalar rho(const typename std::vector<Scalar>::size_type i) const
    {
        return (i < this->perf_densities_.size())
            ? this->perf_densities_[i]
            : 0.0;
    }
 
    Scalar pressure_diff(const unsigned perf) const
    { return perf_pressure_diffs_[perf]; }
 
    using Eval = typename WellInterfaceIndices<FluidSystem,Indices>::Eval;
    using EvalWell = typename StandardWellPrimaryVariables<FluidSystem,Indices>::EvalWell;
 
    Eval connectionRateBrine(double& rate,
                             const double vap_wat_rate,
                             const std::vector<EvalWell>& cq_s,
                             const std::variant<Scalar,EvalWell>& saltConcentration) const;
 
    Eval connectionRateFoam(const std::vector<EvalWell>& cq_s,
                            const std::variant<Scalar,EvalWell>& foamConcentration,
                            const Phase transportPhase,
                            DeferredLogger& deferred_logger) const;
 
    std::tuple<Eval,EvalWell>
    connectionRatePolymer(double& rate,
                          const std::vector<EvalWell>& cq_s,
                          const std::variant<Scalar,EvalWell>& polymerConcentration) const;
 
    std::tuple<Eval,Eval,Eval>
    connectionRatesMICP(const std::vector<EvalWell>& cq_s,
                        const std::variant<Scalar,EvalWell>& microbialConcentration,
                        const std::variant<Scalar,EvalWell>& oxygenConcentration,
                        const std::variant<Scalar,EvalWell>& ureaConcentration) const;
 
    std::tuple<Eval,EvalWell>
    connectionRatezFraction(double& rate,
                            const double dis_gas_rate,
                            const std::vector<EvalWell>& cq_s,
                            const std::variant<Scalar, std::array<EvalWell,2>>& solventConcentration) const;
 
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<Scalar>& perfComponentRates,
                          const Properties& props,
                          DeferredLogger& deferred_logger);
 
    const WellInterfaceIndices<FluidSystem,Indices>& well_; 
 
    std::vector<Scalar> perf_densities_; 
    std::vector<Scalar> perf_pressure_diffs_; 
};
 
}
 
#endif // OPM_STANDARDWELL_CONNECTIONS_HEADER_INCLUDED