AssembleFEGeneralizedNewtonian_decl.hpp

#ifndef AssembleFEGeneralizedNewtonian_DECL_hpp
#define AssembleFEGeneralizedNewtonian_DECL_hpp
 
#include "feddlib/core/AceFemAssembly/AssembleFE.hpp"
#include "feddlib/core/AceFemAssembly/specific/AssembleFENavierStokes.hpp"
#include "feddlib/core/FE/Helper.hpp"
#include "feddlib/core/FEDDCore.hpp"
#include "feddlib/core/LinearAlgebra/Matrix.hpp"
#include "feddlib/core/LinearAlgebra/MultiVector.hpp"
#include "feddlib/core/General/DifferentiableFuncClass.hpp"
// Add the generalized Newtonian Fluid models for the viscosity
#include "feddlib/core/AceFemAssembly/specific/GeneralizedNewtonianModels/CarreauYasuda.hpp"
#include "feddlib/core/AceFemAssembly/specific/GeneralizedNewtonianModels/PowerLaw.hpp"
#include "feddlib/core/AceFemAssembly/specific/GeneralizedNewtonianModels/Dimless_Carreau.hpp"
 
namespace FEDD
{
 
    template <class SC = default_sc, class LO = default_lo, class GO = default_go, class NO = default_no>
    class AssembleFEGeneralizedNewtonian : public AssembleFENavierStokes<SC, LO, GO, NO>
    {
    public:
        typedef Matrix<SC, LO, GO, NO> Matrix_Type;
        typedef Teuchos::RCP<Matrix_Type> MatrixPtr_Type;
 
        typedef SmallMatrix<SC> SmallMatrix_Type;
        typedef Teuchos::RCP<SmallMatrix_Type> SmallMatrixPtr_Type;
 
        typedef MultiVector<SC, LO, GO, NO> MultiVector_Type;
        typedef Teuchos::RCP<MultiVector_Type> MultiVectorPtr_Type;
 
        typedef AssembleFE<SC, LO, GO, NO> AssembleFE_Type;
 
        typedef DifferentiableFuncClass<SC, LO, GO, NO> DifferentiableFuncClass_Type;
        typedef Teuchos::RCP<DifferentiableFuncClass_Type> DifferentiableFuncClassPtr_Type;
 
        typedef InputToOutputMappingClass<SC, LO, GO, NO> InputToOutputMappingClass_Type;
        typedef Teuchos::RCP<InputToOutputMappingClass_Type> InputToOutputMappingClassPtr_Type;
        // smart pointer inside we need a type

        void assembleJacobian() override;

        void assembleRHS() override;

        void assembleJacobianBlock(LO i){TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error, "No implementation");};

        void computeLocalconstOutputField() override;
 
        /*
            \brief Assembly of FixedPoint- Matrix (System Matrix K with current u) 
         */
        void assembleFixedPoint();
 
       SmallMatrixPtr_Type getFixedPointMatrix(){return this->ANB_;};
 
    protected:
        std::string shearThinningModel;
        int dofsElementViscosity_;

        AssembleFEGeneralizedNewtonian(int flag, vec2D_dbl_Type nodesRefConfig, ParameterListPtr_Type parameters, tuple_disk_vec_ptr_Type tuple);

        void assemblyStress(SmallMatrixPtr_Type &elementMatrix);

        void assemblyStressDev(SmallMatrixPtr_Type &elementMatrix);
        void assemblyOutflowNeumannBoundaryTerm(SmallMatrixPtr_Type &elementMatrix){TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error, "Will be added when normal evaluation on the finite element is finished");} ;

        void assemblyOutflowNeumannBoundaryTermDev(SmallMatrixPtr_Type &elementMatrix){TEUCHOS_TEST_FOR_EXCEPTION(true, std::logic_error, "Will be added when normal evaluation on the finite element is finished");};

        void computeShearRate(vec3D_dbl_Type dPhiTrans, vec_dbl_ptr_Type &gammaDot, int dim);
 
        friend class AssembleFEFactory<SC, LO, GO, NO>; // Must have for specfic classes
 
        InputToOutputMappingClassPtr_Type viscosityModel; // viscosity Model can be in theory any Input to output mapping
 
    private:
    };
 
}
#endif