gmpMaterialCubicLaw.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_CUBICLAW_H_
#define _GEMA_PLUGIN_HYDRAULIC_MATERIAL_CUBICLAW_H_
 
#include <gmTrace.h>
#include <assert.h>
 
#include "gmpMaterialInterfaceFlow.h"
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpMaterialCubicLaw : public GmpMaterialInterfaceFlow
{
protected:
  enum CubicLawPropIds
  {
    //GAP_ID = GmpMaterialInterfaceFlow::NUM_PROPERTY_IDS, //!< Id for retrieving the parameter that controls the shape of the curve Sw x pc
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPERTY_IDS,  
  };
 
public:
 
  GmpMaterialCubicLaw(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMaterialInterfaceFlow(typeIndex, typeName, logger) {}
 
  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialCubicLaw(typeIndex, typeName, logger);
  }
 
 
  virtual double cubicLaw(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
 
    double kl = 0.0;
    double d_n = gaussAttrAc(AnOLD_GA_ID)->scalarValueAt(e, ip, coord);
    double uf = fluidViscosity(e, coord, ip);
 
    // computes the cubic law
    if (uf > 0.0)
    {
      if (d_n > 0.0)
      {
        kl = pow(d_n, 3.0) / (12.0 * uf);
      }
 
    }
    else
    {
      assert(uf != 0 && "Ufw Dynamic viscosity must be different to zero");
    }
 
    return kl;
  }
 
  virtual void hydraulicConstitutiveModel(const GmElement* e, GmVector& flowParam, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    flowParam.fill(0.0);
 
    //Compute the longitudinal conductivity according to cubic's law.
    flowParam(0) = cubicLaw(e, coord, ip); 
    //computes the leakoff coefficients
    flowParam(1) = leakoffTop(e, coord, ip);
    flowParam(2) = leakoffBottom(e, coord, ip);
  }
 
};
 
#endif