dummyflash.cpp

Go to the documentation of this file.

#include "dummyflash.h"






DummyFlash::DummyFlash(unsigned nPhases,unsigned nComps,VecDouble &viscs,VecDouble &molarMass) 
   :FlashCompositional(nPhases,nComps)
{
  assert(molarMass.size() == nComps);
  assert(viscs.size() == nPhases);
  m_compMass=molarMass;
  m_viscs=viscs;
}

 const VecDouble& DummyFlash::getComponentsMolarMass() const 
{
  return m_compMass;
}




 void DummyFlash::getPhasesVolume(double P, const FlashData &data, VecDouble &phasesVol) 
{
  data.getPhasesTotalMoles(phasesVol);
  //  assert(phasesVol(0)==data.getMoles(0,0));
  //assert(phasesVol(1)==data.getMoles(1,1));

}




 void DummyFlash::getPhasesViscosities(double P,const FlashData &data,VecDouble &visc) 
{
  visc=m_viscs;
}

void DummyFlash::flash(double P, const VecDouble &compTotalMoles, FlashData &data) 
{
  assert(data.getNPhases() == this->numPhases());
  assert(data.getNComponents() == this->numComponents());
  for (unsigned i=0;i<numPhases();i++)
    for (unsigned j=0;j<numComponents();j++)
      data.setMoles(i,j,0.0);

  for (unsigned i=0;i<numPhases();i++)
  {
    data.setMoles(i,i,compTotalMoles(i));
  }
}


 DummyFlash::~DummyFlash() 
{

}


 double DummyFlash::getFluidCompressibility(double P, FlashData &data) 
{
  return 0.0;
}


 void DummyFlash::getTotalVolumeDerivatives(double P,FlashData &data,VecDouble &dv_dm) 
{
  dv_dm=1;
}

---