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

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

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

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

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

激波/湍流边界层干扰物面剪切应力统计特性

为了揭示激波/湍流边界层干扰区内物面剪切应力统计特性的演化规律,采用直接数值模拟方法对来流马赫数2.9、12°激波角的入射激波与平板湍流边界层相互作用问题进行了研究。通过与风洞试验数据的比较分析,验证了计算结果的可靠性。系统地探究了干扰区内物面剪切应力的典型特征,如预乘谱、概率密度分布和相干结构等。研究结果表明,分离激波的低频振荡运动对流向及展向分量的预乘谱均没有实质影响,其脉动仍以高频特征为主,低频能量变化较小。干扰区内流向剪切应力概率密度函数分布变化剧烈,分离泡内对数正态分布规律不再适用,而展向剪切应力在干扰区内与正态分布较为接近。相较于上游湍流边界层,分离泡内物面剪切应力矢量夹角与幅值的联合概率密度变化显著,峰值概率降低,峰值范围增大。此外,流向剪切应力脉动场的本征正交分解分析指出,主能量模态与分离激波的低频振荡以及下游再附边界层内的Görtler-like流向涡结构密切相关。     

欠膨胀火箭多股喷流干扰流场的数值模拟

本文利用高精度的TVD有限体积差分格式,对欠膨胀火箭多股自由喷流的干扰流场进行了数值研究,给出了两股自由射流的干扰流场的复杂激波系结构,包括内、外层反射激波、相交激波等波系结构。数值结果与实验结果较为吻合。     

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

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

弹体超声速绕流与底喷流干扰流场的数值模拟

利用ＴＶＤ有限体积格式数值模拟了弹体超声速绕流与底喷注相互干扰的复杂流场。采用四步Ｒｕｎｇ－Ｋｕｔｔａ时间推进多步法及“冻结粘性”技术,加速了迭代收敛,节省了计算机时。     

出口边界条件对跨音扩压器流场的影响

本文对音速及等压两种不同的出口边界下跨音速扩压器的流场做了实验研究.扩压器的扩张角为6°,激波波前马赫数均为1.335.实验结果表明:在等压出口边界条件时,激波自恃振荡的幅度及分离包的尺度比音速出口条件时要大得多.     

STUDY OF SWEPT SHOCK WAVE AND BOUNDARY LAYER INTERACTIONS

Ｔｈｅｉｎｔｅｎｓｉｖｅｒｅｓｅａｒｃｈｏｎｔｈｅｓｕｂｊｅｃｔｏｆｓｗｅｐｔｓｈｏｃｋｗａｖｅ／ｔｕｒｂｕｌｅｎｔｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎｓｈａｓｂｅｅｎｃｏｎｄｕｃｔｅｄａｎｄｓｕｐｐｏｒｔｅｄｆｏｒｍａｎｙｙｅａｒｓｂ...     

激波与氦气泡相互作用的实验与数值研究

提出了一种在激波管中实现运动激波与球形气一气界面相互作用的研究方案。以肥皂膜作为气体交界面,高压火花为光源,采用阴影法,在马赫数M=1．2的情况下,在矩形激波管中开展了一系列的实验研究,获得了激波作用下球形氦气泡界面发展变化的相关实验结果。还通过求解二维轴对称的Euler方程,结合网格自适应策略,对实验过程进行了数值模拟。通过与实验结果的对比分析,证实了涡量的产生和分布对于界面的变形和发展有着重要意义。     

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

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

激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰,平板边界层上形成分离区,又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游,引入小扰动的TS波后,扰动波传播通过带有这些复杂流动现象的流场时,扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析,研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现,激波的出现不同于一般的压缩波,在亚声速区与超声速区对扰动波演化的影响是不同的,此外,分离对扰动有稳定作用。     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

超声速流动中横向射流流场的影响参数

采用高精度的Weighted ENO格式，结合两方程湍流模型，准确模拟了二次射流形成的干扰流场，详细地描述了平板上单股射流干扰流场和喷管扩张段二次射流干扰流场中的激波、流动分离和旋涡运动，同实验结果进行了比较。探讨了射流／主流总压比、射流宽度。以及射流与来流夹角对射流穿透深度、分离距离等影响，揭示了二次射流推力短量控制干扰流场的控制机理。     

间断Galerkin方法计算叶型转捩流动

使用间断Galerkin方法研究叶型转捩流动,进行了大涡模拟(LES)和转捩模型的求解,对T106低压透平和Zierke-Deutsch压气机叶栅内的流动进行了计算.通过T106低压透平的计算对LES和转捩模型进行了比较,结果表明两种方法得到的压力分布和分离泡位置均与实验值吻合良好.LES得到的分离泡的轴向位置为0.145~0.165m,转捩模型得到的分离泡的轴向位置为0.150~0.156m.LES可以再现复杂的瞬时流动细节,对于深入研究流动机理很有意义,而转捩模型尽管无法获得足够的流动细节,但是也能预测边界层的分离和转捩现象,并且结果与LES时均结果相差不大,对于工程应用很有价值.使用转捩模型计算Zierke-Deutsch压气机叶栅内的流动也获得了与实验值符合的结果.      

有/无尾迹作用下低压涡轮叶栅分离边界层转捩的大涡模拟

采用经过大量算例验证的可压缩大涡模拟求解器对雷诺数为60154、马赫数为0.402的低压涡轮叶栅T106D-EIZ进行了细致模拟,计算了定常来流和周期性尾迹来流两种工况.对计算结果的分析表明:定常来流工况下,叶片吸力面后部出现大尺寸的层流分离泡,分离剪切层的转捩过程受Kelvin-Helmholtz （K-H）不稳定性控制;尾迹来流工况下,由于来流尾迹的周期性扫掠,时均分离泡尺寸变小,叶栅总压损失降低.对相位平均和瞬态流场的分析表明,尾迹引起的逆射流使分离点后移,形成卷起涡结构,逆射流掠过卷起涡的过程中与其发生强烈的相互作用,产生大量气动损失,而后卷起涡破碎,流动转捩为湍流.      

尾迹对低压涡轮边界层稳定性的影响

在低压涡轮内部上游周期性扫过的尾迹是下游叶片排内最重要的非定常扰动源,会对边界层的失稳和转捩过程产生决定性的影响.在施加低压涡轮负荷分布的平板上,采用热线对雷诺数为130000工况下,有/无尾迹作用时分离剪切层中的速度波动进行了测量.通过对比两种工况下分离剪切层中的扰动发展,分析了尾迹对分离泡稳定性的影响.结果表明:在无尾迹来流工况下,扰动在分离泡前部的增长符合线性不稳定机制,扰动增长至饱和后非线性机制开始起作用.主导分离剪切层的失稳过程的为K-H(Kelvin-Helmholtz)无黏不稳定机制.尾迹与分离剪切层的相互作用产生的扰动的增长在分离泡前部同样满足线性不稳定机制.尾迹加速了分离剪切层的失稳转捩进程,从而抑制了分离.      

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。     

航空涡轮发动机燃烧室内流场的PIV测量

设计了一种基于旋风分离原理的高压粒子发生器,并成功应用于高压状态下的航空涡轮发动机燃烧室内流场的PIV（粒子图像测速法）测量.在氢氧燃烧加热来流温度为813K、燃烧室压力为2.78MPa条件下,应用PIV技术开展了航空涡轮发动机单头部燃烧室复杂内流场测量研究,实现了高温高压条件下强旋流、强扰流、宽速域流场的PIV测量,获得了接近燃烧室工作压力工况下的流场速度和流场精细结构.结果表明:该型燃烧室内流场存在多处旋涡结构,形成回流区;流场旋流作用强,横截面流场存在顺时针大涡;主燃孔射流和掺混孔射流作用明显,射流穿透深度较大,对流场结构影响显著;高温高压状态下,流场结构与常温中压状态类似.