内燃机学报

参数	数值
喷射气体	H₂
环境气体	N₂
喷射压力p_a/MPa	0.4，0.8，1.2，2.0，4.0，6.0
喷射对比压力p_ar	0.31，0.62，0.93，1.54，3.08，4.63
环境压力p_b/MPa	0.1，0.5，1.0，1.5，2.0
环境对比压力p_br	0.08，0.39，0.77，1.16，1.54
试验温度/K	298
孔径/mm	0.35
喷射脉宽/ms	4
分辨率	992×744

参数	数值
喷射气体	H₂
环境气体	N₂
喷射压力p_a/MPa	0.4，0.8，1.2，2.0，4.0，6.0
喷射对比压力p_ar	0.31，0.62，0.93，1.54，3.08，4.63
环境压力p_b/MPa	0.1，0.5，1.0，1.5，2.0
环境对比压力p_br	0.08，0.39，0.77，1.16，1.54
试验温度/K	298
孔径/mm	0.35
喷射脉宽/ms	4
分辨率	992×744

参数	数值
p₀/MPa	0.4	0.8	1.2	2.0	4.0	6.0
ρ_a/(kg·m^-3)	0.235	0.648	0.969	1.608	3.179	4.713
m/(g·s^-1)	0.024	0.047	0.071	0.118	0.235	0.350
P/(kg·m·s^-2)	0.028	0.057	0.085	0.142	0.283	0.422
m^*/(g·s^-1)	0.021	0.040	0.060	0.101	0.199	0.298
P^*/(kg·m·s^-2)	0.024	0.049	0.073	0.121	0.240	0.358

参数	数值
p₀/MPa	0.4	0.8	1.2	2.0	4.0	6.0
ρ_a/(kg·m^-3)	0.235	0.648	0.969	1.608	3.179	4.713
m/(g·s^-1)	0.024	0.047	0.071	0.118	0.235	0.350
P/(kg·m·s^-2)	0.028	0.057	0.085	0.142	0.283	0.422
m^*/(g·s^-1)	0.021	0.040	0.060	0.101	0.199	0.298
P^*/(kg·m·s^-2)	0.024	0.049	0.073	0.121	0.240	0.358

参数	数值
环境压力/MPa	0.1	0.5	1.0	1.5	2.0
氮气密度/(kg·m^-3)	1.131	5.658	11.326	17.001	22.682

参数	数值
环境压力/MPa	0.1	0.5	1.0	1.5	2.0
氮气密度/(kg·m^-3)	1.131	5.658	11.326	17.001	22.682

不同临界状态下氢气射流与混合特性

PDF下载

李欣 , 陈皓 , 罗林雷 , 马志豪 , 王鑫

内燃机学报 | 2025,43(1): 38-46

收起

内燃机学报 | 2025 , 43 (1) : 38 -46

不同临界状态下氢气射流与混合特性

全屏

[Author(id=1325528793892537011, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, orderNo=0, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=210321030265@stu.haust.edu.cn, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325528793997394613, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528793892537011, language=EN, stringName=Xin Li, firstName=Xin, middleName=null, lastName=Li, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325528794244858550, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528793892537011, language=CN, stringName=李欣, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=河南科技大学车辆与交通工程学院，洛阳 471003, bio={"content":"

李欣，硕士研究生，210321030265@stu.haust.edu.cn．

"}, bioImg=null, bioContent=

李欣，硕士研究生，210321030265@stu.haust.edu.cn．

, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325528793787679407, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, xref=null, ext=[AuthorCompanyExt(id=1325528793800262320, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China), AuthorCompanyExt(id=1325528793812845233, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=河南科技大学车辆与交通工程学院，洛阳 471003)])]), Author(id=1325528794307773112, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, orderNo=1, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325528794391659194, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528794307773112, language=EN, stringName=Hao Chen, firstName=Hao, middleName=null, lastName=Chen, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325528794458768059, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528794307773112, language=CN, stringName=陈皓, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=河南科技大学车辆与交通工程学院，洛阳 471003, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325528793787679407, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, xref=null, ext=[AuthorCompanyExt(id=1325528793800262320, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China), AuthorCompanyExt(id=1325528793812845233, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=河南科技大学车辆与交通工程学院，洛阳 471003)])]), Author(id=1325528794538459837, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, orderNo=2, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325528794626540223, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528794538459837, language=EN, stringName=Linlei Luo, firstName=Linlei, middleName=null, lastName=Luo, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325528798195892928, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528794538459837, language=CN, stringName=罗林雷, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=河南科技大学车辆与交通工程学院，洛阳 471003, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325528793787679407, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, xref=null, ext=[AuthorCompanyExt(id=1325528793800262320, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China), AuthorCompanyExt(id=1325528793812845233, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=河南科技大学车辆与交通工程学院，洛阳 471003)])]), Author(id=1325528798292361922, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, orderNo=3, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=mazhihao@haust.edu.cn, emailSecond=null, emailThird=null, correspondingAuthor=1, authorType=1, ext={EN=AuthorExt(id=1325528798380442308, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528798292361922, language=EN, stringName=Zhihao Ma, firstName=Zhihao, middleName=null, lastName=Ma, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325528798451745477, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528798292361922, language=CN, stringName=马志豪, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=河南科技大学车辆与交通工程学院，洛阳 471003, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325528793787679407, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, xref=null, ext=[AuthorCompanyExt(id=1325528793800262320, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China), AuthorCompanyExt(id=1325528793812845233, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=河南科技大学车辆与交通工程学院，洛阳 471003)])]), Author(id=1325528798514660039, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, orderNo=4, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1325528798623711945, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528798514660039, language=EN, stringName=Xin Wang, firstName=Xin, middleName=null, lastName=Wang, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1325528798699209418, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, authorId=1325528798514660039, language=CN, stringName=王鑫, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=河南科技大学车辆与交通工程学院，洛阳 471003, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1325528793787679407, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, xref=null, ext=[AuthorCompanyExt(id=1325528793800262320, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China), AuthorCompanyExt(id=1325528793812845233, tenantId=1146029695717560320, journalId=1205116786560364546, articleId=1325528789853422232, companyId=1325528793787679407, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=河南科技大学车辆与交通工程学院，洛阳 471003)])])]

李欣, 陈皓, 罗林雷, 马志豪, 王鑫

作者信息

河南科技大学车辆与交通工程学院，洛阳 471003

通讯作者:

马志豪，博士，教授，mazhihao@haust.edu.cn．

作者简介:

李欣，硕士研究生，210321030265@stu.haust.edu.cn．

Hydrogen jet and mixing characteristics under different critical states

Xin Li, Hao Chen, Linlei Luo, Zhihao Ma, Xin Wang

Affiliations

Vehicle and Transportation Engineering College，Henan University of Science and Technology，Luoyang 471003，China

出版时间: 2025-01-25 doi: 10.16236/j.cnki.nrjxb.202501005

文章导航

摘要

收起

在内燃机气道和缸内环境下喷射氢气时，喷射压力和环境压力之间会出现高超临界向超临界以及超临界向亚临界的转变，采用高速纹影成像技术，在定容燃烧弹中研究了不同条件下的喷射特性．结果表明：不同条件下氢气射流的贯穿距、射流表面积以及射流体积随时间的增长趋势不同，但射流锥角在1ms后随时间变化不大，差值稳定在2°以内；在超临界喷射进入超临界环境时射流与周围气体相互作用使湍流强度增强，这有利于夹带更多的周围气体．最后提出的氢气射流贯穿距经验公式具有较高的精度，可以利用这些结果优化氢内燃机混合气的形成和燃烧，为喷氢策略提供参考．

关键词

亚/超临界 / 氢气射流 / 混合特性 / 卷吸比

Abstract

收起

When injecting hydrogen gas into the intake air path and in-cylinder environment of an internal combustion engine, there is a transition from supercritical to subcritical between the injection pressure and ambient pressure. High-speed schlieren imaging technology was used to study the injection characteristics under different conditions in a constant volume combustion bomb. The results show that the penetration distance, jet surface area, and jet volume of the hydrogen jet vary with time under different conditions. However, the change in the jet cone angle after 1ms is not significant, and the difference remains within 2°. When the supercritical injection enters the supercritical environment, the interaction between the jet and surrounding gases enhances turbulence intensity, which facilitates the entrainment of more surrounding gases. An empirical formula for the jet penetration distance of the hydrogen jet was proposed, which has proved to be high accuracy. These results can be utilized to optimize the mixture formation and combustion of hydrogen internal combustion engines, providing references for hydrogen injection strategies.

Key words

sub/supercritical / hydrogen jet / mixing characteristics / entrainment ratio

引用本文

李欣, 陈皓, 罗林雷, 马志豪, 王鑫. 不同临界状态下氢气射流与混合特性. 内燃机学报, 2025 , 43 (1) : 38 -46 . DOI: 10.16236/j.cnki.nrjxb.202501005

Xin Li, Hao Chen, Linlei Luo, Zhihao Ma, Xin Wang. Hydrogen jet and mixing characteristics under different critical states[J]. Transactions of CSICE, 2025 , 43 (1) : 38 -46 . DOI: 10.16236/j.cnki.nrjxb.202501005

正文

收起

目前，社会经济的飞速发展加速了化石燃料消耗，导致CO₂排放大幅度增加，全球气候变暖现象加剧，而交通运输产生的温室气体排放总量稳居第一^[1]．由于低碳和可持续发展的要求，高效、清洁的可再生替代燃料受到了各国广泛关注．H₂作为无碳能量载体，燃烧不会产生任何碳排放，此外H₂来源广泛，能够利用化石燃料和可再生能源(生物质、风能和太阳能)等方法制取^[2]，其中利用光伏系统和风电系统获取电力，再通过电解水获取的“绿氢”正成为全球氢能发展的焦点^[3]，为目前全球环境问题和燃料消耗提供了解决方案．H₂的能量密度(120MJ/kg)约是汽油(43.5MJ/kg)的3倍^[4]，被认为是未来最重要的清洁替代燃料方案之一．

氢内燃机热效率高、排放污染较小，但混合气的均匀程度对其性能影响显著^[5]，因而研究氢气射流与混合特性显得尤为重要．刘政等^[6]基于纹影法研究了氩-氧氛围中不同喷射压力和环境压力对氢气射流贯穿距、射流锥角及卷吸率发展的影响，发现对于一定的喷射压力，随着环境压力的升高卷吸率增加，但环境压力升高到一定值后卷吸率反而会减小．Shahsavan等^[7]模拟了氢气在不同比热比的环境气体下射流的发展过程和混合特性表明，相比于氢气与氮气或氙气，氢气在氩气中具有更快的混合速率．Deng等^[8]针对氩循环氢发动机，在定容燃烧弹中对比氩气氛围下氢气和氧气的射流发展时发现，氢气直喷策略优于氧气直喷策略．邓俊等^[9]研究了不同喷射压力和喷射温度对氢气射流特性和卷吸特性的影响表明，氢气对周围气体的卷吸能力会随压力比的增加呈先增后减的趋势，卷吸比随温度的升高而显著下降．

以上研究大多为氢气在不同喷射压力和环境压力情况下的射流发展及混合特性，而在内燃机环境下会出现超临界氢气向亚临界环境以及超临界氢气向超临界环境下的喷射，氢气在亚超临界的特性不同，必然对其喷射特性产生影响，从而影响一定时间内氢气与环境气体的混合，而大部分研究并未考虑亚超临界状态对喷射特性及环境气体混合的影响．Ma等^[10]利用定容燃烧弹研究了液体燃料在亚/超临界环境条件下喷射过程表明，超临界环境与亚临界环境下的喷射在形态学特征上存在较大差异．张纪宝^[11]开展了正庚烷超临界喷射混合特性的模拟研究，得出超临界压力下正庚烷的热力参数随温度和压力呈连续非线性变化的规律．Gao等^[12]将亚/超临界条件下的航空煤油喷入静止的氮气中，观察发现超临界压力下出现激波结构且随着喷射压力的增加，马赫盘的尺寸及喷射锥角都会增加．Hou等^[13]利用模拟进行了超临界和亚临界环境压力下氢射流的分子动力学研究表明，超临界条件下射流径向密度变化率很小，超临界环境压力下氢射流中分子应力是亚临界的88.9%．

当前，超临界现象研究大多集中在亚/超临界喷雾和超临界液滴蒸发与燃烧的模拟，以及部分关于超临界或超超临界锅炉的排放特性研究，超临界气体射流的相关试验鲜见报道．此外，氢气射流研究中的喷射压力范围多为超临界，环境压力范围多为亚临界．因此，笔者开展氢气亚临界喷射进入亚临界环境(简称“亚-亚”)、超临界喷射进入亚临界环境(简称“超-亚”)及超临界喷射进入超临界环境(简称“超-超”)下射流现象学研究是必要的，可以为氢气在氢内燃机上的应用提供理论指导，同时为模拟研究提供试验依据．

1　试验方案

收起

1.1　试验系统设计

试验采用如图1所示的定容燃烧弹系统，主要包括定容燃烧弹、高速摄像机、光学纹影系统、喷射系统和喷射控制系统．

定容燃烧弹两侧有直径为100mm的光学石英玻璃观察窗，可以承受7MPa的环境压力．光学纹影系统由两面凹面镜和两面平面镜组成‘Z’字形光路，纹影法的原理是光强差异的偏转角与折射率梯度呈正比，通过流场的折射率反映密度变化．在试验系统调试中发现，周围气体扰动引起的密度变化会对射流本体的采集产生影响，使射流边缘难以识别，因而在试验开始前花费了大量的时间去调整背景的亮度与刀口的开度，在保证射流边缘可见的前提下，适当降低对比度，减小气流扰动造成的明暗变化，最终选择刀口的切割量为60%进行试验．选用NAC Memrecam GX-8高速摄相机，均在992×744的分辨率和5000幅/s的帧率下进行拍摄．喷射系统由气瓶、管路、压力表和喷嘴组成，分别利用安装在氢气高压管路和容弹主体上的高精度压力表控制喷嘴的喷射压力和容弹内氮气的压力．喷嘴在容弹内伸出的距离设置为6mm，以保证射流不受康达效应的影响．喷射控制系统利用法国EFS公司生产的8426和8427模块控制喷射规律和喷射脉宽．通过电脑控制端可以同时触发高速摄像机与单次喷射仪．

1.2　试验工况

氢内燃机有进气道喷射和缸内直喷两种喷射方式，其喷射压力和环境压力范围广、跨度大，因而笔者选择了表1的试验工况，喷射压力p_a与氢气临界压力p_cr之比定义为喷射对比压力p_ar，环境压力p_b与氢气临界压力p_cr之比定义为环境对比压力p_br，p_ar和p_br大于1时表示氢气处于超临界压力状态，使用的汽油直喷喷射器孔径为0.35mm．试验时固定喷射脉宽为4ms，所有试验工况点都满足环境压力小于喷射压力的0.528倍，因而出口流速始终保持在当地声速，具有最大的理论质量流量且与环境压力无关，仅受喷射压力和喷孔流通截面积的影响，见式(2)．

假设H₂是理想气体，氢气在喷嘴内的流动为绝热等熵过程，当环境压力p_b与喷射压力p_a之比小于0.528时，氢气流动在喷嘴内阻塞．在可压缩流体中，阻塞流条件下的临界压力比ε可以表示为

ε = p * p a = (2 k + 1) k k − 1

(1)

式中：p^*为临界条件下的环境压力；k为氢气的比热比，k＝1.41．

在临界条件下，喷嘴出口截面上气体流动速度达到当地声速，质量流量最大．此时喷嘴出口的质量流量m、动量P分别为

m = ρ u A = A (k ρ a p a) 1 / 2 (2 k + 1) k + 1 2 (k − 1)

(2)

P = m u = A k p a (2 k + 1) k k − 1

(3)

式中：A为喷嘴出口截面积；u为喷嘴出口流速；ρ_a为喷射压力下氢气的密度；

ρ

为气体密度．

将氢气的密度代入式(2)和式(3)得到不同喷射状态下氢气的质量流量m和动量P．考虑喷嘴内部的压力损失，Tsujimura等^[14]研究了不同喷射条件对喷嘴流量系数的影响表明，喷射脉宽大于5ms时流量系数与喷射脉宽无关，稳定在0.90附近；喷射脉宽小于5ms时由于针阀运动对流量影响较大，使得流量系数略微减小．因此，结合试验工况条件取喷嘴流量系数C_d为0.85，将m和P分别乘以C_d得到实际气体喷射下的质量流量m^*和动量P^*．氢气在不同状态下的喷射参数如表2所示．表3显示了在不同环境压力下氮气的密度．

1.3　图像处理

通过纹影法获得的氢气射流图像通过matlab程序进行批量处理，使用一系列图像形态学matlab函数提取射流轮廓，如图2所示．首先喷射序列中的每一张图像(图2(a))减去背景图像以突出显示射流轮廓，如图2(b)所示．随后利用滤波函数和二值化处理消除背景噪声，提高对比度，如图2(c)所示．填充射流二值化图像的空洞，如图2(d)所示．最后使用matlab开运算和闭运算重构射流图像，并利用canny算子获取射流边界轮廓，如图2(e)和图2(f)所示．

二值化处理中阈值R的选择会影响射流的边缘轮廓，进而会对试验结果产生一定的影响．图3显示在3种不同阈值下二值化后的射流轮廓．随着阈值增加会导致边缘轮廓的收缩，同时阈值过大也会导致边缘特征的丢失，为了保证试验结果的可靠性，选择二值化的阈值为0.40．

1.4　特征参数

通过对射流序列图像的后处理获得了射流贯穿距、射流锥角、射流表面积与射流体积等特征参数．

射流贯穿距定义为喷嘴到喷嘴所在的铅锤线与射流前锋面交点的距离，即喷嘴所在的像素点与喷嘴铅锤线与射流前锋面相交的像素点间的距离，通过像素比例尺换算得到射流贯穿距L，射流表面积为射流轮廓边缘所包含的所有像素点的面积S，如图4(a)所示．由于射流前锋面存在涡流的影响，射流形状变化大，射流锥角为射流外部边界的夹角，为了减小射流外围褶皱的影响，对射流主体小于L/2的部分采用了边界拟合的算法，即在纵向上每隔一定的像素点获取射流边界，在横向上的左边界与右边界点的像素点坐标，对获得的左边界点和右边界点分别使用polyfit函数拟合，这两条拟合线的夹角即为射流锥角，如图4(b)所示．射流体积定义为射流主体所占有的体积，假设射流在图片上每行像素点上是直径与宽度相等的微小圆柱体，对不同行像素点上的微小圆柱体的体积V_d叠加即为射流主体的体积V，如图4(c)所示．

为评估不同喷射状态下氢气对于周围气体的卷吸能力，用卷吸比ϕ评估氢气射流的混合特性^[15]．卷吸比定义为

ϕ = (V − V 0) / V 0

(4)

式中：V₀为氢气射流中氢气的体积，通过质量流量计算；V为射流主体的体积，通过纹影图片处理得到．

2　结果及分析

收起

2.1　氢气射流发展

图5为3种状态下射流发展随时间的变化．由于3种状态下喷射参数的不同，喷射时氢气的密度不同，周围气体密度也不同．超-超状态下氢气射流图像更加清晰，亚-亚与超-亚状态下，射流前锋区域较为模糊．不同状态下射流贯穿距均随时间增加，超-超状态下的射流锥角比其他两种喷射状态略微增大，射流近场区域边界较为平滑，而远场区域射流边界由于动量和能量交换以及热质传递产生褶皱．图5中亚-亚、超-亚和超-超状态下雷诺数分别为36280、173019和99813，雷诺数Re由式(5)计算而得．雷诺数越大，射流中从大涡到小涡之间的过渡越复杂，涡系结构越丰富，同时环境压力的增加会使环境密度增大，对射流的阻力也更大，褶皱更加明显，因而超-超状态下射流在0.8ms时射流前锋区域出现了非常明显的褶皱，其射流前锋的宽度更大；在0.8ms之后，超-超状态下射流前锋褶皱依旧明显．亚-亚状态下射流随射流贯穿距的增加，射流前锋轮廓更加不清晰，这主要是由于亚-亚状态下氢气密度小，且随着扩散与环境密度的差异更小，是光的折射减弱造成的．

R e = (ρ a u D e) / μ

(5)

D e = D ρ a / ρ b

(6)

式中：ρ_a为射流流体密度；u为射流速度；μ为流体动力黏度；D_e为当量直径；D为喷嘴孔径；ρ_b为环境密度．

2.2　亚-亚与超-亚对比分析

图6所示对比压力为0.31、0.62和0.93表示氢气亚-亚状态下射流宏观特性随时间的变化，对比压力为1.54、3.08和4.63表示氢气超-亚状态下射流宏观特性随时间的变化．

图6(a)为两种状态下射流贯穿距随时间的变化趋势．由于1ms以前喷射出的氢气量较少，后续从喷嘴喷出的氢气不断地给射流提供动量，推动前端射流向前发展，射流贯穿距离随时间增加而迅速增长，超-亚状态下氢气射流贯穿距的斜率大于亚-亚状态下射流贯穿距斜率，这是由于氢气进入相同的环境，超临界氢气相较亚临界具有更大压差且密度大，氢气射流膨胀程度更大，超临界氢气射流贯穿距略微大于亚临界射流．1ms后贯穿距随射流发展趋势变缓，主要有两方面原因造成：(1)随着时间的推移，射流前端氢气发生扩散并与氮气混合，氢气浓度逐渐减小，前端氢气携带的动量几乎全部耗散，此时射流前端整体质量增加，后续射流携带的动量需要推动更多的气体，进而减弱了射流的强度和贯穿能力；(2)低密度氢气射流在穿过高密度氮气环境的过程中会与氮气分子相互碰撞，从而转移动量和能量，导致射流的速度减小和能量损失，使得射流贯穿能力减弱．超-亚状态下氢气射流的密度较亚-亚状态下的大，质量流量随之增加，出口动量增大，后续射流可以对前端射流提供更多的能量传递，从而导致射流贯穿距的增加，3ms时超-亚状态下氢气射流贯穿距较亚-亚状态平均增加31.5%．

图6(b)为射流锥角随时间的变化趋势．射流锥角的变化趋势以1ms时为界可以分为两段：1ms前的射流初期由于针阀存在节流作用，喷射的氢气量较小时，氢气作为低密度气体，气体分子间距与液体相比要大的多，分子间作用力非常微小，分子可以自由运动，极易变形，因而射流锥角变化较大；当喷嘴流量稳定后，射流形态变得稳定，1ms后射流锥角均在各自的平均值±1°内浮动．在0.2ms时尽管由于存在针阀的运动导致了质量流量不稳定，但超-亚状态下的氢气密度和质量流量仍大于亚-亚状态，而两种状态下的射流轴向贯穿距相差3mm左右．因此，两种状态下质量流量的差异表现在射流锥角上．亚-亚状态下，喷射对比压力从0.31升高到0.93，射流锥角增加3°左右；而在超-亚状态下喷射对比压力从1.54增加到4.63时，射流锥角增加了1.5°左右．超-亚状态下平均射流锥角均大于亚-亚状态下的，这是由于超-亚状态下氢气射流密度大，后续射流对前端射流能量传递的更多，射流边界层与氮气的湍流程度也更大，卷吸周围气体的量随之增加，使得径向贯穿能力增加，导致射流锥角增大．

图6(c)为射流表面积随时间的变化趋势．不同状态下氢气射流的表面积均随时间增长．超-亚状态下氢气射流表面积比亚-亚状态下氢气射流表面积大，3ms时超-亚状态下氢气射流面积较亚-亚状态平均增加33.6%，是因为超-亚状态下氢气射流贯穿距和射流锥角均大于亚-亚状态．在亚-亚状态下，不同喷射压力下射流表面积曲线的斜率相差不大，而超-亚状态下氢气射流表面积曲线的斜率差异较大，可见超-亚状态下氢气射流表面积对喷射压力的改变较敏感，而亚-亚状态下射流表面积对喷射压力的改变不太敏感．

图6(d)为射流体积随着时间的变化趋势．在不同状态下氢气的射流体积随时间而增加．在亚-亚状态下，氢气射流体积增加缓慢，而超-亚状态下氢气射流体积迅速增加．超-亚状态下，超临界氢气会被压缩为更高密度的状态，当氢气射流进入氮气环境中时，由于压差会发生膨胀过程，一些试验和数值模拟^[16-17]已经证明了随着压差的增大，射流的膨胀程度更大，因而超临界氢气膨胀程度更大．同时超临界氢气的密度和黏性较大，与氮气接触的边界湍流强度更强，卷吸周围气体的量增加．因此，随着射流的发展，其射流体积迅速增加，3ms时超-亚状态下氢气射流体积较亚-亚状态平均增加200.7%．

2.3　超-亚与超-超对比分析

图7为超-亚与超-超状态下射流特征参数对比．其中：对比压力为0.08、0.39和0.77表示氢气超-亚状态下射流宏观特性随时间的变化；对比压力为1.16、1.54表示氢气超-超状态下射流宏观特性随时间的变化．

图7(a)为射流贯穿距随时间的变化趋势．射流贯穿距的发展同样呈现为两个阶段，在1ms以前射流贯穿距随时间迅速增加，随后贯穿距随时间的增长趋势变缓，此时超-超状态下氢气射流贯穿距的斜率小于超-亚状态下的氢气射流贯穿距斜率，这是因为超-超状态下，虽然喷嘴出口处氢气状态与超-亚状态相同，即后续射流对前端射流提供的动量相同，但由于容弹内环境状态不同，环境处于氢气的超临界条件，环境密度大，背压和分子间作用力加剧了氢气射流的阻力，氢气射流的动量耗散加快，后续射流更难推动前端射流向前发展．超-超状态下，超临界环境使得其贯穿能力明显小于超-亚状态下的贯穿能力．3ms时超-亚状态下氢气射流贯穿距较超-超状态平均增加47.6%．

图7(b)为射流锥角随时间的变化趋势．与图6(b)所述的亚-亚与超-亚状态下的对比相似，1ms前射流初期锥角由于针阀的运动射流锥角变化较大；1ms后射流后期锥角都在各自的平均值±1°内浮动．由于氢气射流进入超临界环境时环境密度增加导致受到的阻力增加，其轴向穿透能力减弱，更多向径向发展，超-超状态下氢气射流锥角大于超-亚状态下的锥角．超-亚状态下，环境对比压力从0.08升高到0.77，射流锥角增加5°左右；而在超-超状态下环境压力从1.16增加到1.54时，射流锥角增加了1°左右．可以发现，超-超状态下氢气射流锥角对环境变化不敏感，而超-亚状态下氢气射流锥角对环境变化较为敏感．

图7(c)为射流表面积随时间的变化趋势．超-亚与超-超状态下氢气射流的表面积均随时间增大．但超-亚状态下氢气射流表面积比超-超状态下氢气射流表面积大，是因为超-亚状态与超-超状态相比，氢气喷射时的热力学性质相同，然而环境从亚临界过渡到超临界，氢气在高密度介质中穿透能力减弱，限制了氢气射流的扩散．超-亚状态下射流贯穿距大、射流锥角小，其射流面积增加，说明射流面积受贯穿距的发展影响最大．3ms时超-亚状态下氢气射流表面积较超-超状态平均增加51.5%．

图7(d)为射流体积随着时间的变化趋势．在不同状态下氢气射流体积随时间快速增加．超-亚与超-超状态相比，超-亚状态下射流体积大于超-超状态下的射流体积，其原因类似于两种状态下氢气射流表面积的发展．3ms时超-亚状态下氢气射流体积较超-超状态平均增加130.6%．

2.4　卷吸比

图8为氢气射流在3种状态下的卷吸比．实际喷射过程中，在针阀开启初期流通截面积较小，导致射流初期实际的质量流量比计算所得的结果小，因而使得卷吸比出现了负值．但实际针阀运动的时间极短，不影响对3种状态下卷吸比的比较分析．

图8(a)显示了卷吸比随时间的变化趋势，在1ms以前，3种状态下卷吸比均随射流发展增加较快，随后缓慢增加．与贯穿距随时间的变化趋势相同，氢气刚流出喷嘴时具有较大的动量，与周围气体相互作用引起的湍流强度较大，因而氢气射流可以夹带较多的周围气体．随着射流与周围气体间的动量交换，射流与周围气体之间的速度差减小使得湍流强度减小，夹带能力减弱．喷射条件从亚临界到超临界时，氢气射流的卷吸比增大，是因为超临界氢气具有更大的密度，当喷射进入环境中时，会产生更大的膨胀和加速效果，与环境气体的接触面积较大，可以携带更多的环境气体．环境条件从亚临界到超临界时，氢气射流的卷吸比同样增加则是因为环境气体的压力升高会使其密度增加，射流中的高速气流会与高压环境气体相互作用，由于密度差异更大，在相互作用区域会产生更强的湍流和涡旋结构，将周围环境气体混合到射流中，从而增加了环境气体的携带量．在3ms时，超-超状态下氢气射流卷吸比相较亚-亚状态增加了349%，较超-亚状态增加了109%．化学计量比下氢气与空气的体积比为2.38，则理论空燃比下氢气射流的卷吸比为1.38．虽然亚-亚状态下卷吸比较低，但仍能够携带足够的空气．在缸内直喷氢内燃机中，于压缩上止点附近喷入氢气，其喷射条件能满足超-亚或超-超状态，有利于实现氢气的稀薄分层燃烧．

图8(b)为3种状态下卷吸比随贯穿距的变化趋势．卷吸比随贯穿距呈线性增长，贯穿距越大，卷吸比越大．随着喷射条件和环境条件从亚临界过渡到超临界，卷吸比随贯穿距增长的斜率更大．随着时间的增加卷吸比增长趋势变缓，但随着贯穿距的增长依然保持相同的斜率．在环境条件从亚临界到超临界时，尽管超-超状态下贯穿距减小了，但仍具有较好的卷吸能力，因而缸内直喷氢内燃机适当提高缸内压力，有利于氢气射流卷吸更多的空气．

2.5　氢气射流贯穿距经验公式

Turner^[18]在1962年首次提出了欠膨胀气体射流的数学模型，该模型把气体射流的几何形状近似为球形涡头和传递动量的准稳态射流区域．后来Hill等^[19]在Turner^[18]的基础上，通过圆孔喷嘴天然气射流的试验，提出气体射流的贯穿距是孔径、时间、出口总动量流率和周围气体密度的函数，并建立了射流贯穿距的数学公式，如式(7)所示．

y = Γ (M n / ρ 0) 0.25 t 0.5

(7)

式中：y为射流贯穿距；

Γ

为贯穿常数，由试验确定，Hill等^[19]的试验中贯穿常数

Γ

为3.0±0.1；

M n

为动量流率，即为表2中的P^*；

ρ 0

为周围气体密度；t为时间．

Chitsaz等^[20]研究了孔径为0.5～1.0mm和压比(喷射压力/环境压力)在1.5～36.0下的氦气射流贯穿距，得出

Γ

在1.89～2.95内．通过以上研究发现贯穿常数

Γ

由于喷射条件、喷嘴结构和孔径的不同并不是固定值，且大部分工作都是针对天然气进行的．因此，基于现有的试验数据验证式(7)对氢气射流的适用性．

图9给出了y/

(M n / ρ 0) 0.25

随t^0.5的变化规律，可见量纲为1的贯穿距曲线均呈直线，直线的斜率即为

Γ

．但是将

Γ

代入式(7)中得到贯穿距曲线整体上比试验值偏大，笔者根据试验数据提出了氢气喷射的贯穿距经验公式为

y = Γ (M n / ρ 0) 0.25 t 0.5 + a

(8)

式中：a为修正常数，-0.004<a<-0.002．

图10显示射流贯穿距试验值与拟合值的对比．在环境温度为298K、喷嘴孔径为0.35的试验条件下提出的氢气射流贯穿距经验公式具有较好的精度，贯穿常数

Γ

为2.55±0.2．对氢气射流贯穿距经验式(8)中的

Γ

和a是否与氢气浮力有关以及对式(8)的优化将在以后的工作中深入讨论．

3　结论

收起

(1) 不同喷射状态下，射流贯穿距均随时间增加，但射流锥角的变化不大，射流近场区域边界较为平滑，而远场区域射流边界由于与周围气体的动量交换变得更加皱褶．

(2) 超-亚与亚-亚状态相比，3ms时超-亚状态下氢气射流贯穿距较亚-亚状态增加31.5%，射流面积增加33.6%，射流体积增加200.7%；射流锥角的差值在误差带范围内，两种状态下相差不大；氢气超临界喷射有利于增大射流贯穿距．

(3) 超-亚与超-超状态相比，3ms时超-亚状态下氢气射流贯穿距较超-超状态增加47.6%，射流面积增加51.5%，射流体积增加130.6%，但超-超状态下射流锥角较超-亚状态增大23.8%；当环境从亚临界到超临界状态时，射流在径向上的发展较为明显．

(4) 对于3种状态下卷吸比的对比，发现超-超状态下其射流卷吸周围气体的能力更强，在3ms时卷吸比较超-亚状态增加了109%，较亚-亚状态增加了349%，由此可见超临界环境有利于氢气射流混合；在试验条件下提出了关于氢气射流贯穿距经验公式，引入了修正常数使其具有较高的精度．

基金

收起

国家自然科学基金资助项目(51906061)
龙门实验室资助项目(LMQYTSKT038)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

［1］

Verhelst

，Wallner

. Hydrogen-fueled internal combustion engines[J]. Progress in Energy and Combustion Science，2009，35(6)：490-527.

［2］

Zamfirescu

，Dincer

. Ammonia as a green fuel and hydrogen source for vehicular applications[J]. Fuel Processing Technology，2009，90(5)：729-737.

［3］

Mazzeo

，Herdem

M S

，Matera

，et al. Green hydrogen production：Analysis for different single or combined large-scale photovoltaic and wind renewable systems[J]. Renewable Energy，2022，200：360-378.

［4］

Salvi

B L

，Subramanian

K A

. Sustainable development of road transportation sector using hydrogen energy reviews[J]. Renewable and Sustainable Energy Review，2015，51：1132-1155.

［5］

Faizal

，Chuah

L S

，Lee

，et al. Review of hydrogen fuel for internal combustion engines[J]. J Mech Eng Res Dev，2019，42(3)：36-46.

［6］

刘政，钟慧平，龚寅春，等. 氩-氧氛围中氢气射流喷射与混合特性的研究[J]. 内燃机学报，2020，38(2)：110-118.

［7］

Shahsavan

，Morovatiyan

，Mack

J H

. A numerical investigation of hydrogen injection into noble gas working fluids[J]. International Journal of Hydrogen Energy，2018，43(29)：13575-13582.

［8］

Deng

，Zhong

，Gong

，et al. Studies on injection and mixing characteristics of high pressure hydrogen and oxygen jet in argon atmosphere[J]. Fuel，2018，226：454-461.

［9］

邓俊，钟慧平，龚学海，等. 环境气体对氢气喷射与卷吸特性的影响[J]. 内燃机学报，2019，37(1)：68-75.

［10］

，Li

，et al. Evaporation and combustion characteristics of hydrocarbon fuel droplet in sub- and super-critical environments[J]. Fuel，2018，220：763-768.

［11］

张纪宝. 发动机燃油超临界喷射及流动特性研究[D]. 北京：北京交通大学机械与电子控制工程学院，2021.

［12］

Gao

，Lin

，Hui

，et al. Injection characteristics of near critical and supercritical kerosene into quiescent atmospheric environment[J]. Fuel，2019，235：775-781.

［13］

Hou

，Ding

，Yu

，et al. Comparative molecular dynamics study of liquid hydrogen annular jets in subcritical and supercritical environments[J]. International Journal of Hydrogen Energy，2023，51：977-989.

［14］

Tsujimura

，Mikami

，Achiha

，et al. A study of direct injection diesel engine fueled with hydrogen[C]// SAE 2003 World Congress and Exbition，Warrendale，PA，USA，2003：390-405.

［15］

Wang

，Sun

，Luo

，et al. Visualization research on hydrogen jet characteristics of an outward-opening injector for direct injection hydrogen engines[J]. Fuel，2020，280：118710.

［16］

Yip

H L

，Srna

，Liu

，et al. Visualization of hydrogen jet evolution and combustion under simulated direct injection compression ignition engine conditions[J]. International Journal of Hydrogen Energy，2020，45(56)：32562-32578.

［17］

Franquet

，Perrier

，Gibout

，et al. Free underexpanded jets in a quiescent medium：A review[J]. Prog Aero Sci，201，77：25-53.

［18］

Turner

J S

. The starting plume in neutral surroundings[J]. J Fluid Mech，1962，13：356-68.

［19］

Hill

P G

，Ouellette

. Transient turbulent gaseous fuel jets for diesel engines[J]. J Fluids Eng，1999，121：93-101.

［20］

Chitsaz

，Saidi

M H

，Mozafar

A A

，et al. Experimental and numerical investigation on the jet characteristics of spark ignition direct injection gaseous injector[J]. Appl Energy，2013，105：8-16.

2025年第43卷第1期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.16236/j.cnki.nrjxb.202501005

接收时间：2024-03-24
首发时间：2026-11-06
出版时间：2025-01-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2024-03-24
修回日期：2024-09-18

基金

国家自然科学基金资助项目(51906061)

龙门实验室资助项目(LMQYTSKT038)

作者信息

河南科技大学车辆与交通工程学院，洛阳 471003

通讯作者:

马志豪，博士，教授，mazhihao@haust.edu.cn．

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/nrjxb/CN/10.16236/j.cnki.nrjxb.202501005

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT