gmpDGFemPhysics.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/
 
#ifndef _GEMA_PLUGIN_DG_FEM_PHYSICS_H_
#define _GEMA_PLUGIN_DG_FEM_PHYSICS_H_
 
#include "gmpDGConfig.h"
#include <gmpFemPhysicsCommon.h>
 
#include <gmCellAccessor.h>
#include <gmBoundaryConditionAccessor.h>
#include <gmGaussAccessor.h>
#include <gmIntegrationRule.h>
#include <gmElement.h>
#include <gmCellGroupSet.h>
#include <gmElementDof.h>
#include <gmShapeH.h>
 
#include <luaTable.h>
#include <luaStackBalancer.h>
#include <gmModelData.h>
 
#include <gmTrace.h>
 
class GMP_DG_FEM_PHYSICS_API_EXPORT GmpDGFemPhysics : public GmpFemPhysicsCommon
{
 
public:
  GmpDGFemPhysics(GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger);
  virtual ~GmpDGFemPhysics();
 
  // Comments on the base class
  virtual const char * pluginName() const { return "DGFemPhysics"; }
  virtual const char * pluginType() const { return "Galerkin"; };
 
  enum StateVarIds
  {
          C_ID,                       
  };
 
 
  enum PhysicsAttributeIds
  {
          NodsX_ID,                                   
          NodsY_ID,                                   
          NodsZ_ID,                               
          Upwind_ID,            
          Implicit_ID,          
          Limiter_ID,           
          Regular_Mesh_ID,      
          Periodic_Bnd_ID,                      
          Couple_ID,                              
  };
 
  enum ElementPropertyIds {
          PHI_ID,               
          V_ID,                 
          R_ID,                 
          F_ID,                 
          IC_ID,                
          A_ID,                 
  };
 
  enum GaussAttributeIds
  {
          V_D_ID,                           
          C_D_ID,
  };
 
  enum BoundaryContidionsIds {
          BC_C_Inflow_1D_ID,                            
          BC_C_Outflow_1D_ID,                           
          BC_C_Periodic_1D_ID,
          BC_C_Inflow_ID,                   
          BC_C_Outflow_ID,                  
          BC_Periodic_ID,                                   
          BC_C_NODE_ID
  };
 
  enum BoundaryConditionsValueIds {
          BC_CELL_1D_ID,
          BC_C_CELL_ID,
          BC_CELL_FREE_1D_ID,
          BC_C_FREE_CELL_ID,
          BC_C_PERIODIC_1D_ID,
          BC_C_PERIODIC_CELL_ID,
          BC_C_ID,
          BC_C_FREE_ID,
          BC_C_PERIODIC_ID,
          BC_C_NODE_V_ID,
  };
 
  virtual bool loadPrivateData(LuaTable & table);
  virtual const GmElementDof * dofMapping(const GmElement * e) const;
  virtual const QVariantMap* physicsMetaDataMap();
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  virtual FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
          int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  virtual bool fixedNodalDofsBc(QList<int>& nodes, QList<int>&dof, QList<double>& values, bool* constantValues) const;
  bool calcDerivedResults(bool);
  bool calcElementDerivedResults();
 
protected:
        int  BndIndex = -1;                     
        double TOLERANCE = +1.e-3;              
 
private:
        int  numShapeFunctions;                 
        double detAffineTransformationMAT;      
        GmMatrix AffineTransformationMAT;       
        GmMatrix InvAffineTransformationMAT;    
        GmMatrix MatFaceIntPoint;
 
        QList<GmVector> FNeighborhood;          
        QList<GmVector> SNeighborhood;          
        QList<GmVector> TNeighborhood;          
 
        bool  setNumFunctions();                                              
        
        //<! Get of number of shape functions for each element
        int   getNumFunctions()const { return numShapeFunctions;};
        //<! Get of Affine Transform Matrix
        GmMatrix getAffineTransformationMAT()    { return AffineTransformationMAT; };
        //<! Get of Inverse Affine Transform Matrix
        GmMatrix getInvAffineTransformationMAT() { return InvAffineTransformationMAT; };
        //<! Get of Determinant of Affine Transform Matrix
        double getDetAffineTransformMatrix() { return detAffineTransformationMAT; };   
        
        bool TopologyMeshSetup();               
        bool TopologyMeshSetup1D();             
        bool TopologyMeshSetup2D();                     
        bool TopologyMeshSetup3D();                     
 
        bool AffineTransformationSetup();       
        bool InitialConditionProjection();      
        
        /*<! Set up Numerical Flux approximation used in Discontinuous Galerkin Methods apply to Advective Boundary-Value Problems
           User can choose two classical numerical flux, Upwind Numerical or Lax-Wendrolf schemes.      */
        bool NumericalFlux(const GmElement * e, int borderE, GmMatrix &MFluxBorderE, GmMatrix &MFluxBorderN);
 
        bool LinearLimeter();                                      //<! Linear Reconstruction function. 
        GmMatrix MatAvrgCalc(const GmElement *e);
        
        int      TClassicElementToHigher(const GmElement * e, int)const;      //<! Geometric Transformation between Classical Geometric Elemento to Hierarchical.
        int  ElemPosition(int i, int j, int k);                                                   //<! Auxiliar Function to 3D Topology Mesh. (i,j,k) Index.    
        int  Edge_Face_Ngh(int border);                                                                   //<! Return of Edge/Face Neighbor to border.          
        GmVector DofEdge_Face(const GmElement * e,int border)const;               //<! Auxiliar function for identify to edge/face functions defined on boundary element.       
        bool PrintResult();                                                                                               //<! Print Result function.   
        bool ErrorAnalysis();                                             //<! Error Analysis.
        bool PeriodicBND();
        bool Projection(void);
        int  EdgeVd(int i, int border);
};
 
#endif