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;i 3;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;i 3;i++ do f 0.0; printf s %f ,s ; for i