总第 120 期 2006 年 12 月 工 程 设 计 与 研 究 石灰石湿法脱硫反应的动力学数学模型与求解杨运华 〔摘  要〕 石灰石/ 石灰 —石膏湿法烟气脱硫是控制二氧化硫排放的常用方法。以单个石灰石颗 粒为研究对象 ,根据双膜理论 ,建立碱性溶液脱硫的传质模型 ,并利用 Matlab 软件进行模拟。通过对气相中的 SO2向液相(相界面)中的传质扩散和石灰石的溶解扩散进行理论分析 ,提出一系列合理的假设条件 ,进而推导出传质模型 ,并根据该模型得出了传质速率的必要参数。   〔关键词〕 湿法烟气脱硫 ;传质 —反应过程 ;数学模型 ;石灰石   以往人们主要通过经验公式和定性分析来进行烟气脱硫除尘装置的设计 ,这种设计一般可满足工程要求 ,但也不乏失败的例子。对于烟气脱硫除尘 ,每项工程的烟气条件都不尽相同 ,况且对每一烟气脱硫除尘的过程进行试验几乎是不可能的 ,因为试验对设备的要求比较高 ,诸如高温、高压等试验条件难以实现 ,不容易进行实验。这就需要按实际情况建立系统模型来进行研究。 在湿法脱硫技术领域 ,人们从 1969 年起就一直致力于建立系统数学模型 ,期望能够包含所有速率的控制步骤 ,但实际上大多数模型仅适用于某一工艺过程[1] 。本文旨在从反应动力学的角度 ,建立石灰石湿法脱硫的传质模型 ,讨论传质与反应对过程总速率的影响 ,为湿法脱硫技术提供理论依据 ,从而改进湿法烟气脱硫工艺装置的设计和操作。 1  石灰石湿法脱硫过程机理分析 1. 1  反应过程 烟气脱硫过程包括气相中的 SO2向液相(相界面) 中的传质扩散、气 —液相的化学反应 ,以及石灰石的溶解扩散等过程。 由于 SO2 最先被吸收到水中并转化为SO3 ,所以脱硫率由液相中 H2 SO3 与气相中SO2的平衡决定 , SO2与水的连续反应[2] : SO2 + H2 O = H2 SO3 H2 SO3 = H + + HSO - 3 HSO - 3 = H + + SO2 - 3 H2 SO3 + 1 2 O2 = H2 SO4 石灰石被液膜包围 ,在液相中的分解速率由表面反应机理控制 ,在液相中的连续反应为[3] : CaCO3Ca2 + + CO2 - 3 CO2 - 3 + H +HCO - 3 HCO - 3 + H +CO2 (aq) + H2 OCO2 (aq)ΖCO2 (g) 在液膜中溶解的 Ca2 + 、 SO2 - 3 、部分氧化 的 SO2 - 4 离子之间发生复杂的化学反应 : CO2 + + SO2 - 3 + 1 2 H2 O = CaSO3 ·1 2 H2 O CaSO3 + H2 SO4 + H2 O = CaSO4 · 2H2 O + CO2 这是一个复杂的涉及传质、反应及结晶过程的气液固多相体系。对于这种气固液多 52 工  程  设  计  与  研  究 总第 120 期 相反应体系 ,传质过程往往是过程的控制步骤 ,在液相中的反应相对较快 ,并假定为瞬间反应。 1. 2  双膜理论 1969 年 ,Ramachandran 和 Sharma 借助于双膜理论 ,首次导出了伴有固体颗粒溶解的气体吸收反应过程的数学模型[4] ,随后许多研究者在此基础上进行了扩展和修正工作[5~7] 。 根据双膜理论 ,SO2 的传质反应过程如图 1 所示。考虑石灰石固体颗粒反应物在其周围的液膜内发生溶解传质 ,并与溶解于液相中的 SO2气体反应物发生瞬间反应。 图 1  SO2 传质反应过程示意 2  数学模型的设定 2. 1  模型简化 完全真实

matlab气相分解反应动力学,石灰石湿法脱硫反应的动力学数学模型与求解相关推荐

  1. matlab气相分解反应动力学,第十二章 化学动力学.ppt

    快反应 慢反应 速控法 稳态法 (1) 一级反应 (2) 二级反应 高压时 低压时 RRKM理论 §12.7 光化学反应 * 物理化学(下) PHYSICAL CHEMISTRY (11) 碰撞理论 ...

  2. MATLAB机器人机械臂运动学正逆解、动力学建模仿真与轨迹规划

    MATLAB机器人机械臂运动学正逆解.动力学建模仿真与轨迹规划,雅克比矩阵求解.蒙特卡洛采样画出末端执行器工作空间 基于时间最优的改进粒子群优化算法机械臂轨迹规划设计 ID:4610679190520 ...

  3. matlab分析能控条件,一级倒立摆MATLAB仿真、能控能观性分析、数学模型、极点配置.doc...

    一级倒立摆MATLAB仿真.能控能观性分析.数学模型.极点配置 题目一: 考虑如图所示的倒立摆系统.图中,倒立摆安装在一个小车上.这里仅考虑倒立摆在图面内运动的二维问题.倒立摆系统的参数包括:摆杆的质 ...

  4. Matlab:Matlab编程语言应用之数学计算(求极限/渐近线求导数常微分方程求解求微分方程组的解求临界阻尼系数的解)的简介、案例实现之详细攻略

    Matlab:Matlab编程语言应用之数学计算(求极限/渐近线&求导数&常微分方程求解&求微分方程组的解&求临界阻尼系数的解)的简介.案例实现之详细攻略 目录 三.极 ...

  5. Matlab气液相界面,MATLAB,气液相平衡程序,求帮忙改一下。

    function obj=lsqfun3(t12,t21,t13,t31,t23,t32) %t12,t21,t13,t31,t23,t32,是所需要关联的NRTL方程参数: %a1,a2,a3,是组 ...

  6. 气液平衡matlab,气液平衡 实验报告.doc

    气液平衡 实验报告 化工专业实验报告 实验名称: 二元气液平衡数据的测定 实验人员: 同组人 实验地点:天大化工技术实验中心 606 室 实验时间: 2015年4月20日下午14:00 年级: 201 ...

  7. Matlab Robotic Toolbox V9.10工具箱(六):puma560 动力学建模与仿真

    puma560,是经典的机器人模型,很多教材上都用它作为例子. 其D-H矩阵为: j theta d a alpha offset 1 q1 0 0 1.571 0 2 q2 0 0.4318 0 0 ...

  8. MATLAB实战系列(十四)-如何通过YALMIP和CPLEX求解小规模(CVRP)路径调度问题(附MATLAB代码)

    前言 有约束条件的车辆路径问题(CVRP),可以看作是TSP(Traveling Salesman Problem,旅行商问题)的拓展.由于TSP已经被证明是NP难问题,所以CVRP也是NP难问题. ...

  9. MATLAB实战系列(二)- 如何使用YALMIP检验数学模型的正确性?

    参加数学建模竞赛的小伙伴们,为了能让大家在比赛的第一个步骤不出错,即在模型建立这个步骤不出错,因此,今天准备讲一讲如何检验自己构建的数学模型是否合理. 本次推文分为三部分: 1)YALMIP工具箱的安 ...

最新文章

  1. windos 2008 vista 下的端口范围改变
  2. nova-scheduler详解 openstack-ice版
  3. AD5272数字变阻器
  4. 05构建之法阅读笔记之三
  5. ue查找文件中每行第二个单词_UI设计和UE/UX设计有什么区别?它们的晋升路径是什么?...
  6. linux静态网络ip dns怎么设置,Linux下如何配置静态IP设置DNS和主机名?
  7. 监听在微信中打开页面时的自带返回按钮事件
  8. 线段树相关(研究总结,线段树)
  9. springMVC实现的crud操作
  10. cmd代码玩贪吃蛇_关于N行贪吃蛇回答的补充
  11. 2022年最全Java面试题库基础篇
  12. linux设置双屏拼接_Linux 与Windows(A卡、N卡)下折腾双屏、3屏拼接
  13. Visual Assist X 破解
  14. 哔哩哔哩视频下载神器
  15. 计算机怎么打开网络共享,windows电脑如何开启wifi网络共享呢
  16. 软考中项第三章 信息系统集成专业知识
  17. CSS中常用的选择器都有那些?
  18. 实现支付功能并生成二维码
  19. MyBatis(九):MyBatis类型处理器(TypeHandler)详解
  20. 如何使用 Python 构建自己的 AI 个人助理

热门文章

  1. Arduino篇(二)火焰、烟雾、CO等多个传感器组合实验设计
  2. java suggest_lucene的suggest(搜索提示功能的实现)
  3. Sphinx的一个应用实例
  4. 手机拍照光线不够白平衡自动改变_干货 | 手机摄影专业模式设置,您知道吗?...
  5. CS224n课程Assignment3参考答案
  6. C语言标准输入scanf()函数
  7. php larval框架运行环境,Laravel框架的运行环境配置(一)
  8. 计算机专业考研卷子是统一的吗,考研试卷是不是全国统一?考研的试题相同吗?...
  9. 软件需求工程 高校教学平台 需求变更控制文档
  10. 杰里之ESD 测试方式【篇】