gmpMaterialInterfaceMultiPhaseFlow.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_INTERFACE_MULTIPHASE_FLOW_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_INTERFACE_MULTIPHASE_FLOW_H_
 
#include <gmpFemPhysicsCommonMaterial.h>
#include "gmpMaterialInterfaceFlow.h"
#include "gmpHydraulicConfig.h"
#include "gmpHydraulicInterfaceTwoPhaseFlow.h"
//#include "gmpHydraulicMaterialMultiPhase.h"
#include "gmpHydraulicMaterialTwoPhase.h"
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpMaterialInterfaceMultiPhaseFlow: public GmpHydraulicMaterialTwoPhase
{
protected:
    
  enum TwoPhaseInterfacePropertyIds
  { 
      GAP_ID = GmpHydraulicMaterialTwoPhase::NUM_PROPERTY_IDS,  
      KTTOP_ID,                                                 
      KTBOT_ID,                                                 
      IOPEN_ID,                                                 
 
      // --- NO ADDING BELOW THIS LINE
      NUM_PROPERTY_IDS,  
  };
 
public:
 
    enum InterfaceGaussAttributeIds
    {
        An_GA_ID,             
        AnOLD_GA_ID,
 
        // --- NO ADDING BELOW THIS LINE
        NUM_GA_IDS,  
    };
 
  GmpMaterialInterfaceMultiPhaseFlow(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpHydraulicMaterialTwoPhase(typeIndex, typeName, logger) {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialInterfaceMultiPhaseFlow(typeIndex, typeName, logger);
  }
 
  // set interface flow material map
  virtual const QVariantMap* materialMetaDataMap();
 
  // See comments on the base class
  virtual bool addStateItemsToGroup(GmStateDump* state, GmElementMesh* mesh, bool fixedHint, int groupId);
 
  virtual double fractureOpening(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE();
      assert(e);
      return propertyAc(GAP_ID)->scalarValueAt(e, coord, ip);
  }
  
  virtual double oldFractureAperture(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    return gaussAttrAc(AnOLD_GA_ID)->scalarValueAt(e, ip, coord);
  }
 
  virtual double transversalPermeabilityTopFace(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    if (propertyAc(KTTOP_ID) == NULL)
    {
      return 0;
    }
    return propertyAc(KTTOP_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double transversalPermeabilityBottomFace(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE();
      assert(e);
      if (propertyAc(KTBOT_ID) == NULL)
      {
          return 0;
      }
      return propertyAc(KTBOT_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double initiallyOpenElem(const GmElement* e, const GmVector* coord, int ip) const
  {
      S_TRACE();
      assert(e);
      if (propertyAc(IOPEN_ID) == NULL)
      {
          return 1;
      }
      return propertyAc(IOPEN_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double cubicLaw(const GmElement* e, const GmVector* coord, int ip) const;
  
  virtual void hydraulicConstitutiveModel(const GmElement* e, GmVector& flowParam, const GmVector* coord, int ip) const;
  
  virtual bool initialCondition(const GmElement* e, const GmVector* coord, int ip) const;
};
 
#endif