Dimless_Carreau_decl.hpp

#ifndef DIMLESS_CARREAU_DECL_hpp
#define DIMLESS_CARREAU_DECL_hpp
 
#include "feddlib/core/General/DifferentiableFuncClass.hpp"
// #include "feddlib/core/AceFemAssembly/Helper.hpp"
#include "feddlib/core/FEDDCore.hpp"
 
namespace FEDD
{

 
   template <class SC = default_sc, class LO = default_lo, class GO = default_go, class NO = default_no>
   class Dimless_Carreau : public DifferentiableFuncClass<SC, LO, GO, NO>
   {
   public:
      typedef MultiVector<SC, LO, GO, NO> MultiVector_Type;
      typedef Teuchos::RCP<MultiVector_Type> MultiVectorPtr_Type;
      typedef Teuchos::RCP<const MultiVector_Type> MultiVectorConstPtr_Type;
 
      typedef DifferentiableFuncClass<SC, LO, GO, NO> DifferentiableFuncClass_Type;
 
      // Inherited Function from base abstract class

      void setParams(ParameterListPtr_Type params) override;
      void evaluateMapping(ParameterListPtr_Type params, double shearRate, double &viscosity) override;
      void evaluateDerivative(ParameterListPtr_Type params, double shearRate, double &res) override;

      void echoInformationMapping() override;
 
      // New Added Functions
      double getViscosity() { return viscosity_; };

      Dimless_Carreau(ParameterListPtr_Type parameters);
 
   private:
      double viscosity_;
      std::string shearThinningModel_; // for printing out which model is actually used
      double characteristicTime; // corresponds to \lambda in the formulas in the literature here dimensionless
      double fluid_index_n;      // corresponds to n in the formulas being the power-law index
      double nu_0;               // is the zero shear-rate viscosity here dimensionless
      double nu_infty;           // is the infnite shear-rate viscosity here dimensionless
      double shear_rate_limitZero;
      double reference_viscosity; // to obtain actual viscosity we have to multiply with reference viscosity
   };
 
}
#endif