第4章过程系统的动态模拟解说.pptVIP

第4章 过程系统的动态模拟 4.1 动态模拟的意义与应用 4.1.1 化工过程系统的动态特性 动态特性是化工过程系统最基本的特性之一。 间歇过程：变压吸附、变温吸附、化学反应器周期操作等； 连续过程的开停工； 连续过程中原料性质改变、催化剂活性变化等； 意外事故，认为误操作等； …… 集中参数模型：状态变量在系统中呈空间均匀分布 （强烈搅拌的反应罐） 分布参数模型：状态变量在系统内呈非均匀，但一般是连续的空间分布 （管式反应器） 多级集中参数模型：一般用于描述多级串连、级内状态变量均匀分布的过程（板式塔内的传质分离过程） 化工过程系统确定性动态模型的数学表达形式 4.1.5 确定性动态模型的数学处理 正问题—模型方程组的求解 逆问题—模型参数的估计 过程系统的定性分析 定常态多重性、定常态稳定性、参数敏感性、自激振荡，甚至更复杂的时间序列结构。 4.2 连续搅拌罐反应器的动态模拟 选择理由： 通常采用集中参数模型，典型性； 在模型的数学处理方法方面，与其它类型的化工过程系统集中参数模型也有相似性； 常常涉及到非线性系统的定性分析问题，也具有典型性，所运用的分析方法有普遍意义。 4.2.1 动态数学模型 例1：敞口连续操作搅拌罐的流量计算。 进料量为Fi, 原有料液高度为H0，试求取自开工后排料量的变化关系。 设搅拌罐的横截面积为A，排液量与罐中料液的高度成正比关系，即： Fo＝k·H 将初始化条件：t=0时，H=H0代入式，并化简可得： 排液量与时间的变化关系为： 例3-2：搅拌槽内含盐量的动态模型 作盐组分的物料平衡，有： 表明有两项累积量，第一项是因浓度变化而引起的，第二项是由体积变化所引起的，这两项皆与求解有重要关系。 积分，并利用初始条件 t=0时，V＝V0，可以得出： 其中，B为积分常数。 将初期条件：t=0时，C=0代入式，可以解出B，于是可以化简为： 小结 以上例子通过一些理想化的假设，削减了过程的复杂性，使得该过程可以通过数学方式精确求解 对于一般的连续搅拌罐式反应器，除总物料衡算和组分物料衡算外，还存在着伴随化学反应的热效应以及反应罐本身的热衡算。 对于这种复杂的过程，是不太可能通过数学方法精确求解的，一般要通过数值方法进行积分运算，方可求得过程的解。 普遍性的CSTR问题 通常假定反应罐内处于分子级理想混合，且为液相均相反应，因此可以认为反应混合物的温度和组成在反应区里是均匀的， 进一步假定反应区的容积不随时间变化，则加料与排料的流量也可以认为是近似相等的，即Fin ＝Fout=F。 其中，V、F分别代表反应区容积和加料容积流量； Ci 、Ci,f分别代表反应器内和加料中第i组分的浓度； t表示时间； 其中，T、Tf分别代表反应区内和加料混合物的温度； U表示反应液体与冷却剂之间热交换的总传热系数； A表示反应液体与冷却剂之间的总传热面； Tc表示冷却剂平均温度； ?、Cp分别代表反应混合物的平均密度与比热容； （－Hj）表示第j个反应的热效应； Rj表示第j个反应的速率； Ri表示因化学反应引起的第i个组分浓度的变化速率 其中，μi，j表示第j反应计量式中i组分的系数。 初始条件的约束 4.2.2 模型的数学处理与应用 上述动态数学模型的求解在计算数学上是典型的常微分方程组的初值问题，通常可以利用龙格－库塔法（R-K），基尔（Gear）法等通用程序来求数值解。 * 四阶龙格-库塔公式 4.2.2.1 应用1―开工过程分析 计算开工过程所需要的时间：从给定的初始条件出发，求模型的数值解，求取直至状态变量的每一个分量Ci、T接近定常值所需要的时间，就是近似的开工时间。 研究初始条件对开工过程的影响：改变不同的初始条件，通过数值分析考察初始条件（开工条件）的不同对开工时间的影响，了解在开工过程中系统状态变化的经历与初始条件的相互关系，从而可以帮助制订适当的开工方案，达到既缩短开工时间，又不致使开工过程出现某些工艺上不允许的温度和浓度。 4.2.2.2 应用2―动态响应的数字仿真 在控制系统合成过程中，了解被控制对象的输入输出关系是最基本的需要。 传统的方法是在对象上进行实验测试，既耗费人力物力，还可能会干扰系统的正常操作. 利用数字仿真技术来了解对象的动态响应特性，即输入输出关系，就要简单得多。 4.2.2.3 状态空间分析 动态模型方程组的解是非线性代数方程组，若此非线性方程组有多重根，就意味着系统