天津大学化工数学设计作业.docVIP

化工数学设计作业 ----乙烷裂解炉中的衡算 院 系 名 称：化工学院 专 业 名 称：化学工程与工艺 组 内 成 员：组长 二○一一年五月二十八日 目录 一、 化工数学背景 3 二、 我们的课题 4 1. 课题的来源 4 2. 存在的问题 4 3. 报告内容 4 4. 组内成员分工 5 三、 选题详解 5 四、 符号说明 22 五、 致谢 22 六、 参考文献 23 化工数学背景 从化工数学发展过程来看，数学对化学工业起了很大的作用。最初的化学工程规模小，技术落后，经验操作，无数学模型，随着生产规模的扩大，提出了很多数学问题，但因技术落后，无法求解，在二十世纪七十年代，化工过程模型化及化工模型的求解技术得到了迅速发展，很多辅助应用软件应运而生，数学模型在化工界得到了普遍重视。 对于数学应用而言，首先是数学建模，所谓数学建模，就是指当需要从定量的角度分析和研究一个实际问题时，人们就要在深入调查研究、了解对象信息、作出简化假设、分析内在规律等工作的基础上，用数学的符号和语言，把它表述为数学式子，也就是数学模型，然后用通过计算得到的模型结果来解释实际问题，并接受实际的检验。其建模方法主要有两种，一是经验建模，二是根据物理规律建模。所研究的对象主要是动态模型，一般用微分方程或差分方程来描述，建模的过程中常常遇到模型辨识问题。从本质上讲，动力系统建模以及模型辨识问题都是要解决数学问题。 数学建模的目的是实现优化，从而达到合理控制的目的，近几年数学建模在解决控制与优化问题上起到了很好的效果，被广泛应用到控制过程的预测与软测量当中，为解决工业控制过程的优化提供了崭新的途径。我组设计作业正是基于数学模型的建立，从而解决化工生成所面临的问题。 我们的课题 课题的来源 化工流体流动课程中，我们正在学习化工单元操作。乙烷裂解炉是化学工程中常见的传质单元。常常涉及大量的物料衡算和热量衡算，并要研究温度，反应物的转化率，反应速率等随管长的变化。方程常以微分形式存在，难以直接求解。我们将在化工数学课程第四章中学到的用龙格库塔法求解常微分方程组应用于其中，编写成程序，简化了计算，可以得到任意管长下所需的数据，且精度较高。 存在的问题 在反应器中，需求解乙烷的转化率，反应温度随管长的变化。用反应速率方程式得到两个常微分方程，需求解初值问题。并将结果用函数表示。、 报告内容 化工流体流动中计算题 建立模型 得到方程组 用龙格库塔法求解 编写程序 程序检验 运行 得出结果 用最小二乘法拟合函数 得到曲线 组内成员分工 资料查找：王苗 问题分析，计算，编写程序（最小二乘法）：孙瑜 张馨钰 编写程序，计算（龙格库塔法）：张诗阳 排版，核对：杨洋 选题详解 乙烷裂解炉 裂解炉辐射室内炉管分4组，双面辐射。每组炉管10程，均为单管。每程有效长度约为9.15m，管径为Φ133×8.5。原料乙烷流量为384.01kmol/h。摩尔组成为1％，98％，1％。稀释蒸汽流量为193.40kmol/h，出料流量为755.01kmol/h。进口温度为677℃，压力为0.3637MPa；出口温度为847℃，压力为0.2097MPa。乙烷单程转化率为50％，有效热流量为33.0489×kJ/h。 管长核算 核算1组裂解管管长，核算时采用下列简化条件： 不考虑副反应，只计算乙烷转化为乙烯的反应。 压力降与管长成正比。 热强度沿管长均匀分布。 加入的硫化物量极少，忽略不计。 基本计算式为 反应速度方程 考虑到逆反应，根据质量作用定律，乙烷裂解反应速度方程为 （1） 将反应物，生成物用公式C＝n／V＝P／RT换算成压力表达式，又，则上述公式（1）可写成 （2） 式中 T——反应温度，K； K——反应速度常数，1／h，由式（2）表示； ——以分压表示的化学反应平衡常数，由式（3）表示 ,和——乙烷，乙烯和氢的分压，Pa，由（4）表示。 （3） （4） 设乙烷在管长为ι处的转化率为XA，裂解管进口处（l＝0）和当管长为ι时各组的摩尔流量（kmol／h）如下： 管长 C?H6 C?H? H? C?H6 H?O 共计 （A）