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