flash.h

这是本人两年前兼职为某个公司做的石油钻进设计软件

#if !defined(AFX_FLASH_H)
#define AFX_FLASH_H
  
#if _MSC_VER >= 1000
#pragma once
#endif // _MSC_VER >= 1000

#define R         82.06
#define PI        3.1415926
#define iComNum   17
#define TOTALCOMNUM                     18

class  EOSFlashClass
{
     public:
		  EOSFlashClass( );
		  ~EOSFlashClass( );

	 public:
		  int    m_iNum;           //组分数
		  float  *m_pfMi;          //各组分的分子量
		  float  **m_pfKij;        //各组分的相互作用系数
		  float  *m_pfTci;         //各组分的临界温度
		  float  *m_pfPci;         //各组分的临界压力
		  float  *m_pfWi;          //各组分的偏心因子
		  float  *m_pfVci;         //各组分的临界体积
		  
		  float  *m_pfZci;

		  float  m_fAm;            //状态方程系数
		  float  m_fBm;            //状态方程系数
		  float  m_fCm;            //状态方程系数
		  float  m_fwm;            //状态方程系数
		  float  m_fam;            //状态方程系数
		  float  m_fbm;            //状态方程系数
		  float  m_fdamdT;
		  float  m_f2damdT;
		  float  m_fdZdT;
		  float  m_fdZdP;

		  float  *m_pffiL;         //油相逸度
		  float  *m_pffiV;         //气相逸度
		  float  *m_pffi;          //临时存储相逸度
		  float  *m_pfKi;          //平衡常数
		  float  *m_pfAlfai;       //状态方程系数
		  float  *m_pfai;          //状态方程系数
		  float  *m_pfbi;          //状态方程系数
		  float  *m_pfmi;

		  float   *m_pfZi;          //混合物组成
		  float   *m_pfXi;          //油相组成
		  float   *m_pfYi;          //气相组成

		  float   *X1;            //临时存储油相组成
		  float   *X2;            //临时存储油相组成
		  float   *X3;            //临时存储油相组成
		  float   fP1;            //对应X1的压力
		  float   fP2;            //对应X2的压力
		  float   fP3;            //对应X3的压力

		  float   m_fTemp;          //温度
		  float   m_fDewPressure;
		  float   m_fPress;         //压力
          float   m_fPStep;
		  float   m_fV;             //气相摩尔分数
		  float   m_fSg;            //气相饱和度
		  float   m_fDo;            //油相密度
		  float   m_fDg;            //气相密度
		  float   m_fZo;            //油相偏差因子
		  float   m_fZg;            //气相偏差因子
		  float   m_fUo;            //油相粘度
		  float   m_fUg;            //气相粘度
		  float   m_fMo;            //油相分子量
		  float   m_fMg;            //气相分子量

		  float   m_fdZodP;         //Z因子对压力的偏导
		  float   m_fdZgdP;         //

		  int     DS;
		  
		  float*  m_pfPMP_Press;
		  float*  m_pfPMP_MP;
		  int     m_iPMP_Num;

		  float*  m_pfV;             //气相摩尔分数
		  float*  m_pfSg;            //气相饱和度
		  float*  m_pfDo;            //油相密度
		  float*  m_pfDg;            //气相密度
		  float*  m_pfZo;            //油相偏差因子
		  float*  m_pfZg;            //气相偏差因子
		  float*  m_pfUo;            //油相粘度
		  float*  m_pfUg;            //气相粘度
		  float*  m_pfMo;            //油相分子量
		  float*  m_pfMg;            //气相分子量

	 public:
		  int InputCompFixData( float *pfZi );
		  int InitKi( float fT, float fP );
		  int CalEOSCoeff_mi_ai_bi( );
		  int CalEOSCoeff_Alfa_am_bm_Am_Bm( float fP, float fT, float *pfEi );
		  int Calfi( float fZ, float fP, float *pfEi, float *pffi );
		  int CaldamdT( float fT, float *pfEi );
		  int Cald2amdT( float fT, float *pfEi );
		  int CaldZdT( float fP, float fT, float fZ, float *pfEi );
		  int CaldZdP( float fT, float fZ );

		  void    SetComposition( float *_pfZi );
		  int     InitV( float *pfKi );
		  int     InitXi_Yi( float *pfKi );
		  int     CalZ_Factor( float fAm, float fBm, float &_fInitZ );

		  void    CalMolecularWeight(  );
		  void    CalQualityDensity(  );
		  void    CalMoleDensity(  );
		  void    CalVolumeFraction(  );
		  void    CalGasSatuation(  );
		  int     CalViscosity( float *pfEi, float fDen, float &fViscosity );
		  int     CalZV( );
		  int     CalDewPressure(  );
		  int     ObjectEquation(  );
		  int     Step( float *PP, float **H );
		  int     BroyDen(  );
		  int     InitNextFlashSpot(  );

		  void    SetPMP_Pressure( );
		  int     FlashCalculate( float* pfComposite, float fTemp, float fPress );
};

/*
	 public:
		  void    SetNewTimeStepValue( );
		  int     CalThermalDynamicFun(  );

	 public:
		  int     LL;
		  int     CL;
		  int     KK;
		  int     iRes;

		  int     iIterNum;       //进行的迭代次数

		  int     &iNum;          //组分数
		  int     iMaxIterNum;    //控制迭代次数
		  int     iResFlag;       //返回标志


		  float   fOldPress;
		  float   fOldV;             //气相摩尔分数
		  float   fOldSg;            //气相饱和度
		  float   fOldDo;            //油相密度
		  float   fOldDg;            //气相密度
		  float   fOldZo;            //油相偏差因子
		  float   fOldZg;            //气相偏差因子
		  float   fOldUo;            //油相粘度
		  float   fOldUg;            //气相粘度

		  float   fOldKro;           //油相相对渗透率
		  float   fOldKrg;           //气相相对渗透率

		  float   fOldMo;            //油相分子量
		  float   fOldMg;            //气相分子量
		  float   fOldCpo;           //油相热容
		  float   fOldCpg;           //气相热容
		  float   fOldCjo;           //油相焦-汤系数
		  float   fOldCjg;           //气相焦-汤系数
		  float   fOldHo;            //油相焓
		  float   fOldHg;            //气相焓

		  float   fOldCp;            //混合物热容
		  float   fOldCj;            //混合物焦-汤系数
		  float   fOldH;             //混合物焓
		  float   fOldU;             //混合物粘度

};*/
#endif