C语言在闪蒸过程计算中的应用.docVIP

C语言在闪蒸过程计算中的应用 常兴 青海大学盐湖系（青海西宁 810016） Application of Flash Process Calculation by C Language Chang Xing Department of Salt Lake, QingHai University（QingHai XiNing 810016） Abstract：This Paper analysis the calculation of flash process through material balance ．By a example to discuss the basic application of Newton iterative method and C language program in flash process calculation ． Keywords：Flash，Chemical Calculation，C Language 摘要：本文从物料衡算出发，分析了闪蒸过程计算，通过一个计算实例探讨了Newton 迭代求解法以及C语言编程方法在闪蒸过程计算的基本应用。 关键词：闪蒸 化工计算 C语言 中图分类号： 文献标识码： 1问题的提出 闪蒸是连续的单级蒸馏过程。该过程主要用于化工过程混合物的初步分离[1]。 等温闪蒸过程如图1所示，流量F（摩尔流量，下同）、组成为zi（摩尔分数，下同）的混合物经换热器进入 气液分离罐，在罐中经闪蒸分离出互成平衡的 气液两相。汽相流量为V，组成为yi；液相流 量为L，组成为xi。分离罐压力为P，温度为T。 对于等温闪蒸过程的问题是：已知混合 物的总组成zi闪蒸温度T和压力P，要求计算 闪蒸后汽相所占分率e，及平衡的汽液相组成 yi和xi。这个问题必须通过迭代法求解。如果采用 计算机求解则比较方便，程序的通用性强。可以把 同学们从繁重的试算工作中解脱出来。 图1 等温闪蒸 2问题的求解 2．1工作方程： 对于如图1示的闪蒸过程，用物料衡算和相平衡方程得到： …………………………………………（1） 令汽化分率e为 …………………………（2） 由（1）得到： …………………………………………（3） 由及可写出： ………………………………………（4） 当温度、压力及混合物组成给定之后，方程（4）中只有一个未知数e，对于二元物系可直接求解，但对于组分数C≥3的体系用一般手算的方法已显得无能为力，随着计算机技术的进步以及算法语言的不断出现，求解这个问题显得非常容易。 2．2求解过程 对于（4）式采用Newton 迭代法求解。令： …………………………………………（5） 若目标函数为。此时求解得到的e即为问题的解。计算框图如图2 输入 设e，初值 FOR i=1 TO C ; 圆整， 调整e 直至收敛 图2 等温闪蒸计算框图 3计算算例 3.1题目： 组成为60%（mol）苯，25%甲苯和15%对二甲苯的100Kmol液体混合物，在102.3Kpa和100℃下闪蒸。试计算液体和汽体产物的量及组成[2]。 解：此条件下物系可看成是理想体系，因此可以用安托因公式计算饱和蒸汽压。 3．2求解程序[3] #includestdarg.h #includemath.h #define T (100+273.15) static float Z[3]={0.60,0.25,0.15},K[3]; main() {int i; double p[3]; float f,f1,x,y,a,b,s; float V,F,L,PO; static float l[3]={20.7936,20.9065,20.9891}, m[3]={2788.51,3096.52,3346.65}, n[3]={52.36,53.67,57.84}; clrscr();/* K */ PO=101.3; for(i=0;i3;i++) { p[i]=exp(l[i]-(m[i]/(T-n[i]))); K[i]=(p[i]/PO)/1000; printf(k=%f\n,K[i]); } /* s,f,f1*/ s=0.1; for(i=0;i3;i++) {do {f=0.0; printf(s=%f ,s); fo