Research on microphone phase array design based on surrogate model
-
摘要: 气动声学风洞试验过程中,针对目标声源特性进行麦克风相位阵列改进时,为保证风洞试验效率,必须在较短的时间内完成阵列优化设计工作。为满足这一试验需求,最大限度提升麦克风相位阵列设计效率,引入了基于Kriging代理模型的优化设计方法。采用基于点扩散函数的计算程序进行阵列性能分析,获取阵列最大旁瓣水平和分辨率。通过对样本点响应值进行计算,建立Kriging代理模型,进而以计算速度极高的Kriging代理模型作为阵列性能分析方法开展优化搜索,避免了大量调用阵列性能计算程序导致计算耗时过高的问题,显著提升了麦克风相位阵列设计效率。该阵列设计方法能够有效改善麦克风阵列的测量性能,满足声学风洞试验的特殊应用需求。Abstract: In aero-acoustic wind tunnel tests, the improvement of the microphone phase array according to the noise source features must be conducted as fast as possible to ensure the required efficiency of the tests. An optimization method based on the Kriging surrogate model is proposed in this study. The maximum sidelobe level and the resolution of the microphone phase array are obtained by analyzing the results computed from the point spread function. A Kringing surrogate model is constructed based on the response values computed on sample points. Then the surrogate model was utilized to fast assess the maximum sidelobe level and the resolution of the microphone phase array, which avoided the high computational cost caused by computing the point spread function. This method can improve the performance of the microphone phase array effectively for aero-acoustic wind tunnel tests.
-
Keywords:
- microphone phase array /
- optimization /
- surrogate model /
- resolution /
- maximum sidelobe level
-
0 引言
研究声场与火焰的作用有着重要的工程价值和理论意义。一方面,声与火焰的耦合,将有效地组织燃烧过程,不仅强化了燃烧和传热,而且有效降低污染物的排放量,脉动燃烧技术是其典型代表;另一方面,燃烧过程中放热波动与压力脉动满足一定相位关系(瑞利准则),将诱发热声振荡,也称为燃烧不稳定,是航空发动机、燃气轮机等重要动力装置设计与运行中所遇到的技术难题,目前为止,一直没有得到有效的解决。此外,基于流动与燃烧的关系,火焰在声场速度脉动作用下的响应行为,如放热、火焰动力学等,对于认识湍流火焰燃烧特性具有重要意义。
学术界当前主要关注声场作用下贫燃预混火焰的响应行为,代表性的思想是构建描述速度脉动与放热量变化之间关系的火焰传递函数。Durox等[1]测量了不同几何形状的火焰在声场下的火焰传递函数,根据传递函数增益和相位差2方面的数据,比较了形状不同带来的非线性响应差别,分析了几何形状对火焰动力学的影响。Balachandran等[2]实验研究了速度脉动作用下湍流预混火焰的非线性响应行为,重点关注贫燃预混火焰在大振幅声场作用下的饱和机理。Lieuwen [3]以本生型锥焰为建模对象,通过求解恒定火焰传播速度的锋面动力学方程(G-方程),研究了火焰结构特征参量对声场作用下火焰动力学特性的影响,分析了火焰传递函数从线性响应发展为非线性响应的条件。
这些工作存在1个共同的特点:研究对象是化学当量比为1左右的预混火焰,其原因可能有2点,一是工程中燃气轮机燃烧室热声振荡的条件是贫燃预混燃烧,二是理论上描述火焰锋面动力学特性的G-方程是基于一维火焰的假设。值得注意的是,作为燃烧方式中同样重要的部分预混火焰(Partially premixed flame),其在声场下的响应行为则很少被人关注。最近,Kim等[4]以燃气轮机燃烧室内的部分预混旋流燃烧器火焰为对象,研究了火焰在受到声场机理下的响应行为,作者认为声场作用下部分预混火焰的响应行为目前还没有被很好地认识,原因是其间的作用关系更为复杂,涉及化学、火焰动力学等方面的难点更为突出。
本文以Rijke型驻波声场为声环境,以本生型层流部分预混火焰为关注对象,借助多种测试手段,通过分析声波振荡周期内火焰的时间序列红外图像和纹影图像,揭示火焰锋面的运动过程及由此带来的热量演化行为。
1 实验系统与方法
1.1 声场发生系统
激励火焰的声环境由Rijke管实现,管内要实现激励的声场需满足以下条件:(1) Rijke管两端开口,保证管内存在一定的平均流速;(2)距管下端1/4管长处存在热源,且热源需形成一定的温度梯度;(3)保证Rijke管为良好的声保持系统,即管侧壁需完全密封,不发生声泄漏;(4) Rijke管长径比大于10。
综合考虑上述条件,设计Rijke管为两端开口的直管,管径为150mm,管长为1.6m。如图 1所示,在管内距离下端口400mm (管长的1/4)处设置一通透性较好的近似平面热源,在热源加热过程中,热对流产生的管内气体流动与热源放热过程互相作用,使声压与放热满足一定的相位关系,最终激发管内的热声振荡。这种声环境是典型的一维驻波声场[5],图 1右侧给出了压力脉动量和速度脉动量的管内分布特征,管内中点位置压力振幅最大,速度振幅最小,管口两端则与此相反,在距管下端1/4管长位置,该点的压力振幅和速度振幅的乘积最大,最易于激发热声振荡。实验中的Rijke管采用内径150mm、壁厚4mm的碳钢管,热源由螺旋状电炉丝盘绕而成。为确保管内的流动通畅,热源支架由5根宽度为10mm的硅条按72°等角度均布排列。火焰受激励位置选择压力振幅和速度振幅相对较大的上1/4位置。声压特征通过在Rijke管侧壁上设置压力传感器提取,位置与火焰基部(燃烧器出口)对应。
理论上n阶的声压和速度波动表示为:
式中:An为压力振幅;kn=ωn/c为波数;ωn=nπc/L为n阶声场的频率;ρ0为气体密度;c是声速;L是管长。由上述公式可知,管内声场频率由本地声速和管长决定,振幅由热源的输出功率决定。实际测得的声场参数为:火焰位置的声压级(SPL)为144dB,管内声频率为111Hz,对应周期为9ms。
1.2 燃烧系统
燃烧系统的实验装置如图 2所示,甲烷和空气分别由各自供气系统送入混合腔4,流量用热式流量计测量,预混后的气体经混气管2进入烧嘴3,烧嘴由内径13mm、厚度2.5mm的不锈钢管加工而成,在没有声场作用的情况下火焰为本生型。Rijke管上方管段设置了一对平行放置的石英视窗,尺寸为90mm×400mm,供纹影法和红外摄像观测火焰。火焰的特征参量为甲烷流量1000mL/min,一次空气当量比为1.57,燃烧器内气流Re数为790,这种火焰在没有声场作用的条件下表现为双火焰(Double flame)结构。
1.3 纹影和红外拍摄系统
纹影系统的布置如图 3所示,光源处发出的光线,经过多组光学元件后,以平行光穿过Rijke管内的火焰区域,而后经一组球面反射镜汇聚于主反射镜(二)焦点位置的光刀处,纹影图像由光刀后方的相机记录, 刀口布置的方向与狭缝方向设置相同。
拍摄火焰红外图像的CCD光谱如图 4所示,红外CCD的平均响应主要集中在近红外波段。拍摄中对整体火焰结构进行拍摄,相机至于视窗前方一定位置(物距约22cm),另一视窗的内侧加设黑色云母板,目的是为拍摄过程提供慢反射背景,避免由于石英玻璃表面光线发生镜面反射对图像记录构成影响。
1.4 相位锁定同步拍摄系统
通过压力相位锁定的记录方式,对压力信号和红外图像进行同步触发,其系统如图 5所示。
信号发生器1发出连续的TTL电平触发信号,经过信号开关2后被分为2路,一路与红外相机3的触发接口连接,另一路与数据采集器4的触发通道连接,压力传感器6与数据采集器4的其它采集通道连接。相机和数据采集器设置为单次触发模式,对应的作业方式为:由TTL电平信号触发后,仪器以自身的采样频率连续采集信号。具体的试验过程为:断开触发信号开关,分别激活红外相机和数据采集器的采集状态,当二者均为等待触发信号的状态时,闭合信号开关,TTL信号分别触发红外相机和数据采集器,实现了红外时序图像和声场压力信号的同步连续采集。为确保触发时刻足够短,红外相机和数据采集器均设置为上升沿(Rising)触发,信号发生器TTL电平信号的脉冲频率为8kHz,对应的信号周期远小于声周期。此外,红外相机的曝光时间和快门等待时间分别为0.1和10ms,数据采集器的采样频率与相机的快门时间相匹配,设置为10 kHz。
2 结果与讨论
2.1 红外和纹影技术的优势
火焰的红外图像用于描述声周期下火焰锋面的运动特性和热量演化,具有以下优势:
(1)在红外成像的物理过程上,火焰自发光红外光谱波段的光辐射具有较高的能量。火焰红外图像是由燃烧区域能够辐射红外光线的物质在CCD上成像而得到的,对于甲烷-空气的燃烧反应而言,二氧化碳和水是主要的燃烧产物,二者主要集中在红外波段[6-7],即便是缺氧条件下形成碳烟,其辐射强度最大值也处于红外区[8],因此,形成甲烷火焰红外图像的光线具有较大的辐射强度,图像中的平均灰度值更高,火焰区域与背景区域的对比度更大。
(2)甲烷火焰的红外图像,在物理意义上,描述了燃烧反应后产生高温产物的分布及运动。这些高温产物包围在燃烧反应区域附近,近似地表征了反应区域的运动过程。顾璠等[9]利用湍流火焰的红外图像表征大涡逆序结构的做法可说明这一点。
图 6对可见图像和红外图像进行了对比,对火焰整体分别采用单反和红外相机进行拍摄。其中,可见图像通过单反相机获得,快门时间为0.25ms,红外图像的快门时间为0.1ms。外轮廓上,二者都具有对称的尖头(Cusp)结构,数量和空间分布几乎相同,但红外图像比单反图像中的尖头结构尺寸更大,对应总体的火焰幅面更宽。因此,对于描述锋面运动过程中尖头结构的发展与变化而言,红外图像更能够突显特征结构的发展。
火焰的纹影图像则可用于描述径向的热量传递,根据Dunn-Rankin等[10]的观点,纹影信息描述了燃烧区域温度梯度较大的位置,那么火焰纹影轮廓将反映火焰前锋附近的温度梯度,由此可推导对应的热量传递。最终,火焰红外和纹影2类图像的综合信息将有效揭示甲烷火焰的热量演化过程及特征。
2.2 部分预混火焰的锋面结构
图 7给出了采用相位锁定方法拍摄的不同时刻的红外火焰图像,以此描述声场作用下部分预混火焰的锋面运动过程,其中红外相机的曝光时间为0.1ms。靠近燃烧器出口附近的火焰区域为火焰的上游,远离燃烧器出口的火焰顶部区域处为火焰的下游,对不同高度的尖头结构进行编号,在45°相位图像中由上游至下游的4层尖头依次编号为1~4号,1号和2号定义为上游,3号和4号定义为下游,以此描述锋面运动过程和结构特征。
由图 7可知,对于部分预混火焰锋面而言,各个相位下的火焰具有相似的空间结构特征,尖头结构仍是其在声激励下的典型动力学特征。火焰结构呈现层状分布,锋面上尖头的数量相对固定,连续的火焰区域存在约4对;火焰下游的幅面更宽,上游的幅面相对较窄,同时距离下游火焰顶部一定距离处均存在弧线状孤立区域。
对预混火焰的研究表明[3],预混锋面在声场作用下存在2种运动行为,包括沿锋面法向的传播运动和沿锋面切向的对流运动。对流波是2种运动形式产生的主要锋面几何特征,对应着预混锋面上轴对称的褶皱和尖头结构。对流波的波长决定锋面切向相邻2个尖头结构之间的距离,由锋面对流运动形成;对流波的幅值决定尖头结构的明显程度,由锋面传播运动形成。
而根据扩散锋面的研究结果,声激励的作用下,锋面运动由单方向对流转为往复式的变速运动。Hardalupas Y等[11]提出的唯象模型(Phenomenological model)描述,这种往复式的变速运动将在扩散锋面的外围区域形成一系列的涡街,相邻2个涡团的卷吸作用产生了锋面上的尖头结构。因此,尖头结构的明显程度应主要与涡团的运动行为有关。
综上所述,对于部分预混火焰而言,其在声作用下的结构应为对流波和涡团2种行为的综合结果。根据对流波的运动特征,声周期的不同时刻,锋面形成脉动位移最大值的位置不同,在火焰的下游和上游区域分别形成火焰基部和顶部的变形。因此,将火焰锋面上有尖头出现的区域进行划分,4个区域的出现描述了沿锋面法向的传播运动。
2.3 火焰锋面的运动行为
如图 7所示,上半周期:前1/4周期内(0°~90°相位),上游的1号尖头包围燃烧器外壁的面积逐渐增大,红外轮廓的幅面增加,最低点略有下降;2号尖头的尖角上翘,在空间高度略有降低;3号尖头和4号尖头之间连接部分不断减小,4号尖头在空间上绝对高度保持不变,3号尖头处于下降过程,二者之间呈现了明显的分离特征,从整体上看,火焰在该时间段内呈现下降的运动形式。随后(90°~180°相位),90°相位下的4号尖头从火焰本体上脱离,形成了空间上孤立的红外区域,如180°相位的图像所示,这一过程描述了高温产物随火焰锋面运动最终脱离火焰本体的过程;同时,火焰基部红外轮廓的幅面逐渐减小,最低点不断上升,包围燃烧器外壁的特征消失,基部轮廓的宽度在135°后几乎接近燃烧器内径,尖头结构的尖角收缩,整个火焰呈拉伸状,运动形式表现为上升。
下半周期:锋面由向上运动转为向下运动,向上拉伸过程中逐渐减小的尖头尖角在下降过程中不断扩大,在火焰下游区域,顶部逐渐展开,幅面变宽,表现为压缩的锋面变形特征。此后(270°~360°相位),火焰锋面再次上升,呈现出90°~180°相位时间段内的锋面运动行为。
2.4 火焰的热量演化
整个周期内火焰区域的高温产物经历了2次下降和上升的过程,同时扩散锋面上尖头结构周期性地发展。结合锋面的运动行为和尖头结构的发展过程,认为声波作用下部分预混火焰沿轴向的热量演化路径为:火焰基部处的反应热,以高温产物为载体,由上游至下游运动,沿程汇聚不同位置(各层各处的尖头结构)产生的反应热,最终从火焰顶部脱离反应区域,以较大的表面积向下游的气体传递。
图 8为拍摄的火焰时间序列纹影图像,使用的纹影系统光刀沿水平方向移动,因此图中亮暗区域及其对比度反映的是火焰区域沿水平方向的温度梯度,目的是说明径向的热量演化。总体结构上,包括内锥和外轮廓2部分,外轮廓的空间分布方式与红外图像相似,以一层一层的形式存在,但红外图像的尖头结构在纹影图像上表现为“手指状”结构。周期运动过程中,纹影图像的各层结构也存在下降和上升2个阶段,下面就上述存在的异同及反映的实际物理现象进行讨论。
由于红外图像表征的是高温产物的空间分布,而纹影图像描述的是水平方向温度梯度的分布,从热量与温度梯度的关系考虑,尖头结构作为“热源”,所形成的径向温度梯度轮廓线表现“手指状”的弧线是合理的,而“热源”一层一层分布也决定了温度梯度的分布形式。为验证上述分析的正确性,拍摄了纹影和红外信息一体化的图像。具体的做法是,以红外相机作为纹影系统的记录装置,将图 3中的狭缝2尽量关小,使穿过火焰区域平行光的亮度与火焰自身发出的红外光线亮度接近,同时设置红外相机的增益程度,提高CCD的感光度(ISO),获得如图 9所示的图像。图 9中,明亮的区域与红外信息对应,依然表现为尖头状的特征,亮暗相间的区域为纹影信息,外轮廓线为“手指状”。“手指状”的纹影轮廓与“尖状体”位于同一高度,二者在“分层”和“配合”2个方面都比较吻合。
由此可知,声波作用下的部分预混火焰区域中,由于高温产物与周围二次空气之间存在温度差和流速差,在火焰锋面的运动过程中存在径向的热量传递,结构上分层的尖头结构形成分层的“手指状”温度梯度轮廓,升降式的运动方式使不同时刻的径向热量传递过程差别较大。
3 结论
(1)声波作用下的部分预混火焰,在结构上表现为多层对称的尖头结构,尖头结构与对流波和涡团运动有关;锋面在声波振荡的周期内经历2次上升和下降的过程,热量演化由这些结构发展过程和锋面运动方式决定。
(2)轴向热量传递:火焰基部处的反应热,以高温产物为载体,由上游至下游运动,沿程汇聚不同位置(各层尖头结构)产生的反应热,最终从火焰顶部脱离反应区域,以较大的表面积向下游的气体传递。
(3)径向热量传递:形成分层的“手指状”温度梯度轮廓,发生径向传热过程的轴向位置随相位的改变周期性变化。
-
表 1 Kriging模型预测误差
Table 1 Prediction error of Kriging model
最大旁瓣水平/dB 分辨率/m 阵列性能分析结果 -14.696 0.1467 代理模型预测结果 -14.680 0.1460 相对误差 0.07% 0.45% 表 2 分辨率优化阵列设计变量
Table 2 Design variables of resolution optimized array
设计变量 n k a q 变量范围 8~10 0.60~1.40 0.050~0.200 0.600~1.400 初始阵列变量值 9 1.00 0.150 1.000 优化阵列变量值 9 0.89 0.147 0.716 表 3 初始阵列与分辨率优化麦克风阵列性能
Table 3 Performance of initial array and resolution optimized microphone array
初始阵列 优化阵列 频率/
kHz最大旁瓣水平/dB 分辨率
/m最大旁瓣水平/dB 分辨率
/m1 -17.412 0.388 -14.237 0.341 2 -16.882 0.187 -14.330 0.164 3 -15.535 0.123 -11.868 0.108 4 -15.252 0.089 -11.934 0.080 5 -14.771 0.070 -11.992 0.061 6 -14.188 0.057 -10.235 0.048 表 4 初始阵列与最大旁瓣水平优化麦克风阵列性能
Table 4 Performance of initial array and MSL optimized microphone array
初始阵列 优化阵列 频率
/kHz最大旁瓣水平/dB 分辨率
/m最大旁瓣水平/dB 分辨率
/m1 -17.412 0.388 -20.537 0.412 2 -16.882 0.187 -21.922 0.198 3 -15.535 0.123 -19.346 0.129 4 -15.252 0.089 -16.105 0.094 5 -14.771 0.070 -15.698 0.074 6 -14.188 0.057 -14.009 0.061 表 5 初始阵列与优化矩形麦克风阵列性能
Table 5 Performance of initial array and optimized rectangular microphone array
初始阵列 优化阵列 频率
/kHz最大旁瓣水平/dB 分辨率
/m最大旁瓣水平/dB 分辨率
/m1 -12.545 0.302 -16.101 0.320 2 -11.786 0.146 -16.012 0.155 3 -11.798 0.094 -15.450 0.101 4 -11.800 0.068 -15.252 0.073 5 -11.806 0.053 -15.147 0.056 6 -6.523 0.042 -11.271 0.046 -
[1] Brooks T F, Marcolini M A, Pope D S. A directional array approach for the measurement of rotor noise source distributions with controlled spatial resolution[J]. Journal of Sound and Vibration, 1987, 112(1):192-197. DOI: 10.1016/S0022-460X(87)80105-0
[2] Ahmad D V, Radhakrishnan S. Review of acoustic measurement techniques in wind tunnels[R]. AIAA-2000-2204, 2000.
[3] Wei Q K. Compressive sensing based beamforming and its application in aeroacoustic experiment[R]. AIAA-2014-2918, 2014.
[4] 周家检, 郝璇, 张卫民, 等.相阵列技术在民机机体气动噪声研究中的应用[J].空气动力学学报, 2016, 34(1):91-97. http://d.old.wanfangdata.com.cn/Periodical/kqdlxxb201601016 Zhou J J, Hao X, Zhang W M, et al. Application of phased array technique in the research of civil airplane airframe noise[J]. Acta Aerodynamic Sinica, 2016, 34(1):91-97. http://d.old.wanfangdata.com.cn/Periodical/kqdlxxb201601016
[5] 陈正武, 王勋年, 李征初, 等.基于声学风洞的麦克风阵列测试技术应用研究[J].实验流体力学, 2012, 26(3):84-90. DOI: 10.3969/j.issn.1672-9897.2012.03.016 Chen Z W, Wang X N, Li Z C, et al. Application investigation of microphone array measuring and testing technique in anechoic wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2012, 26(3):84-90. DOI: 10.3969/j.issn.1672-9897.2012.03.016
[6] Gade S, Hald J. Array designs optimized for both low-frequency NAH and high-frequency Beamforming[R]. SAE Technical Paper 2005-01-4014, 2015.
[7] Underbrink J R. Practical considerations in focused array design for passive broadband source mapping applications[D]. Pennsylvania: Pennsylvania State University, 1995.
[8] 黄奔.气动噪声源的麦克风阵列识别定位技术研究[D].绵阳: 中国空气动力研究与发展中心, 2014. Huang B. Investigation of aerodynamic noise sources identification technique based on microphone arrays[D]. Mianyang: China Aerodynamics Research and Development Center, 2014.
[9] William M, Humphreys J R. Design and use of microphone direction arrays for aeroacoustic measurement[R]. AIAA-98-0471, 1998.
[10] Thomas F B, William M, Humphreys J R. Effect of directional array size on the measurement of airframe noise components[R]. AIAA-99-1958, 1999.
[11] Takeshi Ito. Aeroacoustic noise measurements in aerodynamic low-speed wind tunnels[C]//Proc of the 26th International Congress of the Aeronautical Sciences. 2008.
[12] Sijtsma P. Phased array beamforming in wind tunnels[R]. NLR-TP-2006-732, 2006.
[13] 陈宝, 李周复, 谭啸, 等.声衬试验段环境下航空声学定位试验技术研究[J].实验流体力学, 2015, 22(5):78-83. http://www.syltlx.com/CN/abstract/abstract10879.shtml Chen B, Li Z F, Tan X, et al. Investigation of aeroacoustic localization technique in lining test section[J]. Journal of Experiments in Fluid Mechanics, 2015, 29(5):78-83. http://www.syltlx.com/CN/abstract/abstract10879.shtml
[14] 郭旺柳, 宋文萍, 许建华, 等.旋翼桨尖气动/降噪综合优化设计研究[J].西北工业大学报, 2012, 30(1):73-79. Guo W L, Song W P, Xu J H, et al. An effective aerodynamics/acoustic optimization of blade tip planform for helicopter rotors[J]. Journal of Northwestern Polytechnical University, 2012, 30(1):73-79.
[15] 聂雪媛, 刘中玉, 杨国伟.基于Kriging代理模型的飞行器结构刚度气动优化设计[J].气体物理, 2017, 2(2):8-16. http://d.old.wanfangdata.com.cn/Periodical/qtwl201702002 Nie X Y, Liu Z Y, Yang G W. Aircraft structure stiffness and aerodynamics optimization design based on Kriging surrogate model[J]. Physics of Gases, 2017, 2(2):8-16. http://d.old.wanfangdata.com.cn/Periodical/qtwl201702002
[16] Kanazaki M, Tanaka K, Jeong S, et al. Multi-objective aerodynamic optimization of elements' setting for high-lift airfoil using Kriging model[R]. AIAA-2006-1471, 2006.
[17] Han Z H. Improving adjoint-based aerodynamic optimization via gradient-enhanced Kriging[R]. AIAA-2012-0670, 2012.
[18] 刘俊.基于代理模型的高效气动优化设计方法及应用[D].西安: 西北工业大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10699-1016909342.htm Liu J. Efficient surrogate-based optimization method and its application in aerodynamic design[D]. Xi'an: Northwestern Polytechnical University, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10699-1016909342.htm
[19] 韩忠华. Kriging模型及代理优化算法研究进展[J].航空学报, 2016, 37(11):3197-3225. http://d.old.wanfangdata.com.cn/Periodical/hkxb201611001 Han Z H. Kriging surrogate model and its application to design optimization:a review of recent progress[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(11):3197-3225. http://d.old.wanfangdata.com.cn/Periodical/hkxb201611001
[20] 夏露, 王丹.基于Kriging自适应代理模型的气动优化方法[J].航空计算技术, 2013, 43(1):13-17. DOI: 10.3969/j.issn.1671-654X.2013.01.004 Xia L, Wang D. Aerodynamic optimization method based on Kriging adaptive surrogate model[J]. Aeronautical Computing Technique, 2013, 43(1):13-17. DOI: 10.3969/j.issn.1671-654X.2013.01.004
[21] Zhang Y, Han Z H, Shi L X, et al. Multi-round surrogate-based optimization for benchmark aerodynamic design problems[R]. AIAA-2016-1545, 2016.
[22] Thomas J M. Aeroacoustic measurements[M]. Berlin:Springer, 2002.
-
期刊类型引用(2)
1. 孟晟,周昊,陶成飞,刘子华. 入口段长度变化对于非预混火焰声场响应的影响及其被动控制. 中国电机工程学报. 2019(01): 212-218+335 . 百度学术
2. 钟英杰,王根娟,王明晓,王炜豪,邓凯. 基于卡塞格林原理的火焰自由基测量系统. 浙江大学学报(工学版). 2017(05): 1044-1050 . 百度学术
其他类型引用(2)