CN113588204B - 一种进气道激波边界层干扰特性测量方法 - Google Patents
一种进气道激波边界层干扰特性测量方法 Download PDFInfo
- Publication number
- CN113588204B CN113588204B CN202110738673.7A CN202110738673A CN113588204B CN 113588204 B CN113588204 B CN 113588204B CN 202110738673 A CN202110738673 A CN 202110738673A CN 113588204 B CN113588204 B CN 113588204B
- Authority
- CN
- China
- Prior art keywords
- pressure
- boundary layer
- image
- air inlet
- frame type
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000035939 shock Effects 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012360 testing method Methods 0.000 claims abstract description 94
- 238000005259 measurement Methods 0.000 claims abstract description 38
- 239000003973 paint Substances 0.000 claims abstract description 32
- 238000009826 distribution Methods 0.000 claims abstract description 30
- 238000005507 spraying Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 14
- 238000011161 development Methods 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 230000003287 optical effect Effects 0.000 claims description 27
- 239000002245 particle Substances 0.000 claims description 26
- 230000010349 pulsation Effects 0.000 claims description 12
- 238000000889 atomisation Methods 0.000 claims description 11
- 230000004907 flux Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 7
- 239000000700 radioactive tracer Substances 0.000 claims description 6
- 238000011426 transformation method Methods 0.000 claims description 6
- 239000011248 coating agent Substances 0.000 claims description 5
- 238000000576 coating method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000012800 visualization Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 239000007850 fluorescent dye Substances 0.000 claims description 3
- 230000035515 penetration Effects 0.000 claims description 3
- 229920000642 polymer Polymers 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000002834 transmittance Methods 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 2
- 238000000926 separation method Methods 0.000 abstract description 8
- 230000000007 visual effect Effects 0.000 abstract description 8
- 230000007547 defect Effects 0.000 abstract description 5
- 238000000917 particle-image velocimetry Methods 0.000 description 15
- 239000007921 spray Substances 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 8
- 238000000354 decomposition reaction Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000004140 cleaning Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- CSCPPACGZOOCGX-UHFFFAOYSA-N Acetone Chemical compound CC(C)=O CSCPPACGZOOCGX-UHFFFAOYSA-N 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000000428 dust Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 229920000742 Cotton Polymers 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000003517 fume Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012113 quantitative test Methods 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000007665 sagging Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
- G01M9/04—Details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/065—Measuring arrangements specially adapted for aerodynamic testing dealing with flow
- G01M9/067—Measuring arrangements specially adapted for aerodynamic testing dealing with flow visualisation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/08—Aerodynamic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种进气道激波边界层干扰特性测量方法，首先根据风洞试验模拟相似性准则为发动机进气道设计缩比框架式试验模型；在缩比框架式试验模型的待测量区域喷涂压敏漆，构建PSP测量系统，用于测量进气道内表面压力分布特性、激波的大小、位置、激波边界层的流态发展和表面涡流近似摩擦力线分布特性；搭建PIV测量系统，用于测量进气道内流场的速度场信息和波系结构；搭建纹影测量系统，用于测量进气道内层流湍流流态变化和边界层厚度。本发明能够有效弥补现有进气道实验项目参数测量的不足，对典型激波边界层干扰流动结构、分离特征、激波串流场提供全局、直观辨识方法。
Description
技术领域
本发明属于航空航天实验技术领域，涉及一种进气道激波边界层干扰特性测量方法。
背景技术
激波边界层干扰是超/高超声速飞行器中普遍存在的一种复杂流动现象，这一现象往往会对飞行器的气动力/热分布产生显著的影响，从而影响整个飞行器的气动性能和防热特性。尤其是对于复杂来流条件下跨声速机翼和叶栅、超燃冲压发动机进气道和隔离段等内流场激波边界层干扰特性分析中，缺乏对复杂流场的波系结构、流态变化、全场动态压力变化等多物理场参数非定常特性的全方位观测和分析，而现有的研究方法对典型激波边界层干扰流动结构、分离特征、激波串流场缺乏全局、直观辨识方法。传统测量方法一般采用单点传感器测量发动机内部压力特性来预测流动分离情况，由于布点不可能太多，对于分离位置判断精度很低，而且对大面积压力脉动无能为力，结果非常不直观而且获得信息量少。进气道实验项目大部分仅限于总压恢复系数、流量系数、进气道出口截面流场均匀性三项参数测量，对于进气道内表面全局压力分布、内流场速度场信息、波系结构、层流湍流流态变化、边界层厚度、激波边界层三维壁面效应、表面涡流近似摩擦力线分布等内流特性无法有效观测。
发明内容
本发明解决的技术问题是：克服现有技术的不足，提出一种进气道激波边界层干扰特性测量方法。
本发明解决技术的方案是：
一种进气道激波边界层干扰特性测量方法，该方法的步骤包括：
步骤一，根据风洞试验模拟相似性准则为发动机进气道设计缩比框架式试验模型，在缩比框架式试验模型上设计光学可视化窗口，用于观测缩比框架式试验模型中的内流场；
步骤二，在缩比框架式试验模型的待测量区域喷涂压敏漆，所述压敏漆为聚合物形成的多孔性涂层结构，孔中包含有荧光探针分子；
步骤三、选择PSP测量使用的第一图像采集设备和第一光源，对压敏漆涂料的动态响应特性进行标定，确定第一图像采集设备的拍摄距离、曝光时间，第一光源的照射距离和光强度、压敏漆涂料的发光-压力关系式；
步骤四、搭建PSP测量系统，所述PSP测量系统包括风洞、高超声速喷管、缩比框架式试验模型、第一光源和第一图像采集设备；
将缩比框架式试验模型和高超声速喷管放置在风洞试验段，且缩比框架式试验模型放置于高超声速喷管的均匀区域内，第一光源和第一图像采集设备放置于风洞光学顶窗外侧，透过风洞光学顶窗，能够观察到缩比框架式试验模型的光学可视化窗口；依据步骤三确定的第一图像采集设备参数、第一光源参数对图像采集设备和光源进行设置；
步骤五、开启PSP测量系统各设备，进行PSP测量，获得进气道内表面压力分布特性、激波的大小、位置、激波边界层的流态发展和表面涡流近似摩擦力线分布特性；
步骤六、搭建PIV测量系统，所述PIV测量系统包括风洞、PIV雾化示踪粒子发生器、高超声速喷管、缩比框架式试验模型、第二光源、第二图像采集设备以及同步触发器；
将缩比框架式试验模型和高超声速喷管放置在风洞试验段，且缩比框架式试验模型放置于高超声速喷管的均匀区域内，第二光源布置在风洞试验段顶部，第二图像采集设备位于风洞试验段一侧观察窗；PIV雾化示踪粒子发生器放在高超声速喷管的前端；同步触发器用于控制第二光源和第二图像采集设备同步工作；
步骤七、先运行第二图像采集设备，拍摄缩比框架式试验模型内流场；然后开启PIV雾化示踪粒子发生器，使其发出的粒子随喷管流场一起进入风洞，第二图像采集设备捕捉粒子从无到有全过程；根据拍摄的图像，利用自相关算法获取进气道内流场的速度场信息和波系结构；
步骤八、搭建纹影测量系统，所述纹影测量系统包括风洞、高超声速喷管、缩比框架式试验模型以及双镜平行光式纹影仪；
将缩比框架式试验模型和高超声速喷管放置在风洞试验段，且缩比框架式试验模型放置于高超声速喷管的均匀区域内，双镜平行光式纹影仪放置于风洞光学顶窗外侧；
步骤九、双镜平行光式纹影仪以平行光入射于风洞试验段，双镜平行光式纹影仪内的单反相机捕捉缩比框架式试验模型内部的密度场图像，获得进气道内层流湍流流态变化和边界层厚度。
所述步骤一中，光学窗口对于激发光源的波长和受激辐射光波长具有良好的透过率，且光学窗口各处等厚。
所述步骤五中，所述进行PSP测量，获得进气道内表面压力分布特性、激波的大小、位置、激波边界层的流态发展和表面涡流近似摩擦力线分布特性的方法如下：
(5.1)开启PSP测量系统各设备，来流通过高超声速喷形成高超声速流场，连续采集N张吹风状态下的图像；
关闭高超声速喷管，采集M张无风状态下的参考图像，N>>M；
(5.2)通过图像比获取流场的壁面表面基于时间序列的压力图像序列，实现方法如下：
计算M张参考图像上同一个像素点的平均灰度值，得到一张参考图像平均值图像；
利用如下公式计算M张参考图像上像素点(xi,yi)的平均灰度值g；
Gm(xi,yi)为第m张参考图像中像素点(xi,yi)的灰度值，i∈[1,I]，I为图像中的像素点总数；
第j张吹风状态下图像上每个像素点的灰度值除以平均值图像上对应像素点的灰度值，得到第j张吹风状态下的压力图像；
j遍历[1,N]，得到N个相对于时间序列的压力图像序列；
(5.3)利用正交变换方法将压力图像序列转换为模态序列，截取前N阶模态，从前N阶模态中剔除与压力脉动无关的噪声和无效模态，将剩余与压力脉动相关的模态通过反正交变换重构得到与原来同时间分布的压力图像序列；
(5.4)根据(5.3)重构得到的压力图像序列，获得每一点的压力值，从而得到压力分布特性、激波的大小、位置、激波边界层的流态发展。
所述步骤(5.3)中，利用正交变换方法将压力图像序列转换为模态序列的方法如下：
寻找一组最佳的标准正交基；
使压力图像序列中每个压力图像数据在该标准正交基上的投影依次递减，每个压力图像数据在该标准正交基上的投影即为模态序列。
所述步骤(5.3)中，前N阶模态的能量之和超过总能量的90％。
所述步骤五中，所述进行PSP测量，获得表面涡流近似摩擦力线分布特性的方法如下：
采用表面近似摩擦力线提取方法，对流场的NS方程进行处理，得到压力梯度-摩擦力方程；
当边界层涡通量未知时，通过压力梯度解算压力梯度-摩擦力方程，求得摩擦力矢量的相对大小和方向。
对流场的NS方程进行处理，得到压力梯度-摩擦力方程的处理方式如下：
流场的NS方程如下：
在近壁面底层，NS方程简化为：
两端同时乘以摩擦力线方程τ＝μω×n得：
引入涡通量Ω＝|ω|2/2，有：
方程中第二项为相对小量，引入边界层涡通量fΩ，得到压力梯度-摩擦力方程为：
其中，ρ为密度，u为速度，t为时间，为导数，/>为梯度，p为压力，μ为粘性系数，为拉普拉斯算子，θ为速度的旋度，F为单位质量彻体力，z为边界层厚度方向，τ为粘性摩擦力，ω为涡量，n为单位外法向矢量，kω为涡量系数。
第一光源为波长400nm的紫外光源，第一图像采集设备为普通相机。
第二光源为YAG双曝光激光器，激光器出光孔处装有导光臂和集成片光源；第二图像采集设备为跨帧数字相机。
本发明与现有技术相比的有益效果是：
(1)本发明方法获得的全局表面压力图像空间分辨率高，压力信息全面，流场可视化更直观，结合粒子图像测速和纹影等方法获取的空间脉动压力场、空间速度场和密度场信息，对进气道内部复杂流场的非定常特性和流态变化及精细结构进行全方位观测和分析，同时为激波边界层干扰、空间拓扑结构分析提供了强有力的分析手段。传统流场测量方法中，只能通过测压传感器或测温传感器获取关心测量区域的单点压力或单点温度，采用压敏漆测量方法获得的流场表面压力是全局测量的，大面积的定量测试方法，现有技术不存在针对连续大面积的表面压力分析和处理方法，目前没有其他方法能够获得这一观测结果。
(2)基于模态进行表面全局压力信息分解和重构能更清晰、高效的获取流场的波系结构和能量分布，是传统测量方法和处理方法无法获取的。本发明中采用高效的降阶方法和psp非接触测量手段相结合方式，对激波边界层干扰引起的表面压力和全局分布特性进行分解再重构，能够清晰地捕捉到激波边界层干扰的区域和位置，帮助判定流场内部结构和特点，目前未有针对进气道内激波边界层和激波串结构分析等相关工作开展。
(3)本发明能够有效观测进气道内表面全局压力分布、内流场速度场信息、波系结构、层流湍流流态变化、边界层厚度、激波边界层三维壁面效应、表面涡流近似摩擦力线分布等内流特性，弥补现有技术的不足，对典型激波边界层干扰流动结构、分离特征、激波串流场提供全局、直观辨识方法。
附图说明
图1为发明提供的方法原理图；
图2为本发明流程图；
图3为PSP测量系统示意图；
图4为PSP获取的表面结构和幅值与纹影对比结果；
图5为粒子图像拍摄显示结果；
图6为进气道内流的波系结构及激波边界层干扰；
图7为POD模态分解重构得到主要压力脉动位置；
图8为表面近似摩擦力大小相对分布。
具体实施方式
下面结合附图对本发明作进一步阐述。
本发明针对复杂流场激波边界层干扰特性分析测量需求，通过对模型进行光学可视化设计加工结合先进流动测试技术，开展对复杂流场激波边界层干扰的波系结构、流态变化、全场动态压力变化等非定常特性的全方位观测和分析，弥补数值模拟和常规风洞试验的不足，形成模型复杂流场可视化测试能力，为翼型、叶栅和进气道的设计、流动控制及改进提供原理性试验依据。图1为本发明的原理图。
为了分析激波边界层干扰特性，本发明提出基于压敏漆和全场时间序列压力分布模态分解技术的方法，该方法是一种对大面积压力分布的高效的降阶方法，其核心思想是寻找一组最佳的标准正交基，使得样本数据在该标准正交基上的投影依次迅速递减，通过分析激波边界层干扰结构特性和能量特征，对压力相关分解的能量模态进行分析和截取，剔除流场噪声、随机脉动等不相干模态，从而可以用较少的基展开获得较高阶数据的近似描述，尤其对激波入射结构判定、影响区域、三维效应等气动特征的描述更直观和定量。用模态分解方法对压力场图像进行分解后，得到了发生压力脉动的位置，波系结构特性和强度，不仅能够获得进气道气动特性研究方法，而且还能掌握流场非定常气动载荷流场显示技术，补充典型激波边界层干扰特性非定常气动载荷数据库。
本发明方法，需要如下准备条件：首先准备三面带可视透明窗口的进气道缩比框架式试验模型，将进气道的主体结构分割为若干个块，对每个块建立缩比试验模型，若干个缩比实验模型拼接后形成缩比框架式试验模型。
将进气道缩比框架式试验模型待测量底部区域喷涂有压敏漆，喷涂好的缩比框架式试验模型安装在风洞试验段内，通过测试设备和风洞进行流场测量，通过侧窗拍摄流场粒子图像和纹影结果，通过顶部测量底面表面发光图像，最终获取进气道内表面全局压力分布、内流场速度场信息、波系结构、层流湍流流态变化、边界层厚度、激波边界层三维壁面效应、表面涡流近似摩擦力线分布等内流特性。
如图2所示，具体步骤如下：
(1)缩比框架式试验模型设计：根据风洞试验模拟相似性准则为发动机进气道设计缩比框架式试验模型，合理简化试验模型，满足PSP、PIV、纹影三种非接触流动显示技术的测量需求。在缩比框架式试验模型待测量区域三面设计光学窗口。该光学窗口对于激发光源的波长和受激辐射光波长须具有良好的透过率，且光学窗口尽量为等厚，以减少光学畸变造成的影响。
(2)在缩比框架式试验模型的底部待测量区域喷涂压敏漆
1.取步骤(1)设计后的模型进行压敏漆涂料喷涂作业，压敏漆涂料采用PtTFPP作为发光基团，其稳定性较强。涂料的激发光源为波长400nm紫外光源，激发出的荧光发射波长区域为600-700nm。压敏漆涂料为由聚合物形成的多孔性的涂层结构，荧光探针分子(PtTFPP)包含在这些微孔之中，增大了空气接触面积，使得响应扩散率增加即降低了反应时间。喷涂前清洁缩比框架式试验模型表面，用酒精、丙酮去除油污灰尘，后用质量好的清洁布擦除，避免使用棉花或纱布。喷涂前在通风柜开启条件下打扫实验室，保证喷涂室内清洁无灰尘，保证喷涂室内不会因为喷枪喷射气流过程中引起灰尘或颗粒混杂进涂料中。喷涂气源需要干净气源，可用2级过滤器清洁气源。喷涂时，选择高压喷枪喷涂气压3～4Bar，喷涂距离18～23cm。如操作人员不够娴熟，可考虑加工固定距离的移动装置辅助喷涂。喷枪与待喷件保持一定垂直距离，匀速往复运行。喷涂过程中，喷枪与被喷件保持垂直90°角，以保证涂料完全均匀的喷涂在模型表面。喷枪要求喷幅重叠50％以上。喷枪移动保持平稳匀速，喷涂速度为30～50cm/s。喷涂速度太快会令油漆太干，表面粗糙，容易产生橘皮；而喷涂太慢较容易产生流挂及下坠。因此将喷枪的喷涂移动速度控制在30～50cm/s可有效防止褶皱、表面粗糙和橘皮现象。喷涂结束后静置10分钟左右，然后放进高温烘箱在60℃下加热5小时，待漆膜完全固化后，停止加热，自然冷却，待用。
(3)对步骤(2)喷涂后的缩比框架式试验模型进行动态标定。动态标定试验应用专门设计的动态标定装置进行标定，可通过激波管、驻波管、高速喷流或振荡射流发生器产生压力阶跃或高频射流的方法对涂料的动态压力响应进行标定。涂料的动态响应特性需满足风洞试验流场特征频率要求。对动态标定后的模型进行发光-压力对应关系校准，校准压力和温度范围大于整个试验测量范围，校准后得到的发光-压力关系式可用于表面压力的定量测量。得到第一图像采集设备3的拍摄距离、曝光时间，第一光源的照射距离和光强度、涂料的发光-压力关系式。
通过动态标定装置产生高频射流或阶跃压力的方法对压敏漆涂层的压力响应进行标定，响应时间在微秒量级的动态压敏漆才可用于本方法的高超声速流场进气道内表面的动态压力获取。
(4)按照图3搭建PSP测量系统，PSP测量系统包括喷管1、缩比框架式试验模型2、第一图像采集设备3和第一光源4，搭建PSP测量系统时，将缩比框架式试验模型和高超声速喷管放置在风洞试验段中，将缩比框架式试验模型放置于高超声速喷管的均匀区域内，第一光源和第一图像采集设备位于风洞光学顶窗5外侧，且透过风洞光学顶窗5，能够观察到缩比框架式试验模型的光学窗口7及缩比框架式试验模型的内流场。图3中，6为压敏漆涂层，光学窗口5材料为石英玻璃。每次试验都连续采集1000张吹风状态下的图像，然后采集100张无风状态下的参考图像。第一图像采集设备3采用35mm或50mm焦距的镜头，在风洞运行时保证采集图像不过曝前提下将光圈调到最大以最大程度的接受PSP辐射出的荧光。在镜头前安装650±10nm的带通滤光片，防止其他波段光线的干扰。依照步骤三校准测量中确定的第一图像采集设备和第一光源参数调节设备进行风洞试验，并通过图像比最终获取流场的壁面表面压力分布数据。风洞试验得到所有图像后进行数据处理。将风洞停止运行后采集的所有参考图像做平均。风洞运行采集到的图像保持时间连续性。然后用光强-压力转换公式计算基于时间序列的压力数据。
进行PSP测量，获得进气道内表面压力分布特性、激波的大小、位置、激波边界层的流态发展的方法如下：
(4.1)开启PSP测量系统各设备，来流通过高超声速喷形成高超声速流场，连续采集N张吹风状态下的图像；
关闭高超声速喷管，采集M张无风状态下的参考图像，N>>M；
(4.2)通过图像比获取流场的壁面表面基于时间序列的压力图像序列，实现方法如下：
计算M张参考图像上同一个像素点的平均灰度值，得到一张参考图像平均值图像；
利用如下公式计算M张参考图像上像素点(xi,yi)的平均灰度值g；
Gm(xi,yi)为第m张参考图像中像素点(xi,yi)的灰度值，i∈[1,I]，I为图像中的像素点总数；
第j张吹风状态下图像上每个像素点的灰度值除以平均值图像上对应像素点的灰度值，得到第j张吹风状态下的压力图像；
j遍历[1,N]，得到N个相对于时间序列的压力图像序列；
(4.3)利用正交变换方法将压力图像序列转换为模态序列，截取前N阶模态，从前N阶模态中剔除与压力脉动无关的噪声和无效模态，将剩余与压力脉动相关的模态通过反正交变换重构得到与原来同时间分布的压力图像序列；
利用正交变换方法将压力图像序列转换为模态序列的方法如下：
寻找一组最佳的标准正交基；
使压力图像序列中每个压力图像数据在该标准正交基上的投影依次递减，每个压力图像数据在该标准正交基上的投影即为模态序列。
前N阶模态的能量之和超过总能量的90％。
图7为POD模态分解重构得到主要压力脉动位置。
(4.4)根据(4.3)重构得到的压力图像序列，获得每一点的压力值，从而得到压力分布特性、激波结构的大小、位置和激波边界层的流态发展。
对于进气道构型的典型激波边界层干扰，存在唇口激波和侧壁面形成扫略激波/边界层干扰诱发侧壁面低能扫略，肩部分离形成旋涡，沿流向形成强剪切流，即角涡。采用基于PSP技术的表面近似摩擦力线提取方法，对流场的NS方程进行简化，得到流场近似的压力梯度-摩擦力方程。
由NS方程：
在近壁面底层，NS方程简化为：
两端同时乘以摩擦力线方程τ＝μω×n得：
引入涡通量Ω＝|ω|2/2，有：
方程中第二项为相对小量，引入边界层涡通量fΩ，得到近似的压力梯度-摩擦力方程为：
其中，ρ为密度，u为速度，t为时间，为导数，/>为梯度，p为压力，μ为粘性系数，为拉普拉斯算子，θ为速度的旋度，F为单位质量彻体力，z为边界层厚度方向，τ为粘性摩擦力，ω为涡量，n为单位外法向矢量，kω为涡量系数。
该方程可通过类光流方程方法进行求解，通过压力梯度解算摩擦力矢量，当边界层涡通量未知时，能求得摩擦力矢量的相对大小和方向。
图8为表面近似摩擦力大小相对分布。
(5)搭建PIV测量系统，所述PIV测量系统包括风洞、PIV雾化示踪粒子发生器、高超声速喷管、缩比框架式试验模型、第二光源、第二图像采集设备以及同步触发器。
第二光源选用YAG双曝光激光器作为照明光源，采用调Q技术，最大脉冲能量350mJ，工作频率1～10Hz，脉冲宽度6-8ns，激光器出光孔处装有导光臂和集成片光源；第二图像采集设备主要由跨帧数字相机、图像采集板和计算机组成，跨帧数字相机图像分辨率2048×2048；同步触发器为一台延时信号发生器作为同步控制器，可同时输出6路延时信号，控制激光器和跨帧数字相机同步工作，同步控制器采用标准TTL信号格式，延时精度0.25ns，由软件控制。激光器、导光臂和集成片光源布置在激光光源布置在风洞试验段顶部，片光通过风洞试验段顶部的光学窗口和模型顶部的光学窗口，照明缩比框架式试验模型内流场对称中心截面，跨帧数字相机从风洞试验段一侧观察窗拍摄模型内流场粒子图像。选择多点标定对拍摄平面进行校准，通过校正算法修正拍摄图像的畸变，还原真实拍摄面图像；通过PIV雾化示踪粒子发生器(放在1的前端，随喷管流场一起进入风洞)降低颗粒直径提高粒子在分离区的跟随性。风洞运行时，先运行PIV图像采集系统记录图像，再开启粒子发生器，捕捉粒子从无到有全过程，确保有效图像的抓取，避免图像过曝光导致的无法计算。根据抓取的图像，利用自相关算法获取流场的速度场信息和波系结构。图5为粒子图像拍摄显示结果。
纹影试验采用双镜平行光式纹影仪(放置于窗口外侧)，以平行光入射于风洞试验段，开展纹影试验，双镜平行光式纹影仪内的单反相机捕捉缩比框架式试验模型内部的密度场图像，获得进气道内层流湍流流态变化和边界层厚度。
图6为进气道内流的波系结构及激波边界层干扰。
图4为PSP获取的表面结构和幅值与纹影对比结果。
本发明针对进气道激波边界层干扰特性，提出一种新型的流场波系结构测量方法，通过对压缩面压敏漆表面压力测量，并对获得的光场图像转换为压力值图像，通过模态分解方法实现流场能量的重构和降噪，最终获得清晰地流场激波形状。
采用压敏漆(Pressure Sensitive Paint，PSP)、粒子图像测速(Particle ImageVelocimetry，PIV)、纹影等光学非接触测试技术能够系统的开展对超燃冲压进气道内部复杂流场的波系结构、流态变化、速度场、全场动态压力变化等多物理场参数非定常特性的全方位观测和分析，弥补数值模拟和常规风洞试验的不足，揭示内流流动机理，形成模型复杂内流场可视化测试能力，为准确判别典型激波/边界层干扰分离类别、位置、干扰区形状、流向连续特性分布及三维效应提供测量参数，为型号研制过程中进气道的设计、流动控制及改进提供原理性试验依据。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (9)
1.一种进气道激波边界层干扰特性测量方法，其特征在于该方法的步骤包括：
步骤一，根据风洞试验模拟相似性准则为发动机进气道设计缩比框架式试验模型，在缩比框架式试验模型上设计光学可视化窗口，用于观测缩比框架式试验模型中的内流场；
步骤二，在缩比框架式试验模型的待测量区域喷涂压敏漆，所述压敏漆为聚合物形成的多孔性涂层结构，孔中包含有荧光探针分子；
步骤三、选择PSP测量使用的第一图像采集设备和第一光源，对压敏漆涂料的动态响应特性进行标定，确定第一图像采集设备的拍摄距离、曝光时间，第一光源的照射距离和光强度、压敏漆涂料的发光-压力关系式；
步骤四、搭建PSP测量系统，所述PSP测量系统包括风洞、高超声速喷管、缩比框架式试验模型、第一光源和第一图像采集设备；
将缩比框架式试验模型和高超声速喷管放置在风洞试验段，且缩比框架式试验模型放置于高超声速喷管的均匀区域内，第一光源和第一图像采集设备放置于风洞光学顶窗外侧，透过风洞光学顶窗，能够观察到缩比框架式试验模型的光学可视化窗口；依据步骤三确定的第一图像采集设备参数、第一光源参数对图像采集设备和光源进行设置；
步骤五、开启PSP测量系统各设备，进行PSP测量，获得进气道内表面压力分布特性、激波的大小、位置、激波边界层的流态发展和表面涡流近似摩擦力线分布特性；
步骤六、搭建PIV测量系统，所述PIV测量系统包括风洞、PIV雾化示踪粒子发生器、高超声速喷管、缩比框架式试验模型、第二光源、第二图像采集设备以及同步触发器；
将缩比框架式试验模型和高超声速喷管放置在风洞试验段，且缩比框架式试验模型放置于高超声速喷管的均匀区域内，第二光源布置在风洞试验段顶部，第二图像采集设备位于风洞试验段一侧观察窗；PIV雾化示踪粒子发生器放在高超声速喷管的前端；同步触发器用于控制第二光源和第二图像采集设备同步工作；
步骤七、先运行第二图像采集设备，拍摄缩比框架式试验模型内流场；然后开启PIV雾化示踪粒子发生器，使其发出的粒子随喷管流场一起进入风洞，第二图像采集设备捕捉粒子从无到有全过程；根据拍摄的图像，利用自相关算法获取进气道内流场的速度场信息和波系结构；
步骤八、搭建纹影测量系统，所述纹影测量系统包括风洞、高超声速喷管、缩比框架式试验模型以及双镜平行光式纹影仪；
将缩比框架式试验模型和高超声速喷管放置在风洞试验段，且缩比框架式试验模型放置于高超声速喷管的均匀区域内，双镜平行光式纹影仪放置于风洞光学顶窗外侧；
步骤九、双镜平行光式纹影仪以平行光入射于风洞试验段，双镜平行光式纹影仪内的单反相机捕捉缩比框架式试验模型内部的密度场图像，获得进气道内层流湍流流态变化和边界层厚度。
2.根据权利要求1所述的一种进气道激波边界层干扰特性测量方法，其特征在于，所述步骤一中，光学窗口对于激发光源的波长和受激辐射光波长具有良好的透过率，且光学窗口各处等厚。
3.根据权利要求1所述的一种进气道激波边界层干扰特性测量方法，其特征在于，所述步骤五中，所述进行PSP测量，获得进气道内表面压力分布特性、激波的大小、位置、激波边界层的流态发展和表面涡流近似摩擦力线分布特性的方法如下：
(5.1)开启PSP测量系统各设备，来流通过高超声速喷形成高超声速流场，连续采集N张吹风状态下的图像；
关闭高超声速喷管，采集M张无风状态下的参考图像，N>>M；
(5.2)通过图像比获取流场的壁面表面基于时间序列的压力图像序列，实现方法如下：
计算M张参考图像上同一个像素点的平均灰度值，得到一张参考图像平均值图像；
利用如下公式计算M张参考图像上像素点(xi,yi)的平均灰度值g；
Gm(xi,yi)为第m张参考图像中像素点(xi,yi)的灰度值，i∈[1,I]，I为图像中的像素点总数；
第j张吹风状态下图像上每个像素点的灰度值除以平均值图像上对应像素点的灰度值，得到第j张吹风状态下的压力图像；
j遍历[1,N]，得到N个相对于时间序列的压力图像序列；
(5.3)利用正交变换方法将压力图像序列转换为模态序列，截取前N阶模态，从前N阶模态中剔除与压力脉动无关的噪声和无效模态，将剩余与压力脉动相关的模态通过反正交变换重构得到与原来同时间分布的压力图像序列；
(5.4)根据(5.3)重构得到的压力图像序列，获得每一点的压力值，从而得到压力分布特性、激波的大小、位置、激波边界层的流态发展。
4.根据权利要求3所述的一种进气道激波边界层干扰特性测量方法，其特征在于，所述步骤(5.3)中，利用正交变换方法将压力图像序列转换为模态序列的方法如下：
寻找一组最佳的标准正交基；
使压力图像序列中每个压力图像数据在该标准正交基上的投影依次递减，每个压力图像数据在该标准正交基上的投影即为模态序列。
5.根据权利要求3所述的一种进气道激波边界层干扰特性测量方法，其特征在于，所述步骤(5.3)中，前N阶模态的能量之和超过总能量的90％。
6.根据权利要求1所述的一种进气道激波边界层干扰特性测量方法，其特征在于，所述步骤五中，所述进行PSP测量，获得表面涡流近似摩擦力线分布特性的方法如下：
采用表面近似摩擦力线提取方法，对流场的NS方程进行处理，得到压力梯度-摩擦力方程；
当边界层涡通量未知时，通过压力梯度解算压力梯度-摩擦力方程，求得摩擦力矢量的相对大小和方向。
7.根据权利要求6所述的一种进气道激波边界层干扰特性测量方法，其特征在于，对流场的NS方程进行处理，得到压力梯度-摩擦力方程的处理方式如下：
流场的NS方程如下：
在近壁面底层，u→0,z→0，NS方程简化为：
两端同时乘以摩擦力线方程τ＝μω×n得：
引入涡通量Ω＝|ω|2/2，有：
方程中第二项为相对小量，引入边界层涡通量fΩ，得到压力梯度-摩擦力方程为：
其中，ρ为密度，u为速度，t为时间，为导数，/>为梯度，p为压力，μ为粘性系数，/>为拉普拉斯算子，θ为速度的旋度，F为单位质量彻体力，z为边界层厚度方向，τ为粘性摩擦力，ω为涡量，n为单位外法向矢量，kω为涡量系数。
8.根据权利要求1所述的一种进气道激波边界层干扰特性测量方法，其特征在于，第一光源为波长400nm的紫外光源，第一图像采集设备为普通相机。
9.根据权利要求1所述的一种进气道激波边界层干扰特性测量方法，其特征在于，第二光源为YAG双曝光激光器，激光器出光孔处装有导光臂和集成片光源；第二图像采集设备为跨帧数字相机。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110738673.7A CN113588204B (zh) | 2021-06-30 | 2021-06-30 | 一种进气道激波边界层干扰特性测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110738673.7A CN113588204B (zh) | 2021-06-30 | 2021-06-30 | 一种进气道激波边界层干扰特性测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113588204A CN113588204A (zh) | 2021-11-02 |
| CN113588204B true CN113588204B (zh) | 2023-12-12 |
Family
ID=78245287
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110738673.7A Active CN113588204B (zh) | 2021-06-30 | 2021-06-30 | 一种进气道激波边界层干扰特性测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113588204B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114119681A (zh) * | 2021-11-12 | 2022-03-01 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种空腔模型脉动压力与噪声测量方法、装置及电子设备 |
| CN113884267B (zh) * | 2021-12-07 | 2022-02-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于脉冲风洞的瞬态喷流试验装置 |
| CN113945354B (zh) * | 2021-12-14 | 2023-04-07 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于辨识膨胀风洞加速段流动分区特性的试验方法 |
| CN114838906B (zh) * | 2022-03-24 | 2024-07-19 | 中国人民解放军国防科技大学 | 对风洞来流能谱进行测量的方法及系统、设备、存储介质 |
| CN115306559B (zh) * | 2022-06-30 | 2024-07-30 | 哈尔滨工业大学 | 一种抑制隔离段激波串突跳的分布式自适应流动控制方法及其控制装置 |
| CN116341106B (zh) * | 2023-03-14 | 2024-06-07 | 南京航空航天大学 | 一种基于流场相似变换的强展向压力梯度压缩面设计方法 |
| CN116202791B (zh) * | 2023-04-27 | 2023-07-25 | 中南大学 | 一种使用低频piv测试高频扫掠射流瞬时特性的方法及装置 |
| CN116735069B (zh) * | 2023-08-15 | 2023-10-31 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种新型流场压力光学测量方法 |
| CN117232395B (zh) * | 2023-11-14 | 2024-01-23 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种压敏漆图像激波位置自动识别方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004156968A (ja) * | 2002-11-05 | 2004-06-03 | Railway Technical Res Inst | 空力音源探査システム及び空力音源探査方法 |
| RU67258U1 (ru) * | 2006-12-28 | 2007-10-10 | Закрытое акционерное общество "Гражданские самолеты Сухого" | Установка для аэродинамических испытаний модели воздухозаборника двигателя летательного аппарата (варианты) |
| JP2013079878A (ja) * | 2011-10-04 | 2013-05-02 | Japan Aerospace Exploration Agency | 流体測定方法および装置 |
| CN103149010A (zh) * | 2013-02-22 | 2013-06-12 | 中国人民解放军国防科学技术大学 | 基于高超声速脉冲风洞的进气道流场npls测量系统及方法 |
| CN106596037A (zh) * | 2016-12-16 | 2017-04-26 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞试验模型绕流密度投影场的视频测量方法 |
| CN108120583A (zh) * | 2018-01-20 | 2018-06-05 | 南京航空航天大学 | 用于测量高超声速风洞试验动态压力和空间流场的装置 |
| CN110057535A (zh) * | 2019-05-27 | 2019-07-26 | 中国空气动力研究与发展中心 | 激波风洞判断进气道起动的壁面丝线流动显示方法及系统 |
| CN110333044A (zh) * | 2019-05-27 | 2019-10-15 | 中国空气动力研究与发展中心 | 激波风洞进气道自起动试验出口气动旋转门装置及方法 |
| CN111792061A (zh) * | 2020-06-30 | 2020-10-20 | 中国航天空气动力技术研究院 | 一种利用激波边界层干扰的气动热试验装置及方法 |
| CN112945501A (zh) * | 2021-02-03 | 2021-06-11 | 中国空气动力研究与发展中心高速空气动力研究所 | 层流机翼转捩位置测量试验方法 |
-
2021
- 2021-06-30 CN CN202110738673.7A patent/CN113588204B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004156968A (ja) * | 2002-11-05 | 2004-06-03 | Railway Technical Res Inst | 空力音源探査システム及び空力音源探査方法 |
| RU67258U1 (ru) * | 2006-12-28 | 2007-10-10 | Закрытое акционерное общество "Гражданские самолеты Сухого" | Установка для аэродинамических испытаний модели воздухозаборника двигателя летательного аппарата (варианты) |
| JP2013079878A (ja) * | 2011-10-04 | 2013-05-02 | Japan Aerospace Exploration Agency | 流体測定方法および装置 |
| CN103149010A (zh) * | 2013-02-22 | 2013-06-12 | 中国人民解放军国防科学技术大学 | 基于高超声速脉冲风洞的进气道流场npls测量系统及方法 |
| CN106596037A (zh) * | 2016-12-16 | 2017-04-26 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞试验模型绕流密度投影场的视频测量方法 |
| CN108120583A (zh) * | 2018-01-20 | 2018-06-05 | 南京航空航天大学 | 用于测量高超声速风洞试验动态压力和空间流场的装置 |
| CN110057535A (zh) * | 2019-05-27 | 2019-07-26 | 中国空气动力研究与发展中心 | 激波风洞判断进气道起动的壁面丝线流动显示方法及系统 |
| CN110333044A (zh) * | 2019-05-27 | 2019-10-15 | 中国空气动力研究与发展中心 | 激波风洞进气道自起动试验出口气动旋转门装置及方法 |
| CN111792061A (zh) * | 2020-06-30 | 2020-10-20 | 中国航天空气动力技术研究院 | 一种利用激波边界层干扰的气动热试验装置及方法 |
| CN112945501A (zh) * | 2021-02-03 | 2021-06-11 | 中国空气动力研究与发展中心高速空气动力研究所 | 层流机翼转捩位置测量试验方法 |
Non-Patent Citations (1)
| Title |
|---|
| 超燃进气道激波/湍流边界层干扰;潘宏禄;李俊红;沈清;;推进技术(09);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113588204A (zh) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113588204B (zh) | 一种进气道激波边界层干扰特性测量方法 | |
| CN103698554B (zh) | 一种流场实时精确测量系统及方法 | |
| CN113155399B (zh) | 高速飞行器表面压力与变形三维连续分布同步测量方法 | |
| CN110441021A (zh) | 一种基于piv的平板湍流边界层流场结构测量装置及方法 | |
| CN113324727A (zh) | 一种针对压缩拐角超声速流场结构的纹影图像处理方法 | |
| Duncan et al. | Effects of step excrescences on a swept wing in a low-disturbance wind tunnel | |
| CN109870288B (zh) | 一种基于激光诱发粒子技术的空气流动显示方法 | |
| Balamurugan et al. | Experiments on localized secondary instability in bypass boundary layer transition | |
| Hsu et al. | Turbulent kinetic energy transport in a corner formed by a solid wall and a free surface | |
| CN108256166A (zh) | 一种用于热图技术的数据处理方法 | |
| CN114459727A (zh) | 用于swtbli非定常特性研究的实验系统及方法 | |
| Simmons et al. | Benchmark smooth body flow separation experiments | |
| CN113390602B (zh) | 一种基于内嵌式温敏漆的气动热全场测量系统及方法 | |
| Winski et al. | Computational and experimental study of plume and shock interaction effects on sonic boom in the NASA Ames 9x7 supersonic wind tunnel | |
| Neel et al. | Characterization of environmental disturbances on hypersonic crossflow instability on the HIFiRE-5 Elliptic Cone | |
| Pehrson et al. | Wall-normal fleet velocimetry in a canonical hypersonic inlet | |
| CN112903236A (zh) | 一种基于焦平面扫描的气动热力参数光学测试装置及方法 | |
| Weisberger et al. | Self-Aligned Focusing Schlieren at the 0.3-M Transonic Cryogenic Tunnel and the National Transonic Facility | |
| Lorenzoni et al. | Aeroacoustic analysis of a rod-airfoil flow by means of time-resolved PIV | |
| Lu et al. | Experimental study on time evolution of shock wave and turbulent boundary layer interactions | |
| Reising et al. | Development and Assessment of a New Particle Image Velocimetry System in the NASA GRC 225 cm2 Wind Tunnel | |
| Cramer et al. | Boiler tube corrosion characterization with a scanning thermal line | |
| Eppink | Effects of random micron-sized roughness on swept-wing transition | |
| Zoppini et al. | Experimental investigation on receptivity of crossflow instability to discrete roughness amplitude and location | |
| CN114689898A (zh) | 一种观测飞秒激光成丝冲击云雾滴的装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |