化工过程系统建模.ppt

《化工过程系统建模.ppt》由会员分享，可在线阅读，更多相关《化工过程系统建模.ppt（59页珍藏版）》请在三一办公上搜索。

1、2023/10/27,1,第11章 化工过程系统建模,2023/10/27,2,工业过程的基本定律与建模方法,11.1 连续性方程11.2 能量方程11.3 传递方程11.4 化学反应动力学11.5 过程微观动力学研究11.6 过程宏观动力学的研究11.7 流体力学过程的动态方程11.8 过程数学模型的确定,2023/10/27,3,11.1 连续性方程 总的连续性方程 对于一个动态系统的质量守恒定律（质量平衡）可以表示如下：在系统工程中，质量的变化率=进入该系统的质量流量-离开该系统的质量流量。,上面的关系式是单位时间质量。对于一个系统只有一个总的连续方程。上式左侧的变化率以数学形式表达时，

2、可以用导数 或偏导数 来表示。,2023/10/27,4,例11-1 下图11-1所示是一个均匀混合的储槽的流量为 其密度为，储槽滞料量为，它的密度因均匀混合与流出的密度相同，该系统的物料平衡方程。,（11-1）,2023/10/27,5,例11-2 当流体流经一个直径不变的圆管时，如图11-2所示。,当流体处于湍流流动状态时，取管道长 一个单元，并假设无径向速度及密度梯度，而只有轴向梯度存在。因此，当流体沿着轴向即Z方向流动时，速度与密度将有变化，因此出现了时间t及距离Z两个变量，所以分别 以及 表示。,2023/10/27,6,通过Z+dz 截面离开系统的质量流量=+,现取管截面为A，管段

3、长为，建立质量平衡方程：,系统内质量的变化率=,通过Z截面进入系统的质量流量=,2023/10/27,7,组份的连续性方程 在化学反应过程中与单纯物料进出时不同，不同组份的质量并不是守恒的，在一个系统中发生了化学反应，如果对生成物来说，则它的摩尔数将增加；而对于反应物来说，则反而减少。,系统中第j个组份摩尔数的变化率=（进入该系统的第j个组份的摩尔数流量）-（离开该系统的第j个组份的摩尔数流量）+（由于化学反应第j个组份的摩尔数生成速率）（11-3）,液体在流动过程中，具有体积流动和分子扩散的对流作用。对于一个系统只有n个独立的连续方程式（对于具有n个组份而言），通常选用总的质量平衡和（n-1

4、）个组份的连续方程式用于求解。,2023/10/27,8,例11-3 如图11-3所示，在连续搅拌的储槽内进行化学反应，通常称为连续搅拌槽式反应器。组份A以不可逆的特定速率k形成一个产品组份B。,2023/10/27,9,解：设进入反应槽组份A的浓度为CA0，以及反应器内和出口的浓度为CA。若是简单的一级反应，在单位体积内组份A的消耗率将与槽内A的瞬间浓度成正比，列出有关项，并代入（11-3）式。,流入系统的A=F0CA0；(F0为进料量，V为反应体积)流出系统的A=FCA；由化学反应使A的生成率=-VkCA在反应器内A的变化率=d(VCA)/dt,2023/10/27,10,上式为常微分方程

5、，这是因为对于集中参数表示的系统，t是仅有的独立变量。,上式的左侧为动态项，右侧一、二项是对流项，右侧的第三项为化学反应生成项。对A与B两组份形成的双元系数，同样可以列出对B组份的连续方程。,（11-5）,或者可以采用另一个总的连续方程来描述，因为，、及 有单值关系。,2023/10/27,11,例11-4 如图11-3所示，假设这一化学反应是连续进行的连串反应,设由A生成B的速率为，由B生成C的速率为，即,对于组份A，B，C的连续方程按一级反应可表示如下：,2023/10/27,12,于是，我们可以用（11-7）式中的三组方程式，或者取（11-7）式中的两组份平衡式和一个总的质量平衡式来表示

6、。,2023/10/27,13,假设其反应为 的 一级反应，在 一段内物料由A转化为B，所以 将下降，密度、速度 及浓度将随时间 及轴向位置 而变。在此，仍假设在径向上无速度和密度梯度。,例11-5 如果流体通过一个管道反应器，如图11-4所示。,2023/10/27,14,并设反应器入口处（=0时）A的浓度为，在 处反应器流出物中A的浓度为，以 段作组份A的连续方程式，流入项可以分为体积流量与扩散两项。由于轴向有浓度差，因此就会产生扩散。虽然在一般情况下对实际系统影响很小，但是把它列写出来可以观察扩散对模型的影响。对于A的扩散流通是可以用 表示：,2023/10/27,15,以 微元作物料平

7、衡：,在 处进入系统A的摩尔流量=体积流量+扩散项=,在 处离开系统的摩尔流量=,在系统内A的生成率=,在系统内A随时间的变化=,由以上得总的方程为：,2023/10/27,16,11.2 能量方程,2023/10/27,17,例11-6 设连续搅拌槽式反应器中，加冷却盘管以移走反应热（kcal/mol),放热时为负值，如图11-5所示。,图11-5 移走热量的连续搅拌反应槽,由于反应所放出的热量与参与反应的A的mol数成正比，即 QG=-VCAK（Kcal/h)，设被盘管移走的热量为-Q，,2023/10/27,18,设被盘管移走的热量为-Q，进入系统有物料温度为T0，反应器温度为T，根据（

8、11-11）式可得：,式中：U内能（kcal/kg);K动能(kcal/kg)位能(kcal/kg);W系统对外作的功(kg.m/h)P系统内的压力(kg/m2);P0进料压力(kg/m2);J换算系数；密度。,（11-12）,2023/10/27,19,11.3 传递方程,有关传递方程均以通量的形式（单位面积的传递率）来表示，且与推动力（温度、浓度或速度梯度）成比例，比例常数是系统的一个物理特性（如热导率、扩散系数或粘度）。,动量传递方程式 动量传递方程式主要描述运动的规律。这里仅研究由于外部原因（重力、压力差、摩擦力等作用）而引起的流体流动。（牛顿定律）,（11-13）,2023/10/2

9、7,20,上式中-单位面积上一层流体对其相邻的一层作相对运动时，分子间相互作用所产生的内摩擦力，也可以看成单位时间、单位面积内由速度高的一层流体传递给速度低的一层流体的动量。,-y方向上的速度梯度；-粘性系数(简称粘度)。,2023/10/27,21,上式中：-瞬时的热流量；-温度梯度；-导热系数。当温度沿着单位法线长度下降1度时，单位时间内传过单位面积的热量。,2023/10/27,22,上式中：-组分A的瞬时传递量；-浓度梯度，作为质量传递过程中的推动力；-分子的扩散系数。,2023/10/27,23,11.4 化学反应动力学,化学反应动力学是研究化学反应速度及影响反应速度的条件的科学，影

10、响化学反应速度的主要条件是：温度、浓度、催化剂等。作为提供反应器设计参数，通常保持催化剂不变，考察温度和浓度对反应速度的影响，求取该反应的宏观动力学方程式。,2023/10/27,24,上式K为速度常数，它与反应物浓度无关,是温度、反应物性质、溶剂等的函数。,可逆反应：,2023/10/27,25,（11-18）,串行反应：,（11-19）,2023/10/27,26,然而，工业上进行的化学反应大部份是机制不够明确的复杂反应，只能通过试验求得反应速度和浓度的关系，如：,（11-20）,2023/10/27,27,式中：K-为某一反应速率；-比例系数；E-活化能；T-绝对温度；R-理想气体常数。

11、,注：这一指数温度关系在化工系统中为最严重的非线性关系之一,2023/10/27,28,11.5 过程微观动力学研究,在化学反应过程中，相互作用是在反应物的分子水平上实现的，可认为这个过程是在微观水平上进行的，所以，对过程的研究属于过程微观动力学的研究。也就是在实验装置上获取有关化学信息的过程。,封闭系统的动力学数据的测定：浓度C的变化速度是时间t的函数，即,（11-22）,2023/10/27,29,浓度C的变化速度是浓度的函数，即,浓度C是时间t的函数，即,注：动力学实验数据一般分两步：先使温度固定不变，确定原 始物质浓度的变化；再求出反应速度常数随温度的变化。,（11-23）,（11-2

12、4）,2023/10/27,30,11.6 过程宏观动力学的研究,从实验装置得到过程微观动力学的信息（中间试验装置、较大的中间试验装置上）获取过程宏观动力学信息。,宏观动力学是研究大的分子团、物流宏观质点的相互作用。其主要任务是：求出和修正过程数学模型中的各个系数并确定模型的实用性。,2023/10/27,31,（1）用阶跃法、脉冲法和伪随机信号造成扰动（注入示踪剂）来确定被研究对象的流体力学模型的结构：是理想混合还是理想置换，采用扩散模型、槽列模型或组合模型。（2）研究各个传热工况参数的影响，修正传热条件和传热面积，分析稳定性。（3）研究对象对于不同扰动通道的动态性质，分析各参数随时间变化的

13、情况，提出控制对象的方案。（4）建立较完整的对象（过程）数学模型。（5）实现数学模型的最优化，再现和修正最优工况。,中间试验装置的研究：,2023/10/27,32,y2,非解析模型（包括：专家知识集、模糊逻辑、人工神经网络）,生产过程(实际现场),解析模型（包括：微分方程、传递函数、差分方程、脉冲传递函数、线性方程以及非线性方程等）,u,y,y1,e1,e2,图11-6 建模示意图,2023/10/27,33,2023/10/27,34,2023/10/27,35,11.7 流体力学过程的动态方程(流体力学的基本方程式)流体力学过程动态方程的列写，是从三个基本方程出发，它们是：连续性方程、运

14、动方程和能量方程。,1）连续性方程 连续性方程表明流体流动过程中的物料的平衡关系。设在流体流束中取介于截面1-1与11-2之间，其长度为无限小dz的微元，如下图所示。,2023/10/27,36,依据物料平衡关系，单位时间内流入微元的质量与离开微元的质量差等于微元内蓄存量的变化率：,上式中：,2023/10/27,37,上面(11-25)式为一维非平稳渐变流的连续性方程。若用平均流速V来表示，则流速与流量的关系是：,（11-26）,式(11-26)和(11-27)是一般情况下的连续性方程。对于各种特定情况，可化为其它形式，如：,2023/10/27,38,2023/10/27,39,流体不可压

15、缩的情况,此时，由(11-27)可得,（11-29）,稳态情况,对于平稳流来说，,由(11-25)式可得,（11-30）,（11-31）,2023/10/27,40,也就是说，在稳定流动情况下，沿着长度z没有变化，即,如果流体不可压缩，则,2023/10/27,41,2）运动方程 运动方程是流体力学过程力平衡关系的微分表达形式。对于下面图中所示的微元，讨论压力、重力、摩擦力和惯性力之间的平衡关系。,2023/10/27,42,应用于牛顿第二定律，可以得到,上式中：,2023/10/27,43,由于流速v是距离z和时间t的函数，即v=v(z,t)，因此上式可写成,（11-33）,将(11-32)

16、式代入(11-33)式，可得,（11-34）,由图11-10可知：,2023/10/27,44,式(11-34)和(11-35)都是一般情况下的运动方程。当忽略摩擦力的影响，即FR=0，则有,（11-36）,对于稳定情况下，所以上式可写成,（11-37）,2023/10/27,45,又因在稳定情况下，P和v都仅是z的函数，没有其它的自变量，所以 可写成，即,3）能量方程 能量方程表达了流体过程中的能量平衡关系。能量方程可以按能量的平衡关系直接推导，也可以由运动方程积分得出。对于没有摩擦、不可压缩的平稳流，把式(11-38)积分后得到：,（11-39）,2023/10/27,46,每项除以g，并

17、引入重度,（11-40）,上式是著名的柏努利方程，式中每一项都具有明确的物理意义。各个截面上速度头、压头和位头之和为定值。,2023/10/27,47,11.8 过程数学模型的确定,（1）分布函数 任何有液体和气体移动的化工过程的数学模型，其结构首先是由流体力学参数决定的，并表现出物流质点在所研究的系统中停留时间的分布性质。,这种分布的性质服从统计规律，可根据通过系统的信号形式求得。在系统的输入端以阶跃、脉冲或谐振扰动的形式注入示踪剂。,2023/10/27,48,引入脉冲扰动时，示踪剂的统计分布函数（C曲线）可写成：,(11-41),停留时间分布函数 表示在输出物流中，在小于 时间内示踪剂在

18、系统中出现的分布速率。,2023/10/27,49,分布函数 可表示为：,上式中：,-入口处的初始浓度。,(11-43),2023/10/27,50,阶跃扰动和脉冲扰动时的分布函数之间的联系为：,(11-46),由于 值表示输出物流中年龄为 物质的分布率，将此数乘以停留在年龄为 的物流单元中的物质浓度，即可确定离开实际设备的物质平均浓度为：,(11-47),2023/10/27,51,例11-7 确定设备的数学模型时，由于引入函数形式的扰动（脉冲式注入示踪剂），在设备输出端得出下列示踪剂浓度的数值：,表11-1,试作示踪剂分布曲线。,2023/10/27,52,解：由方程（11-43）确定分布

19、函数。为此，先求出采样 时间间隔,分布函数 随时间的值列于表11-2。,表11-2,2023/10/27,53,为了求出函数C，将时间变成无量纲形式，把C()变成C的形式。为此，按方程（11-42）、（11-43）求出在设备的平均停留时间。,由方程（11-44）求出无量纲时间,由此，根据方程（11-45）可得,2023/10/27,54,图11-11 示踪剂分布曲线,2023/10/27,55,上式表示，如果函数 描述某个系统随时间的行为，则它可用复数变量 的相应函数F(s)代替。,（2）频率特性 频率特性是通过以一定的频率（通常为正弦信号）往设备中周期性地注入示踪剂而获得的，此时，在输入和输

20、出端上得到振幅和相位都不同的正弦信号。,应用拉普拉斯算子法，频率分析可将时间函数变换成频率函数。拉普拉斯算子定义为：,(11-48),2023/10/27,56,称为原函数，称为象函数，其关系为：,(11-49),拉普拉斯逆变换可表示为：,(11-50),2023/10/27,57,若将输出时间函数进行拉普拉斯变换并除以输入函数的拉普拉斯变换，则得到传递函数：,(11-52),上式中：是输出、输入的振幅比；是输出信号相对于输入信号的相位差。,在较宽的频率范围内作传递函数，可得出频率特性图，然后与典型的反应器模型的已知频率特性进行比较，不同的模型结合起来就可以对实际设备的传递函数作出估计。,2023/10/27,58,例11-8 以一级反应的带搅拌器的流通设备为例讨论传递函数的求法。,解：设 是体积流量；为输入；为输出；k为反应速度 常数；V为设备的容积；为时间。,根据物料平衡，有：,(11-53),进行拉氏变换后，得,(11-54),2023/10/27,59,输出和输入之比为：,(11-55),引入时间常数，则有,或写成：,(11-56),(11-57),