gmpAdvDiffusionFemPhysics.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_ADV_DIFFUSION_FEM_PHYSICS_H_
#define _GEMA_PLUGIN_ADV_DIFFUSION_FEM_PHYSICS_H_
 
#include "gmpAdvDiffusionConfig.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_ADV_DIFFUSION_PHYSICS_API_EXPORT GmpAdvDiffusionFemPhysics : public GmpFemPhysicsCommon
{
public:
  GmpAdvDiffusionFemPhysics(GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger);
  virtual ~GmpAdvDiffusionFemPhysics();
 
  // Comments on the base class
  virtual const char* pluginName() const { return "AdvDiffusionFemPhysics"; }
 
  // Comments on the base class
  virtual const char* pluginType() const { return "advDiffusion"; }
 
  virtual bool loadPrivateData(LuaTable & table);
  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;
 
  virtual bool calcDerivedResults(bool nonLinearSolver);
 
protected:
        int  BndIndex = -1;                     
private:
  enum PhysicsAttributeIds
  {
          NodsX_ID,                                   
          NodsY_ID,                                   
          NodsZ_ID,                               
          Transient_ID,               
          Error_ID,                   
          staticCondensation_ID,      
          Couple_ID,                              
  };
 
  enum StateVarIds
  {
    P_ID,   
  };
 
  enum GaussAttributeIds
  {
         P_A_ID,
  };
 
  enum ElementPropertyIds
  {
          PHI_ID,                     
          K_ID,                       
          F_ID,                       
          PIC_ID,                     
          QIC_ID,                     
          PT_ID,                                          
          QT_ID,                                          
          GPT_ID,                                         
          DQT_ID                                          
  };
 
  enum BoundaryConditionIds
  {
          BC_P_ID,                   
          BC_Q_ID,                   
          BC_P_NODE_ID,
  };
 
  enum BoundaryConditionValueIds
  {
          BC_P_V_ID,                   
          BC_Q_V_ID,                   
          BC_P_NODE_V_ID,                   
  };
 
  /********************************************************************/
  QList<int> BCNodesGhostP;
  QList<int> BCNodesGhostQ;
  QList<int> BCEdge;
  QList<GmVector> BCElements;
  QList<GmVector> FNeighborhood;          
 
  int PressureDoF;
  int FluxDoF;
  int LSDoF;
  /********************************************************************/
  virtual bool loadPrivateData2D(LuaTable & table);
  virtual bool loadPrivateData3D(LuaTable & table);
 
  virtual bool ICProjection2D(const GmElement *e, QList<int>&  nodes, QList<double> & values)
  {
          ICPProjection2D(e, nodes, values);    
//        ICQProjection2D(e, nodes, values);  
          return true;
  };
 
  virtual bool ICPProjection2D(const GmElement *e, QList<int>&  nodes, QList<double> & values);
  virtual bool ICQProjection2D(const GmElement *e, QList<int>&  nodes, QList<double> & values);
 
  virtual bool ICProjection3D(const GmElement *e, QList<int>&  nodes, QList<double> & values)
  {
          ICPProjection3D(e, nodes, values);
          ICQProjection3D(e, nodes, values);
          return true;
  };
 
  virtual bool ICPProjection3D(const GmElement *e, QList<int>&  nodes, QList<double> & values);
  virtual bool ICQProjection3D(const GmElement *e, QList<int>&  nodes, QList<double> & values);
 
  virtual const GmElementDof* dofMapping(GmCellType type) const;
  virtual const GmElementDof* dofMapping2D(GmCellType type) const;
  virtual const GmElementDof* dofMapping3D(GmCellType type) const;
 
  virtual bool PrintResult(void);
  virtual bool ErrorAnalysis(void);
 
  bool Projection(void);
 
  int  ElemPosition(int i, int j, int k);                                                         //<! Auxiliar Function to 3D Topology Mesh. (i,j,k) Index.    
  int  TClassicElementToHigher(const GmElement * e, int)const;        //<! Geometric Transformation between Classical Geometric Elemento to Hierarchical.
 //<! Auxiliar function for identify to edge/face functions defined on boundary element.        
  GmVector DofEdge_Face(const GmElement * e, int border, const GmBoundaryCondition* bc)const;
  
};
 
#endif