FlashCO2H2O Class Reference

#include <flashco2h2o.h>

Inheritance diagram for FlashCO2H2O:

[legend]

Collaboration diagram for FlashCO2H2O:

[legend]

List of all members.

Public Member Functions
	FlashCO2H2O (OrthoMesh &mesh, const VecDouble &por, double referenceT)
	~FlashCO2H2O ()
virtual void	updateDynamicModuleData ()
virtual void	execute ()
Private Attributes
const VecDouble &	m_por
double	m_T
double	m_volCell
double	m_volMolarCO2
double	m_volMolarH2O
Static Private Attributes
static bool	bFirstTime = true

---

Detailed Description

Definition at line 11 of file flashco2h2o.h.

---

Constructor & Destructor Documentation

FlashCO2H2O::FlashCO2H2O	(	OrthoMesh &	mesh,
		const VecDouble &	por,
		double	referenceT
	)

Definition at line 9 of file flashco2h2o.cpp.

                                                                               :FlashCompositional(2,2),m_por(por),m_T(referenceT)  
{
  m_volCell = mesh.cellVolume();
  m_volMolarH2O= 18.16E-6; // m³/mol
}

FlashCO2H2O::~FlashCO2H2O

(

)

[inline]

Definition at line 23 of file flashco2h2o.h.

{}

---

Member Function Documentation

void FlashCO2H2O::execute

(

)

[virtual]

Reimplemented from FlashCompositional.

Definition at line 24 of file flashco2h2o.cpp.

{
  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)   
  }
}

void FlashCO2H2O::updateDynamicModuleData

(

)

[virtual]

Definition at line 19 of file flashco2h2o.cpp.

{
}

---

Member Data Documentation

bool FlashCO2H2O::bFirstTime = true [static, private]

Definition at line 18 of file flashco2h2o.h.

const VecDouble& FlashCO2H2O::m_por [private]

Definition at line 14 of file flashco2h2o.h.

double FlashCO2H2O::m_T [private]

Definition at line 15 of file flashco2h2o.h.

double FlashCO2H2O::m_volCell [private]

Definition at line 16 of file flashco2h2o.h.

double FlashCO2H2O::m_volMolarCO2 [private]

Definition at line 17 of file flashco2h2o.h.

double FlashCO2H2O::m_volMolarH2O [private]

Definition at line 17 of file flashco2h2o.h.

---

The documentation for this class was generated from the following files:

• flashco2h2o.h
• flashco2h2o.cpp