资源描述：
《复杂精馏塔的数学模型建立和模拟计算》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、复杂精馏塔的数学模型建立和模拟计算应用化学0311马迪卓200349032前言简单蒸馏与平衡蒸馏只能将混合物进行初步的分离。为了获得较高纯度的产品，应通过多级的蒸馏过程，使混合物的气、液两相经过多次混合接触和分离。进行质量和热量的传递，使混合物中的组分达到更高程度的分离，这一目标可采用精馏的方法予以实现。精馏塔是一圆形筒体，塔内装有多层塔板或填料，塔顶设有冷凝器，塔底装有再沸器，塔中部适宜位置设有进料板。精馏操作是分离液体混合物最为广泛采用的重要分离单元之一。在精馏过程中所采用的塔通常是一股进料分离出两种产品，即塔顶馏出轻组分产
2、品D,塔底获得重组分产品W,称该塔为常规精馏塔或简单精馏塔。当分离多股组分相同而组分各异的混合物时，常采用多股进料塔。当分离多组分混合物时，其中某产品纯度要求不是很高，其挥发能力仅次于塔顶产品，故可采用抽侧线产品塔，以降低投资费。如果进料中含有少量轻组分，若要将其冷凝下来需要较高品位的冷剂。为此，可采用带部分冷凝器的精馏塔，将该少量轻组分以气象排出，然后单独处理。如果精馏过程允许带水，且要分离比水轻的组分，则可采用直接蒸汽加热方式，这样既提高了塔釜传热效率，又省出一台再沸器。如果混合物中待分离的两组分相对挥发度接近1或形成恒沸物
3、则应选择加入质量分离剂进行分离，如恒沸精馏或萃取精馏。对于那些批量小，组分经常变化的混合物分离通常选择间歇精馏。为了对上述精馏塔进行模拟计算，则需分别建立数学模型。本文下面就介绍用三对角线矩阵法来实现对复杂精馏塔的数学模拟。数学建模1．数学模型复杂精馏塔模型如图。设塔含N个理论板或平衡级。理论版序号从上之下编排，塔顶冷凝器序号为1，塔底再沸器序号为N。塔内理论板数为(N-2)该塔的每一个理论板上均有一个进料Fj,气相侧线Uj和一个中间冷却（或加热）器热流量为Qj。塔内的任意理论塔板如图。平衡级的严格计算应满足4组方程式物料衡算方
4、程（M方程），相平衡方程（E方程），摩尔分数加和方程即归一方程（S方程），热衡算方程（H方程）A．物料衡算方程（M方程）对塔内任意一个平衡级j上的组分i进行物料恒算可得B．相平衡方程（E方程）C．摩尔分数加和方程即归一方程（S方程）D.热衡算方程（H方程）以上公式中，i=1,2,3,……mj=1,2,3,……Nx(i,j)——平衡级j上的组分i的液相摩尔分数；y(i,j)——平衡级j上的组分i的气相摩尔分数；z(i,j)——平衡级j上的组分i的进料摩尔组成；K(i,j)——平衡级j上的组分i的平衡常数；L(j)——平衡级j流出进
5、入下一级的液相流量，kmol/h;V(j)——平衡级j流出进入上一级的气相流量，kmol/h;hj(j)——平衡级j上液相摩尔焓,kJ/kmol;Hj(j)——平衡级j上气相摩尔焓,kJ/kmol;Hfj(j)——平衡级j上进料摩尔焓,kJ/kmol;以上M、E、S、H方程分别含mN、mN、N、及N个方程，共含2N（m+1）个方程，该方程组共含变量为（4mN+9N）个。为求解M、E、S、H方程，应根据系统自由度的概念，给定其中2mN+7N个设计变量，使之具有唯一解。在需要给定的设计变量中，部分是由已知条件给定，如Fj、zfi,j
6、、Wj、Vj、Pj等。而另一部分则由补充方程，例如物流焓与组成、温度压力的关系等，确定、、。选择热力学状态方程确定。在求解M、E、S、H中2N（m+1）个方程时由于方程多、变量多，且变量间存在非线性关系，采用解析方法求解较难，通常采用数值方法，按顺序进行求解。2．数学模型求解先将M方程、E方程联立，有E方程消去M方程中气相组成并整理得并将上式表示为以下简单形式：式中组分i在塔内各理论级上液相组成的分布关系式可表示为以下非线性方程的形式。也可表示为以下矩阵形式或上面的矩阵即三对角线矩阵方程，是一非线性方程组，其求解过程是：首先根据
7、初值条件求得系数阵A即列矩阵D，使上面的非线性方程组线性化，而线性方程组可采用多种方法求解。计算步骤1．给定塔模拟计算的条件进料液的给定条件：压力：常压；烃化液流量：F=150kmol/h;烃化液温度：T=60度烃化液组成（摩尔分数）：苯＝0.5982; 甲苯＝0.00602; 邻二甲苯＝0.3091; 乙苯＝0.08668;2．假定各级的液相组成及温度分布初值设第一块板的液相组成为x11=0.98x21=0.001 x31=0 x41=0设第NT板的液相组成为x1NT
8、=0.002;x2NT=0.2035;x3NT=0.3125;x4NT=0.482设j表示理论板的块数，j=(2:NT-1)。建立各组分在各板上的组成的线性方程：x1j=0.98-(0.98-0.002).*(j'-1)/(NT-1)x2j=0.001+(0.2
