RelpermDiagnostics.hpp

Go to the documentation of this file.

/*
  Copyright 2015 Statoil ASA.
 
  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_RELPERMDIAGNOSTICS_HEADER_INCLUDED
#define OPM_RELPERMDIAGNOSTICS_HEADER_INCLUDED
 
#include <vector>
#include <utility>
 
#if HAVE_CONFIG_H
#include "config.h"
#endif // HAVE_CONFIG_H
 
#include <opm/core/grid.h>
#include <opm/core/grid/GridManager.hpp>
#include <opm/core/grid/GridHelpers.hpp>
#include <opm/core/utility/linearInterpolation.hpp>
#include <opm/parser/eclipse/EclipseState/EclipseState.hpp>
#include <opm/parser/eclipse/Deck/Deck.hpp>
#include <opm/common/OpmLog/OpmLog.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/SsfnTable.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/MiscTable.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/MsfnTable.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/SgcwmisTable.hpp>
#include <opm/parser/eclipse/EclipseState/Tables/SorwmisTable.hpp>
#include <opm/material/fluidmatrixinteractions/EclEpsScalingPoints.hpp>
 
namespace Opm {
 
    class Sof2Table;
    class SgwfnTable;
 
    class RelpermDiagnostics 
    {
    public:
        template <class GridT>
        void diagnosis(const EclipseState& eclState,
                       const Deck& deck,
                       const GridT& grid);
 
    private:
        enum FluidSystem {
            OilWater,
            OilGas,
            WaterGas,
            BlackOil,
            Solvent
        };
        
        FluidSystem fluidSystem_;
 
        enum SaturationFunctionFamily {
            FamilyI,
            FamilyII,
            NoFamily
        };
  
        SaturationFunctionFamily satFamily_;
 
        std::vector<Opm::EclEpsScalingPointsInfo<double> > unscaledEpsInfo_;
        std::vector<Opm::EclEpsScalingPointsInfo<double> > scaledEpsInfo_;
 
 
        void phaseCheck_(const EclipseState& es);
 
        void satFamilyCheck_(const EclipseState& eclState);
 
        void tableCheck_(const EclipseState& eclState);
 
        void unscaledEndPointsCheck_(const Deck& deck,
                                     const EclipseState& eclState);
 
        template <class GridT>
        void scaledEndPointsCheck_(const Deck& deck,
                                   const EclipseState& eclState,
                                   const GridT& grid);
 
        void swofTableCheck_(const Opm::SwofTable& swofTables,
                             const int satnumIdx);
        void sgofTableCheck_(const Opm::SgofTable& sgofTables,
                             const int satnumIdx);
        void slgofTableCheck_(const Opm::SlgofTable& slgofTables,
                              const int satnumIdx);
        void swfnTableCheck_(const Opm::SwfnTable& swfnTables,
                             const int satnumIdx);
        void sgfnTableCheck_(const Opm::SgfnTable& sgfnTables,
                             const int satnumIdx);
        void sof3TableCheck_(const Opm::Sof3Table& sof3Tables,
                             const int satnumIdx);
        void sof2TableCheck_(const Opm::Sof2Table& sof2Tables,
                             const int satnumIdx);
        void sgwfnTableCheck_(const Opm::SgwfnTable& sgwfnTables,
                              const int satnumIdx);
        void sgcwmisTableCheck_(const Opm::SgcwmisTable& sgcwmisTables,
                                const int satnumIdx);
        void sorwmisTableCheck_(const Opm::SorwmisTable& sorwmisTables,
                                const int satnumIdx);
        void ssfnTableCheck_(const Opm::SsfnTable& ssfnTables,
                             const int satnumIdx);
        void miscTableCheck_(const Opm::MiscTable& miscTables,
                             const int miscnumIdx);
        void msfnTableCheck_(const Opm::MsfnTable& msfnTables,
                             const int satnumIdx);
    };
 
} //namespace Opm
 
#include <opm/core/props/satfunc/RelpermDiagnostics_impl.hpp>
 
#endif // OPM_RELPERMDIAGNOSTICS_HEADER_INCLUDED