BlackoilPVT.hpp

Go to the documentation of this file.

/*
  Copyright 2010 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_BLACKOILPVT_HEADER_INCLUDED
#define OPM_BLACKOILPVT_HEADER_INCLUDED


#include "MiscibilityProps.hpp"
#include "BlackoilDefs.hpp"
#include <opm/parser/eclipse/Deck/Deck.hpp>
#include <boost/scoped_ptr.hpp>
#include <string>


namespace Opm
{
    class BlackoilPVT : public BlackoilDefs
    {
    public:
        void init(Opm::DeckConstPtr deck);

        double getViscosity(double press,
                            const CompVec& surfvol,
                PhaseIndex phase) const;
    CompVec surfaceDensities() const;
        double B   (double press,
                    const CompVec& surfvol,
            PhaseIndex phase) const;
        double dBdp(double press,
                    const CompVec& surfvol,
            PhaseIndex phase) const;
        double R   (double press,
                    const CompVec& surfvol,
            PhaseIndex phase) const;
        double dRdp(double press,
                    const CompVec& surfvol,
            PhaseIndex phase) const;

        void getViscosity(const std::vector<PhaseVec>& pressures,
                          const std::vector<CompVec>& surfvol,
                          std::vector<PhaseVec>& output) const;
        void B(const std::vector<PhaseVec>& pressures,
               const std::vector<CompVec>& surfvol,
               std::vector<PhaseVec>& output) const;
        void dBdp(const std::vector<PhaseVec>& pressures,
                  const std::vector<CompVec>& surfvol,
                  std::vector<PhaseVec>& output_B,
                  std::vector<PhaseVec>& output_dBdp) const;
        void R(const std::vector<PhaseVec>& pressures,
               const std::vector<CompVec>& surfvol,
               std::vector<PhaseVec>& output) const;
        void dRdp(const std::vector<PhaseVec>& pressures,
                  const std::vector<CompVec>& surfvol,
                  std::vector<PhaseVec>& output_R,
                  std::vector<PhaseVec>& output_dRdp) const;

    private:
    int region_number_;
        const MiscibilityProps& propsForPhase(PhaseIndex phase) const;

    boost::scoped_ptr<MiscibilityProps> water_props_;
    boost::scoped_ptr<MiscibilityProps> oil_props_;
    boost::scoped_ptr<MiscibilityProps> gas_props_;
    CompVec densities_;
        mutable std::vector<double> data1_;
        mutable std::vector<double> data2_;
    };

}


#endif // OPM_BLACKOILPVT_HEADER_INCLUDED