facefluxco2brine.h

Go to the documentation of this file.

#ifndef  _MY_FaceFluxCO2Brine_
#define  _MY_FaceFluxCO2Brine_
#include "facefluxfunction.h"
#include "flashcompositional.h"
#include "fchaventmobility.h"
#include "fbucleyleverettgravitymob.h"
#include "dfchaventmobility.h"
#include "dfbucleyleverettgravitymob.h"


class FaceFluxCO2Brine : public FaceFluxFunction
{
   private:
  FlashCompositional &m_flash;
  Function1D &m_mob,&m_gravMob,&m_Dmob,&m_DgravMob;
  /* FChaventMobility  m_mob; */
  /* FBucleyLeverettGravityMob m_gravMob; */
  /* DFChaventMobility  m_Dmob; */
  /* DFBucleyLeverettGravityMob m_DgravMob; */

  VecDouble m_vNormalVel,;
  const VecDouble *m_pP;
  const VecDouble &m_K;
  VecDouble m_densities1,m_densities2;
  double m_Sw1,m_Sw2,m_cell1,m_cell2;
  double m_grav;
  OrthoMesh &m_mesh;
 protected:
    virtual void fluxAtFace(VecDouble &vFlux, const FaceInfo &face,const VecDouble &Q1,const VecDouble &Q2);
    virtual double maxLocalCharVelocity(const FaceInfo &face,const VecDouble &Q1, const VecDouble &Q2);

 public:
    FaceFluxCO2Brine(OrthoMesh &mesh,FlashCompositional &flash,const VecDouble &K,GeneralFunctionInterface &fMob,GeneralFunctionInterface &fDMob,GeneralFunctionInterface &fGrav,GeneralFunctionInterface &fDGrav,double grav); 

  virtual ~FaceFluxCO2Brine(){}
  


  

  virtual const VecDouble& getPressure(){assert(m_pP);return *m_pP;}
  virtual const VecDouble& getK(){return m_K;}
  virtual void  updateDynamicData(DynamicBase &dynMod);
  virtual void  maxGlobalCharVelocity(double vel[3]);
  virtual void transformToComponentMoles(VecDouble &v1);
};

#endif

---