膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

压缩拐角激波与旁路转捩边界层干扰数值研究

为了研究激波与旁路转捩边界层的干扰机理,采用直接数值模拟(DNS)方法对来流马赫数Ma∞=2.9,24°压缩拐角内激波与转捩边界层的相互作用进行了系统的研究。考察了旁路转捩干扰下压缩拐角内分离区形态和激波波系结构的典型特征。比较了转捩干扰与湍流干扰流动结构的差异,并分析了造成差异的原因。研究了拐角内转捩边界层的演化特性,探讨了转捩干扰下脉动峰值压力和峰值摩阻的分布规律及形成机制。研究结果表明:相较于湍流干扰,两侧发卡涡串的展向挤压使得分离区起始点以V字型分布,且分离激波沿展向以破碎状态为主,激波脚呈现多层结构;拐角内的干扰作用急剧加速了边界层的转捩过程;转捩干扰下的拐角内峰值脉动压力以单峰结构出现在分离区的下游,同时干扰区内的强湍动能和高雷诺剪切应力使得其局部峰值摩阻系数要高于湍流干扰。     

马赫数6柱-裙激波/边界层干扰直接模拟

高超声速激波/边界层干扰比超声速工况下具有更强的可压缩效应,再附之后会形成极高的局部力/热载荷,严重影响飞行器飞行安全。而激波/湍流边界层干扰区附近流动的三维特性使得流动更加复杂而难以预测。采用直接数值模拟对高超声速条件下的柱-裙激波/湍流边界层干扰进行了详细研究,特别是对Görtler涡结构对分离泡、物面压力和热流造成的展向差异开展了定性和定量分析。研究发现,干扰区附近的分离泡结构呈现出明显的三维效应,且Görtler涡展向分离位置截面的分离泡要明显小于再附位置,而这两个截面上分离泡的运动基本同步,没有明显的延迟或超前现象。物面压力和热流在展向出现显著的不均匀性,展向再附位置的平均压力和热流要比展向分离位置分别高13%和16.2%,脉动压力和热流比展向分离位置分别高28%和20%。截面流向速度特征正交分解结果显示两个位置上的能量都集中在剪切层附近,并且展向再附位置上低频模态占有更高的能量。在采用模态重构流场分析分离区面积发现,展向分离位置重构误差更小,而展向再附位置上的重构误差收敛更快。     

马赫6柱-裙构型激波/湍流边界层干扰摩阻统计特性

为探究激波/边界层干扰表面摩阻统计特性,对马赫6柱-裙构型的激波/湍流边界层流场进行了直接数值模拟,推导了柱坐标形式的时均摩阻分解公式,将其与对流效应、流向不均匀性效应、分子黏性耗散效应、曲率效应和湍动能耗散关联起来,在此基础上对摩阻脉动和时均摩阻的统计特性进行了研究。概率密度分布结果表明干扰区内摩阻脉动概率密度明显偏离正态分布,摩阻脉动在该区域的间歇性较强,这与湍流脉动在激波干扰后存在显著内在可压缩性效应有关;谱功率密度结果表明摩阻脉动的能量在干扰前中频部分的能量占据主导,峰值对应频率约为0.14,干扰区后附近区域峰值位置转移到高频区域。采用摩阻分解公式对激波干扰前后的时均摩阻进行分解,结果表明时均摩阻的主要贡献来自分子黏性耗散效应和湍动能耗散效应;在激波干扰前分子黏性应力效应占据主导位置,干扰后湍动能耗散效应占据主导位置;激波干扰后对流效应明显增强,曲率效应的贡献在激波干扰前后基本保持不变;流向不均匀项贡献由正贡献变为负贡献,且所占比例有所上升,这主要是由压力梯度项的变化导致的。对分解后各项积分内函数沿法向分布的分析表明,分解后各项在边界层内分布呈现较大的差异,特别地,湍动能耗散项...     

高超声速有攻角锥裙直接数值模拟

采用直接数值模拟方法对有攻角的高超声速7°～34°锥裙开展了数值研究,通过对比0°、90°、180°周向子午面,评估了三维横流效应对激波/边界层干扰的影响规律和作用机制,包括壁面压力、摩阻、热流分布,分离泡非定常运动,再附边界层演化等。研究发现,不同周向站位均出现流动分离,横流区、迎风区内复杂激波结构与边界层相互作用导致壁面压力、摩阻、热流显著升高。热流与压力的比值在干扰区上升后由于再附降低,而热流与摩阻的雷诺比拟关系在分离区则完全失效。分离泡面积脉动的功率谱结果表明,分离泡非定常膨胀/收缩运动呈低频特征,且分离泡呼吸与激波低频振荡在横流区密切相关,在迎风区存在迟滞,而在背风区不相关。速度脉动场的本征正交分解结果表明,分离区的低频特征与低阶模态相应的剪切层附近大尺度结构相关。对下游再附边界层演化分析指出,攻角的存在导致雷诺应力在再附点附近大幅增强,其流向分量的恢复最为迅速,雷诺应力分量的峰值位置在背风区沿流向持续外移,而在迎风区、横流区已迅速向内层恢复。此外,雷诺应力各向异性不变量分布进一步表明干扰下游的近壁区湍流各向异性峰值在背风区弱于迎风区。     

侧壁约束下压缩激波/湍流边界层相互作用低频不稳定性实验研究

侧壁面约束条件下压缩激波/湍流边界层相互作用（Shock-Wave/Turbulent Boundary-Layer Interactions, SWTBLIs）呈现出区别于经典二维SWTBLIs的强三维特性，尤其是低频不稳定性问题自发现以来就受到了广泛关注，而其与下游压缩环境之间的关联性尚不明确。本文基于带侧壁约束的超声速压缩拐角构建三维受限的SWTBLIs流动，并通过自由射流风洞实验进行了系统研究，其中来流Ma∞=2.5、湍流边界层发展厚度对矩形流道宽度的占比约0.08、压缩拐角在12°~24°的宽范围内变化。通过不同压缩拐角下流场结构演化及壁面压力脉动信号谱特性的对比，揭示了侧壁面约束条件下压缩拐角对低频大尺度振荡时间-空间演化行为的影响规律和作用机制。研究发现：压缩拐角角度较小时，侧壁面约束导致的“拖尾效应”有助于延缓压缩拐角附近流动分离模式的过渡，并抑制压缩激波的不稳定性；压缩拐角达到20°及以上时，压缩拐角前沿逐渐发展为大尺度流动分离模式，侧壁面约束的角区将首先发展出能量集中的低频脉动，并逐步演化为两种峰值频率分别为约50Hz、200Hz量级的不同大尺度间歇性低频振荡叠加行为，且伴随着压缩激波形态的显著改变；两种低频振荡行为都属于SWTBLIs流场的全局性不稳定振荡，但沿着展向振荡强度差异显著，角区内的振荡强度远远高于中心主流区域。     

高分辨率激波/边界层干扰时间演化过程分析

在超声速风洞中开展了湍流边界层与圆柱相互作用流场研究，试验马赫数为3.4和3.8。圆柱安装在试验段底板上，安装位置的边界层为充分发展的湍流边界层，研究了圆柱直径和高度对流场结构和压力脉动的影响。采用基于纳米示踪的平面激光散射(NPLS)技术获取了流向和展向流场精细结构，激波系和马蹄涡结构均可清晰分辨。通过展向流场图像可以发现干扰区内激波与湍流结构的相互作用具有明显的非定常性。采用动态压力传感器测量了圆柱前方相互作用区域的压力脉动特性，在激波足区域压力呈现11～38 kHz的宽频分布，推测主要由激波足与涡结构相互作用及滞止区涡结构的破碎引起。随着圆柱高度的增加，激波足附近测点对应的特征频率有所降低；上游测点则发现了0～3 kHz低频区能量的增强，这主要是由分离区引起的，表明在一定高度范围内高度的增加增强了流动分离。     

超声速压缩拐角激波/边界层干扰动力学模态分解

压缩拐角激波与边界层干扰问题广泛存在于高速飞行器的外部和内部流动中,其非定常复杂流场结构对飞行器气动性能影响显著。动力学模态分析将有助于进一步加深理解激波与边界层干扰流场不同特征频率对应的流动结构及动力学特性,为揭示其复杂流动机理提供参考。本文采用动态模态分解(DMD)方法对来流马赫数为2.9、24°压缩拐角内激波与超声速边界层干扰下的非定常流动进行了模态分析。评估了稀疏改进动态模态分解方法在压缩拐角流动中的适用性,研究了湍流干扰和转捩干扰下典型特征频率对应的动力学模态空间结构差异及其原因,分析了转捩边界层展向非均匀性对低频/高频模态动力学机制的影响规律。研究发现,湍流干扰与转捩干扰下拐角干扰区内均存在两类截然不同的动力学模态:低频模态和高频模态。低频模态结构集中在分离激波及分离泡剪切层的根部,表征为分离泡的大尺度膨胀和收缩运动;高频模态空间分布则以平均声速线附近正负交替结构为主,对应为边界层内不稳定波沿剪切层往下游的传播。转捩边界层的展向结构对低频模态运动特性影响明显,而对高频模态的影响则相对较小。     

Direct numerical simulation of impinging shock wave and turbulent boundary layer interaction over a wavy-wall

The interaction of an impinging oblique shock wave with an angle of 30° and a supersonic turbulent boundary layer at Ma_∞=2.9 and Re_θ = 2400 over a wavy-wall is investigated through direct numerical simulation and compared with the interaction on a flat-plate under the same flow conditions. A sinusoidal wave with amplitude to wavelength ratio of 0.26 moves in the streamwise direction and is uniformly distributed across the spanwise direction. The influences of the wavy-wall on the interaction, including the characterization of the flow field, the skin-friction, pressure and the budget of turbulence kinetic energy, are systematically studied. The region of separation grows slightly and decomposes into four bubbles. Local peaks of skin-friction are observed at the rear part of the interaction region. The low-frequency shock motion can be seen in the wall pressure spectra. Analyses of the turbulence kinetic energy budget indicate that both diffusion and transport significantly increase near the crests, balanced by an amplified dissipation in the near-wall region. Proper orthogonal decomposition analyses show that the most energetic structures are associated with the separated shock and the shear layer over the bubbles. Only the bubbles in the first two troughs are dominated by a low-frequency enlargement or shrinkage.     

转捩对压缩拐角激波/边界层干扰分离泡的影响

为了研究转捩对压缩拐角内分离泡结构的影响,进行了来流马赫数2.9,24°压缩拐角激波/转捩边界层干扰的直接数值模拟（DNS）。通过在拐角上游平板的不同流向位置处添加周期性吹吸扰动激发流动转捩,使得转捩不同阶段进入拐角入口,从而在拐角内产生激波/转捩边界层的相互干扰。计算得到的平均速度剖面、壁面压力分布以及分离泡大小与风洞试验及以往直接数值模拟结果吻合较好,验证了计算结果的可靠性。研究了转捩过程对角部干扰区内分离泡结构的影响规律,分析比较了不同转捩阶段下角部分离区内湍动能的生成、耗散和分配机制。研究结果表明：转捩初期的拟序涡结构对分离泡尺度及形状影响最大,发卡涡包在角部拐点附近发生展向融合,并在角部区域形成湍流斑,此时分离泡尺度最小,形状呈现中间高两边低的山峰型。随着转捩的发展,分离区内湍动能生成和近壁区的耗散逐步降低,此时输运项起到了主要的平衡作用。     

Characteristics of reattached boundary layer in shock wave and turbulent boundary layer interaction

The reattached boundary layer in the interaction of an oblique shock wave with a flatplate turbulent boundary layer at Mach number 2.25 is studied by means of Direct Numerical Simulation(DNS). The numerical results are carefully compared with available experimental and DNS data in terms of turbulence statistics, wall pressure and skin friction. The coherent vortex structures are significantly enhanced due to the shock interaction, and the reattached boundary layer is characterized by large-scale...     

MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

近壁湍流边界层壁面局部扰动诱导形成相干结构的动力学模型

在壁面局部脉冲扰动作用下,采用直接数值模拟的方法,数值研究了湍流边界层近壁区两个相干结构形成的理论机制和非线性作用问题.然后,通过数值计算的手段,进一步探讨了近壁湍流边界层流中两个相干结构之间的非线性作用机理,分别确定它们在流向、展向方向上排列的不同结构类型分布的两个相干结构之间相互作用的效果有何区别以及它们之间间距的大小对相干结构的产生和发展有什么影响,从而可能合理地解释平面剪切湍流中低速条纹结构出现的原因.