氨滑动弧等离子体激励的数值建模_赵凝秋.pdf

第 卷 第期空军工程大学学报 年月 收稿日期：基金项目：国家自然科学基金青年基金（，）作者简介：赵凝秋（），男，河南商丘人，硕士生，研究方向为等离子体点火助燃。：引用格式：赵凝秋，朱益飞，李振阳，等氨滑动弧等离子体激励的数值建模空军工程大学学报，（）：，（）：氨滑动弧等离子体激励的数值建模赵凝秋，朱益飞，李振阳，尹博（西安交通大学机械学院，西安，；空军工程大学航空动力系统与等离子体技术全国重点实验室，西安，）摘要在双碳背景下，氨的使用及其燃烧难题的解决受到国内外学者的重视。针对改善氨燃烧问题，研究了二维流体动力学零维反应动力学燃烧动力学思路，分析了滑动弧等离子体条件下氨的裂解和燃烧特性。利用 计算软件建立了考虑电磁的二维流体计算模型，研究了滑动弧在反应器中的演化；建立了适用于氨的零维滑动弧模型，研究了滑动弧中全电离转捩非平衡个阶段的温度和组分演化规律；利用 计算软件建立了燃烧模型，评估了等离子体对 混合气燃烧特性的影响，此外对在混合气中添加进行了分析。结果表明，等离子体能够促进氨的裂解，降低氨的点火延迟时间，并且提高其层流燃烧速度，但是在混合气中将部分等离子体条件下的 替换为后，整体燃烧效果有所下降。关键词氨裂解；滑动弧；等离子体调控燃烧；数值模拟 中图分类号 文献标志码文章编号 （），（，；，），；，；，；全球气候变化成为越来越多人关注的议题，减少碳排 放，改善 温室 效应 成为人 们 共 同的 祈愿。年我国发布的 新时代的中国能源发展 白皮书中明确提到，中国要大力推动化石能源到非化石能源的转化。在此背景下，可替代的低污染燃料被众多学者提出，如乙醇、二甲醚（）、生物柴油和氢气等，与此同时新型无碳燃料氨（）的使用发展同样吸引了人们的兴趣。氨（）是一种氮氢化合物，是重要的氢基燃料，燃烧产物仅有水、氮气和少量的氮氧化合物是氨气的巨大优势。由于氨产业基础较强，并且 的条件下就可以压缩成液体，氨气可大规模地用于地方能源储存以及持续运输。此外氨气相比于氢气储存成本更低，相同条件下纯氢气的储存成本大约是氨气的 倍。氨中氢的质量分数达到 ，并且能量密度为 ，可与化石能源媲美。氨气由于刺激性气味的特性，发生泄露时很容易被察觉，能够在危机发生时极大的减少损失。但是想要将氨气发展为运输燃料，仍然要解决个主要问题：一是燃烧质量差、层流燃烧速度慢、可燃性范围窄和点火极限较高，导致燃烧时循环变异性高；二是排放量较高，超过了制定的标准排放量。但这些问题都可以借助等离子激励燃烧进行解决。等离子体是区别于物质固态、液态、气态种状态之外的第种状态，其在流动控制 和点火助燃 领域都有广泛的研究和应用。作为等离子体点火助燃应用的形式之一，滑动弧是一种由于气流作用，在个或者多个电极之间发生的具有周期性自振荡的放电过程，因具有的气体快速加热能力和丰富的反应动力学特性，有望成为激励燃烧的主流方案，越来越受到研究团队的青睐。近年来，等人建立了低电子能的等离子体动力学模型来进行非平衡等离子体对合成氨的研究。等人利用 分子与 分子的相似性，将经过实验验证的 等离子体机理通过简单地取代转化为 机理，并在此机理中加入目前已知仅有的 个包含电子碰撞截面的数据库：（）网站上的 和 数据库，最终模拟了等离子体辅助裂解和燃烧的纳秒放电过程，得出等离子体具有降低点火延迟时间的功效，而后其研究团队也对反应路径进行了分析。等人首先组装建立了 机理，并且运用 代码和 代码相结合的零维求解器进行了测试，运用不同的脉冲放电频率和脉冲数，在常压和中高温（）的条件下分析了其对点火延迟时间的影响，结果分析较高数量的脉冲数可以显著减少点火延迟时间，较高的脉冲频率可以允许较低的初始点火温度。等人研究了在预混式旋流燃烧器与滑动弧放电系统中不同氨含量与空气流量对预混 火焰稳定性和排放特性的影响，得到放电增加了自由基浓度，提高了火焰稳定性，在氨含量较高时，滑 动 弧 可 以 有 效 的 减 少 排 放 最 高 达到 。但是对于等离子体激励氨气点火助燃的数值模拟仍然困难重重，主要是由于等离子体反应动力学机理以及计算模型的缺乏，并且等离子体激励氨燃料裂解的研究少之又少，等离子体条件下氨的燃烧特性更是鲜有报道。基于此，本文研究了氨在等离子体条件下的裂解和燃烧特性，省略了等离子体通道导通并初步形成平衡电弧的阶段，重点关注电弧形成后其演化过程、对燃料的裂解作用和对燃烧的影响。模型构建 二维流体动力学模型使用基于有限元法的商用软件 建立了考虑电磁的二维流体计算模型，此模型计算滑动弧的时空演化，以此得到零维滑动弧模型的输入条件。反应器三维几何结构如图所示，反 应 器由 高 压电 极 和地电极 构成，入 口半径，高压电级半径，出口半径，入口流速 ，腔体总深度。结合实验估算平衡电弧的温度和体积作为此模型的初始条件。在局域温度平衡假设下，电导率可以看成是温度和压力的函数 ，所以滑动弧的温度可以根据实验截取到的平均功率和电导率估算获得：（）式中：为滑动弧电源功率；为滑动弧内部的电场强度；为滑动弧体积。其中电导率可从数据库获得（如 ，见西安交通大学电气工程学院气体放 电 等 离 子 体 数 据 库。：）。第期赵凝秋，等：氨滑动弧等离子体激励的数值建模滑动弧体积则需要依托于实验图像计算，也可以进行估算得到（见式）。（）式中：为滑动弧直径，单位为；为环境压力，单位为；参数取决于电极结构和气体混合物。图反应器模型三维视图 多物理场模型计算的是平衡滑动弧，虽然不能准确的计算出滑动弧的温度，但可以模拟出滑动弧的形态和长度，此模型需要求解斯托克斯方程组：（）（）（）（）（）（）（）（）式中：为密度；为流体的速度；为压力；为流体黏度；为单位矩阵，为动量守恒方程中的洛伦兹力项；为气体的比热容；为气体温度；为热导率。除此之外，还需在模型中添加麦克斯韦方程组以描述电磁场的变化。（）（）（）（）（）式中：为电导率；为磁感应，是磁矢量势。气体物理参数（、）取自气体放电等离子体数据库（）。模型使用的边界条件如表所示。入口流场为流量 ，平均功率分布在阴极上，速度在壁面处具有无滑移条件，在出口处压力为 （）。表 模型的边界条件位置流场热传递电场磁场阳极阴极 入口 出口 壁面 纯氨气等离子体机理构建整合建立的氨气等离子体机理主要包括电子碰撞反应、带电粒子反应中的电荷转移反应和解离复合反应、激发态反应。其中电子碰撞反应来自 网站上的 （见 ：，）和 （见 ：，）数据库，带电粒子和激发态反应主要来源文献整合，并且根据一些计算方法验证，激发态反应均采用 和 估算。零维反应动力学模型结合零维等离子体动力学求解器 建立了适用于氨气的零维滑动弧等离子体计算模型，此模型忽略掉滑动弧内部温度和粒子的不均匀性，以一个点来表示整段弧，考虑了反应加热、滑动弧与周围空气对流换热、滑动弧的增长对温度的影响等。使用此模型可以得到滑动弧个重燃周期随时间的温度演化和组分演化规律。此处建立的滑动弧的等离子体反应动力学模型，方程如下：（），（）（）式中：（）是第个反应对第个粒子的源项；为气体温度；为等离子体总密度；为比热比；为反应热源项；为其他热源，在不同的阶段有不同的源项。此处使用实验和 流体仿真结合的形式确定模型的输入条件。零维滑动弧计算模型需要以下输入条件：）有效输入功率；）滑动弧随时间变化的长度；）滑动弧随位置变化的速度；）滑动弧平均温度。当电压足够时，反应器中高压电极和地电极在距离最短处附近击穿，进而形成电弧放电。在形成电弧后，随着纯氨气流的推动使滑动弧拉长以及对流换热和辐射等产生热损失，当输入功率增加的能量不足以抵消热损失时，气体温度开始下降。但是空军工程大学学报 年电场依然存在，导致电弧中电子温度远高于气体温度，当电弧被拉长到一定长度时，电弧熄灭并开始在电极最短处重燃。结合这个物理过程对零维模型做了简化和设计，图给出了模型的阶段所对应的滑动弧物理过程。模型在第阶段赋值较大的电子温度和根据式（）估算到的气体温度，使电极间的物质在较短时间内发展成为平衡电弧的物理状态。在第阶段时电弧（局部热平衡等离子体）已经形成，此时电子温度等于气体温度，气体温度来源考虑了不同的热增加和热损失，此时功率输入的能量还能维持电弧稳定。随着电弧长度继续增长，在达到一定的转捩条件后，滑动弧进入第阶段非平衡阶段。此时电弧中热损失远大于热输入，气体温度迅速下降，但是由于依然存在电场可以给电子输入能量，导致电子温度远大于气体温度。最后到达个时间节点后电弧熄灭，完成个周期的循环。其中前个阶段电子能量基于麦克斯韦分布，在第阶段电子温度则是基于玻尔兹曼分布。图滑动弧物理过程和模型阶段对应图滑动弧一个周期的持续时间可以从实验功率中得到，部分实验功率如图所示，图中有较多功率增长而后降低的峰值，表示着滑动弧发生多次重燃熄灭周期。在每个周期中功率最低时，标志着滑动弧开始重燃，随着弧长增加，功率也随之增大。为了使用求解器进行计算，此处需要截取一段滑动弧重燃到熄灭的功率，本次研究截取的是图中粉红色背景下的功率，红色较深处代表滑动弧开始重燃，颜色变淡后表示滑动弧长度增加弧内温度下降，在功率最高时滑动弧熄灭。从图中可以看出，滑动弧在个重燃周期内时间范围为 ，从重燃到熄灭持续了 。图滑动弧反应器的实验功率模型计算结果与讨论 二维流体模型计算结果在二维流体计算模型中得到的滑动弧的演化情况，在 时刻的位置和温度如图所示，此时刻对应着滑动弧长度为最大值，即将熄灭，温度在 附近。图滑动弧温度在整个周期内长度和速度演化如图、图所示。在此重燃周期内的起始点 时刻时滑动弧重燃，此时滑动弧在整个周期内长度最短，平均温度最高。随着时间的推移，滑动弧长度开始增加伴随着温度下降，速度则是主要由结构和入口流量决定，由于出口相较于入口较小，所以在位置接近出口距离处速度骤增。滑动弧经过 的演化，在 时再次熄灭。图滑动弧长度演化图滑动弧速度演化第期赵凝秋，等：氨滑动弧等离子体激励的数值建模 零维滑动弧模型计算结果在零维滑动弧模型中得到该段滑动弧温度演化、粒子密度演化等数据。在截取的重燃周期内，得到的滑动弧气体温度如图所示，前设置的电子温度和气体温度都较高且保持不变，滑动弧内为高度电离状态，对应滑动弧的重燃时刻，此时模型处于全电离阶段，在本文中此阶段并不是关注的重点。随后滑动弧进入平衡转捩阶段，滑动弧会受到滑长度增大等影响产生热损失，气体温度开始降低，但是此阶段电子温度仍等于气体温度。在 时滑动弧熄灭，滑动弧熄灭之前气体温度维持在 的高温并且持续下降，滑动弧熄灭之后温度迅速下降至 ，与同功率相似结构的实验 对比（见图）良好。图一个周期内滑动弧气体温度演化图反应器出口处气体温度实验与仿真对比由于入口处通入的是纯氨气体，所以主要聚焦于等离子体对氨的裂解作用。气体从入口到出口大约需要 ，此时间范围内主要分子密度的演化如图所示。其中 在等离子体的作用下分解成了大量的燃料碎片 和，而它们中的一部分又会在等离子体作用下继续分解，最终形成大量的 和。约化电场的大小在全电离阶段为 左右，在此阶段较高的电子温度使总粒子数密度快速增加，因此约化电场大小急剧降低；在平衡转捩阶段总粒子数密度变化较小，约化电场变化同样较小；时滑动弧熄灭等离子体通道断开，此后作用在粒子上的电场强度为。后粒子到达反应器出口处，与空气结合进行燃烧。图滑动弧中主要分子演化 燃烧模型计算结果在燃烧模型中采用 等人组建的以氨气为主要成分的新型氨气机理 ，此机理包括 合成气子机理、反应子机理、以及含碳氢化合物反应的子机理。将纯氨气体经滑动弧裂解后的产物（反应器出口的粒子密度）代入到燃烧模型中进行计算，通过在模型中设置燃料和空气当量比的不同，研究有无等离子体混合气火焰传播速度和点火延迟时间变化情况。层流火焰燃烧速度在反应器出口处，经滑动弧等离子体裂解后的氨气与氧化物空气结合燃烧，此研究设置未燃气体温度值为 。在不同当量比下的层流传播速度如图 所示，横坐标为当量比，纵坐标为层流火焰传播速度。随着当量比从 增加到 ，火焰传播速度整体上先增大后减小。未施加等离子体时，当量比为 左右其火焰传播速度最大，施加等离子体后此当量比向着