gmpMaterialEnhancedDualPorosity.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_MATERIAL_ENHANCEDDUALPOROSITY_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_ENHANCEDDUALPOROSITY_H_
 
#include "gmpHydraulicMaterial.h"
#include <gmpFemPhysicsData.h>
#include "gmpMaterialGeneralizedDualPorosity.h"
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMaterialEnhancedDualPorosity : public GmpHydraulicMaterialGeneralizedDualPorosity
{
protected:
  GmpFemPhysicsData* _matFracData;   
  //GmValueAccessor*   _anAcc;          //!< fracture aperture accessor
  GmElementMesh*        _mesh;       
 
  enum DualPorosityElementPropertyIds
  {
    GAP_ID = GmpHydraulicMaterialGeneralizedDualPorosity::NUM_PROPERTY_IDS, 
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpHydraulicMaterialEnhancedDualPorosity(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpHydraulicMaterialGeneralizedDualPorosity(typeIndex, typeName, logger) 
  {
    _mesh = NULL;
        _matFracData = NULL;
  }
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpHydraulicMaterialEnhancedDualPorosity(typeIndex, typeName, logger);
  }
 
  enum fracDir
  {
    horizontal,     
    vertical,       
    oriented,       
  };
 
  // set dual porosity material map
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual bool setFractureData(GmpFemPhysicsData* dataF) { _matFracData = dataF; return true; };
  
  //virtual void setFractureApertureAcc(GmValueAccessor* anAcc) { _anAcc = anAcc; };
 
  virtual void setMeshAccessor(GmElementMesh* meshAc) { _mesh = meshAc; };
 
  virtual bool fillFracturePermeability(const GmElement* e, GmMatrix& Kper, const GmVector* coord, int ip) const;
  
  virtual bool fillShapeFactor(const GmElement* e, double& ws, GmMatrix Km, const GmVector* coord, int ip)const;
  
  virtual GmMatrix fillFractureLocalPermeability(int dim, double bi, double si) const
  {
          GmMatrix Kfrl(dim, dim);
          Kfrl.fill(0);
 
    //if ((bi == 0) || (si == 0))
    //{
    //  return Kfrl;
    //}
    // fills fracture permeability
          if (dim == 2)
          {
                  // fills bidimensional rotation matrix
                  Kfrl(0, 0) = pow(bi, 3) / (12 * si);
          }
          if (dim == 3)
          {
                  // fills threedimensional rotation matrix
                  Kfrl(0, 0) = pow(bi, 3) / (12 * si);
                  Kfrl(1, 1) = pow(bi, 3) / (12 * si);
          }
          return Kfrl;
  };  
 
  virtual double meanDistance(double dx, double dy) const;
 
  virtual int checkFractureOrientation(GmVector& data, int d) const;
 
  virtual double fillShapeFactorOrtogonalFracture(GmMatrix& Xe, GmVector& data) const;
 
  virtual double fillShapeFactorOrientedFracture(GmMatrix& Xe, GmVector& data) const;
 
  virtual double fractureAperture(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    return propertyAc(GAP_ID)->scalarValueAt(e, coord, ip);
  }
 
private:
  GmMatrix fillRotationMatrix(const GmMatrix& Xe) const;
};
 
#endif