gmpHydraulicTwoPhaseUtils.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_TWOPHASEUTILS_H_
#define _GEMA_PLUGIN_HYDRAULIC_TWOPHASEUTILS_H_
 
#include<gmElement.h>
#include <gmpHydraulicMaterialTwoPhase.h>
#include <gmVector.h> 
#include <gmMatrix.h>
 
namespace GmpHydraulicTwoPhaseUtils
{
  // structure for nodal variables
  // Saturation terms
  struct saturationTerms
  {
    double   Slm;        // Mean liquid saturation
    double   Sgm;        // Mean gas saturation
    GmVector Sl;         // Liquid saturation degree vector at the current time step
    GmVector SlOld;      // Liquid saturation degree vector at the previous time step
    GmVector Sg;         // Gas saturation degree vector
    GmVector SgOld;      // Gas saturation degree vector at previous time step
    GmVector dSldPc;     // Derivative of liquid saturation degree in terms of capillary pressure
    GmVector dSlmdPl;    // Derivative of mean liquid saturation in terms of nodal liquid pressure
    GmVector dSlmdPg;    // Derivative of mean liquid saturation in terms of nodal gas pressure
 
    explicit saturationTerms(size_t np)
      :Sl(np), SlOld(np), Sg(np), SgOld(np),
      dSldPc(np), dSlmdPl(np), dSlmdPg(np) {}
  };
 
  // Density terms
  struct densityTerms
  {
    double   rholm;      // Mean liquid density
    double   rhogm;      // Mean gas density
    GmVector rhol;       // Liquid density at the nodes vector
    GmVector rhog;       // Gas density at the nodes vector
    GmVector rholOld;    // Liquid density vector at previous time step
    GmVector rhogOld;    // Gas density vector at previous time step
    GmVector drholmdPl;  // Derivative of mean liquid density in terms of nodal liquid pressure
    GmVector drhogmdPg;  // Derivative of mean gas density in terms of nodal gas pressure
    explicit densityTerms(size_t np)
      : rhol(np), rhog(np), rholOld(np),
      rhogOld(np), drholmdPl(np), drhogmdPg(np) {}
  };
 
  // Viscosity terms
  struct viscosityTerms
  {
    double mul = 0.0;    // Dynamic viscosity of fluid
    double mug = 0.0;    // Dynamic viscosity of gas
    GmVector dmulmdPl;   // Derivative of mean liquid viscosity in terms of nodal liquid pressure
    GmVector dmugmdPg;   // Derivative of mean gas viscosity in terms of nodal gas pressure
 
    explicit viscosityTerms(size_t np)
      : dmulmdPl(np), dmugmdPg(np) {}
  };
 
  // Mass terms
  struct massTerms
  {
    GmVector ml;         // Liquid mass vector at the nodes: ml = rhol * Sl
    GmVector mg;         // Gas mass vector at the nodes: mg = rhog * Sg
    GmVector Dml;        // Liquid mass increment vector
    GmVector Dmg;        // Gas mass increment vector
    GmVector dmldPl;     // Derivative of liquid mass (rhol*Sl) in terms of nodal liquid pressure at the nodes
    GmVector dmldPg;     // Derivative of liquid mass (rhol*Sl) in terms of nodal gas pressure at the nodes  
    GmVector dmgdPl;     // Derivative of gas mass (rhog*Sg) in terms of nodal liquid pressure at the nodes
    GmVector dmgdPg;     // Derivative of gas mass (rhog*Sg) in terms of nodal gas pressure at the nodes
 
    explicit massTerms(size_t np)
      : ml(np), mg(np), Dml(np), Dmg(np), dmldPl(np),
      dmldPg(np), dmgdPl(np), dmgdPg(np) {}
  };
 
  // Pressure terms
  struct pressureTerms
  {
    GmVector Pl;                        // Liquid phase pore pressure
    GmVector Pg;                        // Gas phase pore pressure
    GmVector PlOld;                     // Liquid phase pore pressure at previous time
    GmVector PgOld;                     // Gas phase pore pressure at previous time
    GmVector Ps;                        // Nodal effective pore pressure vector: Ps = Pl * Sl + Pg * Sg
    GmVector PsOld;                     // Nodal effective pore pressure vector at previous time: PsOld = PlOld * SlOld + PgOld * SgOld
    GmVector Pc;                        // Nodal capillarity pressure: Pc = Pg - Pl
    GmVector PcOld;                     // Nodal capillarity pressure at previous time: PcOld = PgOld - PlOld
    const double dPcdPl = - 1.0;        // Derivative of capillary pressure in terms of liquid pressure
    const double dPcdPg = 1.0;              // Derivative of capillary pressure in terms of gas pressure
 
    explicit pressureTerms(size_t np)
      : Pl(np), Pg(np), PlOld(np), PgOld(np),
      Ps(np), PsOld(np), Pc(np), PcOld(np) {}
  };
 
  // Relative permeability terms
  struct permeabilityTerms
  {
    double   krl;        // Liquid relative permeability
    double   krg;        // Gas relative permeability
    double   dkrlmdSl;   // Derivative of liquid relative permeability in terms of the mean liquid saturation
    double   dkrgmdSl;   // Derivative of gas relative permeability in terms of the mean liquid saturation
    GmVector dkrlmdPl;   // Derivative of liquid relative permeability in terms of nodal liquid pressure
    GmVector dkrlmdPg;   // Derivative of liquid relative permeability in terms of nodal gas pressure
    GmVector dkrgmdPl;   // Derivative of gas relative permeability in terms of nodal liquid pressure
    GmVector dkrgmdPg;   // Derivative of gas relative permeability in terms of nodal gas pressure
 
    explicit permeabilityTerms(size_t np)
      : dkrlmdPl(np), dkrlmdPg(np), dkrgmdPl(np), dkrgmdPg(np) {}
  };
  // Mass mobility terms
  struct massMobilityTerms
  {
    double massMobilityL; // mass mobility of liquid phase
    double massMobilityG; // mass mobility of gas phase
    // Mass mobilities derivatives. The densities and the fluid viscosities are considered a function of the corresponding phase only. 
    GmVector dmassMobilityLdPl; // Derivative of liquid mass mobility in terms of nodal liquid paressure
    GmVector dmassMobilityLdPg; // Derivative of liquid mass mobility in terms of nodal gas paressure
    GmVector dmassMobilityGdPl; // Derivative of gas mass mobility in terms of nodal liquid paressure
    GmVector dmassMobilityGdPg; // Derivative of gas mass mobility in terms of nodal gas paressure
 
    explicit massMobilityTerms(size_t np)
      : dmassMobilityLdPl(np), dmassMobilityLdPg(np), dmassMobilityGdPl(np), dmassMobilityGdPg(np) {}
  };
 
  // Jacobian matrix terms
  struct jacobianTerms
  {
    GmMatrix dfildpl;       // Derivative of the liquid phase pressure internal force vector wrt to the liquid pore-pressures
    GmMatrix dfildpg;       // Derivative of the liquid phase pressure internal force vector wrt to the gas pore-pressures
    GmMatrix dfigdpl;       // Derivative of the gas phase pressure internal force vector wrt to the liquid pore-pressures
    GmMatrix dfigdpg;       // Derivative of the gas phase pressure internal force vector wrt to the gas pore-pressures
 
    explicit jacobianTerms(size_t np)
      : dfildpl(np, np), dfildpg(np, np), dfigdpl(np, np), dfigdpg(np, np) {}
  };
 
  // Methods used by two phase formulations
  // Evaluates nodal variables
  //bool evaluateNodalVariables(const GmElement* e, saturationTerms& Sat, densityTerms& Den,
  //  viscosityTerms& Mu, massTerms& Ms, pressureTerms& Pres);
 
  // Evaluates relative permeabilities
  GMP_HYDRAULIC_PHYSICS_API_EXPORT bool evaluateRelativePermeability(const GmElement* e, const GmpHydraulicMaterialTwoPhase* hmat, saturationTerms& Sat, permeabilityTerms& kr, int np);
 
  // Evaluates the hydraulic gravity vector for each ip
  GMP_HYDRAULIC_PHYSICS_API_EXPORT void evaluateHydraulicGravityVector(const GmElement* e, GmMatrix& Kh, GmVector& fgrav, GmMatrix& Bp,
    const GmVector& grav, double c);
 
  GMP_HYDRAULIC_PHYSICS_API_EXPORT bool evaluateAdvectiveForceAndJacobian(const GmElement* e, GmVector& Fil, GmVector& Fig, jacobianTerms& K,
    GmMatrix& H, massMobilityTerms& massMob, densityTerms& den, viscosityTerms& mu, permeabilityTerms& kr, pressureTerms& pres);
 
  GMP_HYDRAULIC_PHYSICS_API_EXPORT bool evaluateMassStorageForceAndJacobian(const GmElement* e, GmVector& Fil, GmVector& Fig, jacobianTerms& K,
    saturationTerms& sat, massTerms& mass, pressureTerms& pres, GmVector& iN, double Cp, double phi, double dt, int np);
 
  GMP_HYDRAULIC_PHYSICS_API_EXPORT bool evaluateGravityForceAndJacobian(const GmElement* e, GmVector& Fil, GmVector& Fig, jacobianTerms& K,
    GmVector& fgravp, massMobilityTerms& mMb, saturationTerms& sat, densityTerms& den, double phi);
 
}
#endif