gmpHydraulicBrinkmanFlow.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_HYDRAULIC_BRINKMANFLOW_H_
#define _GEMA_PLUGIN_HYDRAULIC_BRINKMANFLOW_H_
 
#include "gmpHydraulic.h"
#include <gmpFemPhysicsCommon.h>
#include <gmTrace.h>
#include <math.h>
 
 
class GmGaussAccessor;
 
class GmpHydraulicBrinkmanFlow : public GmpHydraulic
{
public:
  // constructor
  GmpHydraulicBrinkmanFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                           const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpHydraulicBrinkmanFlow();
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
 
  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
 
  // fill elemental matrix and vectors (K, Fi)
  virtual GmpFemPhysics::FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  virtual bool fixedNodalForcesBc(QVector<int>& nodes, QVector<int>&dof, QVector<double>& values) const;
  
  virtual bool fixedNodalDofsBc  (QVector<int>& nodes, QVector<int>&dof, QVector<double>& values, bool* constantValues) const;
  
  virtual bool fillNeMatrix(GmMatrix& Nec, GmMatrix& Np, const GmVector& Ne, const GmVector& N) const;
  
  virtual void fillFluidMomentMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Ac, const GmMatrix& Bp, const GmMatrix& N, double c);
  
  virtual void fillStorageCompMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Sc, const GmVector& N, double c);
 
  virtual void fillFluidConservMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Cm, const GmMatrix& Bp, const GmMatrix& Np, double c)const;
  
  virtual void fillStiffnessMatrix(const GmElement* e, GmMatrix& K, const GmMatrix& Ac, const GmMatrix& Cm, const QVector<int>& iVe, const QVector<int>& iPe) const;
  
  virtual bool calcDerivedResults(bool nonLinearSolver);
 
  virtual void fillElementVelocities(const GmElement* e, GmVector& Vfe);
  
  virtual void fillDampingMatrix(const GmElement* e, GmMatrix& C, const GmMatrix& Sc, const QVector<int>& iVe) const;
  
  virtual void fillInternalForceVector(const GmElement* e, GmVector& Fi, const GmMatrix& Ac, const GmMatrix& Cm, const GmVector& pe, const GmVector& vfe,  const QVector<int>& iVe, const QVector<int>& iPe) const;
  
  virtual void FillDofIndexPosition(const GmElement* e, QVector<int>& iPe, QVector<int>& iPfe);
  
  virtual GmpFemPhysics::FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
          int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  virtual GmpFemPhysics::FemResultType fillBcElementData(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
          int border, GmMatrix& elemMat, GmVector& elemVec);
  
  virtual bool fillFixedNodalPressureBcLists(int dofId, QVector<int>& nodes, QVector<int>&dof, QVector<double>& values) const;
 
  virtual int numBCData() const { return _bcData.size(); };
 
  void addBCData(int id, double value) { _bcData.insert(id, value); };
 
  virtual double bcData(int id) const { return _bcData.value(id); };
 
  virtual QList<int> nodeIds() const { return _bcData.uniqueKeys(); };
 
  void clearBcData() { _bcData.clear(); };
 
  //virtual GmpFemPhysics::FemResultType fillPressureValues(const GmElement* e, GmVector& elemVec, int border);
  virtual GmpFemPhysics::FemResultType fillPressureValues(const GmElement* e, double press, int border);
protected:
        // See comments on the base class
        virtual const GmElementDof* dofMapping(GmCellType type) const;
 
  enum BrinkmanStateVarIds
  {
    VF_ID = GmpHydraulic::NUM_STATEVAR_IDS, 
    PL_BC_ID,                               
 
    // --- NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };
 
 enum BrinkmanBoundaryConditionIds
 {
   FIXED_VF_BC_ID = GmpHydraulic::NUM_BC_IDS,  
   
   // --- NO ADDING BELOW THIS LINE
   NUM_BC_IDS,  
 };
 
 enum BrinkmanBoundaryConditionValueIds
 {
   BC_VX_ID = GmpHydraulic::NUM_BCV_IDS,      
   BC_VY_ID,                                  
   BC_VZ_ID,                                  
   BC_PL_ID,                                  
 
   // --- NO ADDING BELOW THIS LINE
   NUM_BCV_IDS,  
 };
 
//private:
  QMap< int, double >  _bcData;     
 
};
 
#endif