gmpHydraulicDualFlow.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_DUALFLOW_H_
#define _GEMA_PLUGIN_HYDRAULIC_DUALFLOW_H_
 
#include "gmpHydraulic.h"
#include <gmpFemPhysicsCommon.h>
#include <gmpFemPhysicsData.h>
#include <gmSimulationData.h>
#include <gmPropertySet.h>
#include <gmTrace.h>
#include <math.h>
 
 
class GmGaussAccessor;
 
class GmpHydraulicDualFlow : public GmpHydraulic
{
public:
  // constructor
  GmpHydraulicDualFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                       const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpHydraulicDualFlow();
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
 
  // 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 void fillFractureStorageCompMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Sfc, const GmVector& N, double c) const;
  virtual void fillTransferMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Wfm, const GmVector& N, double c) const;
  virtual void fillFluidFlowMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Kfc, const GmMatrix& Bp, double c) const;
  virtual void fillElementFluidPorePressure(const GmElement* e, GmVector& Pfe);
  virtual void fillDampingMatrix(const GmElement* e, GmMatrix& C, const GmMatrix& Sc, const GmMatrix& Sfc, const QVector<int>& iPe, const QVector<int>& iPfe) const;
  virtual void fillStiffnessMatrix(const GmElement* e, GmMatrix& K, const GmMatrix& Kc, const GmMatrix& Kfc, const GmMatrix& Wfm, const QVector<int>& iPe, const QVector<int>& iPfe) const;
  virtual void fillInternalForceVector(const GmElement* e, GmVector& Fi, const GmMatrix& Kc, const GmMatrix& Kfc,
      const GmMatrix& Wfm, const GmVector& pe, const GmVector& pfe, const QVector<int>& iPe, const QVector<int>& iPfe) const;
  virtual void FillDofIndexPosition(const GmElement* e, QVector<int>& iPe, QVector<int>& iPfe) const;
  virtual bool calcDerivedResults(bool nonLinearSolver);
  // Returns the fracture property set
  virtual GmPropertySet* fracturePropertySet(int id) const;
 
protected:
        // See comments on the base class
        virtual const GmElementDof* dofMapping(GmCellType type) const;
  
  // See comments on the base class
  virtual bool checkAndLoadPrivateData(LuaTable& table);
 
  enum DualFlowStateVarIds
  {
    PF_ID = GmpHydraulic::NUM_STATEVAR_IDS,    
  
    // --- NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };
 
  enum DualFlowNodeAttributeIds
  {
    VFF_NA_ID = GmpHydraulic::NUM_NA_IDS,    
    VEF_NA_ID,                               
 
    // ------    NO ADDING BELOW THIS LINE
    NUM_NA_IDS,  
  };
 
  enum DualFlowGaussAttributeIds
  {
          VEF_GA_ID = GmpHydraulic::NUM_GA_IDS,  
 
          // --- NO ADDING BELOW THIS LINE
          NUM_GA_IDS,  
  };
 
 enum DualFlowBoundaryConditionIds
 {
   FIXED_PF_BC_ID = GmpHydraulic::NUM_BC_IDS,  
   FIXED_QF_BC_ID,                             
 
   // --- NO ADDING BELOW THIS LINE
   NUM_BC_IDS,  
 };
 
 enum DualFlowBoundaryConditionValueIds
 {
   BC_PF_ID = GmpHydraulic::NUM_BCV_IDS,       
   BC_QWF_ID,                                  
 
   // --- NO ADDING BELOW THIS LINE
   NUM_BCV_IDS,  
 };
 
 enum DualFlowPhysicsAttributeIds
 {
   FRACPROPSET_ID = GmpHydraulic::NUM_ATTRI_IDS, 
 
   // ------         NO ADDING BELOW THIS LINE
   NUM_ATTRI_IDS,    
 };
 
 
 virtual void calcElementDerivedResultAtPoints(const GmElement* e, int resultId, const GmMatrix& evalPoints, bool evalAtIp, GmMatrix& result);
 
 void calcElementVelocityAtPoints(const GmElement* e, const GmMatrix& evalPoints, GmMatrix& v);
 void calcElementFlowAtPoints(const GmElement* e, const GmMatrix& evalPoints, bool evalAtIp, GmMatrix& q);
 
//private:
  GmpFemPhysicsData _fractureData;   
  GmValueAccessor*  _anAcc;          
};
 
#endif