gmpCoupledThermoHydraulic.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_THERMOHYDRAULIC_H_
#define _GEMA_PLUGIN_THERMOHYDRAULIC_H_
 
 
#include "gmpCoupledTHFemPhysics.h"
#include <gmpFemPhysicsCoupled.h>
#include <gmTrace.h>
 
class GmGaussAccessor;
 
class GMP_COUPLED_TH_PHYSICS_API_EXPORT GmpCoupledThermoHydraulic: public GmpCoupledTHFemPhysics
{
public:
  // constructor
  GmpCoupledThermoHydraulic(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
    const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpCoupledThermoHydraulic();
 
  enum GaussAttributeIds
  {
    T_GA_ID , //= GmpCoupledTHFemPhysics::NUM_GA_IDS,            //!< Base Id for Gauss attribute(s) used to store the calculated temperature at the current state
    TOLD_GA_ID,         
 
    // ------    NO ADDING BELOW THIS LINE
    NUM_GA_IDS,  
  };
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
  
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  
  virtual bool calcDerivedResults(bool nonLinearSolver);
  
  virtual void gravityForce(const GmElement* e, const GmVector* coord, int ip, GmVector& GFc, GmVector& N, double c);
  
  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
 
protected:
  // See comments on the base class
  virtual bool checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& cellTable, LuaTable& gaussTable);
 
  bool   _setInitialCond; 
 
  // See comments on the base class
  virtual const GmElementDof* dofMapping(GmCellType type) const;
  virtual void fillDofIndexPosition(const GmElement* e, QVector<int>& iPe, QVector<int>& iTe);
 
  virtual void fillFluidCompressibilityMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Sc, GmVector& Np, double c);
  virtual void fillThermalHydraulicCouplingMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Ze, const GmVector& Np, double a, double af, double c) const;
  virtual bool assemblyElemDampingMatrix(const GmElement* e, GmMatrix& elemC, const GmMatrix& Se, const GmMatrix& Ze,
      const GmMatrix& Xe, const  QVector<int>& iPe, const QVector<int>& iTe) const;
  virtual bool assemblyElemStiffnessMatrix(const GmElement* e, GmMatrix& elemK, const GmMatrix& He, const GmMatrix& Ie,
      const GmMatrix& Ge, const QVector<int>& iPe, const QVector<int>& iTe) const;
  virtual bool getPorePressureGaussPoint(const GmElement* e, const GmVector* coord, int ip, GmVector& Pg, const GmVector& Np, const GmVector& Pe) const;
  virtual void assemblyInternalForceVector(const GmElement* e, GmMatrix& elemFi, const GmMatrix& He, const GmMatrix& Ge, const GmMatrix& Ie,  const GmVector& Pe, const GmVector& Te, const QVector<int>& iPe, const QVector<int>& iTe) const;
 
  // Sets initial conditions
  virtual bool setInitialConditions();
 
private:
  bool  _hasElementGeneration;  
};
 
#endif