PolymerPropsAd.hpp

Go to the documentation of this file.

/*
  Copyright 2014 SINTEF ICT, Applied Mathematics.
  Copyright 2014 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_POLYMERPROPSAD_HEADED_INLCUDED
#define OPM_POLYMERPROPSAD_HEADED_INLCUDED

#include <cmath>
#include <vector>
#include <opm/autodiff/AutoDiffBlock.hpp>
#include <opm/autodiff/AutoDiffHelpers.hpp>
#include <opm/polymer/PolymerProperties.hpp>

namespace Opm {
    
    class PolymerPropsAd 
    {
    public:
        double rockDensity() const;
        
        double deadPoreVol() const;

        double cMax() const; 

        const std::vector<double>& shearWaterVelocity() const;

        const std::vector<double>& shearViscosityReductionFactor() const;

        double plyshlogRefConc() const;

        bool hasPlyshlogRefSalinity() const;

        bool hasPlyshlogRefTemp() const;

        double plyshlogRefSalinity() const;

        double plyshlogRefTemp() const;

        double shrate() const;

        double viscMult(double c) const; // multipler interpolated from PLYVISC table

        typedef AutoDiffBlock<double> ADB;
        typedef ADB::V V;

        V viscMult(const V& c) const;

        PolymerPropsAd(const PolymerProperties& polymer_props);

        ~PolymerPropsAd();
        
        V 
        effectiveInvWaterVisc(const V& c,const double* visc) const;

        ADB 
        effectiveInvWaterVisc(const ADB& c,const double* visc) const;
        
        V 
        polymerWaterVelocityRatio(const V& c) const;

        ADB
        polymerWaterVelocityRatio(const ADB& c) const;

        V
        adsorption(const V& c, const V& cmax_cells) const;

        ADB
        adsorption(const ADB& c, const ADB& cmax_cells) const;

        V
        effectiveRelPerm(const V& c, const V& cmax_cells, const V& relperm) const;


        ADB
        effectiveRelPerm(const ADB& c, const ADB& cmax_cells, const ADB& krw) const;

        bool computeShearMultLog(std::vector<double>& water_vel, std::vector<double>& visc_mult, std::vector<double>& shear_mult) const;


    private:
        const PolymerProperties& polymer_props_;
    };
    
} //namespace Opm

#endif// OPM_POLYMERPROPSAD_HEADED_INLCUDED