FaceFluxSimpleBlackOilMass Class Reference

#include <facefluxsimpleblackoilmass.h>

Inheritance diagram for FaceFluxSimpleBlackOilMass:

[legend]

Collaboration diagram for FaceFluxSimpleBlackOilMass:

[legend]

List of all members.

Public Member Functions
	FaceFluxSimpleBlackOilMass (OrthoMesh &mesh, FlashCompositional &flash)
	~FaceFluxSimpleBlackOilMass ()
virtual void	fluxAtFace (VecDouble &vFlux, OrthoMesh::Face_It &faceIt, int face, int cell1, int cell2, const VecDouble &Q1, const VecDouble &Q2)
virtual void	updateDynamicData (DynamicBase &dynMod)
virtual void	maxGlobalCharVelocity (double vel[3])
Private Attributes
VecDouble	m_vNormalMassFlux

---

Detailed Description

Definition at line 9 of file facefluxsimpleblackoilmass.h.

---

Constructor & Destructor Documentation

FaceFluxSimpleBlackOilMass::FaceFluxSimpleBlackOilMass	(	OrthoMesh &	mesh,
		FlashCompositional &	flash
	)

Definition at line 10 of file facefluxsimpleblackoilmass.cpp.

   :FaceFluxSimpleBlackOil(mesh,flash),m_vNormalMassFlux(mesh.numFaces())
{
}

FaceFluxSimpleBlackOilMass::~FaceFluxSimpleBlackOilMass

(

)

Definition at line 5 of file facefluxsimpleblackoilmass.cpp.

{

}

---

Member Function Documentation

void FaceFluxSimpleBlackOilMass::fluxAtFace	(	VecDouble &	vFlux,
		OrthoMesh::Face_It &	faceIt,
		int	face,
		int	cell1,
		int	cell2,
		const VecDouble &	Q1,
		const VecDouble &	Q2
	)			[virtual]

(phasesVolume(0)+phasesVolume(1));

(phasesVolume(0)+phasesVolume(1));

Definition at line 24 of file facefluxsimpleblackoilmass.cpp.

{
  assert(vFlux.size()==2);
  assert((unsigned) face < m_vNormalVel.size());
  static FlashData data1(2,2),data2(2,2);
  static VecDouble phasesVolume1(2);
  static VecDouble phasesVolume2(2);
  static VecDouble vFlux1(2);
  static VecDouble vFlux2(2);
  const VecDouble &P = getPressure();
  m_bUnsaturated=false;

  m_flash.flash(P(cell1),Q1,data1);
  m_flash.flash(P(cell2),Q2,data2);

  m_flash.getPhasesViscosities(P(cell1),data1,m_vViscosities1);
  m_flash.getPhasesViscosities(P(cell2),data2,m_vViscosities2);
  
  m_flash.getPhasesVolume(P(cell1),data1,phasesVolume1);
  m_Sw1=phasesVolume1(0);

  
  m_mob.setParameters(m_vViscosities1);
  double mobO=m_mob(m_Sw1);

  vFlux(0)=data1.getMoles(0,0)*mobO/(data1.getMoles(0,0)*mobO + data1.getMoles(0,1)*mobO + _div(data1.getMoles(1,1),phasesVolume1(1))*(1-mobO)*phasesVolume1(0));
  vFlux(1)=1-vFlux(0);
  


  m_flash.getPhasesVolume(P(cell2),data2,phasesVolume2);
  m_Sw2=phasesVolume2(0);
  mobO=m_mob(m_Sw2);
  m_mob.setParameters(m_vViscosities2);


  vFlux2(0)=data2.getMoles(0,0)*mobO/(data2.getMoles(0,0)*mobO + data2.getMoles(0,1)*mobO + _div(data2.getMoles(1,1),phasesVolume2(1))*(1-mobO)*phasesVolume2(0));
  vFlux2(1)=1-vFlux2(0);

  
  vFlux+=vFlux2;
  vFlux*= m_vNormalMassFlux(face)/2.0;

  if (data1.getMoles(1,1) == 0.0 || data2.getMoles(1,1)==0.0)
  {
    //printf("Unsaturated at %g\n",faceIt->barycenter()[0]);
    m_bUnsaturated=true;
  }


//Store the cells just to check at getLocalCharVelocity
  m_cell1=cell1;
  m_cell2=cell2;

  /*
  VecDouble vF(2);
  FaceFluxSimpleBlackOil::fluxAtFace(vF, faceIt, face, cell1, cell2,Q1,Q2);
  printf("vFlux(0)=%g %g\n",vFlux(0),vF(0));
  ********/
}

void FaceFluxSimpleBlackOilMass::maxGlobalCharVelocity

(

double

vel[3]

)

[virtual]

Reimplemented from FaceFluxSimpleBlackOil.

Definition at line 86 of file facefluxsimpleblackoilmass.cpp.

{
  double maxVel=NumericMethods::vectorMaxAbsValue(m_vNormalMassFlux);
  double maxVel2=NumericMethods::vectorMaxAbsValue(m_vNormalVel);
  printf("Max Vel: %g %g",maxVel,maxVel2);
  double P=NumericMethods::vectorMinAbsValue(getPressure());
  static FlashData data(2,2);
  data.allocateOwnMemory();
  double fact=pow((P/10000 + 1.0)/(0.05*P+ 1),2.0)*(0.06*P + 1);

  m_flash.getPhasesViscosities(P,data,m_vViscosities1);
  m_Dmob.setParameters(m_vViscosities1);
  m_Dmob.updateInflectionPoint();
  vel[X]= fact*maxVel*m_Dmob.getMaxNorm(0,1);
  vel[Y]=0;
  vel[Z]=0;
  return;

}

void FaceFluxSimpleBlackOilMass::updateDynamicData

(

DynamicBase &

dynMod

)

[virtual]

Reimplemented from FaceFluxSimpleBlackOil.

Definition at line 15 of file facefluxsimpleblackoilmass.cpp.

{
  FaceFluxSimpleBlackOil::updateDynamicData(dynMod);
  dynMod.getNormalMassFluxAtFaces(m_vNormalMassFlux);
  //  m_vNormalMassFlux.print(std::cout,8);
}

---

Member Data Documentation

VecDouble FaceFluxSimpleBlackOilMass::m_vNormalMassFlux [private]

Definition at line 12 of file facefluxsimpleblackoilmass.h.

---

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

• facefluxsimpleblackoilmass.h
• facefluxsimpleblackoilmass.cpp