小型聚焦纹影实验台的设计与搭建_王定奇.pdf

书书书测控技术2023 年第 42 卷第 2 期试验与测试收稿日期:2022 04 17引用格式:王定奇 小型聚焦纹影实验台的设计与搭建 J 测控技术，2023，42(2):62 67WANG D Q Design and Construction of Small Focusing Schlieren Experiment BenchJ Measurement Control Technology，2023，42(2):62 67小型聚焦纹影实验台的设计与搭建王定奇(中国飞行试验研究院 发动机所，陕西 西安710089)摘要:聚焦纹影系统能对垂直于光轴的特定平面内的流场细节进行显示，避免了视场范围内其余非均匀流场对背景成像的干扰。为实现对收缩喷管出口超声速流场的观测，设计并搭建了小型聚焦纹影实验台和超声速流动实验台。超声速流动实验台由 2 组收缩喷管组成，相距 15 cm，垂直摆放，利用高压气源供气，使喷管出口继续膨胀产生超声速气流。首次利用 CAD 设计胶片上格栅间距尺寸，验证了在胶片纸上进行源光栅、刀口栅改造的正确性。综合考虑超声速流动区域的大小，选取了锐利聚焦深度和非锐利聚焦深度等设计指标。完成了对小型超声速流动实验台流场的拍摄，获得了其流场结构图。通过对比聚焦纹影、计算流体力学(CFD)仿真和传统纹影的流场图，说明此聚焦纹影实验台原理正确，可以达到预期的显示效果。同时，针对聚焦纹影系统高速流场的显示，提出合理的建议。关键词:聚焦纹影;数值计算;光栅;菲涅耳光源;流场显示中图分类号:V211 72文献标志码:A文章编号:1000 8829(2023)02 0062 06doi:10 19708/j ckjs 2023 02 011Design and Construction of Small Focusing Schlieren Experiment BenchWANG Ding-Qi(Engine Department，Chinese Flight Test Establishment，Xi an 710089，China)Abstract:Focusing schlieren system can display the details of the specific plane perpendicular to the optic axisin the flow field，so the interference of other non-uniform flow fields in the field of view on the background ima-ging is avoided In order to realize the observation of the supersonic flow field at the exit of the shrinking noz-zle，a small focusing schlieren experiment bench and a supersonic flow experiment bench are designed andbuilt The supersonic flow experiment bench is composed of two sets of shrinking nozzles，15 cm apart andplaced vertically The high-pressure air source is used to supply air to make the nozzle outlet continue to ex-pand to generate supersonic airflow It is the first time to use CAD to design the grid spacing size on the film，which verifies the correctness of the source grating and knife-edge grid transformation on the film paper Con-sidering the size of the supersonic flow area，the design indicators such as sharp focus depth and non-sharp fo-cus depth are selected The filming of the flow field of the small supersonic flow experiment bench is comple-ted，and the flow field structure diagram is obtained By comparing the flow field diagram of focusing schlierenand computational fluid dynamics(CFD)simulation and traditional schlieren，it shows that the principle of fo-cusing schlieren experiment bench is correct，and the expected display effect can be achieved At the sametime，reasonable suggestions are put forward for the display of the high-speed flow field of the focusing schlierensystemKey words:focusing schlieren;numerical calculation;grating;Fresnel lens;flow visualization26纹影法是一种通过光线的偏折情况来显示流场的非接触方法1 2，被广泛应用在空气动力学、燃烧学等领域的研究工作中。随着流体力学的发展，对流动分离、漩涡和激波干扰等复杂流动现象的研究和探索变得更加迫切，流动显示技术在帮助人们认识和了解这些复杂流动现象方面具有明显的优点，具有其他方法无法比拟的作用。传统的流场测量需要借助于压力和温度耙来测量流场参数，但是将受感部件深入流场，在获得参数的同时也对流场产生了扰动。尤其对于超声速流动，若直接采用测量耙，其前方会产生正激波，无法得到流场真实的信息，因而实验中使用传统的探针测量流场很难满足多样化的测试需求3。纹影技术的发展不仅带动了流体力学的研究，而且也推动了其工程应用。现代光学和计算机图像处理技术的迅速发展，大大提高了流动显示技术的使用程度。传统纹影方法沿着光路的流场积分获取整个流场信息，因此无法满足对特定区域流场结构的观测，而基于传统纹影发展而来的聚焦纹影系统，通过聚焦透镜，可以针对特定平面进行聚焦，在像平面上得到聚焦区域内的流场信息，而将非聚焦区域流场信息当作背景信息模糊处理，从而使纹影系统能够准确观测流动细节4 6。20 世纪 90 年代开始，ice 等7 开始利用聚焦纹影技术研究啸音激励对超声速流场影响;Cook 等8 通过聚焦纹影图像，提取出流场的密度信息，并开展了数值仿真对比研究;Hargather 等9 研究了聚焦纹影技术与粒子图像测速(Particle Image Velocimetry，PIV)相结合的湍流边界层速度分布测量可行性，当聚焦景深在40 mm 时，能看出涡密度图像。近年来，我国相关学者在流动显示技术方面取得了显著成绩，吴文堂等10 针对高超声速流场的波系结构采用彩色滤光片和高速摄影技术，研究了斜激波相互作用产生的正规反射和马赫反射激波结构，得到了彩色纹影照片;中国空气动力研究与发展中心将聚焦纹影显示技术和图像处理技术应用在激波风洞试验中，对流场密度进行了定量测量，证明了纹影法的可行性11。谢爱民等12 提出了激波风洞流场密度测量的聚焦纹影技术和图像密度场处理技术，并将试验测量的密度值与数值计算结果进行对比，试验获取的密度变化规律与数值模拟吻合度较好。本文通过搭建小型聚焦纹影实验台对超声速流动实验台产生的高速气流的波系结构进行观测，并将传统纹影、聚焦纹影和数值仿真获取的流场图像进行对比分析。1纹影显示原理光线的传播速度为v=c0n(1)式中:c0为光速;n 光的折射率。由式(1)可以得到，光在密度大的媒介中传播速度慢，在密度小的介质中速度较快。光通过非均匀的流场时，方向发生改变，如图 1 所示。图 1光线在非均匀流场的偏折可推导出折射角为d=drd(2)dr=dc(3)=d/c(4)联立式(2)式(4)得到:d=dccdd(5)式中:d 为光线在 时间内走过的距离;dr 为 r1比 r2多走的距离;dc 为 r1与 r2的速度差。由式(2)可得到 dc/c=dn/n，将此公式带入可得:d=1ndndd(6)由式(6)可以看出，光线在单位长度上的偏折角正比于其法线方向的折射率梯度，偏折方向在密度增加的梯度方向。光线经过整个扰动区后，偏折角应为式(6)沿光路上的积分。=211ndndd(7)式中:1和 2为扰动区沿光路的边界。在多数情况下 的绝对值很小，光线在偏离扰动区时 x、y 坐标的变化可以忽略不计，仅有一个小偏角，所以光路积分可用沿 z 轴的积分公式代替。=z2z11ndnddz(8)实际使用中仪器检测的偏折方向不一定与 重合，光线在 x、y 方向的偏折角为=z2z11ndnddz(9)式中:z1和 z2分别为扰动区在 z 方向的边界。如果在扰动区垂直方向放置一个像平面，光线穿过扰动区后在屏上移动距离在 x、y 方向的投影为36小型聚焦纹影实验台的设计与搭建x=l211nnxdy=l211nnyd(10)式中:l 为像平面到扰动区中心距离。传统纹影中光线偏折在整个成像平面上所引起的光强相对变化13 15 为I/I0=Kfaz2z1xdz(11)式中:I=I I0;f 为聚焦透镜的焦距;K 为格拉斯通 戴尔常数;为气流密度;a 为未被刀口切割的光源宽度。而聚焦纹影光强相对的变化为I/I0=Kfa(L f)z2z1x(L f)(z)dz(12)式中:L 为源光栅与凸透镜距离;(z)为传递函数16，如图 2 所示。该传递函数可简化为用一系列厚度为“d”的薄片来代替曲线的积分函数。I/I0=Kfa(L f)dj(L lj)(x)dz(13)式中:lj为凸透镜与成像平面距离;dj为不同流场区域厚度。图 2纹影传递函数图2纹影系统总体结构2 1传统纹影系统结构传统的纹影系统由 1 个点光源和 2 个相同焦距的凸透镜、刀口、像平面和非均匀流场区域组成。当纹影系统的测试区域流场均匀时，位于焦点上的点光源发出光线穿过凸透镜后形成一束平行光，平行光线在穿过第 2 个透镜时又汇聚到一点，在像平面上形成倒立的像。若在第 2 个透镜的焦点处用刀口遮挡部分光线，当光线穿过实验区域的非均匀流场，由于密度梯度的不同，使得平行光束中的一部分光线发生偏折，在刀口平面上，光源狭缝像相对于刀口有一个位移量。刀口处的光线一部分被遮挡，一部分透过，从而在背景中形成暗场和亮场，像平面上的亮暗区域与实验区域的密度梯度对应起来。传统纹影系统示意图如图 3 所示，光通过非均匀流场发生偏折，折射线(由蓝色实线表示)在 2 个透镜间的密度梯度将会在平面上的镜头上显示出来。图 3传统纹影系统示意图2 2聚焦纹影系统结构聚焦纹影系统由一个扩展光源、菲涅耳透镜、源光栅、透镜、刀口栅和像平面组成。扩展光源发出的光通过菲涅耳透镜产生聚焦效果并显著增强光源的亮度。光线到达源光栅时，明暗相间的条纹使光线一部分通过栅格，形成多狭缝光源。聚焦透镜的作用是将非均匀流场中光线偏折反映在背景图像中;源光栅和刀口栅都是由明暗交替的条纹组成的，在光路上形成共轭像17。在共轭像的位置放置刀口栅，充当传统纹影中的刀口。聚焦纹影系统能把观测区域聚焦在狭窄的平面附近，使得实验区域内的三维流动很容易观察，并消除观察域外的流场扰动18。在流场显示区域有扰动引起流场折射率的变化，光线一部分被刀口栅遮挡，观测区域内密度的变化同样会在像平面中显示出来，如图 4 所示。图 4聚焦纹影系统示意图2 3小型聚焦纹影系统搭建小型聚焦纹影实验台主要由 LED 光源(如图 5 所示)、菲涅耳透镜、格栅、聚焦透镜和成像原件(CCD 相机)组成。为了使光源发出的光线在穿过非均匀流场后仍然能得到明亮的拍摄视野，选用功率为 30 W 的 LED灯。菲涅耳透镜多由聚乙烯材料制成，镜片表面一面为光面，另一面刻录了由小到大的同心圆，其面积大，厚度为 2 mm。菲涅耳透镜可以显著提高光源的亮度，实验中根据需求选用了焦距为 220 mm、外形尺寸为210 mm 130 mm 的菲涅耳透镜，如图 6 所示。46测控技术 2023 年第 42 卷第 2 期图 5LED 光源图 6菲涅耳透镜光栅是聚焦纹影系统的关键核心部件，源光栅和刀口栅由一系列明暗相间的条纹组成，光栅应在亮条纹中有较好的透光率，而在暗条纹中阻挡光线的透过。采用胶片打印技术，用 CAD 调试好格栅的间距尺寸后再打印到胶片上，如图 7 所示。实验中为了增大调节的精确度，将源光栅和刀口栅的条纹沿竖直方向摆放，这样左右移动格栅调节刀口切割量时更加方便。聚焦透镜除满足普通透镜的成像外，还能消除色差和球差成像，因此采用了光学性能较好的 K9 平凸透镜，使光束聚焦和准直19，如图 8 所示。根据透镜成像公式可以确定各光学元件的安放位置和测试区域的聚焦位置。根据测试流场和窗口的距离，同时考虑系统的景深影响，设计聚焦透镜的焦距为 210 mm、直径为 50 mm。聚焦纹影系统参数如表 1 所示。图 7源光栅和刀口栅图 8聚焦透镜表 1聚焦纹影系统参数参数名称选取值参数名称选取值源光栅亮纹(a)1mm凸透镜直径(D)50mm源光栅暗纹(b)1mm混合线取值(n)0 459刀口栅亮纹1425mm灵敏度(w)245mm刀口栅暗纹0726mm放大倍率(m)273光线波长()532nm锐利景深(DS)203mm菲涅耳透镜焦距(F)220mm混合线()4 641凸透镜焦距(f)210mm3超声速流动实验台的设计超声速流动实验台由空气压缩机、供气管路、喷管和固定支架等结构组成，如图 9 所示。空气压缩机为喷管提供高压气体，气体压力为 1 2 MPa，其额定功率11 kW，转速 1450 r/min。压缩机出口接高压软管，软管通过三通阀分出两路气管，分别引入垂直喷管和水平喷管。软管的内直径 5 mm，喷管出口直径 2 mm。图 9超声速流动实验台示意图喷管采用收缩喷管，当喷管出口静压与总压之比小于临界值时，在收缩喷管出口处可以达到声速，并且出口处的压力远高于环境压力，所以气流可以在管外继续膨胀达到超声速流。喷管的加工采用 3D 打印。考虑到实验件的安全性，采用能承受较大压力的树脂材料，另一方面加厚喷管的壁面。高压软管一端接通压气机，另一端连接喷管;由于软管的内表面光滑，沿程损失系数很小，可以忽略。高压软管与喷管连接处采用黏性好、强度高的704 硅胶填充。用胶带将喷管固定到正方体的相邻两截面上，喷管出口截面相互垂直，且轴线相距 150 mm。实验台支架可以在垂直和水平方向进行高度和偏转角度的调节。表 2超声速流动实验台参数参数名称选取值参数名称选取值喷管进口总压(P*)10atm三通阀流量(Qv1)0 025m3/s喷管进口总温(T*)350K喉道截面压力(P7)6 8atm环境压力(P0)1atm流量效率()02进口截面(A1)01963cm3出口半径(r)1mm体积流量(Qv)0005m3/s4结果分析4 1聚焦纹影观测结果确定各实验器材和参数后，按照设计依次排列，利用水平仪调节各部件校准光轴。首先将水平仪摆放好，确定系统的光轴，再依次摆放 LED 光源、菲涅耳透镜、源光栅、小型超声速流动验证实验台、聚焦透镜、刀口和单反相机，并按照设计的尺寸调整部件之间的位置关系，如图 10 所示。当小型超声速流动验证实验台中竖直流场处于聚焦纹影系统的锐利景深厚度范围内，背景图像中只显示了竖直方向流场，并将外界流场虚化，如图 11 所示。调整小型超声速流动验证实验台的位置，使横向流场落在聚焦纹影的锐利景深厚度范围内，相机拍摄的流场图像只能显示横向气流的流场，如图 12 所示。喷管出口可以看到清晰的马赫盘;通过分别显示横向和纵向气流，可以看出聚焦纹影系统的景深在较小的厚度范围内。56小型聚焦纹影实验台的设计与搭建图 10聚焦纹影实验台图 11竖直聚焦纹影流场图图 12横向流场聚焦纹影图4 2仿真流场对比验证使用 UG 建立喷管模型，ICEM 生成网格，Fluent开展计算，Tecplot 进行后处理。计算网格总量为 15万，湍流模型采用 k-omega SST，喷管进口设置为压力入口，喷管内表面设置为无滑移壁面，模型网格如图13 所示。进口压力为 911 kPa，总温 300 K，喷管出口环境压力为 101 kPa。图 13模型网格喷管出口流场的截面马赫数云图，如图 14(a)所示。颜色的深度与马赫数相关，随着颜色的加深马赫数变大，可以看出喷管出口的马赫数最高可达 2 9Ma，满足了风洞试验中所需要的超声速气流。与图 14(b)所示的聚焦纹影拍摄的图像进行对比，可以看出两流场都在出口外部形成了马赫盘，流场分布相似，较好地模拟了实验流场。4 3传统纹影与聚焦纹影对比Z 型结构纹影仪的布置中，2 个球面反射镜之间为平行光线，使得球面反射镜之间的流动都可以反映在背景图像中。光源发出扩展光经聚光镜后照射到球图 14流场对比面反射镜上，调整好其与光线的位置，使得经过球面反射镜后为平行光，在第二个球面反射镜后调整刀口的位置，使刀口恰好落在焦点处。Z 型结构纹影仪在调节时，尽可能使光源狭缝像与刀口在同一平面内，且刀口的刃边与狭缝平行，其方向与被观察区域的折射率梯度方向垂直。同时，刀口必须切割在透镜的焦点上，且当切割量为 50%时效果最明显20 21。由于气体密度不同导致光线的折射率不同，使偏折角的大小可在灰度背景图中反映出来，故便于观察。Z 型结构纹影光路图如图 15 所示。图 15Z 型结构纹影光路图传统纹影和聚焦纹影图像喷管超声速气流对比如图 16 所示。传统纹影在成像时，可以显示 2 个方向的超声速流场。以横向的喷管出口喷流为实验观察对象时，竖直方向上的气流为扰动源。在观察横向气流时，竖直方向上的气流也会清晰地显示在背景图像中，这说明传统的纹影系统的检测区域为 2 个凹透镜之间的任何区域，当实验段为其中的一小块区域时，会受到外界流动的干扰并在背景图像上反映出来，给图像分析带来了较大的干扰。而聚焦纹影只把需要观察的流动区域显示出来，将垂直方向气流的扰动虚化，将流动的细节更充分地显示出来。5结束语通过自行搭建的聚焦纹影实验台，开展了超声速气流的流场显示，通过聚焦纹影可以聚焦某个较窄视场范围的流场，并得出以下结论。66测控技术 2023 年第 42 卷第 2 期图 16传统纹影和聚焦纹影喷管超声速气流对比 聚焦纹影系统能够实现特定区域非均匀流场的显示，并消除观测区域外的非均匀流场，提高了对流场流动细节的观测能力。在设计小型超声速流动验证实验台的基础上，设计并搭建了满足小型超声速流动验证实验台的聚焦纹影系统，能够对超声速流场局部进行显示。利用 CAD 调试好格栅的间距尺寸后再打印到胶片上，验证了在胶片纸上进行源光栅和刀口栅改造的正确性，且缩短了实验周期，降低了成本。通过对比聚焦纹影结果与 CFD 计算和传统纹影结果，证明了此系统原理正确，能实现预期显示效果。本项目在聚焦纹影实验台搭建过程中，选用低成本的胶片纸制作源光栅和刀口栅，材料未能完全满足光学实验所需的精度，导致聚焦纹影景深范围内的流场显示不够清楚，在背景图像中出现了光源的干扰。后续可以在前期实验的基础上做以下改进。将光栅与桌面的垂直度精细调节，保证光栅切割的一致性;在刀口光栅平移过程中，采用机械结构助力调节，定量移动刀口栅。改进光栅的材质，强化光栅刀口对光线的切割量，避免黑色条纹透过光线对背景图像产生干扰。增大菲涅耳光源和透镜的尺寸，更好地显示小景深大视场的非均匀流场。参考文献:1WEINSTEIN L M Large-field high-brightness focusing schlier-en system J AIAA Journal，1993，31(7):1250 1255 2SETTLES G S Shadowgraph techniquesM/Schlierenand shadowgraph techniques Berlin，Heidelberg:Springer，2001:143 163 3徐翔，谢爱民，吕治国，等 聚焦纹影显示技术在激波风洞的初步应用 J 实验流体力学，2009，23(3):75 79XU X，XIE A M，LYU Z G，et al Application of focusingschlieren visualization system in shock tunnel experiment J Journal of Experiments in Fluid Mechanics，2009，23(3):75 79(in Chinese)4中国人民解放军总装备部军事训练教材编辑工作委员会 流动显示技术 M 北京:国防工业出版社，2002Editorial Committee of Military Training Textbooks of theGeneral Equipment Department of the Chinese PeoplesLiberation Army Flow display technology M Beijing:Na-tional Defense Industry Press，2002(in Chinese)5李桂春 风洞试验光学测量方法 M 北京:国防工业出版社，2008LI G C Optical instrumentation for wind tunnel testing M Beijing:National Defense Industry Press，2008(inChinese)6何霖，易仕和，赵玉新 基于背景成像纹影技术在流场密度场测量中的应用 C/第十二届全国激波与激波管学术会议论文集 2006:286 289HE L，YI S H，ZHAO Y X Application of background ima-ging schlieren technology in flow field density field meas-urementC/Proceedings of the 12th National AcademicConference on Shock Waves and Shock Tubes 2006:286 289(in Chinese)7ICE E J，AMAN G Enhanced mixing of a rectangularsupersonic jet by natural and induced screechC/AIAAShear Flow Conference 1993 8COOK S，CHOKANI N Quantitative results from the focu-sing schlieren techniqueC/31st Aerospace SciencesMeeting 1993:630 9HAGATHE M J，SETTLES G S Natural-background-ori-ented schlieren imagingJ Experiments in Fluids，2010，48(1):59 68 10 吴文堂，洪延姬，王殿恺，等 斜激波与弓形激波相互作用的彩虹纹影实验研究J 推进技术，2015，36(5):671 677WU W T，HONG Y J，WANG D K，et al Experimental in-vestigation on oblique shock and bow shock interaction withrainbow schlierenJ Journal of Propulsion Technology，2015，36(5):671 677(in Chinese)11 黄训铭，谢爱民，郑蕾，等 聚集纹影图像密度场处理技术 J 四川兵工学报，2015，36(6):77 81HUANG X M，XIE A M，ZHENG L，et al Technique re-search of focusing schlieren images processing for gettingflow densityJ Journal of Sichuan Ordnance，2015，36(6):77 81(in Chinese)12 谢爱民，黄洁，徐翔，等 激波风洞流场密度测量的聚焦纹影技术 J 实验流体力学，2013，27(2):82 86XIE A M，HUANG J，XU X，et al Focusing schlieren tech-nique applied to measure the flow density in shock tunnel J Journal of Experiments in Fluid Mechanics，2013，27(2):82 86(in Chinese)13 MIANDA D A，CISTIANO K L，TIANA D A De-sign，assembly and calibration of an optoelectronic devicethat measures three wavelengths in the near infrared forbreast cancer early detectionC/IEEE 5th ColombianWorkshop on Circuits and Systems 2014:1 4(下转第 93 页)76小型聚焦纹影实验台的设计与搭建frequency underwater acoustic signal using a improved em-pirical mode decompositionJ Journal of Applied Acous-tics，2016，35(4):316 323(in Chinese)9赵龙，高井祥，李增科，等 基于 CEEMD-Wavelet-SavGol模型的 GPS 多路径修正算法 J 测绘通报，2017(11):1 5ZHAO L，GAO J X，LI Z K，et al GPS multipath correctionalgorithm based on CEEMD-wavelet-SavGol model J Bul-letin of Surveying and Mapping，2017(11):1 5(in Chi-nese)10 刘霞，李文 基于变分模态分解的相关能量熵自适应阈值去噪 J 吉林大学学报(信息科学版)，2019，37(5):518 525LIU X，LI W Adaptive threshold denoising of relevant ener-gy entropy based on variational modal decompositionJ Journal of Jilin University(Information Science Edition)，2019，37(5):518 525(in Chinese)11 尹丽，陈富民，张琦，等 采用集合经验模态分解和改进阈值函数的心电自适应去噪方法J 西安交通大学学报，2020，54(1):101 107YIN L，CHEN F M，ZHANG Q，et al ECG adaptive denois-ing method based on EEMD and improved threshold func-tion J Journal of Xian Jiaotong University，2020，54(1):101 107(in Chinese)12 周慧 基于经验模态分解和小波分析的靶船外弹道测量数据去噪方法 J 测控技术，2012，31(7):31 34ZHOU H Target ship exterior ballistic data denoising meth-od based on empinical mode decomposition and wavelettransformationJ Measurement Control Technology，2012，31(7):31 34(in Chinese)13 苏建芳，吴钦木 基于小波包分析的电机滚动轴承故障诊断 J 测控技术，2019，38(4):64 67SU J F，WU Q M Fault diagnosis of motor rolling bearingbased on wavelet packet analysis J Measurement Con-trol Technology，2019，38(4):64 67(in Chinese)14 BEKIOS S，NGUYEN D K，SANDOVAL L J，et al Infor-mation diffusion，cluster formation and entropy-based net-work dynamics in equity and commodity markets J Euro-pean Journal of Operational esearch，2017，256(3):945 961 15 何林泽，韩先平 高速弹丸轨迹精度影响因素研究J 测控技术，2020，39(11):7 12HE L Z，HAN X P Influence factors of trajectory accuracyof high-speed projectile J Measurement Control Tech-nology，2020，39(11):7 12(in Chinese)16 倪淑燕，程乃平，崔建华，等 基于 GNSS 的飞行器测姿系统设计与实现 J 测控技术，2020，39(1):71 75NI S Y，CHENG N P，CUI J H，et al Design and implemen-tation of aircraft attitude measurement system based onGNSSJ Measurement Control Technology，2020，39(1):71 75(in Chinese)(上接第 67 页)14 岳茂雄，王如琴，姚向红，等 高速聚焦纹影改进及应用 J 实验流体力学，2013，27(5):88 93YUE M X，WANG Q，YAO X H，et al Improved high-speed focusing schlieren system and its applicationJ Journal of Experiments in Fluid Mechanics，2013，27(5):88 93(in Chinese)15KASHITANI M，YAMAGUCHI Y，KITANO H PreliminaryStudy for quantitative measurement of flow fields by focu-sing Schlieren methodC/49th AIAA Aerospace Sci-ences Meeting including the New Horizons Forum and Aer-ospace Exposition 2011 16 李炜龙 纹影仪流场密度定量测量技术研究 D 西安:西安工业大学，2017LI W L Study on quantitative schlieren techniques for den-sity measurement of flow field D Xian:Xian Technologi-cal University，2017(in Chinese)17 张敏莉，易仕和，赵玉新 超声速短化喷管的设计与实验研究 J 空气动力学学报，2007，25(4):500 503ZHANG M L，YI S H，ZHAO Y X The design and experi-mental investigations of supersonic length-shorted nozzle J Acta Aerodynamica Sinica，2007，25(4):500 503(in Chinese)18 乐嘉陵 吸气式高超声速技术研究进展J 推进技术，2010，31(6):641 649LE J L Progress in air-breathing hypersonic technology J Journal of Propulsion Technology，2010，31(6):641 649(in Chinese)19 俞刚，范学军 超声速燃烧与高超声速推进 J 力学进展，2013，43(5):449 472YU G，FAN X J Supersonic combustion and hypersonicpropulsionJ Advances in Mechanics，2013，43(5):449 472(in Chinese)20 席锋，秦华锋，胡莉，等 可视化聚焦纹影系统实现流场中的三维扰动成像 J 光学学报，2021，41(9):90 95XI F，QIN H F，HU L，et al Three-dimensional perturbationimaging in fluid field with visualized focusing schlieren sys-tem J Acta Optica Sinica，2021，41(9):90 95(in Chi-nese)21BATHEL B F，WEISBEGE J M Compact，self-alignedfocusing schlieren system J Optics Letters，2021，46(14):3328 333139基于本征模函数能量的外测数据误差分析方法