flashco2h2o.cpp

Go to the documentation of this file.

#include "flashco2h2o.h"
#include "flashessential.h"
#include "flash.h"
#include "fluxforco2inj.h"
#include "arrayofvecdouble.h"
bool FlashCO2H2O::bFirstTime=true;

// reference temperature must be in Kelvins degree
FlashCO2H2O::FlashCO2H2O(OrthoMesh &mesh,const VecDouble &por, double referenceT):FlashCompositional(2,2),m_por(por),m_T(referenceT)  
{
  m_volCell = mesh.cellVolume();
  m_volMolarH2O= 18.16E-6; // m³/mol
}





void FlashCO2H2O::updateDynamicModuleData() 
{
}


 void FlashCO2H2O::execute() 
{
  TransportBase &trans=getTransportModule();
  const VecDouble& vP = getDynamicModule().getPressureAtCells(); 
  if (bFirstTime){
    double referencePressure=vP.mean_value()*1e-5;
    m_volMolarCO2= 1.0E-6 * ComputeVolumeCO2(m_T, referencePressure, Flash::SUPERCRITICAL); 
    printf("AllanFlash Using Co2 vol Molar %g, Pressure: %g \n",m_volMolarCO2,referencePressure);
    bFirstTime=0;
    //    m_volMolarCO2=61.9578e-6;
  }

  //Get the necessary fields from the transport and dynamic module
  
  ArrayOfVecDouble &sol=trans.getSolutionAtCells();

  assert(sol.size() == vP.size());

  //for each cell i in the mesh
  for (unsigned i=0;i<vP.size();i++)
  {
    double Sd = sol(i,FluxForCO2Inj::_Sd); //volume of dissolved CO2 per pore volume
    double Sc = sol(i,FluxForCO2Inj::_Sc); //Saturation of CO2
    double P =  vP(i)*1E-5; //pression in the cell
    double por = m_por(i);
    assert(P > 75);

    double Vco2 = Sc*por*m_volCell;
    double Vd   = Sd*por*m_volCell;
    double Vw   = (1-Sc-Sd)*por*m_volCell;
 
    double nH2O  = Vw/m_volMolarH2O;
    double nCO2  = Vco2/m_volMolarCO2;
    double nCO2d = Vd/m_volMolarCO2;

    double nT = nH2O + nCO2 + nCO2d;

    double zH2O = nH2O/nT;
    double zCO2 = 1.0 - zH2O;
    double zNaCl = 0.0; // assuming no salt in the system.

    // Compute the flash equilibrium.
    Flash::BrineCO2Data brineco2Data;

    Flash::ProgressiveFlashBrineCO2(brineco2Data, m_T, P, zNaCl, zH2O, zCO2);

    // Computing the new saturations in the cell.
    sol(i,FluxForCO2Inj::_Sc) = brineco2Data.yCO2 * brineco2Data.betaC * m_volMolarCO2 * nT / (m_volCell * por); // (yCO2 * betaC * vCO2 * nT)/(Vcell * porosity)
    sol(i,FluxForCO2Inj::_Sd) = brineco2Data.xCO2 * brineco2Data.betaA * m_volMolarCO2 * nT / (m_volCell * por); // (xCO2 * betaA * vCO2 * nT)/(Vcell * porosity)   
  }
}

---