gmpTwoPhaseHydroMechanical.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_TWOPHASE_HYDROMECHANICAL_H_
#define _GEMA_PLUGIN_TWOPHASE_HYDROMECHANICAL_H_
 
 
#include "gmpCoupledHMFemPhysics.h"
#include "gmpHydroMechanical.h"
#include <gmpFemPhysicsCoupled.h>
#include <GmpHydraulicMaterialTwoPhase.h>
#include <gmpHydraulicTwoPhaseFlow.h>
#include <gmpHydraulicTwoPhaseUtils.h>
#include <gmTrace.h>
 
class GmGaussAccessor;
 
namespace tp = GmpHydraulicTwoPhaseUtils;
 
class GmpTwoPhaseHydroMechanical: public GmpHydroMechanical
{
public:
  // constructor
  GmpTwoPhaseHydroMechanical(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                     const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
 
  // destructor
  virtual ~GmpTwoPhaseHydroMechanical();
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
  virtual bool checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& cellTable, LuaTable& gaussTable);
  //virtual bool getFluidDensities(const GmElement* e, double& rhol, double& rhog, const double PlRef, const double PgRef, const double Pl, const double Pg);
  
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  virtual bool calcDerivedResults(bool nonLinearSolver);
  virtual bool updatePorosity();
  virtual double porosityRate(const GmElement* e, const double phiOld, const double dt);
  virtual double volumetricStrainRate(const GmElement* e, const double dt);
  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex, int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  // See comments on the base class
  virtual bool checkLoadedData();
  virtual bool addStateItemsToGroup(GmStateDump* state, bool fixedHint, int groupId);
 
  //enum TwoPhaseHMPhysicsAttributeIds
  //{
  //  //!< Id of the physics attributes
  //  UPDATE_DSET_APERTURE = GmpHydroMechanical::NUM_ATTRI_IDS,   //!< Flag to update the discontinuity set aperture
 
  //  // ------         NO ADDING BELOW THIS LINE
  //  NUM_ATTRI_IDS,    //!< The number of physics attribute ids above
  //};
 
protected:
  // See comments on the base class
  virtual const GmElementDof* dofMapping(GmCellType type) const;
  virtual void FillDofIndexPosition(const GmElement* e, QVector<int>& iUe, QVector<int>& iPe, QVector<int>& iPge);
  virtual bool assembleElementData(const GmElement* e, GmMatrix& J, GmVector& r,
      GmVector& ru, GmVector& rl, GmVector& rg, GmMatrix& Juu, GmMatrix& Jul, GmMatrix& Jug, GmMatrix& Jlu, GmMatrix& Jgu,
      tp::jacobianTerms& K, QVector<int>& iUe, QVector<int>& iPl, QVector<int>& iPg);
 
  enum TwoPhaseHMVarIds
  { 
    PG_ID = GmpHydroMechanical::NUM_STATEVAR_IDS,
 
    // ------     NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };
 
  // IDs for the physics attributes
  enum TwoPhaseHMPhysicsAttributeIds
  { 
      UPDATEPOROSITY = GmpHydroMechanical::NUM_ATTRI_IDS,         
          // ------         NO ADDING BELOW THIS LINE
          NUM_ATTRI_IDS,    
  };
 
  enum TwoPhaseHMNodeAttributeIds
  { 
          RHOl_NA_ID,
          RHOlOLD_NA_ID,
          RHOg_NA_ID,
          RHOgOLD_NA_ID,
          PlRef_NA_ID,
          PgRef_NA_ID,
 
      // ------     NO ADDING BELOW THIS LINE
      NUM_NA_IDS,  
  };
 
  // ! IDs for TwoPhaseHydroMechanical physics element attributes
  enum TwoPhaseHMCellAttributeIds
  { 
          POR_CA_ID,
          POROLD_CA_ID,
    Sl_CA_ID,
          // ------     NO ADDING BELOW THIS LINE
          NUM_CA_IDS,  
  };
 
};
 
#endif