gmpThermoInterface.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_THERMO_INTERFACE_H_
#define _GEMA_PLUGIN_THERMO_INTERFACE_H_
 
 
#include <gmpFemPhysicsCommon.h>
#include <gmGaussAccessor.h>
#include "gmpThermoFemPhysics.h"
#include <gmTrace.h>
 
class GmGaussAccessor;
 
 
class GmpThermoInterface : public GmpThermoFemPhysics
{
public:
  // constructor
  GmpThermoInterface(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                     const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpThermoInterface();
 
  enum InterfacePropertyIds
  {
    IOPEN_ID = GmpThermoFemPhysics::NUM_PROPER_IDS,  
    Kti_ID,    
    Htc_ID,    // !< Id for retrieving the heat transfer coef. at top accessor
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPER_IDS,  
  };
 
  virtual const QVariantMap* physicsMetaDataMap();
  virtual bool checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& cellTable, LuaTable& gaussTable);
  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual GmpFemPhysics::FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  virtual bool calcDerivedResults(bool nonLinearSolver);
 
  // See comments on the base class
  virtual bool addStateItemsToGroup(GmStateDump* state, bool fixedHint, int groupId);
 
  virtual void fillStorageCompMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Sc, const GmVector& N, double c);
 
  virtual void fillFlowMatrix(const GmElement* e, GmMatrix& Kc, const GmMatrix Bp, const GmVector flowParam, const GmVector N, double c);
 
  //virtual bool setInitialGaussAttributes(const GmElement* e, const GmVector* coord, int ip) const;
  virtual bool setInitialConditions() const;
 
  virtual double fillBtMatrix(const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bt) = 0;
 
  virtual void fillElementTemperature(const GmElement* e, GmVector& pe);
 
  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 bool fillThermalConductivityParameters(const GmElement* e, GmVector& flowParam, const GmVector* coord, int i);
  
  virtual double heatTransferCoeficient(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    if (propertyAc(Htc_ID) == NULL)
    {
      return 0;
    }
    return propertyAc(Htc_ID)->scalarValueAt(e, coord, ip);
  }
 
  virtual double fractureHeatConductivity(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    assert(e);
    if (propertyAc(Kti_ID) == NULL)
    {
      return 0;
    }
    return propertyAc(Kti_ID)->scalarValueAt(e, coord, ip);
  }
 
protected:
  virtual const GmElementDof* dofMapping(GmCellType type) const;
 
  void fillMiddleCoordinates(const GmElement* e, GmMatrix Xnode, GmMatrix& MX);
};
 
#endif