CarreauYasuda_decl.hpp

#ifndef CARREAUYASUDA_DECL_hpp
#define CARREAUYASUDA_DECL_hpp
 
#include "feddlib/core/General/DifferentiableFuncClass.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 CarreauYasuda : 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_; }
 
 

    CarreauYasuda(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
    double fluid_index_n;        // corresponds to n in the formulas being the exponent
    double nu_0;                 // is the zero shear-rate viscosity
    double nu_infty;             // is the infinite shear-rate viscosity
    double inflectionPoint;      // corresponds to a in the formulas in the literature
    double shear_rate_limitZero; // In the formulas the shear rate is in the denominator so we have to ensure that it is
  };
 
}
#endif