Infrared molecular tagging velocimetry
-
摘要: 分子标记测速法(MTV)和粒子成像测速法(PIV)常被用于流动显示和流场成像测量,但在示踪粒子跟随性差、示踪粒子分布不均匀时,示踪粒子的引入会给PIV带来速度测量系统误差,而不需引入示踪粒子的MTV因荧光寿命长度限制,主要应用在高速和超声速流动测量中。为了发展无需示踪粒子、可适用于低速气流场的二维流速成像方法,本文介绍了一种基于激光诱导红外荧光的新型分子标记测速法,并在二氧化碳气体轴对称湍流射流中进行了速度测量与验证。在红外分子标记测速法中,通过红外脉冲激光选择性激发气体小分子的共振振动能级跃迁实现分子的标记,随后通过红外相机对不同时刻下跟随流场流动的激发态分子记录其发射荧光分布,进而处理得到流动速度场信息。通过考虑分子振动能量传递过程模型、有限荧光寿命、横向速度分量和分子扩散运动对荧光分布的影响,实现从荧光分布图像定量获取速度场分布。将该方法应用于5~51 m/s速度的二氧化碳湍流射流中,得到了射流轴向速度的径向分布,速度测量的相对不确定度优于8%,径向空间分辨率达到107 μm,且该速度分布与湍流射流理论结果及前人实验测量结果符合较好。利用该方法分辨了射流在不同轴向位置的径向速度分布,观察到了射流从势核区到充分发展区的演变。利用该方法可以获得低速气流场的高分辨瞬时速度成像,后续通过提高红外激光的脉冲功率、激发效率和重复频率,可以进一步提高该方法的测量精度、空间分辨率和时间分辨率,从而有潜力在近壁面流动、微尺度流动和大梯度流动等粒子示踪不适合用场景中提供定量流场测量的有效途径。Abstract: Molecular Tagging Velocimetry (MTV) and Particle Imaging Velocimetry (PIV) are often used for flow visualization and velocity field imaging. However, the requirement for tracer particles can bring systematic errors to the velocity measurement of PIV method when the tracer particles have poor followability and uneven distribution. In the case of MTV, although particle seeding is not required, the finite fluorescence lifetime of the tracer molecules typically constrains its use only in high-speed and supersonic flows. To develop a velocity field imaging method with no requirements of tracer particles and suitable for low-speed flows, a novel MTV method based on infrared (IR) laser-induced fluorescence is developed and verified in axisymmetric turbulent jet of carbon dioxide. Resonant vibrational transition of the small gas molecule is selectively excited by an infrared pulsed laser to achieve molecular tagging, and the fluorescence distributions of the excited molecules at different instants are then imaged by an infrared camera, from which the velocity distributions are deduced. The effects of the molecular vibrational energy transfer process model, finite fluorescence lifetime, lateral velocity component and molecular diffusion motion on the fluorescence distribution are analyzed to improve the accuracy of the velocity measurement. The proposed method has been successfully verified in carbon dioxide turbulent jets with velocities ranging from 5 m/s to 51 m/s, and the radial distribution of the axial velocity in the main region of the jet is measured. The radial spatial resolution can reach 107 microns, and the velocity distribution is consistent with the theoretical calculation of turbulent jets and previous experimental results. The relative uncertainty of velocity measurement is better than 8%. This method can be used to obtain high-resolution instantaneous velocity imaging of low-speed flow field. Subsequently, by improving the pulse power, excitation efficiency and repetition frequency of the infrared laser, the measurement accuracy, spatial resolution and temporal resolution of this method can be further improved. Therefore, the proposed method bears great potential to provide a quantitative velocity field imaging method in the near-wall flow, micro-scale flow and large gradient flow where it is difficult to introduce tracer particles.
-
0 引 言
高超声速飞行器在大气层内飞行,其热环境极其严酷,防热材料及热结构需承受极大的热力载荷,尤其是飞行器前体、进气道压缩面、发动机唇口、舵翼等因局部热流高、温度梯度大、噪声、振动等造成热力耦合失效的风险极大,因此,必须在地面开展相应的考核试验[1-4]。电弧风洞[5-6]是当前能够在地面近似模拟飞行热环境参数(焓值、热流、压力、气体组分等)的试验设备,但电弧风洞的功率有限,即使是50 MW量级的设备也仅能模拟考核飞行器的局部构件。因此,通常截取飞行器典型部位置于电弧风洞中进行试验,如平板类部件一般采用矩形或半椭圆喷管直边与试验件对接,让喷管的超声速气流直接过渡到平板试验件上进行试验考核[7-10],如图1所示。
目前矩形或半椭圆喷管的壁面通常为夹层水冷结构。由于喷管喉道处热流很高,一般采用紫铜或铜合金材料制作并对其进行强制冷却,这会导致近壁面气流焓值降低。如果喷管扩张段采用水冷壁,将会造成试验件表面气流恢复焓降低,使风洞试验考核与真实飞行之间存在显著差距。如式(1)、(2)所示,高温气流对模型壁面产生的热流为:
$$ q = a({h_{\text{r}}} - {h_{\text{w}}}) $$ (1) $$ {h_{\text{r}}} = {h_{\rm{e}}} + r{v^2}/2 $$ (2) 式中:
$ q $ 为热流,$ a $ 为对流换热系数,$ {h_{\text{r}}} $ 为气流恢复焓,${h_{\rm{w}}}$ 为壁面焓,$ {h_{\text{e}}} $ 为气流静焓,$ v $ 为气流速度,$ r $ 为恢复系数。在喷管壁面水冷的条件下,气流的静焓降低,而在喷管壁面隔热的条件下,喷管壁面温度升高,减少了焓降,这与飞行器防热材料表面高温气体流动相似。研究喷管冷壁及热壁对试验件考核的影响,缩小“天地”差距,对防热系统的精确考核和精细设计具有重要的意义。
目前,世界上拥有电弧风洞的机构较少,喷管焓降的相关研究更少。美国NASA Ames中心60 MW电弧风洞(IHF)配套了隔热底面设计的半椭圆喷管,在喷管底边出口950 mm长度范围内铺设了高密度碳化硅陶瓷耐高温材料,以减少高焓试验时的气流焓降,提高模拟的真实程度,但至今尚未开展电弧风洞喷管壁温对平板试验的影响研究[5]。本文通过研制隔热半椭圆喷管,近似模拟飞行条件下的热壁条件和边界层剖面。在中国空气动力研究与发展中心(CARDC)FD–15B电弧风洞不同试验状态下,使喷管壁面达到不同的温度,对测试模型表面冷壁热流、绝热模型平衡温度以及喷管出口近边界层气流温度进行测试分析,以获得喷管不同壁温对平板试验模型的热影响。
1 试验方法
1.1 试验装置
试验在CARDC的FD–15B电弧风洞中进行。采用电弧风洞半椭圆喷管平板试验技术[11],被电弧加热的高温气流,经转接段和半椭圆喷管加速后进入试验段;在半椭圆喷管底边放置平板试验模型,使喷管直边的附面层光滑地过渡到平板模型。隔热半椭圆喉道半轴尺寸为30 mm × 41 mm,喷管出口半轴尺寸为167 mm × 90 mm,名义马赫数为4。如图2所示,喷管分为喉道段和扩张段。为了防止喉道段高热流部分被高温气流损坏,采用水冷壁结构;但在喷管扩张段,热流大幅度下降,若采用常规的冷却结构设计,将导致气流总焓下降过大。对于喷管扩张段,本文采用4层结构,包括防热板、隔热板、水冷板和结构基体。防热板采用耐高温材料铌钨合金制作,长度为334 mm,厚度为10 mm;隔热板采用氧化铝陶瓷纤维板制作,厚度为25 mm。
喷管底板经过防热、隔热处理,隔热板后的温度较低,传到水冷板的热量也较少。因此,水冷板的主要作用是防止基体结构受热变形。
1.2 测试方法
为了准确测试试验的来流总焓,低温(<2000 K)时一般采用总温探针法;温度高于总温探针测试上限时,用驻点热流和驻点压力探针测量驻点热流和压力,采用Fay-Riddell公式[12]间接计算来流总焓(简称“F–R探针法”)。F–R探针法的两种探针均安装在风洞旋转支架上轮流进行扫描测试。热流测试模型、温度测试模型安装在横向送进机构上,如图3所示。
流场核心区的总温探针(带滞止室)如图4所示,探针有2个相距10 mm的分支测点A、B,测点与喷管底边均距离20 mm,感温元件为Ir–IrRh40热电偶。该结构的总温探针由航空工业北京长城计量测试技术研究所热风洞标定[13]。本次试验的总温探针在1500~1900 K的平均标准不确定度为10.3。
为研究不同喷管壁面温度对喷管出口近壁面气流总温的影响,使用近壁面气流温度测试探针。设置#1~#4共4个探针,每个探针测点间隔5 mm,如图4所示。温度传感器为直径1.5 mm的K型接壳型热电偶。测试探针通过风洞横向送进支架送入流场。
为了研究不同喷管壁面温度对喷管出口平板试验模型表面冷壁热流的影响,在冷壁热流平板测试模型上布置小型水冷戈登计(探头直径为4.5 mm),如图5左图所示,可长时间测试壁面热流。F–R探针驻点及冷壁热流平板测试模型上的戈登计均由弧光灯标定系统标定[14],标定后重复性精度在1%以内。冷壁热流平板测试模型尺寸为100 mm × 100 mm,与半椭圆喷管底边平齐连接,迎角为0°。测点布置如图5右图所示,测点1~7为热流测点,
${\rm{P}}1、{\rm{P}}2 $ 为压力测点。在喷管防热板背面采用K型热电偶进行背温测试,分析防热板受热情况,并监测工作温度是否超过许用温度。隔热平板测温模型尺寸为100 mm × 100 mm,厚2 mm,表面带高发射率黑色涂层,以便使用IGAR 12–LO比色高温计测试表面温度,测点位置为图5的测点3,比色高温计的测试精度为读数的0.5% ± 1 ℃。如图6所示,为减少传热,模型四周留有隔热缝,背面采用气凝胶材料隔热,以钢丝拉紧固定。不锈钢板背面焊接K型热电偶,其测试精度为 ± 0.75%T(T为热电偶的实测温度值)。背面温度测点与冷壁热流平板测试模型的2、3、4、6、P2测点对应。试验时,冷壁热流平板测试模型与隔热平板测温模型依次被送至流场的相同位置,上表面与喷管底边平齐。改变电弧风洞试验状态,获得冷壁热流及隔热平板模型背面平衡温度进行对比分析。
为减少不同车次试验运行参数误差对试验结果的影响,在同一车次完成流场总温、驻点压力、冷壁热流、表面压力等参数测试。
1.3 试验状态
表1为电弧风洞试验状态参数,状态1、2的总焓由总温探针测得,状态3的总焓较高,使用F–R探针法获得。表中数据为3次试验的平均值。
表 1 电弧风洞试验状态参数Table 1 Average values of test parameters in arc-heated wind tunnel状态 总压/MPa 总焓/(MJ·kg−1) 驻点压力/kPa 驻点热流/(MW·m−2) 气体流量/(kg·s−1) 总功率/MW 1 2.57 1.04 360.76 1.32 6.79 14.53 2 1.98 2.12 277.40 2.86 3.43 15.84 3 1.66 2.55 236.70 3.75 2.45 15.43 电弧风洞同一车次总温探针测试的温度数据如图7所示,T11、T12为状态1的A、B两个分支测点,T21、T22为状态2的A、B两个分支测点。由图7可知,当探针达到稳定时,流场的总温波动很小,两支热电偶的测试值相差均小于1%。状态3的焓值由驻点热流和驻点压力探针间接测得,根据试验验证,同一次试验测试的波动也小于1%,因此设备运行参数的影响均较小。
2 结果及分析
2.1 试验结果分析
喷管防热板背面特征点随电弧风洞运行时间的温升曲线如图8所示,在试验状态1,100 s后上升斜率变小,200 s后温度趋于平衡,其他测点情况类似。在喷管壁面不同平衡温度条件下,对喷管出口近壁面气流总温、测试模型表面冷壁热流、隔热模型平板平衡温度的影响进行测试。喷管壁面温度对模型表面冷壁热流的影响如图9和表2所示,所有数值均为3次试验的平均值。图中r1–1为状态1第1个测点的冷壁热流,r2–1为状态2第1个测点的冷壁热流,依此类推。在不同时刻将装有戈登计的测试模型送入流场,戈登计单次测试时间为7 s,可获得稳态冷壁热流。在试验150 s后,冷壁热流上升不明显,基本保持平衡。状态3焓值较高,在100 s时温度已超过1000 K,但尚未达到平衡;由于防热板承受温度能力有限,因此仅测试了100 s时的冷壁热流。
表 2 不同总焓及防热板温度下的表面冷壁热流Table 2 Heat flux on cold wall with different total enthalpies and plate temperatures状态 时间/s 特征点背温T6/K 测点1热流/(kW·m−2) 测点2热流/(kW·m−2) 测点3热流/(kW·m−2) 10 309.0 182.0 197.0 180.0 1 300 443.0 204.0 224.0 207.0 增量/% 43.4 12.1 13.7 15.0 10 293.0 390.0 370.0 380.0 2 300 1096.0 420.0 420.0 430.0 增量/% 274.1 7.7 13.5 10.5 10 303.0 405.0 418.0 414.0 3 100 1046.0 424.0 456.0 447.0 增量/% 245.2 4.7 9.1 8.0 从图9中可以看到状态1、2、3下冷壁热流随喷管壁温增加的情况:冷壁热流增加了4.7%~15.0%,随着来流总焓增大,热流增量变小。喷管壁面温度上升到一定程度后,冷壁热流基本保持不变,这是因为壁面达到一定温度后,边界层气流与壁面温差减小,热损失减小,测试板热流不再大幅随壁面温度上升。
喷管壁面温度对近壁气流总温的影响如图10和表3所示,图中T1–0为状态1喷管冷壁温度,T1–1为状态1热壁温度,T2–0为状态2冷壁温度,T2–1为状态2热壁温度。风洞运行不同时间,喷管壁面温度不同,分别在风洞运行7 s和140 s时将近壁面总温排架送入流场,测试近壁气流总温的变化。风洞运行7 s时,喷管壁面温升小,近似为冷壁;140 s时,温升基本达到平衡。由表3可知,试验状态气流总焓越高,冷/热壁条件下气流总温温差相差越大。
隔热平板测温模型的平均温度对比如表4所示,其中,喷管防热板壁面近似热平衡时,隔热平板测温模型正面比色高温计温度为T1,背面5个热电偶测点的平均温度为Tb1;冷壁时(即扩张段采用水冷却壁面)隔热平板测温模型正面比色高温计温度为T2,背面5个测点的平均温度为Tb2。来流焓值越高,温差越大。根据本文的隔热结构设计,焓值在1.00~2.55 MJ/kg范围内,喷管冷却对防热试验模型表面平衡温度的影响在5%以内。
表 3 近壁面气流总温变化Table 3 Variation of total air temperature near wall状态 测点高度/mm 0 5 10 15 1 冷壁温度/K 641 757 873 978 140 s温度/K 714 789 929 985 增量/% 11.4 4.2 6.4 0.7 2 冷壁温度/K 901 1022 1272 1408 140 s温度/K 1058 1214 1359 1459 增量/% 17.4 18.8 6.8 3.6 表 4 隔热平板测温模型平均温度对比Table 4 Comparison of average temperature on stainless-steel plate状态 T1/K T2/K Tb1/K Tb2/K 表面温差率/% 背面温差率/% 1 735 723 648 643 0.67 0.78 2 1127 1084 1034 1004 3.89 2.94 3 1249 1201 1109 1063 3.92 4.23 2.2 数值模拟分析
由于试验状态焓值不高,采用FLUENT软件进行计算,选用
$ k - \omega $ 模型,防热板、隔热板与气流耦合传热。采用结构网格划分计算模型,网格总数为886129。边界层加密,并基于壁面参数的网格雷诺数准则控制第一层网格的高度,边界层网格的雷诺数小于100。考虑了网格数对计算结果的收敛和影响,依据残差值判定收敛性。该算例的残差值均控制在10−6以下且残差不再跳动时,判定为收敛。最终网格如图11所示。数值模拟了状态2下壁温对冷壁热流及附面层的影响尺度范围和百分比,分析了焓降大小。由于试验获得的温度测试结果是单点的,防热板上温度分布不均,因此进行了简化处理,对平均冷壁温度为300和500 K的情况进行分析。
图12为喷管防热板取等温壁边界条件时冷壁热流和近壁面气流总温的计算对比图。从图中可看出:测试模型的热流在喷管防热板热壁时增大7%~9%,气流总温在距离喷管壁面0~10 mm范围内最大差值为6%,且最大差值产生于近壁面5 mm以内。
由于喷管喉道及后部采用水冷铜制造,对高温气流有冷却效果,因此焓值有所下降。喷管扩张段数值计算结果(图13)表明:喷管防热板前端绝热温度受喷管喉道段冷却壁面影响而下降,与冷却部位连接处温度最低,绝热温度沿流向逐渐上升,最后达到绝热温度稳定值。
表5为对试验和数值模拟结果的分析。采用本文设计的防热板结构后,与理想绝热条件相比,隔热平板测温模型表面测试温度和近壁面气流测试温度的焓值下降均小于5%。
表 5 焓值降低量Table 5 Decrease of enthalpy状态 总焓/(MJ·kg−1) 计算恢复焓/(MJ·kg−1) 隔热模型恢复焓/(MJ·kg−1) 近壁气流恢复焓/(MJ·kg−1) 隔热模型焓降/% 近壁气流焓降/% 1 1.04 0.763 0.726 0.76 4.85 0.39 2 2.12 1.216 1.164 1.16 4.28 4.61 3 2.55 1.373 1.306 1.32 4.88 3.86 3 结 论
为了研究喷管壁面温度对平板试验的影响,设计了具有防/隔热壁面的喷管,研究了不同喷管壁面温度对喷管出口近壁面气流总温、模型表面冷壁热流、模型表面温度的影响。研究表明:使用本文研制的防/隔热喷管,气流总焓在1.00~2.55 MJ/kg范围内,喷管防热板受热温度上升后,喷管出口处测试模型的热流比冷壁面增加4.7%~15.0%,且随着来流总焓增大,该热流增量相对变小。
在上述喷管及试验条件下,近壁面气流总温随喷管壁面温度升高而升高0.7%~19.0%,总焓越高,冷、热壁条件下近壁面气流总温温差越大。喷管防热板热壁和冷壁时,喷管出口处的不锈钢隔热平板模型表面及背面平均温差在5%以内。与理想绝热条件相比,隔热平板测温模型表面测试温度和近壁面气流总温的焓值下降均小于5%。
以上结果表明:对于电弧风洞喷管直连平板的材料及热结构考核试验,需要考虑喷管壁面冷却对考核结果带来的影响,并相应提高来流总焓或采取防隔热措施,以减小或消除平板模型表面由于恢复焓降低导致试验状态热流等偏低对考核结果的影响。
-
-
[1] 杨富荣, 陈力, 闫博, 等. 干涉瑞利散射测速技术在跨超声速风洞的湍流度测试应用研究[J]. 实验流体力学, 2018, 32(3): 82–86. DOI: 10.11729/syltlx20170103 YANG F R, CHEN L, YAN B, et al. Measurement of turbulence velocity fluctuations in transonic wind tunnel using Interferometric Rayleigh Scattering diagnostic technique[J]. Journal of Experiments in Fluid Mechanics, 2018, 32(3): 82–86. doi: 10.11729/syltlx20170103
[2] 李晓辉, 王宏伟, 黄湛, 等. 层析粒子图像测速技术研究进展[J]. 实验流体力学, 2021, 35(1): 86–96. DOI: 10.11729/syltlx20190160 LI X H, WANG H W, HUANG Z, et al. Research advances of tomographic particle image velocimetry[J]. Journal of Experiments in Fluid Mechanics, 2021, 35(1): 86–96. doi: 10.11729/syltlx20190160
[3] 吴戈, 李韵, 万明罡, 等. 平面激光诱导荧光技术在超声速燃烧火焰结构可视化中的应用[J]. 实验流体力学, 2020, 34(3): 70–77. DOI: 10.11729/syltlx20190168 WU G, LI Y, WAN M G, et al. Visualization of flame structure in supersonic combustion by Planar Laser Induced Fluorescence technique[J]. Journal of Experiments in Fluid Mechanics, 2020, 34(3): 70–77. doi: 10.11729/syltlx20190168
[4] OJO A O, FOND B, VAN WACHEM B G M, et al. Thermographic laser Doppler velocimetry[J]. Optics Letters, 2015, 40(20): 4759–4762. doi: 10.1364/ol.40.004759
[5] GU C L, ZOU X, ZUO Z, et al. Doppler velocimeter based on dual-comb absorption spectroscopy[J]. Photonics Research, 2020, 8(12): 1895–1903. doi: 10.1364/prj.398876
[6] LI F B, ZHANG H B, BAI B F. A review of molecular tagging measurement technique[J]. Measurement, 2021, 171: 108790. doi: 10.1016/j.measurement.2020.108790
[7] WANG H P, YANG Z X, LI B L, et al. Predicting the near-wall velocity of wall turbulence using a neural network for particle image velocimetry[J]. Physics of Fluids, 2020, 32(11): 115105. doi: 10.1063/5.0023786
[8] 许德辰, 张悦, 刘欣乐, 等. 基于粒子追踪测速的壁面摩擦应力测量[J]. 实验流体力学, 2022, 36(2): 131–138. DOI: 10.11729/syltlx20210156 XU D C, ZHANG Y, LIU X L, et al. Measurement of wall-shear stress via micro-particle tracking velocimetry[J]. Journal of Experiments in Fluid Mechanics, 2022, 36(2): 131–138. doi: 10.11729/syltlx20210156
[9] TIEN W H, DABIRI D, HOVE J R. Color-coded three-dimensional micro particle tracking velocimetry and application to micro backward-facing step flows[J]. Experiments in Fluids, 2014, 55(3): 1684. doi: 10.1007/s00348-014-1684-x
[10] RAFFEL M, WILLERT C E, SCARANO F, et al. Particle image velocimetry: a practical guide[M]. Array Cham: Springer, 2018.
[11] RAGNI D, SCHRIJER F, VAN OUDHEUSDEN B W, et al. Particle tracer response across shocks measured by PIV[J]. Experiments in Fluids, 2011, 50(1): 53–64. doi: 10.1007/s00348-010-0892-2
[12] SAMIMY M, LELE S K. Motion of particles with inertia in a compressible free shear layer[J]. Physics of Fluids A: Fluid Dynamics, 1991, 3(8): 1915–1923. doi: 10.1063/1.857921
[13] MOLEZZI M J, DUTTON J C. Application of particle image velocimetry in high-speed separated flows[J]. AIAA Journal, 1993, 31(3): 438–446. doi: 10.2514/3.11349
[14] PAN F, SÁNCHEZ-GONZÁLEZ R, MCILVOY M H, et al. Simultaneous three-dimensional velocimetry and thermo-metry in gaseous flows using the stereoscopic vibrationally excited nitric oxide monitoring technique[J]. Optics Letters, 2016, 41(7): 1376–1379. doi: 10.1364/ol.41.001376
[15] YE J F, SHI D Y, SONG W Y, et al. Investigation of turbulence flow characteristics in a dual-mode scramjet combustor using hydroxyl tagging velocimetry[J]. Acta Astronautica, 2019, 157: 276–281. doi: 10.1016/j.actaastro.2018.12.040
[16] BATHEL B, JOHANSEN C, DANEHY P, et al. Review of fluorescence-based velocimetry techniques to study high-speed compressible flows[C]//Proceedings of the 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition. 2013. doi: 10.2514/6.2013-339https://doi.org/
[17] CHEN F, LI H X, HU H. Molecular tagging techniques and their applications to the study of complex thermal flow phenomena[J]. Acta Mechanica Sinica, 2015, 31(4): 425–445. doi: 10.1007/s10409-015-0464-z
[18] MIRZAEI M, DAM N J, VAN DE WATER W. Molecular tagging velocimetry in turbulence using biacetyl[J]. Physical Review E, 2012, 86(4): 046318. doi: 10.1103/physreve.86.046318
[19] 杨文斌, 陈力, 闫博, 等. 基于飞秒激光电子激发标记测速技术的剪切流场速度测量[J]. 实验流体力学, 2022, 36(4): 94–102. DOI: 10.11729/syltlx20210060 YANG W B, CHEN L, YAN B, et al. Transient velocity measurement of shear flow using Femtosecond Laser Electronic Excitation Tagging[J]. Journal of Experiments in Fluid Mechanics, 2022, 36(4): 94–102. doi: 10.11729/syltlx20210060
[20] MILES R, COHEN C, CONNORS J, et al. Velocity measurements by vibrational tagging and fluorescent probing of oxygen[J]. Optics Letters, 1987, 12(11): 861–863. doi: 10.1364/OL.12.000861
[21] BALL C G, FELLOUAH H, POLLARD A. The flow field in turbulent round free jets[J]. Progress in Aerospace Sciences, 2012, 50: 1–26. doi: 10.1016/j.paerosci.2011.10.002
[22] KAUSHIK M, KUMAR R, HUMRUTHA G. Review of computational fluid dynamics studies on jets[J]. American Journal of Fluid Dynamics, 2015, 5(3A): 1–11. doi: 10.5923/s.ajfd.201501.01
[23] ABRAMOVICH G N, GIRSHOVICH T A, KRASHE-NINNIKOV S I, et al. The theory of turbulent jets[M]. Moscow: MIT Press, 1984.
[24] RAJARATNAM N. Turbulent jets[M]. New York: Elsevier, 1976.
[25] SONG Z H, WANG W T, ZHU N, et al. Gas velocimetry based on infrared laser-induced fluorescence[J]. Physics of Fluids, 2021, 33(12): 125126. doi: 10.1063/5.0074367