gmpMaterialDualPorosity.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_DUALPOROSITY_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_DUALPOROSITY_H_
 
#include "gmpHydraulicMaterial.h"
#include "gmpMaterialSaturated.h"
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMaterialDualPorosity: public GmpHydraulicMaterialSaturated
{
protected:
  enum DualPorosityElementPropertyIds
  {
    Sx_ID = GmpHydraulicMaterialSaturated::NUM_PROPERTY_IDS,  
    Bx_ID,                                                    
    Sy_ID,                                                    
    By_ID,                                                    
    Sz_ID,                                                    
    Bz_ID,                                                    
    Bpf_ID,                                                   
    Uf_ID,                                                   
    shapeFactor_ID,                                           
    Kf_User_ID,                                               
    PorF_User_ID,                                             
    Wff_User_ID,                                              
 
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpHydraulicMaterialDualPorosity(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpHydraulicMaterialSaturated(typeIndex, typeName, logger) {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpHydraulicMaterialDualPorosity(typeIndex, typeName, logger);
  }
 
  enum shapeFactorType{
    warrenRoot,
    kazemi,
    limAzis,
        coat, 
    user,
          none
  };
 
  // set dual porosity material map
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual double fractureSpacingX(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
          if (propertyAc(Sx_ID) == NULL)
          {
                  return 0;
          }
    return propertyAc(Sx_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureApertureX(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
          if (propertyAc(Bx_ID) == NULL)
          {
                  return 0;
          }
    return propertyAc(Bx_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureSpacingY(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
          if (propertyAc(Sy_ID) == NULL)
          {
                  return 0;
          }
    return propertyAc(Sy_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureApertureY(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
         if (propertyAc(By_ID) == NULL)
         {
                 return 0;
         }
    return propertyAc(By_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureSpacingZ(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
          if (propertyAc(Sz_ID) == NULL)
          {
                  return 0;
          }
    return propertyAc(Sz_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureApertureZ(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
          if (propertyAc(Bz_ID) == NULL)
          {
                  return 0;
          }
    return propertyAc(Bz_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fracturePoreCompressibility(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
        if (propertyAc(Bpf_ID) == NULL)
        {
                return 0;
        }
    return propertyAc(Bpf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureViscosity(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    return propertyAc(Uf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double userDefinedFracturePorosity(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    return propertyAc(PorF_User_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual void userDefinedFracturePermeability(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Kf) const
  {
    S_TRACE();
    assert(e);
    GmpFemPhysicsCommon::fillMatrixFromPropertyAc(propertyAc(Kf_User_ID), e, coord, ip, Kf);
  }
 
  virtual double userDefinedShapeFactor(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    return propertyAc(Wff_User_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual bool shapeFactorModel(const GmElement* e, GmpHydraulicMaterialDualPorosity::shapeFactorType& type, const GmVector* coord, int ip) const;
  virtual void setfrac(double& sx, double& sy, double&sz, int& fx, int& fy, int&fz, int&n, double&sm) const;
  virtual double fillFracturePorosity(const GmElement* e, const GmVector* coord, int ip) const;
  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;
  
};
 
 
#endif