颗粒群平衡模拟的随机模型与燃煤可吸入颗粒物高效脱除的研究

发布时间：2020-07-12 12:31

【摘要】： 可吸入颗粒物(空气动力学直径小于10μm的颗粒物,简称为PM或PM10)已经成为大气环境污染的突出问题.由于PM的微观性和其生成、生长和演变所经历的复杂物理化学过程,无论是对其的采样和物理化学特征分析、在燃烧过程中的形成机理和微尺度动力学演变规律,还是在自然条件和外加条件下的非线性和非稳态运动规律、各种高效捕集策略的设计和优化等方面,人们都缺乏足够的认识,特别是缺乏相关的理论描述和定量描述.本文正是在此背景下,为燃煤PM复杂的动力学演变过程(包括碰撞、凝并、破碎、冷凝/蒸发、成核和沉积等动力学事件)构建颗粒群平衡模拟随机模型的完整数值模拟体系结构,并以此为平台对静电除尘器和湿式除尘器的除尘过程、自然环境中PM的干沉降和湿沉降过程进行颗粒群平衡模拟,对这些除尘机理进行深入分析,最终以此为基础对两种高效PM除尘技术(静电增强湿式除尘器和静电布袋混合除尘器)进行可行性分析和运行优化,为其工业化应用提供理论基础和技术指导.本文的研究成果包括: (1)零维颗粒群平衡模拟的随机算法本文引入“加权虚拟颗粒”概念而发展了一种全新的零维颗粒群平衡模拟随机算法——多重Monte Carlo (Multi-Monte Carlo, MMC)算法.该算法具有时间驱动Monte Carlo (MC)、常数目法和常体积法的特点,由于跟踪数目较少的虚拟颗粒群而具备适用于工程计算的计算代价,由于保持恒定的虚拟颗粒总数目而具有稳定的计算精度,由于保持恒定的计算区域体积和基于时间驱动技术而具有与两相湍流模型无缝耦合的扩展性.特殊工况下MMC算法结果与相应的理论分析解之间的良好吻合,证明该算法能够以较高计算精度和计算效率描述初始单分散性或多分散性颗粒群工况、单动力学事件或复合动力学事件工况等. 本文首次建立了MC算法对颗粒尺度分布时间演变过程描述精度的定量评判方法.以此为基础,定量比较了时间驱动直接模拟MC、阶梯式常体积法、常数目法和MMC算法对于各种动力学事件的描述精度,并定性分析了各种MC数值误差的来源和计算代价的主要贡献源.针对MMC算法所表现的“常数目误差”,通过建立“异数目权值虚拟颗粒群策略”以及相应的“异数目权值虚拟颗粒凝并准则”、在时间驱动MC中引入接受-拒绝法来数值实现标准Markov过程、发展“常数目方案”和“阶梯式常数目方案”来恢复虚拟颗粒数目等手段,对MMC算法进行改进.多种MC算法之间的定量比较表明,改进的MMC算法已经成为时间驱动MC中计算精度和计算效率最高的算法之一. 基于事件驱动MC对解藕误差的免疫及高计算精度和效率的认识,本文发展了另外一种全新的零维颗粒群平衡模拟随机算法——事件驱动常体积(EDCV)法.该算法仍然引入“加权虚拟颗粒”的概念和保持常体积特征,并把“等数目权值虚拟颗粒群策略”和“异数目权值虚拟颗粒群策略”、累积概率法和接受-拒绝法、“常数目方案”和“阶梯式常数目方案”以及“重整方案”等技术统一在一个框架下.通过对多种主流MC的定量比较,证明EDCV法成为目前零维颗粒群平衡模拟中计算精度和计算效率最高的随机算法之一. (2)多维颗粒群平衡模拟的随机算法通过把描述颗粒动力学演变的改进MMC算法与描述两相湍流场的欧拉-拉氏模型的无缝耦合,首次发展了一种多维颗粒群平衡模拟的随机算法——多维MMC算法.该算法引入网格划分技术和设置合适的时间步长,采用随机过程来判断碰撞/凝并等动力学事件的发生、寻找碰撞/凝并伙伴、建立碰撞动力学等.对细微颗粒流和粗重颗粒流中颗粒碰撞和凝并过程进行了数值模拟,模拟结果与直接数值模拟(DNS)结果实现定量上的较好吻合,证明所构建的多维多重Monte Carlo算法为四向耦合的两相湍流模型和多维颗粒群平衡模拟提供一个高精度和高效率的数值描述方案.采用所发展的MC算法考察了轴对称突扩通道有旋气固两相流中颗粒碰撞对于两相流场的影响,发现颗粒之间的相互碰撞使得颗粒速度和雷诺应力重新分配并趋于各向同性,而颗粒湍动能及颗粒-流体脉动速度关联降低. (3)传统除尘装置和自然环境对PM的除尘机理研究为了理解外加条件对PM动力学特性的作用机理,采用事件驱动常体积法对典型的静电和湿式除尘器进行了颗粒群平衡模拟.发现静电力与惯性力的相互竞争导致0.1~1μm的细微颗粒难以被静电除尘器捕集,而布朗扩散和惯性碰撞机制对0.1~1μm区间颗粒物影响最小也导致其难以被湿式除尘器的液滴所捕集.对300MW燃煤锅炉配备的静电除尘器进口和出口烟尘颗粒进行现场采样,并进行实际静电除尘器的颗粒群平衡模拟,模拟结果和是试验结果达到定量的吻合. 为了描述从除尘器逃逸出来的PM在自然环境中的沉降过程,并借鉴大自然的自身净化机制,采用多重Monte Carlo算法对自然环境中PM的干沉降和湿沉降过程进行数值模拟.发现大自然对0.3μm附近的中等尺度颗粒物难以自身净化,降雨量的增加有利于各种尺度颗粒物的湿去除,而雨滴几何平均尺度越小、雨滴几何标准偏差越小,越有利于小尺度和中等尺度颗粒物的湿去除,但稍稍不利于大尺度颗粒物的湿去除. (4)高效PM除尘技术的可行性分析和运行优化通过对自然条件和外加条件下PM动力学演变规律的认识,采用复合外加条件提高PM的除尘效率.对两种高效PM除尘技术(静电增强湿式除尘器和静电布袋混合除尘器)进行了颗粒群平衡模拟,并以此为基础对这些除尘技术进行了可行性分析.模拟结果表明两种除尘技术均可能达到99%以上的整体质量除尘效率和整体数量除尘效率.对静电增强湿式除尘器的优化运行分析表明,增大含尘气流输运速度或减小液滴喷射速度、增大液气比或液滴荷质比、减小液滴几何平均尺度和几何标准偏差(使得喷雾液滴越细和越均匀),可以提高该种除尘器对PM的除尘效率.
【学位授予单位】：华中科技大学
【学位级别】：博士
【学位授予年份】：2007
【分类号】：X701.2
【图文】：

烟尘,模态,三峰分布,元素

超微米(supermicron)颗粒区, 烟尘颗粒群的尺度范围从纳米量级甚至跨越到毫米量级.图1.1为三种美国煤燃烧产生的烟尘颗粒尺度分布[35]. 燃煤产生的烟尘颗粒三峰分布正好与大气环境中 PM 的三峰模态相对应,即爱根核模态 (<0.08μm) 、积聚模态(~0.08μm-2μm)和粗粒子模态(>2μm)[40].美国环保署的科学家 Linak[2]曾经指出,元素周期表中几乎所有元素都存在于煤中.燃煤 PM 是煤中各种元素、各种组分、各种赋存形态混合物在高温放热过程中经历一图 1.1 燃煤产生的颗粒尺度分布

燃煤飞灰,颗粒,凝并

图 1.2 燃煤飞灰颗粒生成机理吸入颗粒物动力学演变过程的数值模拟灰颗粒的动力学演变过程进行定量描述将有助于揭示燃烧源一次核、破碎等)、生长(冷凝、凝并等)、沉积和迁移等机制, 是建立完善与控制理论和技术的基础. 但是燃煤产生的飞灰颗粒存在如下显著特分布复杂, 尺度谱通常满足双峰或三峰分布; ②颗粒尺度谱0μm 跨越 5 个数量级; ③飞灰颗粒生成机理复杂, 飞灰所经历的动力撞、凝并(或团聚、凝聚、聚合等)、破碎、异相冷凝/蒸发(或表面沉、沉积(或捕集、清除等)、化学反应、迁徙等; ④大部分研究以实验成机理限于定性描述甚至猜测, 定量描述 PM 形成机理的理论模型乏, 而且也缺少相关的化学、热力学和动力学参数.的数值模拟方法通过随机性轨道模型、确定性轨道模型、拉氏 La颗粒的轨道经历效应和历史效应, 并通过这些详细的轨道信息来判的发生, 如通过两颗粒的轨道交叉来判断颗粒碰撞、凝并事件的发

示意图,时间驱动,示意图,颗粒数

个标准的 Markov 过程轮流检查每颗颗粒是否将要发生某种类型的事件, 以及选择发生该种事件的颗粒(称为主颗粒), 对于凝并事件, 该 MC 采用累积概率法(thecumulative probabilities method, 也可称为倒数法, the inverse method)实现凝并伙伴的寻找. 该 MC 首次采用阶梯式倍增子系统的方法来恢复模拟颗粒数目. 当子系统内颗粒数目降为原来一半时, 通过拷贝现有颗粒的复制体来恢复子系统内颗粒数目, 虽然这种操作对随机变量的期望值和几何标准偏差不会产生影响, 但是将增大相对误差,Liffman[144]通过数值模拟发现, 相对误差符合如下关系式 δ<2/N(t)1/2, N(t)为子系统内模拟颗粒数目. 时间驱动 MC 由于解耦过程而导致计算效率和计算精度降低, 但由于几乎所有的两相湍流模型的数值求解均是基于时间驱动的模式, 所以它存在考虑颗粒尺度分布空间扩散以及颗粒拉氏轨迹跟踪的友好扩展性.

【引证文献】