HenryLawObiBlunt Class Reference

#include <henrylawobiblunt.h>

Inheritance diagram for HenryLawObiBlunt:

[legend]

Collaboration diagram for HenryLawObiBlunt:

[legend]

List of all members.

Public Member Functions
	HenryLawObiBlunt (OrthoMesh &mesh, double Kh, double vn, double R, double T, double mc, double mv)
virtual	~HenryLawObiBlunt ()
virtual void	execute ()
virtual void	updateDynamicModuleData ()
virtual void	printOutput ()
Private Attributes
double	m_Kh
double	m_c
double	m_mc
double	m_mw

---

Detailed Description

The flash class that implement the HenryLaw

Definition at line 10 of file henrylawobiblunt.h.

---

Constructor & Destructor Documentation

HenryLawObiBlunt::HenryLawObiBlunt	(	OrthoMesh &	mesh,
		double	Kh,
		double	vn,
		double	R,
		double	T,
		double	mc,
		double	mv
	)

Construct the flass module with the constants used in the Henry Law

Parameters:

	Kh	Henry Law Coefficien
	vn	Apparent volume of CO2 at infinite dissolution
	R	Universal constant of gas
	T	Temperature (Kelvin)
	mc	Molar density of the CO2
	mw	Molar density of the water

Returns:

Definition at line 15 of file henrylawobiblunt.cpp.

  :FlashBase(2,2),
   m_Kh(Kh),
    m_c(-vn/(R*T)),
    m_mc(mc),
    m_mw(mv)
{
}

HenryLawObiBlunt::~HenryLawObiBlunt

(

)

[virtual]

Definition at line 85 of file henrylawobiblunt.cpp.

{

}

---

Member Function Documentation

void HenryLawObiBlunt::execute

(

)

[virtual]

Implements FlashBase.

Definition at line 27 of file henrylawobiblunt.cpp.

{
  TransportBase &trans=getTransportModule();
  //Get the necessary fields from the transport and dynamic module
  ArrayOfVecDouble &sol = trans.getSolutionAtCells();
  const VecDouble &vP = getDynamicModule().getPressureAtCells();
  
  assert(sol.size() == vP.size());


  //for each cell i in the mesh
  for (unsigned i=0;i<vP.size();i++)
  {
    double Cd = sol(i,FluxForCO2Inj::_Sd); //volume of dissolved CO2 per pore volume
    double Sc = sol(i,FluxForCO2Inj::_Sc); //Saturation of CO2
    double P = vP(i); //pression in the cell

    double Xd = (P/m_Kh)*exp(m_c*P); //Maximum molar fraction of dissolved co2
    if (Xd > 1.0)
    {
      printf("Houston we have a problem Xd == %g\n",Xd);
      throw new Exception("Finishing");
    }
    if (Sc >= 1.0)
      continue;
    
      
    //with Xd, we calculate the max amount of CO2 that can dissolve in water per pore volume
    double NDco2 = m_mw*(1.0-Sc - Cd)*Xd/(1.0-Xd);   


    double NTco2 = (Cd+Sc)*m_mc; //total moles of CO2 in the mixture
    double CdN,ScN;
    
    //Now compute the dissolution
    if (NTco2 >= NDco2)
    {
      //Ok, we have enough CO2 to dissolve into the water
      //dissolve it all.
      ScN = Sc + Cd - NDco2/m_mc; 
      CdN = NDco2/m_mc;
    }
    else
    {
      ScN = 0.0;
      CdN = (Cd +Sc);
    }
    sol(i,FluxForCO2Inj::_Sd)=CdN; //volume of dissolved CO2 per pore volume
    sol(i,FluxForCO2Inj::_Sc)=ScN; //Saturation of CO2
  } 
}

virtual void HenryLawObiBlunt::printOutput

(

)

[inline, virtual]

Implements FlashBase.

Definition at line 21 of file henrylawobiblunt.h.

{}

void HenryLawObiBlunt::updateDynamicModuleData

(

)

[virtual]

Definition at line 79 of file henrylawobiblunt.cpp.

{
}

---

Member Data Documentation

double HenryLawObiBlunt::m_c [private]

Definition at line 13 of file henrylawobiblunt.h.

double HenryLawObiBlunt::m_Kh [private]

Definition at line 13 of file henrylawobiblunt.h.

double HenryLawObiBlunt::m_mc [private]

Definition at line 13 of file henrylawobiblunt.h.

double HenryLawObiBlunt::m_mw [private]

Definition at line 13 of file henrylawobiblunt.h.

---

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

• henrylawobiblunt.h
• henrylawobiblunt.cpp