华南理工大学 化工系统分析与合成 第3章 化工系统模拟方法及策略1.pptxVIP

第3章

化工系统模拟方法及策略;3.2序贯模块法(SequentialModularMethod);3.1;3.1.1基本物性计算;（1）饱和蒸汽压计算;(2)蒸发潜热计算;也可以采用Riedel方程式:;（3）气体定压比热容计算;以甲烷的定压比热计算为例，这次以AspenPlus模拟软件计算获得的数据作为参考数据，以500K~600K的甲烷比热作为计算比较基准。获得甲烷摩尔热容的数据：;（4）液体定压比热容计算;在高对比温度下（0.8Tr0.99）C?L、CPL和CSL三者的数值并不相等，可利用两个近似关系式进行互相转换计算：;（5）气体黏度计算;如果是极性气体计算粘度碰撞积分时，要考虑偶极距的影响，需对?V进行校正，对于极性气体，计算粘度碰撞积分时，要考虑偶极距的影响，可按下式进行校正：;极性分子的d、?o/k和δ数据如查不到，则可采用下列近似式计算：;【例3-1】求乙醇在150℃时的低压气体粘度（实验值为123.2?P）。已知乙醇的Tc=516.2K,Pc=63.8bar,Vc=167cm3.mol-1,?=0.635,;解：;由于，所以;计算所得的粘度和实验值只有-1.06%误差，表明式（3-17）的粘度计算公式具有较好的精度;（6）液体粘度计算;气体和液体产生粘滞性原因的不同，可以从两者的粘度和温度关系的不同得到证明：对气体来说，温度升高则分子碰撞更加频繁，于是粘度随温度升高而增大。而对液体来说，温度升高动能增大，引力作用变得不那么重要，于是粘度随温度升高而变小。液体粘度估算的Andrade方程如下;对于多原子气体，Eucken提出：;【例3-3】估算2-甲基丁烷蒸气在1bar和373K时的导热系数。文献值为2.2×10-2W.(m.K)-1。已知Tc=460.4K，Pc=33.4bar，Vc=306cm3.mol-1，Zc=0.267,ω=0.227，M=72.151，373K时理想气体热容Cp=144.2J.(mol.K)-1。;解:1）应用Eucken法

由于无计算状态下的粘度数据，需先用前面介绍的方法估算粘度，这里省去，计算得到黏度为85.7μP（1P=1g/cm·s=0.1kg/m·s，1μP=1×10-6×0.1kg/m·s=1×10-7kg/m·s），CV用下式计算;2)Ubbelohde法;3）应用Stiel-Thodos法;（8）液体导热系数计算;纯液体导热系数还可用下列关联式估算：;【例3-4】估算CCl4在293K时液体的导热系数，文献值为0.103W.m-1.K-1。已知CCl4的Tc=556.4K，Tb=349.9K，M=153.823，CPL(273K)=130.7J.mol-1.K-1，?L(273K)=0.0106mol.cm-3，CPL(293K)=132.0J.mol-1.K-1，?L(293K)=0.0103mol.cm-3;3）使用Missenard-Riedel法;3.1.2线性方程组计算;对于下，下面方程线性方程组：;目录;3.1.3非线性方程计算;二分法是求方程近似解的一种简单直观的方法。设函数f(x)在[a,b]上连续，且f(a)f(b)0，则f(x)在[a,b]上至少有一零点，这是微积分中的介值定理，也是使用对分法的前提条件。计算中通过对分区间，逐步缩小区间范围的步骤搜索零点的位置，图3-2给出该方法的计算示意图;在化工设计过程中，常常需要用到各种物性数据，有些物性数据是直接以数字的形式出现，但许多物性数据是以公式的形式出现，并且有些公式是影性函数，无法直接将已知条件代入公式求解，如关于重要热力学数据P-V-T的著名马丁-侯方程：;图3-4excel计算气体摩尔体积;functionqiugen

clearall

clc

x0=-4;%给定初值

x1=fsolve(@f,x0)

x0=1.2;

x2=fzero(@f,x0)

x0=3;

x3=fzero(@f,x0)

%下面的系数向量需和自定义方程中的系数一致，注意系数按降幂排列

c=[11-2815];

x4=roots(c)

c=[101-280015];

x5=roots(c)

functionf=f(x)%读者可以改变下面的方程次数及系数，只要按

范例中的模式书写即可

f=x^3+x^2-28*x+15;;核心代码如下：

forp=1:30

i=i+1;

v(i)=fsolve(@f,x0,[],t,p)

eer(i)=f(v(i),t,p)

end

yp=1:1:30

plot(yp,v,r,yp,e