PowerLaw_decl.hpp

#ifndef POWERLAW_DECL_hpp
#define POWERLAW_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 PowerLaw : 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
      /*
       \brief Get the current viscosity value
       \return the scalar value of viscosity
       */
      double getViscosity() { return viscosity_; }

      PowerLaw(ParameterListPtr_Type parameters);
 
   private:
      double viscosity_;
      std::string shearThinningModel_; // for printing out which model is actually used
      double powerlaw_index_n;    // corresponds to n in the formulas being the power-law index
      double powerlaw_constant_K; // is the zero shear-rate viscosity
      double nu_infty;            // is the artificial infinite shear-rate viscosity needed for making bounds for the viscosity
      double nu_zero;             // is the artificial zero shear-rate viscosity needed for making bounds for the viscosity
      double shear_rate_limitZero;
 
   };
 
}
#endif