MiscibilityDead.hpp

Go to the documentation of this file.

//===========================================================================
//                                                                           
// File: MiscibilityDead.hpp                                                  
//                                                                           
// Created: Wed Feb 10 09:05:47 2010                                         
//                                                                           
// Author: Bjørn Spjelkavik <bsp@sintef.no>
//                                                                           
// Revision: $Id$
//                                                                           
//===========================================================================
/*
  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 SINTEF_MISCIBILITYDEAD_HEADER
#define SINTEF_MISCIBILITYDEAD_HEADER

#include "MiscibilityProps.hpp"
#include <opm/core/utility/UniformTableLinear.hpp>
#include <opm/parser/eclipse/EclipseState/EclipseState.hpp>

namespace Opm
{
    class MiscibilityDead : public MiscibilityProps
    {
    public:

        MiscibilityDead(const PvdoTable& pvdoTable);
        MiscibilityDead(const PvdgTable& pvdgTable);
    virtual ~MiscibilityDead();

        virtual double getViscosity(int region, double press, const surfvol_t& surfvol) const;
        virtual double B(int region, double press, const surfvol_t& surfvol) const;
    virtual double dBdp(int region, double press, const surfvol_t& surfvol) const;
    virtual double R(int region, double press, const surfvol_t& surfvol) const;
    virtual double dRdp(int region, double press, const surfvol_t& surfvol) const;

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

    private:
    // PVT properties of dry gas or dead oil
        Opm::utils::UniformTableLinear<double> one_over_B_;
        Opm::utils::UniformTableLinear<double> viscosity_;
    };

}

#endif // SINTEF_MISCIBILITYDEAD_HEADER