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