计算流体相平衡

计算流体相平衡，仅供大家参考，共同学习，共同进步。

#include
#include
#include
#include 
double antoine（double）;	//运用Antoine Equation方程计算
double yidu（doubledoubledouble）;//运用PR方程计算逸度
double deltayidu（doubledouble）;//计算气相和液相逸度差
double jiexi（int doubledoubledouble）;//运用解析法解关于Z的三次方程

main（）
{
	double tc=562.2;//苯的临近温度
	double *ta=NULL;//动态创建三个数组，分别存储温度和饱和蒸汽压及迭代压力
	double *pa=NULL;
	double *pe=NULL;
	double tb;
	double t=353.35; //起始温度
	int count=0;
	int i;
	int j;
	tb=t;
	while（1）//温度的动态划分
	{
		count=count+1;
		if（tb<=tc-30）
			tb=tb+5;
		else
			tb=tb+0.1;
		if（tb>tc）
			break;
	}
	ta=new double[count];
	pa=new double[count];
	pe=new double[count];
	ta[0]=t;
	pa[0]=antoine（t）;
	pa[0]=pa[0]*0.001315789;
	for（j=1;j<=count-1;j++）//运用Antoine Equation方程计算饱和蒸汽压
	{
		

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4650  2020-04-27 18:28  计算流体相平衡.c