Adding PLYSHLOG related to PolymerPropsAd

opm/polymer/fullyimplicit/PolymerPropsAd.cpp

@@ -60,6 +60,64 @@ namespace Opm {
 		return polymer_props_.cMax();
 	}
 
+    const std::vector<double>&
+    PolymerPropsAd::shearWaterVelocity() const
+    {
+        return polymer_props_.shearWaterVelocity();
+    }
+
+    const std::vector<double>&
+    PolymerPropsAd::shearViscosityReductionFactor() const
+    {
+        return polymer_props_.shearViscosityReductionFactor();
+    }
+
+    double
+    PolymerPropsAd::plyshlogRefConc() const
+    {
+        return polymer_props_.plyshlogRefConc();
+    }
+
+    bool
+    PolymerPropsAd::hasPlyshlogRefSalinity() const
+    {
+        return polymer_props_.hasPlyshlogRefSalinity();
+    }
+
+    bool
+    PolymerPropsAd::hasPlyshlogRefTemp() const
+    {
+        return polymer_props_.hasPlyshlogRefTemp();
+    }
+
+    double
+    PolymerPropsAd::plyshlogRefSalinity() const
+    {
+        return polymer_props_.plyshlogRefSalinity();
+    }
+
+    double
+    PolymerPropsAd::plyshlogRefTemp() const
+    {
+        return polymer_props_.plyshlogRefTemp();
+    }
+
+    double
+    PolymerPropsAd::viscMult(double c) const
+    {
+        return polymer_props_.viscMult(c);
+    }
+
+    V
+    PolymerPropsAd::viscMult(const V& c) const
+    {
+        int nc = c.size();
+        V visc_mult(nc);
+        for (int i = 0; i < nc; ++i) {
+            visc_mult[i] = polymer_props_.viscMult(c[i]);
+        }
+        return visc_mult;
+    }
 
 
 	

opm/polymer/fullyimplicit/PolymerPropsAd.hpp

@@ -41,9 +41,36 @@ namespace Opm {
 		/// \return 	The max concentration injected.
        	double cMax() const; 
 
+        /// \ return    The water velcoity or shear rate in the PLYSHLOG table
+        const std::vector<double>& shearWaterVelocity() const;
+
+        /// \ return    The viscosity reducation factor in the PLYSHLOG table
+        const std::vector<double>& shearViscosityReductionFactor() const;
+
+        /// \ return    The reference polymer concentration for PLYSHLOG table
+        double plyshlogRefConc() const;
+
+        /// \ return    The flag indicating if reference salinity is specified in PLYSHLOG keyword
+        bool hasPlyshlogRefSalinity() const;
+
+        /// \ return    The flag indicating if reference temperature is specified in PLYSHLOG keyword
+        bool hasPlyshlogRefTemp() const;
+
+        /// \ return    The reference salinity in PLYSHLOG keyword
+        double plyshlogRefSalinity() const;
+
+        /// \ return    The reference temperature in PLYSHLOG keyword
+        double plyshlogRefTemp() 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;
+		/// \param[in] c		Array of n polymer concentraion values.
+		/// \return 			Array of n viscosity multiplier from PLVISC table.
+
 		/// Constructor wrapping a polymer props.	
         PolymerPropsAd(const PolymerProperties& polymer_props);