AssembleFENavierStokes_decl.hpp

#ifndef ASSEMBLEFENAVIERSTOKES_DECL_hpp
#define ASSEMBLEFENAVIERSTOKES_DECL_hpp
 
#include "feddlib/core/AceFemAssembly/AssembleFE.hpp"
#include "feddlib/core/FEDDCore.hpp"
#include "feddlib/core/FE/Helper.hpp"
#include "feddlib/core/LinearAlgebra/Matrix.hpp"
#include "feddlib/core/LinearAlgebra/MultiVector.hpp"
 
namespace FEDD {
 
template <class SC = default_sc, class LO = default_lo, class GO = default_go, class NO = default_no>
class AssembleFENavierStokes : public AssembleFE<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;
 

    void assembleJacobian() override;

    void assembleRHS() override;

    void assembleJacobianBlock(LO i) override {};
    
    void setCoeff(SmallMatrix_Type coeff);

    void assembleFixedPoint();
 
    SmallMatrixPtr_Type getFixedPointMatrix(){return ANB_;};
 
   protected:

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

    void assemblyLaplacian(SmallMatrixPtr_Type &elementMatrix);

    void assemblyAdvection(SmallMatrixPtr_Type &elementMatrix);
    
    void assemblyAdvectionInU(SmallMatrixPtr_Type &elementMatrix); 

    void assemblyDivAndDivT(SmallMatrixPtr_Type &elementMatrix);
 
    friend class AssembleFEFactory<SC,LO,GO,NO>; // Must have for specfic classes
 
    void buildTransformation(SmallMatrix<SC>& B);
 
    void applyBTinv(vec3D_dbl_ptr_Type& dPhiIn,
                    vec3D_dbl_Type& dPhiOut,
                    SmallMatrix<SC>& Binv);
 
    //tuple_disk_vec_ptr_Type returnTuple(); /// @todo return tuple in case or check tuple

    int dofsVelocity_;
    int dofsPressure_;
 
    std::string FETypeVelocity_;
    std::string FETypePressure_;
 
    int numNodesVelocity_;
    int numNodesPressure_;
 
    int dofsElementVelocity_ ;
    int dofsElementPressure_ ;
 
    int dofsElement_;
 
    vec_dbl_Type solutionVelocity_;
    vec_dbl_Type solutionPressure_;
 
    SmallMatrixPtr_Type constantMatrix_;
    SmallMatrixPtr_Type ANB_;
 
    SmallMatrix_Type coeff_;
 
    double viscosity_ ;
    double density_ ;
 
    std::string linearization_;
 
   private:
 
    
 };
 
}
#endif