高压液氨喷雾燃烧的数值模拟和试验验证.pdf

第 45 卷（2023）第 3 期Vol.45（2023）No.3柴油机Diesel EngineDOI：10.12374/j.issn.1001-4357.2023.03.001高压液氨喷雾燃烧的数值模拟和试验验证王宁1，李铁1，2，周昕毅1，2，陈润1，2，李世琰1，2，依平1，2（1.上海交通大学海洋工程国家重点实验室，上海 200240；2.上海交通大学动力装置及自动化研究所，上海 200240）摘要：以液氨喷雾混合气形成和燃烧过程为对象进行数值模拟和试验验证。根据试验获得的气液相贯穿距数据对三维模型进行参数标定。结果表明：建立的喷雾模型可以在高温高压环境条件下准确预测液氨喷雾的贯穿发展过程；在喷雾混合气形成过程中液氨气化吸热导致喷雾中心区域降温效果明显，降温高达100 K；氨较高的汽化潜热使之较难蒸发，喷雾边缘处氨的质量分数多在0.1上下。利用特征时间燃烧模型对液氨喷雾燃烧过程进行初步计算，可以较好地模拟整个燃烧过程。关键词：氨；高压直喷；喷雾燃烧；数值模拟中图分类号：TK421.2文献标志码：A文章编号：1001-4357（2023）03-0001-06Numerical Simulation and Experimental Verificationof Spray Combustion of High-Pressure Liquid AmmoniaWANG Ning1，LI Tie1，2，ZHOU Xinyi1，2，CHEN Run1，2，LI Shiyan1，2，YI Ping1，2（1.State Key Laboratory of Ocean Engineering，Shanghai Jiao Tong University，Shanghai 200240，China；2.Institute of Power Plants and Automation，Shanghai Jiao Tong University，Shanghai 200240，China）Abstract：Numerical simulation and experimental verification of the mixture formation and combustion process of liquid ammonia spray were carried out.According to the gas and liquid penetration data obtained in the experiment，the parameters of the three-dimensional model were calibrated.The results show that：The spray model can accurately predict the penetration and development process of liquid ammonia spray under high temperature and high pressure conditions；During the formation of the spray mixture，the vaporization of liquid ammonia absorbs heat，resulting in a significant cooling effect in the center of the spray as high as 100 K；Due to the poor evaporation caused by the high latent heat of vaporization of ammonia，the mass fraction of ammonia at the spray edge is mostly around 0.1.The combustion process of ammonia spray is preliminarily calculated by the characteristic time combustion model，which can simulate the whole combustion process well.Key words：ammonia；high-pressure direct injection；spray combustion；numerical simulation2022年内燃动力碳中和与排放控制学术年会绿色船舶动力系统分会场专栏收稿日期：2022-10-31基金项目：国家自然科学基金（52020105009、51776125）第 45 卷第 3 期 2 柴 油 机0引言当前，降低温室气体排放，加速向低碳经济社会转型是当务之急，引进替代燃料成为交通运输业的绿色环保解决方案之一。在众多的替代燃料中，氨燃烧后不产生二氧化碳，且易于储存和运输，已成为最具前景的替代燃料之一。此外，我国的制氨技术已经非常成熟，产量占世界氨总产量的25%以上。利用煤或者可再生能源（如风能、太阳能等）制取氨使之作为燃料替代石油制品切实可行，并且符合我国国情1。近年来，氨在内燃机和燃气轮机等动力机械中越来越受到关注。LHUILLIER 等 2、FRIGO等3-4、MORCH等5分别研究以氢气为助燃剂的氨燃料火花点火发动机的燃烧和排放特性。他们得出的普遍结论是：与汽油发动机相比，氨燃料发动机具有更低的热效率、更多的未燃氨排放、更高的氮氧化物排放、更低的燃烧效率和更小的输出功率。有研究6-10采用柴油引燃氨和空气的预混气，在燃烧后产生的废气中也发现较多的未燃氨。由于化学反应性低、火焰传播速度慢、可燃性范围窄等特点，在燃烧室的缝隙区和边界区存在氨的藏匿，在预混燃烧模式下难以充分燃尽。因此，提高燃烧效率、减少未燃氨逃逸是一个必要的研究课题。为了发动机的高效、清洁运行，有必要研究和开发合适的氨燃烧方式。通过缸内高压直喷的方式来实现液氨的扩散燃烧在提高燃烧效率、减少氨逃逸方面具有较大的潜力，同时可以弥补氨能量密度低的缺点，提高发动机功率。FRANKL等11对单缸发动机的高压直喷双燃料（high pressure dual fuel，HPDF）模式进行数值研究，将氨作为主燃料直接喷射入发动机气缸，同时直喷柴油作为引燃燃料。他们发现在HPDF模式下氨可以高效稳定地燃烧。通过数值仿真对氨在 HPDF 模式下和预混合低压双燃料（low pressure dual fuel，LPDF）模式下的燃烧状况进行比较12，结果表明：在 HPDF模式下的燃烧在NOx排放、温室气体排放、氨逃逸和燃烧效率等方面均优于LPDF模式。上述发动机整机的模拟工作表明氨燃料发动机在HPDF模式下具有巨大的潜力，但目前关于高压氨燃料喷雾混合气形成和燃烧过程的基础研究比较缺乏，而这对高压氨燃料喷雾燃烧机理的理解至关重要。OKAFOR等13-14研究在旋流燃烧器和微型燃气轮机中液氨喷雾的燃烧特性，发现随着氨质量分数和旋流器入口速度的增加，火焰稳定性下降。由于瞬态蒸发量和氨的蒸发潜热较大，周围环境气体温度骤降约50，大幅度降低氨燃烧的稳定性。该研究提出通过预热旋流空气来提高火焰稳定性。PELE 等15研究汽油直喷（gasoline direct injection，GDI）发动机条件下的液氨喷雾特性，并将结果与乙醇和汽油喷雾进行比较，结果表明：氨喷雾的宏观几何形状与其他燃料不同，氨喷雾对周围环境密度和温度的敏感性高于乙醇和汽油喷雾，这可能是由于氨喷雾在较低的环境压力下沸点较低而产生闪沸所致。如上所述，目前发表的文献集中研究在较低的环境压力和温度以及相对较低的喷射压力下氨燃料喷雾的特性。本课题组16利用自主设计开发的高压液氨直喷系统进行试验，探究在柴油机条件下氨燃料喷雾燃烧的特性。在室温和高温惰性条件下，捕捉氨喷雾的宏观特性。另外，在 900 K和 1 000 K环境温度以及氧体积分数为21%的环境条件下建立高压氨燃料喷雾火焰，试验数据（如气液相贯穿距等）为高压液氨喷雾模型的标定提供有力支撑，以期通过仿真模拟的方式更好地理解氨燃料喷雾混合气形成和燃烧的过程。本研究旨在开发适用于液氨喷雾的仿真模型，分析高压氨喷雾混合气形成和燃烧的特性。1试验与仿真方法1.1试验装置与方法试验在预燃式定容弹中进行。在定容弹中建立类似于柴油机上止点（top dead center，TDC）的高温高压条件，营造出的高温惰性环境以及有氧环境可分别用于液氨喷雾蒸发和燃烧的研究。通过高速相机记录下液氨喷雾形态，并通过后处理提取相关特性参数。喷油器安装在定容弹顶部，喷孔大小为0.22 mm。液氨通过自主设计开发的高压喷射系统加压后经喷油器射入定容弹内。图 1（a）为试验所用的纹影法光路。如图 1（a）所示，LED光源经小孔射入后作为点光源，光经过凹面镜后穿过容弹中喷雾，将气相和液相喷雾形态传递至高速相机（NAC Memrecam HX-6）。另外，同样采用上述Z型光路捕捉喷雾火焰形态，同时在相机前添加中性滤片以避免过度曝光。图1（b）为试验所用的背景光散射法光路。如图1（b）所示，LED光源发出的光束经散射膜产生漫散射2023 年 5 月 3 王宁等：高压液氨喷雾燃烧的数值模拟和试验验证背景，利用高速摄像机（Phantom v2012）捕捉喷雾液相的图像。在上述2种测试方法中测试用的高速相机的曝光时间为1 s，分辨率为512 512，拍摄速度为50 000 帧/s。1.2仿真模型与方法在CONVERGE软件中对喷雾蒸发和燃烧过程进行数值模拟。采用大涡模拟（large eddy simulation，LES）方法17求解湍流，液滴破碎采用KH-RT模型18，液滴蒸发采用Frossling模型19，液滴碰撞采用 NTC 模型20。另外，采用 Shell 着火模型21和特征时间燃烧（characteristic time combustion，CTC）22模型对喷雾燃烧过程进行初步计算。在 仿 真 时 将 计 算 域 简 化 为 48 mm48 mm97 mm 的长方体，如图2所示。整个计算区域的基础网格设置为 1 mm，最小网格尺寸为 0.125 mm。在喷雾发展的区域采用固定网格细分，根据速度梯度和温度梯度采用3级自适应网格细分，最大网格数量为11 704 658。在本研究中，喷雾蒸发和燃烧过程模拟到开始喷射后 2.0 ms。关于数值模型的相关构造和验证的更多细节可以在文献23中找到。仿真与试验条件如表 1 所示。在喷射压力为60 MPa、初始液氨温度为350 K条件下，探究不同环境温度（800 K、900 K、1 000 K）和环境密度（15 kg/m3、18 kg/m3、21 kg/m3）对液氨喷雾蒸发特性的影响。此外，在环境密度为18 kg/m3、环境温度为900 K、环境氧体积分数为21%条件下初步模拟液氨喷雾燃烧的过程。2结果与讨论2.1模型标定与网格尺寸敏感性验证根据气液相贯穿距的试验数据，对上述模型所涉及的参数进行标定，并进行网格尺寸敏感性研究，以确定在喷射压力为60 MPa、环境温度为900 K、环境密度为18 kg/m3、环境氧体积分数为0条件下的最佳网格尺寸。不同网格尺寸模拟结果对比图如图3所示，使用优化后的模型参数得到的模拟结果与试验结果相符。从图3可以看出，当最小网格尺寸小于0.125 mm时，模拟结果几乎与网格尺寸无关。因此，为了在可接受的成本下提高计算精度，以下的计算均采用0.125 mm网格尺寸。2.2喷油速率修正在喷雾模拟时须输入喷油速率的数据，本研究中的喷油速率通过Boach长管法24测得。另外利用课题组开发的喷油速率修正模型23对原始喷油速率进行修正，考虑Bosch长管法与实际喷射时喷嘴内的流动差异，原始和修正后的喷油速率（fuel injection rate，FIR）曲线如图4所示。图5为喷油速率修正前后模拟结果对比，修正后的