gmpHydraulicEnhancedDualFlow.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_ENHANCEDDUALFLOW_H_
#define _GEMA_PLUGIN_HYDRAULIC_ENHANCEDDUALFLOW_H_
 
#include "gmpHydraulicDualFlow.h"
#include <gmpFemPhysicsData.h>
#include <gmpFemPhysicsCommon.h>
#include <gmDiscontinuitySet.h>
#include <gmDiscontinuity.h>
#include "gmDiscontinuityAccessor.h"
#include <gmMathUtils.h>
#include <gmTrace.h>
#include <math.h>
 
 
class GmGaussAccessor;
 
class GmpHydraulicEnhancedDualFlow : public GmpHydraulicDualFlow
{
public:
  // constructor
  GmpHydraulicEnhancedDualFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                       const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpHydraulicEnhancedDualFlow();
 
  enum EmbeddedElemPropertyIds
  {
    UFw_ID = GmpHydraulicDualFlow::NUM_PROPER_IDS, 
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPER_IDS,  
  };
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
 
  // fill elemental matrix and vectors (K, Fi)
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // fill element data for continuous enriched element data
  virtual FemResultType fillContinuousEnrichedElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // Fills fracture fluid flow matrix
  virtual bool assemblyFractureFlowMatrix(GmMatrix& Kf, GmMatrix & Ke, GmVector& iPe);
 
  // Get permeability and shapefractor
  void getPermeabilityAndShapeFactor(const GmElement* e, GmMatrix& Km, double & Gw, double& factor);
 
  // Get equivalent distance
  double getEquivalentDistance(const GmElement* e, const GmElement* ef);
 
  // Sets the initial conditions for hydraulic embedded fracture models.
  virtual bool setInitialConditions();
 
  // Get the number of discontinuities crossing and element
  int getNumberOfDiscontinuities(const GmElement* e);
  int getNumberOfDiscontinuityElements(const GmElement* e);
  void getDiscontinuityGeometry(const GmElement* e, int idiscontinuity, GmVector& Xref, GmMatrix& mn);
  void get2DDiscontinuityGeometry(const GmElement* e, int idiscontinuity, GmVector& Xref, GmMatrix& mn);
  void get3DDiscontinuityGeometry(const GmElement* e, int idiscontinuity, GmVector& Xref, GmMatrix& mn);
  void getIdDiscontinuitySubElementInMesh(const GmElement* e, int k, int& idDiscontinuityMesh);
  bool fillDiscontinuityNodesGlobalToLocalMap(const GmElement* e, QMap<int, int>& globalToLocal);
  FemResultType fillDiscontinuitySetElementData(const GmElement* e, GmMatrix& Hd, GmMatrix& Wf, GmMatrix& Sd, GmMatrix& Tcd, const QMap<int, int>& globalToLocal);
  void getDiscontinuitySubElementLocalData(const GmElement* e, int k, int& idDiscontinuityElem, QVector<int>& dof, const GmElement* ef, const QMap<int, int>& dNodesGlobalToLocal);
  void assembleDiscontinuityElementMatrix(GmMatrix& Kg, const GmMatrix& Kl, const QVector<int>& dof) const;
 
protected:
        // See comments on the base class
  virtual bool checkAndLoadDiscontinuitySetAccessors(LuaTable& table);
 
private:
  const char*       _pluginType;     
 
};
 
#endif