基于GCV-FFT方法的超声速压缩拐角流场气动光学效应计算

doi:10.19596/j.cnki.1001-246x.8065

计算物理 ›› 2020, Vol. 37 ›› Issue (3): 284-298.DOI: 10.19596/j.cnki.1001-246x.8065

基于GCV-FFT方法的超声速压缩拐角流场气动光学效应计算

周志超, 许凌飞, 任天荣, 顾村锋

上海机电工程研究所, 上海 201109

收稿日期:2019-03-26 修回日期:2019-05-27 出版日期:2020-05-25 发布日期:2020-05-25
通讯作者: 顾村锋,E-mail:13816429275@163.com
作者简介:周志超(1984-),男,高级工程师,研究方向:空气动力学、气动光学,E-mail:send2zhou@sina.com
基金资助:
上海市自然科学基金（16ZR1415900）资助项目

Aero-optical Effects in Supersonic Compression Ramp Flow Fields:GCV-FFT Method

ZHOU Zhichao, XU Lingfei, REN Tianrong, GU Cunfeng

Shanghai Electro-Mechanical Engineering Institute, Shanghai 201109, China

Received:2019-03-26 Revised:2019-05-27 Online:2020-05-25 Published:2020-05-25

摘要/Abstract

摘要： 研究超声速压缩拐角流动，基于Rytov近似，采用广义卷积-快速Fourier变换（GCV-FFT），直接求解非均匀介质中的Helmholtz方程，获得激光光束的散射场和衍射场的复振幅分布.给出包括光强、光程差、斯特列尔比、相对光强、光束质心位置和有效光束宽度等气动光学参量.可以看出，激光光束在湍流场区域内传输，光强度分布逐渐出现起伏和偏离.穿过流场区域后，由于衍射作用，光束产生较大幅度的扭曲、偏离和破碎.压缩拐角区的流场对光束质量的影响明显大于充分发展湍流区.从相对光强的分布形态来看，光束的破碎似乎与流场的小尺度相干结构相关.采用传统的OPD均方差或加权均方差值估算获得的斯特列尔比的误差较大，更为准确的方法应该通过计算PSF来进行评估.

关键词: 气动光学, 超声速压缩拐角, 湍流边界层, Helmholtz方程, GCV-FFT方法

Abstract: Supersonic compression ramp flow is studied. Based on Rytov approximation, a generalized convolution-fast Fourier transform (GCV-FFT) method is used to solve directly Helmholtz equation in inhomogeneous medium. Scattering field of the laser beam and complex amplitude distribution of diffracted field are obtained. Aero-optical parameters including light intensity, optical path difference, Strehl ratio, relative light intensity, beam centroid position and effective beam width are given. It can be seen that the laser beam is transmitted in the turbulent flow field, and the light intensity distribution appears gradually undulates and deviates. As passing through flow field, the beam produces a large degree of distortion, deviation and breakage due to the effect of diffraction. Effect of flow field in compressed ramp region on beam quality is significantly greater than that in fully developed turbulent region. From distribution of relative light intensity, the beam breakage seems to be related to small-scale coherent structure of flow field. The error of Strehl ratio obtained by traditional OPD mean square error or weighted mean square error estimation is greater. More accurate method should be evaluated with PSF method.

Key words: aero optics, supersonic compression corner, turbulent boundary layer, Helmholtz equation, GCV-FFT method

中图分类号:

周志超, 许凌飞, 任天荣, 顾村锋. 基于GCV-FFT方法的超声速压缩拐角流场气动光学效应计算[J]. 计算物理, 2020, 37(3): 284-298.

ZHOU Zhichao, XU Lingfei, REN Tianrong, GU Cunfeng. Aero-optical Effects in Supersonic Compression Ramp Flow Fields:GCV-FFT Method[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(3): 284-298.

参考文献

[1] JUMPER E J, FITZGERALD E J. Recent advances in aero-opotics[J]. Progress in Aerospace Sciences 2001,37:299-339.
[2] 殷兴良. 高速飞行器气动光学传输效应的工程计算方法[J]. 中国工程科学, 2006, 8(11):74-79.
[3] 李桂春. 气动光学[M]. 北京, 国防工业出版社, 2007.
[4] 李新亮, 傅德薰, 马延文, 等. 压缩折角激波-湍流边界层干扰直接数值模拟[J]. 中国科学:物理学力学天文学, 2010, 40(6):791-799.
[5] BOOKEY P B, WYCKHAM C, SMITS A J, et al. New experimental data of STBLI at DNS/LES accessible Reynolds numbers[R]. AIAA 2005-309.
[6] 武宇, 易仕和, 陈植, 等. 超声速层流/湍流压缩拐角流动结构的实验研究[J]. 物理学报, 2013,62(18):184702
[7] 童福林, 李欣, 于长平, 等. 高超声速激波湍流边界层干扰直接数值模拟研究[J]. 力学学报, 2018,50(2):197-208.
[8] WHITE M D, VISBAL M R. Investigation of turbulent shock boundary interaction unsteadiness and its effects on aero-optics in a Mach 2 ramp flow[R]. AIAA 2014-2358.
[9] MANI A, WANG M, MOIN P. Resolution requirments for aero-optical simulations[J]. Journal of Computational Physics, 2008, 227:9008-9020.
[10] WANG M, MANI A, GORDEYEV S. Physics and computation of aero-optics[J]. Annu Rev Fluid Mech, 2012, 44:299-321.
[11] GORDEYEV S, JUMPER E, HYDEN T E. Aero-optics of supersonic boundary layers[R]. AIAA 2011-1325.
[12] WHITE M D. High-order paraxial beam approximation for aero-optics[J]. Journal of Computational Physics, 2010, 229:5465-5485.
[13] 陈勇, 郭隆德, 张龙, 等. 带凹窗斜劈高速湍流气动光学效应研究[J]. 空气动力学学报, 2010, 28(6):666-671.
[14] SHI K T, MA H D. Aero-optical effects in compressible mixing layers[J]. Chinese J Comput Phys, 2010, 27(1):65-72.
[15] 甘才俊, 李烺, 马汉东, 等. 可压缩混合层光学传输效应理论分析与实验研究[J]. 物理学报, 2014, 63(5):054703
[16] 李睿劬, 宫建, 毕志献, 等. 高超声速平板边界层转捩的气动光学诊断技术[J]. 空气动力学学报, 2017, 35(1):136-140.
[17] PAN H L, LI J H, CHENG X L, et al. Model coefficient revision on supersonic turbulent pulse for aero-optic effect[J]. Chinese J Comput Phys, 2018, 35(2):194-204.
[18] 丁浩林, 易仕和, 朱杨柱, 等. 不同光线入射角度下超声速湍流边界层气动光学效应的实验研究[J]. 物理学报, 2017, 66(24):244201.
[19] 冯定华. 高速流动精细数值模拟实验研究及其在气动光学中的应用[D]. 长沙:国防科学技术大学, 2010.
[20] 田之丰. 超声速光学头罩流场精细结构及其气动光学效应的机理研究[D]. 长沙:国防科学技术大学, 2011.
[21] 郭广明, 刘洪, 张斌, 等. 混合层流场中涡结构对流速度的特性[J]. 物理学报, 2016, 65(7):074702.
[22] 马科斯·玻恩, 埃米尔·沃耳夫. 光学原理[M]. 第七版. 北京:电子工业出版社, 2009.
[23] ANDREW L C. Laser beam propagation through random media[M]. SPIE, 2005.
[24] CHACKO N, LIEBLING M, BLU T. Discretization of continuous convolution operators for accurate modeling of wave propagation in digital holography[J]. J Opt Soc Am A, 2013, 30:2012-2020.
[25] XU L F, ZHOU Z C, REN T R. Study of a weak scattering model in aero-optic simulations and its computation[J]. J Opt Soc Am A, 2017, 34:594-601.
[26] UNSER M. Sampling:50 years after Shannon[J]. Proc IEEE, 2000, 88:569-587.
[27] UNSER M. Splines:A perfect fit for signal and image processing[J]. IEEE Signal Process Mag, 1999, 16:22-38.

基于GCV-FFT方法的超声速压缩拐角流场气动光学效应计算

Aero-optical Effects in Supersonic Compression Ramp Flow Fields:GCV-FFT Method

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 6

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	吴国正, 王发杰, 程隋福, 张成鑫. 基于物理信息神经网络的内部声场正反问题数值计算[J]. 计算物理, 2022, 39(6): 687-698.
[2]	潘宏禄, 李俊红, 程晓丽, 马汉东. 气动光学效应湍流脉动模型系数修正[J]. 计算物理, 2018, 35(2): 194-204.
[3]	史可天, 马汉东. 可压缩混合层气动光学效应研究[J]. 计算物理, 2010, 27(1): 65-72.
[4]	陆昌根. 湍流边界层近壁区多个相干结构的数值模拟[J]. 计算物理, 2002, 19(5): 383-387.
[5]	邬吉明, 余德浩. 三维Helmholtz方程外问题的自然积分方程及其数值解[J]. 计算物理, 1999, 16(5): 449-456.
[6]	密柯, 章本照. 含尘离心风机叶轮磨损研究的一种二维边界层求解方法[J]. 计算物理, 1992, 9(4): 473-476.