gmpThermalFlow.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_THERMALFLOW_H_
#define _GEMA_PLUGIN_THERMALFLOW_H_
 
#include "gmpThermoFemPhysics.h"
#include <gmTrace.h>
#include <math.h>
 
 
template <class T> class GmpThermalFlow : public T
{
public:
  GmpThermalFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger)
    : T(pluginType, simulation, id, description, logger)
  {}
 
protected:
  virtual double axisymmetricFactor(const GmElement* e, const GmMatrix& X, const GmVector& N) const { Q_UNUSED(e); Q_UNUSED(X); Q_UNUSED(N); return 1.0; }
  
  //virtual double fillBtMatrix(const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bp)
  //{
  //  S_TRACE();
  //  Q_UNSED(N);
  //  assert(shape);
  //  double detJ = 0;
  //  if (issoParametric())
  //  {
  //    detJ = shape->shapeCartesianPartialsFromCoord(ip, X, B, true);
  //  }
  //  else
  //  {
  //    detJ = shape->shapeCartesianPartialsFromJacobian(ncoord, J, Bp, true);
  //  }
 
  //  return detJ;  
  //}
 
};
 
template <class T> class GmpThermalAxFlow: public T
{
public:
    GmpThermalAxFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger)
      : T(pluginType, simulation, id, description, logger)
    {}
 
protected:
  virtual bool isAxisymmetric() { return true; }
 
  // Comments on the base class
  virtual double axisymmetricFactor(const GmElement* e, const GmMatrix& X, const GmVector& N) const
  {
    S_TRACE();
    int n = e->numNodes();
    // evaluate radio at integration points
    double rGauss = 0.0;
    for (int j = 0; j < n; ++j)
      rGauss += N(j)*X(j, 0);
 
    return 2 * arma::datum::pi * rGauss;
  }
};
 
#endif