含激波、旋涡和声波的复杂多尺度流动的直接数值模拟

摘要/Abstract

摘要：

本文主要报道了我们近年来在银河并行机上采用五阶WENO格式所做的一系列直接数值模拟研究，主要包括激波与单旋涡相互作用、激波与旋涡对相互作用、激波与三维纵向旋涡的相互作用，以及可压缩各向同性湍流。研究的主要目的是揭示激波与旋涡间相互作用中的激波动力学特性、旋涡变形、旋涡破裂和声波的产生机理，以及湍流等多尺度复杂流动的流场结构和流动机理。研究表明，高阶WENO格式具有很好的分辨率和稳定性，是研究上述包含强间断与复杂流场结构的流动的理想数值方法。研究发现，激波与强旋涡相互作用具有多级特征，即激波与初始旋涡的相互作用、反射激波与变形旋涡的相互作用、小激波与变形旋涡的相互作用。激波与旋涡对相互作用中产生的声波包含两个区域：线性区和非线性区。在线性区，激波与旋涡对相互作用产生的声波是激波分别与每个旋涡单独作用产生的声波的线性叠加；而在非线性区则与激波和耦合旋涡对的作用有关。在激波与纵向旋涡的相互作用过程中，发现旋涡破裂区存在多螺旋结构。在高初始湍流马赫数的各向同性湍流脉动场中，也发现了广泛报导的“小激波”的存在，这是可压缩湍流有别于不可压缩湍流的显著结构特征。

关键词: 直接数值模拟, 并行计算, 激波旋涡相互作用, 声波, 可压缩各向同性湍流

Abstract:

In this paper,we offer an overview of our recent DNS studies on the interaction between a shock wave and a single planar vortex,a pair of planarvortices or a longitudinal vortex,compressible isotropic turbulence through directly solving the two and three dimensional unsteady compressible NavierStokes equations using a fifth order weighted essentially nonoscillatory(WENO) finite difference scheme based on the YH parallel computer.The main purpose of these studies is to reveal the feature of shock dynamics,vortex deformation or vortex breakdown and the mechanism of sound generation in the interaction between a shock wave and vortices,as well as to explore the flow structure and mechanism of turbulence. These studies have demonstrated the excellent resolution and stability properties of the high order WENO schemes,making such schemes an ideal numerical tool for the study of shock vortex interactions in which both strong discontinuities and complex flow structures coexist.Through these studies,it is found that there is a multistage feature in the interaction between a shock wave and a strong vortex, which contains the interaction of the shock wave and the initial vortex,of the reflected shock wave and the deformed vortex and of the shocklets and the deformed vortex.The sound generated by the interaction between a shock wave and a vortex pair contains two regimes:the linear regime and the nonlinear regime.In the linear regime,the sound wave generated by the interaction of a shock wave and a vortex pair equals to the linear combination of the sound wave generated by the interactions between the same shock wave and each vortex.The second regime corresponds to the shock interaction with a coupled vortex pair.In the interaction between the shock and a longitudinal vortex,we find that there is a multihelix structure in the region of vortex breakdown.Our simulation of the compressible isotropic turbulence also confirms the existence of shocklets at sufficiently high turbulent Mach number,which is the most noticeable influence of compressibility on the structure of turbulence.

Key words: DNS;parallel computing;shockvortex interaction;sound wave;compressible isotropic turbulence

张树海，李虎，刘旭亮. 含激波、旋涡和声波的复杂多尺度流动的直接数值模拟[J]. J4, 2012, 34(8): 99-107.

ZHANG Shuhai,LI Hu,LIU Xuliang. A Direct Numerical Simulation of the Complex MultiScale Flow with Shock,Vortex and Sound Wave[J]. J4, 2012, 34(8): 99-107.

参考文献

［1］Hollingsworth M A, Richards E J. A Schlieren Study of the Interaction between a Vortex and a Shock Wave in a Shock Tube［R］. British Aeronautical Research Council Report 17985,Fluid Motion Subcommittee,17985,FM 2323,1955.
［2］Ribner H S.The Sound Generated by Interaction of a Single Vortex with a Shock Wave［R］.UTIA Report 61,University of Toronto,Institute of Aerospace Studies,1959.
［3］Inoue O, Hattori Y.Sound Generation by ShockVortex Interactions［J］.Journal of FluidMech,1999,380(1):81116.
［4］Dosanjh D S, Weeks T M.Interaction of a Starting Vortex as well as A Vortex Street with a Traveling Shock Wave［J］.AIAA Journal, 1965,3(2):216223.
［5］Naumann A,Hermanns E.On the Interaction Between a Shock Wave and A Vortex Field［J］.AGARD CP131,1973,23:123.
［6］Ribner H S. Cylindrical Sound Wave Generated by ShockVortex Interaction［J］.AIAA Journal, 1985,23:17081715.
［7］Weeks T M, Dosanjh D S.Sound Generation by ShockVortex Interaction［J］.AIAA Journal,1967,5:660669.
［8］Guichard L,Vervisch L, Domingo P.TwoDimensional Weak ShockVortex Interaction in a Mixing Zone［J］.AIAA Journal, 1995,33:17971802.
［9］Ellzey J L, Henneke M R, Picone J M, et al. The Interaction of a Shock with a Vortex:Shock Distortion and the Production of Acoustic Waves［J］.Phys Fluids,1995,7(1):172184.
［10］Ellzey J L, Henneke M R.The ShockVortex Interaction:the Origins of the Acoustic Wave［J］.Fluid Dyn Res, 1997,21(3):171184.
［11］Ellzey J L,Henneke M R.The Acoustic Wave from a ShockVortex Interaction:Comparison Between Theory and Computation［J］.Fluid Dyn Res, 2000,27(1):5364.
［12］Zhang S,Zhang Y T,Shu C W. Multistage Interaction of a Shock Wave and a Strong Vortex［J］.Phys Fluids,2005,17(11):116101.
［13］Zhang S, Zhang Y T, Shu C W. Interaction of an Oblique Shock Wave with a Pair of Parallel Vortices:Shock Dynamics and Mechanism of Sound Generation［J］.Phys Fluids,2006,18(12):126101.
［14］Zhang S, Jiang S, Zhang YT, et al. The Mechanism of Sound Generation in the Interaction Between a Shock Wave and Two CounterRotating Vortices［J］.Phys Fluids, 2009,21(7):076101.
［15］Werle H.Quelques Résultats Expérimentaux Sur Les Ailes en Fleches,Aux Faibles Vitesses［J］. Obtenus en Tunnel Hydrodynamique La Recherche Aéronautique, 1954,41.
［16］Sarpkaya T.On Stationary and Traveling Vortex Breakdown［J］.Journal of Fluid Mech, 1971,45(3):545559.
［17］Faler J H, Leibovich S.Disrupted States of Vortex Flow and Vortex Breakdown［J］.Phys Fluids,1977,20:13851400.
［18］Escudier M P.Observations of the Flow Produced in a Cylindrical Container by a Rotating Endwall［J］.Exp Fluids,1984,2(4):189196.
［19］Brucker C, Althaus W.Study of Vortex Breakdown by Particle Tracking Velocimetry (PTV) Part 3:TimeDependent Structure and Development of Breakdown Modes［J］.Exp Fluids, 1995,18(3):174186.
［20］Spall R E.Transition from Spiral to BubbleType Vortex Breakdown［J］.Phys Fluids, 1996,8(5):13301332.
［21］Zhang S,Zhang H, Zhu G.Numerical Simulation of Vortex Structure in the Leeside of Supersonic Delta Wing［J］.ACTA Aerodynamic Sinica,1996,14(4):430435.
［22］Zhang S,Zhang H, Zhu G.Numerical Simulation of Vortical Flows Over Delta Wings and Analysis of Vortex Motion［J］.ACTA Aerodynamica Sinica,1997,15(1):121129.
［23］Kandil O,Kandil H,Liu C H.Shock/Vortex Interaction and Vortex Breakdown Modes［C］∥IUTAM Symposium of Fluid Dynamics of High Angle of Attack,SpringerVerlag,1992:192212.
［24］Krause E,Thomer O, Schroder W.Strong ShockVortex Interactiona Numerical Study［J］.Comput Fluid Dyn Journal,2003,12:266278.
［25］Benjamin T B.Theory of the Vortex Breakdown Phenomenon［J］.Journal of Fluid Mech, 1962,14:593629.
［26］Bossel H H.Vortex Breakdown Flow Field［J］.Phys Fluids,1969,12(3):498508.
［27］Blackmore D,Brns M, Goullet A.A Coaxial Vortex Ring Model for Vortex Breakdown［J］.Physica D,2008,237(22):28172844.
［28］Brons M,Shen W Z,Sorensen J N,et al.The Influence of Imperfectionson the Flow Structure of Steady Vortex Breakdown Bubbles［J］.Journal of Fluid Mech,2007,578:453466.
［29］Leibovich S.The Structure of Vortex Breakdown［J］.Ann Rev Fluid Mech, 1978,10:221246.
［30］Zhang S,Zhang H,Shu C W.Topological Structure of Shock Induced Vortex Breakdown［J］.Journal of Fluid Mech,2009,639:343372.
［31］Zhang Y T,Shi J,Shu C W,et al.Numerical Viscosity and Resolution of HighOrder Weighted Essentially NonOscillatory Schemes for Compressible Flows with High Reynolds Numbers［J］.Phys Rev,2003,E 68,046709.
［32］Liu XD,Osher S, Chan T. Weighted Essentially NonOscillatory Schemes［J］.Journal of Comput Phys, 1994, 115(1):200212.
［33］Jiang G S, Shu C W. Efficient Implementation of Weighted ENO Schemes［J］.Journal of Comput Phys, 1996, 126(1):202228.
［34］Harten A, Engquist B, Osher S,et al.Uniformly High Order Essentially NonOscillatory Schemes, III［J］. Journal of Comput Phys, 1987, 71(2):231303.
［35］Shu C W, Osher S.Efficient Implementation of Essentially NonOscillatory Shock Capturing Schemes［J］. Journal of Comput Phys, 1988, 77(2):439471.
［36］Shu C W, Osher S.Efficient Implementation of Essentially NonOscillatory Shock Capturing Schemes II［J］. Journal of Comput Phys, 1989, 83(1):3278.
［37］Inoue O,Takahashi T, Hatakeyama N.Separation of Reflected Shock Waves due to Secondary Interaction with Vortices:Another Mechanism of Sound Generation［J］.Phys Fluids, 2002,14(10):37333744.
［38］Melander M V,Zabusky N J, Mcwilliams J C.Symmetric Vortex Merger in Two Dimensions:Causes and Conditions［J］.Journal of Fluid Mech, 1988,195:303340.
［39］Mitchell B E,LeleS K,Moin P.Direct Computation of the Sound from a Compressible Corotating Vortex Pair［J］.Journal of Fluid Mech, 1995,285:181202.
［40］李虎.可压缩均匀各向同性衰减湍流直接数值模拟研究:［硕士学位论文］［D］.绵阳：中国空气动力学研究与发展中心，2011.

[1]	刘屹成, 刘晓燕, 严馨. 并行平衡级联支持向量机[J]. 计算机工程与科学, 2023, 45(07): 1170-1177.
[2]	臧照虎, 李晨, 王耀华, 陈小文, 郭阳. 面向众核系统的层次化栅栏同步机制[J]. 计算机工程与科学, 2022, 44(11): 1901-1908.
[3]	张勇, 张曦, 万云博, 何先耀, 赵钟, 卢宇彤. 非结构有限体积CFD计算的网格重排序优化[J]. 计算机工程与科学, 2022, 44(10): 1721-1729.
[4]	范培勤, 过武宏, 韩梅, 唐帅, 张驰, . 水声环境特征参数并行预报方法研究[J]. 计算机工程与科学, 2021, 43(11): 1920-1925.
[5]	龚昊, 刘莹, 冯建周, 赵仁良, 冷佳旭, . 基于GPU加速的脉冲多普勒雷达信号处理[J]. 计算机工程与科学, 2021, 43(07): 1141-1149.
[6]	焦育威, 王鹏, 辛罡, . 基于采样尺度自适应的多尺度量子谐振子优化算法并行化[J]. 计算机工程与科学, 2021, 43(07): 1200-1209.
[7]	俞茂学, 贾东宁, 魏志强, 许佳立, 马广浩. 一种基于国产异构众核处理器的C++智能源码转换框架[J]. 计算机工程与科学, 2021, 43(06): 997-1005.
[8]	丁哲昭, 储根深, 胡长军, 李扬. 基于申威众核处理器的圣维南求解程序的并行与优化[J]. 计算机工程与科学, 2021, 43(05): 820-829.
[9]	丁峻宏, 苗新强, 李根国. 面向异构超算的结构分析高效并行计算方法[J]. 计算机工程与科学, 2020, 42(12): 2133-2140.
[10]	徐传福, 车永刚, 李大力, 王勇献, 王正华. 天河超级计算机上超大规模高精度计算流体力学并行计算研究进展[J]. 计算机工程与科学, 2020, 42(10高性能专刊): 1815-1826.
[11]	刘夏真, 袁武, 马文鹏, 胡晓东, 陆忠华, 张鉴. CCFD重叠网格并行算法设计和优化[J]. 计算机工程与科学, 2020, 42(10高性能专刊): 1833-1841.
[12]	赵钟，何磊，何先耀. 风雷（PHengLEI）通用CFD软件设计[J]. 计算机工程与科学, 2020, 42(02): 210-219.
[13]	李红军, 郭阳, 贾润. 归一化积相关算法加速方法研究及FPGA实现[J]. 计算机工程与科学, 2019, 41(11): 1905-1910.
[14]	李安民,计卫星,廖心怡，高建花，谈兆年，王一拙，石峰. 一种面向异构计算的结构化并行编程框架[J]. 计算机工程与科学, 2019, 41(03): 424-432.
[15]	聂含伊，杨希，张文喆. 面向多领域的高性能计算机应用综述[J]. 计算机工程与科学, 2018, 40(增刊S1): 145-153.