gmpHydraulicMultiPhaseFlow.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_MULTIPHASEFLOW_H_
#define _GEMA_PLUGIN_HYDRAULIC_MULTIPHASEFLOW_H_
 
#include <gmpFemPhysicsCommon.h>
#include <gmpHydraulic.h>
 
#include <gmMathUtils.h>
#include <gmTrace.h>
 
class GmGaussAccessor;
 
class GMP_HYDRAULIC_PHYSICS_API_EXPORT GmpHydraulicMultiPhaseFlow : public GmpHydraulic
{
public:
  // constructor
  GmpHydraulicMultiPhaseFlow(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                             const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);
  // destructor
  virtual ~GmpHydraulicMultiPhaseFlow();
 
 
  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
 
  virtual GmpFemPhysics::FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  virtual bool fixedNodalForcesBc(QVector<int>& nodes, QVector<int>&dof, QVector<double>& values) const;
 
  virtual bool fixedNodalDofsBc  (QVector<int>& nodes, QVector<int>&dof, QVector<double>& values, bool* constantValues) const;
 
  // See comments on the base class
  virtual bool calcDerivedResults(bool nonLinearSolver);
 
  virtual void fillStorageCompMatrix(const GmElement* e, const GmVector* coord, int ip, 
          GmMatrix& Sww, GmMatrix& Swc, GmMatrix& Scw, GmMatrix& Scc, const GmVector& N, double c) const;
 
  virtual void fillFlowMatrix(const GmElement* e, const GmVector* coord, int ip, GmMatrix& Kc, GmMatrix& Kwc, 
          GmMatrix& Kcc, const GmMatrix& Bp, double c) const;
 
  virtual void fillElementFluidPorePressure(const GmElement* e, GmVector& Pce) const;
 
  virtual void FillDofIndexPosition(const GmElement* e, QVector<int>& iPe, QVector<int>& iPfe) const;
 
  virtual void fillDampingMatrix(const GmElement* e, GmMatrix& C, const GmMatrix& Sww, const GmMatrix& Swc, const GmMatrix& Scw, const GmMatrix& Scc, const QVector<int>& iPw, const QVector<int>& iPc) const;
 
  virtual void fillStiffnessMatrix(const GmElement* e, GmMatrix& K, const GmMatrix& Kww, const GmMatrix& Kwc, const GmMatrix& Kcc, const QVector<int>& iPw, const QVector<int>& iPc) const;
 
  virtual void fillInternalForceVector(const GmElement* e, GmVector& Fi, const GmMatrix& Kww, const GmMatrix& Kwc, const GmMatrix& Kcc, const GmVector& Pw, const GmVector& Pc, const QVector<int>& iPw, const QVector<int>& iPc) const;
  
  virtual void gravitySupplyFluid(const GmElement* e, const GmVector* coord, int ip, GmVector& gFp, GmVector& gFpc, const GmMatrix& Bp, double c);
 
protected:
  // See comments on the base class
  virtual const GmElementDof* dofMapping(GmCellType type) const;
 
  enum MultiPhaseStateVarIds
  {
    PC_ID = GmpHydraulic::NUM_STATEVAR_IDS, 
 
    // --- NO ADDING BELOW THIS LINE
    NUM_STATEVAR_IDS,  
  };
 
  enum MultiPhaseBoundaryConditionIds
  {
    FIXED_PC_BC_ID = GmpHydraulic::NUM_BC_IDS, 
    FIXED_QC_BC_ID,                            
 
    // --- NO ADDING BELOW THIS LINE
    NUM_BC_IDS,  
  };
 
  enum MultiPhaseBoundaryConditionValueIds
  {
    BC_PC_ID = GmpHydraulic::NUM_BCV_IDS,  
    BC_QWC_ID,                             
 
    // --- NO ADDING BELOW THIS LINE
    NUM_BCV_IDS,  
  };
 
};
 
#endif