gmpMaterialKhoei.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_KHOEI_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_KHOEI_H_
 
#include "gmpHydraulicMaterialMultiPhase.h"
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMaterialKhoei : public GmpHydraulicMaterialMultiPhase
{
protected:
    enum KhoeiPropIds
    {
        Bc_ID = GmpHydraulicMaterialMultiPhase::NUM_PROPERTY_IDS, 
 
        // --- NO ADDING BELOW THIS LINE
        NUM_PROPERTY_IDS,  
    };
 
public:
 
  GmpHydraulicMaterialKhoei(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpHydraulicMaterialMultiPhase(typeIndex, typeName, logger) {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpHydraulicMaterialKhoei(typeIndex, typeName, logger);
  }
 
  // Sets multiPhase material map
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual double waterSaturation(const GmElement* e, const GmVector* coord, int ip, double pc) const;
 
  virtual double derivativeWaterSaturation(const GmElement* e, const GmVector* coord, int ip, double pc) const;
 
  virtual double waterRelativePermeability(const GmElement* e, const GmVector* coord, int ip, double Sw) const;
 
  virtual double gasRelativePermeability(const GmElement* e, const GmVector* coord, int ip, double Sw) const;
 
  virtual double linearInterpolate(const GmVector& x, const GmVector& y, double xi) const;
 
  virtual double BcParameterKhoei(const GmElement* e, const GmVector* coord, int ip) const
  {
      if (propertyAc(Bc_ID) == NULL)
      {
          return 0.0;
      }
      return propertyAc(Bc_ID)->scalarValueAt(e, coord, ip);
  }
 
};
 
#endif