36精馏系统模型.docVIP

复杂精馏系统模型 通过微分质量、动量和能量衡算，可以获得控制体内部各局部处物理量，如速度、密度、压力、温 度、组分浓度等随时间和位置变化的普遍规律。 质量衡算： 能量衡算： 精馏稳态数学模型 精馏塔是一种应用极为广泛的单元操作设备。多元物系精馏系统的一种基本处理方法是把精馏塔视为由一系列理想的气液平衡级构成。下图所示是一种普遍化的精馏塔原理示意图。 下图所示精馏塔由0,1,2,…,n,n+1，共n+2个理想平衡级构成，其中，第0级为分凝器，第n+1级为再沸器，见图中塔右侧所列序数。描述这种普遍化精馏塔的数学模型方程包括： (1) 平衡级中组分物料平衡式，共c ( (n+2)个方程，c为组分数，简称M方程组。 第0级 V1y1,i－V0y0,i－(S0+L0)x0,i=0 第j级 Fj zj,i－Lj-1xj-1,i+Vj+1,iyj+1,i－(Lj + Sj)xj,i－(Vj + Gj)yj,i=0 第n+1级 Ln xn,i－Ln+1xn+1,i－(Vn + Gn)yn+1,i=0  (2) 平衡级中组分的相平衡式，共c ( (n+2)个方程，简称E方程组。 yj,i =kj,I ( xj,i (3) 平衡级气液两相中组分的分子分率加和式，共3n+4个方程，简称S方程组。 1－=0 (0≤j≤n+1) 1－=0 (0≤j≤n+1) 1－=0 (0≤j≤n) (4) 平衡级的热平衡式，共(n+2)个方程，简称H方程组。 第0级 V1H1－V0H0－(S0+L0)(h0－Q0=0 第j级 Fj(j+Vj+1Hj+1+Lj-1hj-1+Qj－(Lj +Sj)hj－(Vj +Gj)Hj=0 第n+1级 Vnhn+Qn+1－Ln+1hn+1－(Vn +Gn)Hn+1=0 以上四类合称MESH方程组 共计(2c+1)(n+2)+3n+4=2c(n+2)+4n+6个方程。 大方程组中包含的变量列于下表。 对于直叙类问题，结构变量(理论级数)、输入变量(进料量、进料组成、进料热焓)和操作变量(侧线气相和液相出料量、各级热负荷量、回流比和再沸比)均为指定。其余的为需要计算求值的变量。在计算变量中，由于回流比和再沸比为指定，故各级气相上升量中第n+1级已确定，各级液相下降量中第0级已确定，这两者为非独立变量。而气相组成 、液相组成 和气液平衡常数这三组变量中，只能有两组为独立变量，另一组必定为非独立变量。在下表中取液相组成为非独立变量。 表 普遍化精馏塔方程组涉及的变量 编号 变量名称 变量符号 变量数 直叙类问题的变量分类 指定变量 计算变量 独立 非独立 01 理论级数 N 1 1 02 进料量 进料组成 进料热焓 Fj,j=1,…,n zj,i,j=1,…n (j,j=1,…n n c(n n n c(n n 03 各级气相上升量 气相组成 气相热焓 Vj,j=1,…,n+1 yj,i,j=1,…n+1 Hj,j=1,…n+1 n+2 c(n+2) n+2 n+1 c(n+2) n+2 1 04 各级液相下降量 液相组成 液相热焓 Lj,j=1,…,n+1 xj,i,j=1,…n+1 hj,j=1,…n+1 n+2 c(n+2) n+2 n+1 n+2 1 c(n+2) 05 侧线气相出料量侧线液相出料量 Gj,j=1,…,n+1 Sj,j=1,…,n+1 n+2 n+2 n+2 n+2 06 各级热负荷量 Qj,j=1,…,n+1 n+2 n+2 07 回流比 再沸比 R=V0/L0 B=Vn+1/Ln+1 1 1 1 1 08 气液平衡常数 kj,i,j=1,…n+1 c(n+2) c(n+2) 10 总计 4cn+6c+9n+17 cn+5n+9 2c(n+2)+4n+6 c(n+2)+2 根据以上分析，总共有2c(n+2)+4n+6个独立变量，正好可由MESH方程组所包含的2c(n+2)+4n+6个方程式求解。 精馏动态数学模型 精馏塔稳态模拟过去已提出许多算法，主要有逐板计算法、三对角矩阵法和松弛法。逐板计算法是从塔的一端到另一端逐板求解，但这种方法要求以塔顶和塔底的初始估计值作为计算的初值，计算的稳定性很差；还有一种方法是三对角矩阵法，它克服了逐板计算法的不稳定性，但是计算需要的时间很长；松弛法是多组分精馏的第三种算法，它的基本思路是：将一个系统的动态模型积分到“无限长”的时间，求出精馏的稳态分布曲线，将非线性方程组表示成一个微分方