gmpMaterialBrinkman.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_BRINKMAN_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_BRINKMAN_H_
 
#include "gmpHydraulic.h"  
 
#include <gmpFemPhysicsCommonMaterial.h>
#include <gmTrace.h>
#include <assert.h>
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMaterialBrinkman : public GmpFemPhysicsCommonMaterial
{
protected:
  enum BrinkmanElementPropertyIds
  {
    K_ID,          
    rhof_ID,       
    muf_ID,        
    mufef_ID,      
   
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpHydraulicMaterialBrinkman(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpFemPhysicsCommonMaterial(typeIndex, typeName, logger) {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpHydraulicMaterialBrinkman(typeIndex, typeName, logger);
  }
 
  // set dual porosity material map
  virtual const QVariantMap* materialMetaDataMap();
 
  virtual void permeability(const GmElement* e, const GmVector* coord, int ip, GmMatrix& K) const;
  
  virtual double fluidDensity(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
        //Accessor 
        GmCellAccessor* rhowAcc = propertyAc(rhof_ID);
 
        if (rhowAcc == NULL)
        {
                return 0.0;
        }
    return rhowAcc->scalarValueAt(e, coord, ip);
  }
  
  virtual double fluidViscosity(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
          return propertyAc(muf_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fluidEffectiveViscosity(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE(); assert(e);
    return propertyAc(mufef_ID)->scalarValueAt(e, coord, ip);
  }
};
 
 
#endif