Two-dimensional distribution measurement of direct-connect scramjet combustion flow field based on TDLAS multi-absorption lines
-
摘要:
本文针对超燃冲压发动机燃烧室扩张段非均匀流场温度和水汽浓度二维分布的高分辨率测量需求,发展了先进的可调谐激光吸收光谱(TDLAS)燃烧场分布重建技术,该技术通过增加激光测量光路上扫描获得的水汽吸收谱线数目,实现场分布重建问题求解方程数量的增加;通过联立所有交叉光路下吸收光谱获得的吸光度方程,构建以温度和浓度为未知数的最优化目标函数;并采用全局寻优模拟退火算法对目标函数进行求解,实现温度场和水汽分压场的重建。发动机直连台架试验中,采用正交光路布局,设计共16条测量光路(水平5条、垂直11条)的方形光机结构,集成TDLAS测量系统。对5只DFB激光器采用分时直接吸收探测方式,测量频率4 kHz,每条测量光路下可扫描获得5条水汽吸收谱线(
7467.77 、7444.36 、7185.60 、7179.75 和6807.83 cm),系统在高温炉上开展了多温度台阶标定测试,温度测量偏差在2.7%以内。外场试验中,对16条光路下同步采集到的吸收光谱数据进行离线处理,获得了发动机燃油点火、燃烧、熄火各个状态下的温度场和水汽分压场分布数据。试验结果表明:TDLAS多线吸收测量技术能够实现场分布准确稳定测量,满足发动机复杂燃烧流场诊断和恶劣工况工程应用需求。Abstract:Aiming at the demand of two-dimensional distribution high-resolution measurement of temperature and water vapor concentration in non-uniform scramjet combustion chamber expansion section, advanced Tunable Diode Laser Absorption Spectroscopy (TDLAS) recon-struction method has been developed. By increasing the number of water vapor absorption lines obtained by scanning the laser wavelength, the number of equations for solving the reconstruction problem correspondingly increased, combining the absorbance equations of all absorption spectra under all laser paths, constructing the optimization objective function with temperature and concentration as unknowns, and using the global optimization simulated annealing algorithm to reconstruct the temperature and water vapor concentration distribution. In the direct-connect scramjet combustion test, the orthogonal optical path layout is adopted, and the square optical mechanical structure with 16 measuring optical paths of 5 horizontal and 11 vertical channels is designed. TDLAS measurement system is assembled, and the time division multiplexed direct absorption detection method is adopted for 5 DFB lasers, with the measurement frequency of 4 kHz. Five water vapor absorption spectral lines (
7467.77 、7444.36 、7185.60 、7179.75 and6807.83 cm) can be obtained at each measured optical path, the system has carried out thermometric validation by using high-temperature furnace on the laboratory, and the temperature measurement deviation is within 2.7%. In the test, the absorption spectrum data synchronously collected under 16 optical paths are processed offline, and the distribution data of temperature field and water vapor partial pressure under various states of ignition, combustion and flameout are obtained. The test results show that TDLAS multi-absorption measurement technology can realize accurate and stable reconstruction, and meet the engineering application requirements of complex combustion flow field diagnosis and bad working conditions. -
0 引 言
飞机结冰是指在结冰气象条件下过冷水滴或过冷雨在飞机机体表面某些部位凝结形成冰层的现象(也可能是水蒸气产生的凝华结冰)[1]。为确保结冰条件下的飞机适航性,必须防止或除去影响飞机正常飞行的积冰。积冰形状通常并不规则,识别积冰界线是定性定量识别积冰几何形状、对积冰进行三维重建的必要条件,也是实施防除冰的前提条件。
目前,探测飞机积冰的主要技术有压电技术、超声导波技术[2]、激光技术、红外热波探测技术及基于神经网络的冰形预测技术[3]等。压电技术是通过压电器件组成的结冰传感器来探测冰层的形状[4]。白天等[5]精准地探测出不同结冰状态对压电双晶片悬臂梁结构特征的影响。超声导波技术则是通过传感器检测导波信号的变换程度获取积冰图像,其准确度取决于传感器的种类和布局。与前两者相比,以激光技术进行积冰探测速度较快、精确度较高。国内学者基于激光技术开展了积冰探测的深入研究,康含玉等[6]通过改进中心线提取算法实现了冰形表面激光光带中心线的快速提取。随着人工智能的快速发展,研究者开始利用神经网络对积冰形状参数进行预测,柴聪聪等[7]利用BP神经网络建立了冰形特征参数预测模型,通过输入影响积冰形状的参数得到了冰形特征参数,该方法预测的冰角角度、积冰范围误差较低。
与上述积冰探测方式相比,红外热波探测不受云、雾、雨等气象条件限制,且可在不接触积冰样件的前提下便捷地进行探测。国外一些研究者已尝试将红外热波探测用于风力涡轮叶片结冰探测[8]。在红外热波探测过程中,探测结果易受探测样件和探测环境条件的影响,导致光照不均匀、红外图像和背景融合等问题,对识别目标体造成一定干扰。因此,在利用红外热波探测技术对积冰进行探测的过程中,需要选择合适算法识别积冰界线(包括外部界线和内部界线),从而更为准确地识别冰形。
国内外研究者主要通过传统边缘检测算子、融合改进后的边缘检测算子以及阈值分割和形态学融合的方法解决积冰探测图像处理问题。Wang等[9]采用基于小波变换、浮动阈值以及数字图像形态学和最佳阈值相融合的算法对积冰界线进行探测;陆佳政等[10]采用自适应分割阈值法识别积冰形状。这两种方法可以清晰识别积冰的外部界线,但由于界线区域灰度值差异不明显,无法准确识别积冰的内部界线。
为实现积冰界线(外部界线和内部界线)探测与识别,本文搭建了闪光红外热波探测实验台,将传统算法与高斯–拉普拉斯金字塔和面积滤波相结合的算法进行比较,确定识别积冰界线的优选方法。
1 红外热波积冰探测原理及实验
与被动式红外探测相比,主动式红外无损探测能够更好地呈现被探测样件的细节特征,常被应用于航空领域[11]。本文利用高能闪光氙灯作为主动激励装置,通过红外热成像系统和图像处理控制系统对积冰样件(以制冷设备制取)进行定量识别。
积冰受热后的传热控制方程如下:
$$\begin{split}& \frac{\partial T}{\partial x}\left({\lambda }_{x}\frac{\partial T}{\partial x}\right)+\frac{\partial T}{\partial y}\left({\lambda }_{y}\frac{\partial T}{\partial y}\right)+\\& \frac{\partial T}{\partial z}\left({\lambda }_{z}\frac{\partial T}{\partial z}\right)+{\varphi }_{{\rm{ s}}}=\rho c\frac{\partial T}{\partial x}\end{split} $$ (1) 式中:T为温度;λx、λy和λz分别为材料在x、y和z方向的热传导系数;φs为均匀内热源的生热率;ρ和c分别为材料的密度和比热。
以高能闪光氙灯向积冰样件注入能量,由于积冰样件形状不规则,在热阻效应作用下,当热波传递至内部界线处时,积冰样件的阶跃型内部界线以内部温度变化差异的形式表现出来[12]。采用红外热成像仪连续观测和记录积冰样件内部的温度场变化,生成红外热序图。对红外图像进行分析处理,即可实现积冰界线的快速探测。根据上述原理,将采集到的红外热序图的灰度图进行图像处理,识别出积冰样件的外部界线和内部界线。
在实验中,分别对规则圆形积冰样件、规则方形积冰样件及对应的非规则阶跃型积冰样件进行积冰探测。实验设备包括用于积冰探测的闪光红外热波探测系统和用于积冰样件制作的冷冻冰柜。将纯净水注入模具,在冰柜中冷冻为积冰样件,再将厚度不同的积冰样件堆叠冷冻,以模拟积冰样件内部阶跃状态,如图1所示。
2 积冰界线探测识别的传统算法
2.1 传统边缘检测理论
边缘检测的目的是找到图像中像素值发生急剧变化的像素点集合,通常以图像边缘的形式表示。同时,边缘检测可以有效抑制图像的非主要信息,保留重要的图像特征属性结构。不同边缘检测算子的思路不同,对于不同属性结构的图像,目前还无法使用同一边缘检测算子获取同样的边缘[13]。较为传统的边缘检测是采用一阶或二阶微分算子对图像进行处理。Roberts算子、Sobel算子和Prewitt算子是较为常见的一阶微分算子,Laplacian算子和Canny算子是较为常见的二阶微分算子[14]。Roberts算子是一种通过局部差分计算边缘的算子,其卷积模板为一个2 × 2的矩阵(具有水平和垂直方向):
$$ {\boldsymbol{G}}_{x} = \left[ {\begin{array}{*{20}{c}} {\text{1}}&{\text{0}} \\ {\text{0}}&{{{ - 1}}} \end{array}} \right] ,\; {\boldsymbol{G}}_{y} = \left[ {\begin{array}{*{20}{c}} {\text{0}}&{{{ - 1}}} \\ {\text{1}}&{\text{0}} \end{array}} \right] $$ (2) 计算2 × 2卷积模板对角线像素之差的近似梯度幅值,计算方法如式(3)~(5)所示:
$$ \Delta f_{x}(x, y)= f_{x}(x, y)-f_{x}(x+1, y+1) $$ (3) $$ \Delta f_{y}(x, y)= f(x, y+1)-f(x+1, y) $$ (4) $$ \boldsymbol{G}(x, y)=\sqrt{[\Delta {f}_{x}(x, y)]^2+[\Delta {f}_{y}(x, y)]^2 }$$ (5) Roberts算子对噪声十分敏感,无法有效抑制噪声。Prewitt算子是对Roberts算子的优化,其卷积模板为:
$${\boldsymbol{G}}_x= \left[ {\begin{array}{*{20}{c}} {{{ - 1}}}&{{{ - 1}}}&{{{ - 1}}} \\ {{0}}&{{0}}&{{0}} \\ {{1}}&{{1}}&{{1}} \end{array}} \right],\; {\boldsymbol{G}}_y= \left[ {\begin{array}{*{20}{c}} {{{ - 1}}}&{{0}}&{{1}} \\ {{{ - 1}}}&{{0}}&{{1}} \\ {{{ - 1}}}&{{0}}&{{1}} \end{array}} \right] $$ (6) 采用Prewitt算子对图像相邻像素点的灰度进行水平、垂直方向一阶微分,即可得到边缘。实际上,Prewitt算子是对水平和垂直方向进行非归一化均值差分,因而可以有效地对图像噪声进行平滑,但识别边缘不够准确。
Sobel算子是对Prewitt算子的改进,其卷积模板是一个3 × 3的矩阵,对图像邻接像素点的灰度进行上下左右四领域的灰度值加权差分:
$$ {\boldsymbol{G}}_x=\left[ {\begin{array}{*{20}{c}} {{{ - 1}}}&{{{ - 2}}}&{{{ - 1}}} \\ {{0}}&{{0}}&{{0}} \\ {{1}}&{{2}}&{{1}} \end{array}} \right] ,\; {\boldsymbol{G}}_y= \left[ {\begin{array}{*{20}{c}} {{{ - 1}}}&{{0}}&{{1}} \\ {{{ - 2}}}&{{0}}&{{2}} \\ {{{ - 1}}}&{{0}}&{{1}} \end{array}} \right] $$ (7) 对图像邻接像素点的灰度进行水平和垂直方向加权平均差分;根据卷积模板得到图像的灰度差近似值,当灰度差近似值大于某一阈值时,则判断其为边缘。因此,Sobel算子对图像噪声具有平滑作用,同时还可以准确判断边缘的方向。
与上述一阶微分算子不同,Laplacian算子和Canny算子是二阶微分算子。Laplacian算子是n维欧几里得空间中的一个二阶微分算子,定义为梯度的散度,可以更加精准地识别边缘位置,同时对图片进行锐化。作为较常见的边缘检测算子,Canny算子既能精准地提取边缘,又能有效地抑制噪声。以Canny算子进行边缘检测的主要步骤为:对图像进行高斯滤波降噪;利用一阶偏导计算对角线上4个方向的梯度幅值和方向;对梯度进行非极大值抑制;进行双阈值检测和边缘连接。
除此之外,基于相位一致性对红外图像进行边缘检测[15]比Canny算子检测得到的边缘更加精准。微分算子是针对灰度值进行边缘检测的,而相位一致性是基于频域幅值的比较进行边缘检测,边缘相似的特征在同一阶段出现的频率较高,将这些特征利用K–means进行分类进而得到边缘。
2.2 传统边缘检测对积冰界线的识别
由于红外热成像仪容易受到实验环境和样件自身温度的影响,红外热序图存在复杂噪声较多、图像模糊等问题[16]。因此,在进行积冰界线检测时,需要先对红外热序图进行主成分分析,在保留图像大部分原有特征的前提下降低图像噪声。对规则圆形、规则方形及其对应的非规则阶跃型积冰样件红外热序图的主成分分析结果如图2所示。在此基础上,分别对主成分分析图进行传统边缘检测,检测结果如图3~6所示。
3 基于高斯–拉普拉斯金字塔和面积滤波的边缘检测改进算法
积冰样件红外图像的内部界线不明显,受检测环境的影响相对较大,无法通过传统的边缘检测取得理想的内部界线。另外,传统的红外图像阈值分割及边缘检测技术仅能解决积冰样件红外图像处理中某一类的问题[17]。
针对上述问题,本文对传统边缘检测算子进行改进,以便更准确地识别积冰界线(外部界线和内部界线)。为达到预期效果,将高斯–拉普拉斯金字塔算法和面积滤波算法相结合并加以改进,以解决因阶跃区域红外图像灰度值差异不明显所导致的无法准确检测内部界线的问题[18]。
首先,利用高斯–拉普拉斯金字塔算子突出图像的高频结构,加强图像特征信息;然后,通过面积滤波过滤掉较小的噪声,保存图像边缘信息。算法流程如图7所示。
高斯–拉普拉斯金字塔算子是由高斯滤波函数和拉普拉斯算子组成的[19]。高斯金字塔的作用是对图像向下逐级进行取样分解,并在取样前对原图像进行高斯滤波。式(8)和(9)分别为高斯滤波器公式和高斯金字塔公式:
$$ g(x,y)=\frac{1}{2\mathrm{\pi }{\delta }^{2}}{\mathrm{exp}}\left({-\frac{{x}^{2}+{y}^{2}}{2\mathrm{\pi }{\delta }^{2}}}\right) $$ (8) $$ F(n+1)=F(n)\times g(x,y)\times S(\mathrm{s}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{e}1,\mathrm{s}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{e}2) $$ (9) 式中:F(n)为第n次分解的图像;S(scale1, scale2)为图像在行和列上的取样率;两个方向x和y服从方差为δ2的正态分布。本文中行和列的取样率均为2,因此,每一次取样后的图片尺寸缩小为原图像的1/4,分辨率下降为原图像的1/2[20]。
拉普拉斯金字塔的作用是保存高斯金字塔向下取样分解过程中相邻两级之间的残存差,即图像的高频细节。拉普拉斯金字塔公式如下:
$$ \begin{split}&L(n)=F(n)-F(n+1)\times S(\mathrm{s}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{e}1,\mathrm{s}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{e}2)\times g(x,y)\\\end{split} $$ (10) 本文先采用高斯金字塔将原图像向下取样分解两次,且两次的取样率皆为2;再将向下取样的图像通过拉普拉斯金字塔变换得到提取特征后的图像;最后将高斯–拉普拉斯变换后的图像二值化。以方形阶跃型积冰样件为例,通过高斯–拉普拉斯滤波后得到的图像如图8所示。
经过高斯–拉普拉斯金字塔处理之后,图像的基本特征已经被提取出来,但是红外热序图存在的噪声依然较多。因此,在通过高斯–拉普拉斯金字塔保存图像高频特征的同时,图像噪声也被突出了。
在图8的二值图像中,图像噪声和所需提取的图像边缘可理解为面积值不同的连通域。因此,首先对高斯–拉普拉斯滤波处理后的图像进行连通域面积计算,将相对面积较小的连通域反色(即删除较小连通域);通过图像形态学变换,将已过滤掉较小噪声的图像在原图像基础上进行融合,加深图像外部界线和阶跃型内部界线的结构信息;最后,对加深图像特征的灰度图进行Canny算子边缘检测,得到具有阶跃型界线的积冰样件检测结果,如图9所示。
为进一步检验算法的有效性,选择圆形阶跃型积冰样件的主成分分析图进行测试,得到的边缘检测结果如图10所示。实验表明,改进后的算法可以识别出圆形阶跃型积冰样件的内部界线和方形阶跃型积冰样件的内部界线,适用于识别不规则积冰表面界线。
4 分析与讨论
采用上述传统边缘检测算法和基于高斯–拉普拉斯金字塔与面积滤波的新算法,对积冰样件的红外热序图进行界线探测,对探测结果进行比较分析,可以发现:
传统的边缘检测算法主要是针对图像灰度值的突变来获取图像边缘,检测步骤是先对图像进行降噪和平滑滤波,再计算边缘区域。因此,在使用闪光红外热波探测系统对规则型积冰样件及其对应的阶跃型积冰样件进行积冰内部界线识别时,由于阶跃型积冰样件红外图像内部界线的灰度值变化并不明显,传统的边缘检测算法无法探测出积冰样件阶跃型内部界线。但是,规则型积冰样件外部界线的灰度值突变较为明显,采用传统的边缘检测算法可以检测出外部界线。
相位一致性算法可以识别出内部界线,但是由于图像的相位一致性差值并不明显,所以无法准确清晰地绘制出相位一致性较高的边缘区域。
本文采用的高斯–拉普拉斯金字塔和面积滤波相结合的算法,其主体思想是:首先提取积冰探测红外热序图的主要图像成分及结构信息;而后,考虑到噪声在灰度图中的连通域面积比图像内部界线及外部界线连通域面积小,可以通过面积滤波算法过滤掉目标探测不需要的噪声。因此,采用高斯–拉普拉斯和面积滤波相结合的算法可以更为完整地保留图像特征结构。
高斯–拉普拉斯金字塔和面积滤波相结合的边缘检测算法优先考虑图像特征,而后针对图像整体特征进行降噪,不会将探测结果视为噪声而导致无法识别内部界线的问题;但检测结果的噪声高于传统算法。因此,该算法在识别不规则积冰内部界线方面具有一定优势。
5 结 论
积冰界线识别是积冰几何形状定性与定量识别的基础。通过本文研究发现,传统的边缘检测算法仅能识别出部分积冰样件的外部界线,不能准确探测出积冰样件的阶跃型内部界线。与传统边缘检测算法相比,高斯–拉普拉斯金字塔与面积滤波相结合的边缘检测算法更适用于识别阶跃型内部界线。采用该算法获取积冰样件更为精确的内部界线、进一步降低噪声对结果的影响是后续的研究方向。
-
表 1 筛选的水汽吸收谱线光谱参数
Table 1 Selected water vapor absorption spectral parameters
$ \lambda $/
nm$ {v}_{0} $/
cm−1${E}{{'}{'} }$/
cm−1${S}_{vi,\;{\rm{H}}}$/
(cm−1atm−1)${S}_{vi,\;{\rm{M}}}$/
(cm−1atm−1)1339 7467.77 2551.48 1.27 × 10−5 1.093 × 10−5 1343 7444.36 1790.71 1.16 × 10−3 1.100 × 10−3 1392 7185.60 1045.06 1.97 × 10−2 1.905 × 10−2 1393 7179.75 1216.19 5.97 × 10−3 5.814 × 10−3 1469 6807.83 3319.45 6.17 × 10−7 6.032 × 10−7 -
[1] ČARNOGURSKÁ M, PŘÍHODA M, KOŠKO M, et al. Verification of pollutant creation model at dendromass combustion[J]. Journal of Mechanical Science and Technology, 2012, 26(12): 4161–4169. doi: 10.1007/s12206-012-0877-6
[2] JACKSON K, GRUBER M, BUCCELLATO S. HIFiRE flight 2 project overview and status update 2011[C]//Proc of the 17th AIAA International Space Planes and Hypersonic Systems and Technologies Conference. 2011.doi: 10.2514/6.2011-2202
[3] SUR R, SUN K, JEFFRIES J B, et al. TDLAS-based sensors for in situ measurement of syngas composition in a pressurized, oxygen-blown, entrained flow coal gasifier[J]. Applied Physics B, 2014, 116(1): 33–42. doi: 10.1007/s00340-013-5644-6
[4] BRYNER E, SHARMA M, DISKIN G, et al. Tunable diode laser absorption technique development for determination of spatially resolved water concentration and temperature[C]//Proc of the 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. 2010: 299. doi: 10.2514/6.2010-299
[5] 赵然, 蔡文哲. 航空发动机燃烧室出口温度分布调试与测试技术研究[J]. 战术导弹技术, 2018(5): 84–89. ZHAO R, CAI W Z. Debugging and testing technology for temperature distribution of aero engine outlet[J]. Tactical Missile Technology, 2018(5): 84–89.
[6] BOLSHOV M A, KURITSYN Y A, ROMANOVSKII Y V. Tunable diode laser spectroscopy as a technique for combustion diagnostics[J]. Spectrochimica Acta Part B: Atomic Spectroscopy, 2015, 106: 45–66. doi: 10.1016/j.sab.2015.01.010
[7] CAI W, KAMINSKI C F. Tomographic absorption spectroscopy for the study of gas dynamics and reactive flows[J]. Progress in Energy and Combustion Science, 2017, 59: 1–31. doi: 10.1016/j.pecs.2016.11.002
[8] LIU C, XU L J, CAO Z. Measurement of nonuniform temperature and concentration distributions by combining line-of-sight tunable diode laser absorption spectroscopy with regularization methods[J]. Applied Optics, 2013, 52(20): 4827–4842. doi: 10.1364/AO.52.004827
[9] ZHANG G L, LIU J G, XU Z Y, et al. Characterization of temperature non-uniformity over a premixed CH4-air flame based on line-of-sight TDLAS[J]. Applied Physics B, 2016, 122(1): 3. doi: 10.1007/s00340-015-6289-4
[10] BUSA K, BRYNER E, MCDANIEL J, et al. Demonstration of capability of water flux measurement in a scramjet combustor using tunable diode laser absorption tomography and stereoscopic PIV[C]//Proc of the 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition. 2011: 1294. doi: 10.2514/6.2011-1294
[11] ZHANG L F, WANG F, ZHANG H D, et al. Simultaneous measurement of gas distribution in a premixed flame using adaptive algebraic reconstruction technique based on the absorption spectrum[J]. Chinese Optics Letters, 2016, 14(11): 111201–111205. doi: 10.3788/col201614.111201
[12] 宋俊玲, 洪延姬, 王广宇, 等. 光线分布对基于TDLAS温度场二维重建的影响[J]. 红外与激光工程, 2014, 43(8): 2460–2465. SONG J L, HONG Y J, WANG G Y, et al. Influence of beam distribution on two-dimensional temperature field reconstruction using TDLAS[J]. Infrared and Laser Engineering, 2014, 43(8): 2460–2465.
[13] XIA H, KAN R, XU Z, et al. Two-step tomographic reconstructions of temperature and species concentration in a flame based on laser absorption measurements with a rotation platform[J]. Optics and Lasers in Engineering, 2017, 90: 10–18. doi: 10.1016/j.optlaseng.2016.09.005
[14] CAI W W, KAMINSKI C F. A numerical investigation of high-resolution multispectral absorption tomography for flow thermometry[J]. Applied Physics B, 2015, 119(1): 29–35. doi: 10.1007/s00340-015-6012-5
[15] MA L, LI X S, SANDERS S T, et al. 50-kHz-rate 2D imaging of temperature and H2O concentration at the exhaust plane of a J85 engine using hyperspectral tomography[J]. Optics Express, 2013, 21(1): 1152–1162. doi: 10.1364/OE.21.001152
[16] 辛明原, 金星, 宋俊玲. 基于多波长激光吸收光谱技术的气体温度二维重建[J]. 航天器环境工程, 2016, 33(3): 269–274. DOI: 10.12126/see.2016.03.007 XIN M Y, JIN X, SONG J L. Two-dimensional reconstruction of gas temperature based on multi-wavelength laser absorption spectroscopy[J]. Spacecraft Environment Engineering, 2016, 33(3): 269–274. doi: 10.12126/see.2016.03.007
[17] 林鑫, 李飞, 王宽亮, 等. 旋流火焰燃烧温度场二维吸收光谱诊断[J]. 气体物理, 2019, 4(2): 55–61. LIN X, LI F, WANG K L, et al. Diode laser absorption tomography for swirling flame diagnostics[J]. Physics of Gases, 2019, 4(2): 55–61.
[18] 臧益鹏, 许振宇, 黄安, 等. 基于改进模拟退火算法的非均匀燃烧场分布重建[J]. 物理学报, 2021, 70(13): 111–122. ZANG Y P, XU Z Y, HUANG A, et al. Distribution reconstruction of non-uniform combustion field based on improved simulated annealing algorithm[J]. Acta Physica Sinica, 2021, 70(13): 111–122.
-
期刊类型引用(5)
1. 王海峰,展京霞,陈科,陈翔,陈梓钧. 战斗机大迎角气动特性研究技术的发展与应用. 空气动力学学报. 2022(01): 1-25 . 百度学术
2. 付京来,陈旭立. 6级海况下的某型发射平台回转俯仰耦合动力学研究. 兵工自动化. 2022(07): 68-71 . 百度学术
3. 孙亚东,李涌. 基于ADAMS的某型发射平台回转俯仰耦合动力学分析. 现代制造技术与装备. 2020(07): 96-98+101 . 百度学术
4. 付军泉,史志伟,陈杰,周梦贝,吴大卫,潘立军. 一种翼身融合布局飞行器的偏离特性分析. 实验流体力学. 2020(06): 32-37 . 本站查看
5. 肖中云,刘刚,江雄,王建涛. 飞行器双自由度耦合振荡的数值模拟方法. 航空学报. 2017(11): 123-132 . 百度学术
其他类型引用(2)