带凹腔的支板火焰稳定器三维大涡模拟

为研究支板冷态流场,应用LES(大涡模拟)的方法对带凹腔的支板火焰稳定器在Ma=0.06和 Ma=0.09两种工况下分别进行三维数值模拟.对比试验和仿真结果表明:三维大涡模拟结果具有与试验相似的流场结构,定量对比试验和LES结果的尾涡脱落频率的最大相对误差为9.4%;支板尾部截面流向时均速度的最大相对误差为7.7%.表明数值模拟方法和边界条件设置的正确性和可信度.又对支板的三维流场进行分析:速度越大,支板尾涡的尺寸和强度越大;支板流场存在明显的三维效应,并且随着速度变大,三维效应越明显.      

模型固体火箭发动机内涡脱落过程的大涡模拟

为研究后向台阶结构引起的涡脱落对固体火箭发动机内流动稳定性的影响，采用大涡模拟技术，对模型发动机内的流动进行了数值模拟，获取了发动机燃烧室内压强随时间的变化曲线及其振荡频率，并和已有的计算结果及实验结果进行了对比分析。结果表明采用LES方法的计算结果明显优于采用k-ε湍流模型的计算结果；不考虑燃烧时模型发动机内由涡脱落引起的压强振荡的振幅较小，但频率较高，且随来流速度的增加而增大；在后向台阶这种结构中，涡一般会在台阶的边沿形成，并在其下游周期性地脱落；在保证y^＋较小的情况下，网格大小对计算结果的影响不明显。     

微缝共振器的非线性吸声机理

为模拟高声压级时微缝处涡声相互作用现象并计算微缝声阻抗,发展了一种离散涡模型.联合微缝离散涡模型和一维声传播模型,计算了微缝共振器的声阻抗和吸声系数.数值结果表明,微缝声学性质主要取决于涡脱落和脱落涡层的演化;共振器的吸声系数随着入射声压级的增加而增大;入射声波频率对共振器吸声系数有显著影响.研究结果和先前直接数值模拟...     

基于 LES 方法的增升装置气动噪声特性分析

在气动噪声数值计算中,流场的求解精度对涡流扰动的细节计算以及声学的求解结果有着重要的影响。本文应用 LES 方法对增升装置的流场进行数值模拟,采用可穿透积分面的 Ffcows Wil1iams-Hawkings(FW-H)积分方法进行远场噪声计算。采用圆柱绕流算例对本文的数值计算方法进行了验证,验证结果表明：本文所使用的LES 方法能准确地捕捉到涡脱落、流动分离等非定常流动现象,可为远场气动噪声的计算提供精确的近场流动的数值解;基于 FW-H 的声类比方法能够精确高效求解远场气动噪声。在此基础上,对增升装置噪声产生的流动特性、远场特性、风速影响等进行了数值模拟研究。结果表明：缝翼产生气动噪声的主要原因是,流动在缝翼和主翼之间的凹槽形成的不稳定波以及缝翼钝后缘的小脱落涡;襟翼产生气动噪声的主要原因是,襟翼附近由于流动分离产生的高频的小尺度不稳定涡和低频的大尺度涡。     

三维效应对超声速湍流燃烧流动大涡模拟的影响研究

本文对空间发展Mac=0.61超声速湍流燃烧混合层流动开展了二维/三维大涡模拟（LES）研究,探讨了三维效应在超声速湍流燃烧流动的大涡模拟中的重要性。研究发现,相较于三维LES计算结果,二维LES大幅低估了超声速湍流燃烧混合层的增长率,湍动能的二维LES结果显著偏大而雷诺剪切应力显著偏小。因而,二维LES难以准确模拟超声速湍流燃烧混合层的增长率和湍流脉动特性。二维流动中漩涡拉伸、压缩机制失效,致使二维LES的漩涡结构与三维LES结果显著不同：二维LES计算所得流场呈现清晰的旋涡生成、发展、并对和破碎等演化现象,而三维LES的湍流结构呈不规则状。湍流与燃烧间存在强烈的耦合相互作用,由于二维LES难以准确模拟湍流特性,燃烧场进而在二维LES下呈现一系列特殊的、不同于真实燃烧流动情况的数值现象。漩涡运动带动燃料和氧化剂掺混,掺混界面为火焰前锋。由于漩涡结构不同,致使二维LES的火焰前锋面面积不同于三维LES结果,导致其难以获得准确的燃料效率,最终导致二维LES计算获得的统计平均火焰温度显著偏低。     

开缝圆柱近场尾流的大涡模拟和实验测量

基于Fluent软件平台采用大涡模拟(LES)的方法对开缝圆柱绕流进行了三维数值仿真,计算了五组不同雷诺数(Re=1500、3000、4400、5837、7200)情况下缝宽比s/d=0.15的开缝圆柱绕流流场。通过定性的流动显示和定量的PIV实验,对s/d=0.15情况下开缝圆柱内部缝隙流动以及旋涡脱落情况进行了细致研究,并测量了基准圆柱与开缝圆柱的脱落涡频率。数值模拟与实验测量结果符合得很好。实验和计算的结果表明,在一定雷诺数范围内,开缝圆柱缝隙内部流体呈现周期性的震荡,边界层底层经历"吸入"和"吹出"现象,根本上改变了近区尾流的结构。开缝圆柱尾迹呈现出明显的三维流向涡形态。缝隙内流体的上下震荡对圆柱绕流的周期性产生了显著影响,开缝圆柱的旋涡脱落频率和斯特劳哈尔数均高于相同雷诺数下的基准圆柱。     

空间发展混合流的大涡模拟

采用大涡模拟方法数值模拟了空间发展混合流，通过在平均流中引入小扰动，研究了混合流大涡结构的产生和演化过程，捕捉了展向涡的卷起、配对、合并，以及二次流向涡的出现等大尺度的三维拟序结构。与实验及直接数值模拟的结果相比，吻合程度较好，表明本文所采用的LES方法对湍流的数值模拟是切实可行的。     

低旋流预混燃烧稳燃机理的大涡模拟

采用大涡模拟(LES)方法对比了高低旋流的预混燃烧特性,并分析了低旋流的稳燃机理。结果表明:LES耦合部分搅拌反应燃烧模型(PaSR)可以准确的捕捉旋流预混燃烧的细小火焰结构和流场分布。不同旋流强度引发不同的流场特性,低旋流诱发较弱的回流区,而高旋流会诱发很强的回流区,这种强烈的流场特征往往会引起振动频率极低的进动涡核拟序结构(PVC)。通过涡脱落和燃烧温度分布的瞬时耦合关系,发现火焰剪切层引起的内外涡脱落频率差是低旋流火焰稳定的最主要因素,低旋流燃烧受涡旋结构的影响交替出现W和V型火焰。     

低亚声速和跨声速矩形柱绕流的大涡模拟研究

采用基于非结构网格的有限体/有限元混合格式和大涡模拟的方法求解可压缩Navier-Stokes方程,研究了不同长宽比矩形柱低亚声速和跨声速绕流的流动特性。在雷诺数为22000时,对来流马赫数等于0.1和0.75,截面长宽比分别为1∶1、2∶1、3∶1和4∶1的矩形柱绕流进行了大涡模拟,以研究长宽比和压缩性对矩形柱绕流流场的影响。马赫数为0.1时,Strouhal数随着长宽比的增大先降低再增大然后再降低;长宽比为3∶1和4∶1时会有流动的再附产生;柱体上表面的三维特性在长宽比大时更明显。马赫数为0.75时,Strouhal数随着长宽比的增大逐渐减小;湍流脉动和涡脱落受到抑制;方柱的近尾迹区域,有两种形成机制不同的局部超声速区。     

DNS/LES方法在剪切湍流模拟中的应用

提出采用DNS/LES混合方法进行自由剪切类混合流动湍流研究,该方法兼顾了直接数值模拟(DNS)方法在转捩线性范围内的高精确度模拟和大涡模拟(LES)方法在转捩非线性及全湍流区的高效模拟两方面优点,通过对低对流Mach数流场计算表明该方法可以较好地模拟流场的转捩过程.同时,文中还对比实验分析了湍流场对光学传输效应的影响,给出了与实验一致的结论.     

Influence of thermal inhibitor position and temperature on vortex-shedding-driven pressure oscillations

Vortex-acoustic coupling is one of the most important potential sources of combustion instability in solid rocket motors (SRMs). Based on the Von Karman Institute for Fluid Dynamics (VKI) experimental motor, the influence of the thermal inhibitor position and temperature on vortex-shedding-driven pressure oscillations is numerically studied via the large eddy simulation (LES) method. The simulation results demonstrate that vortex shedding is a periodic process and its accurate frequency can be numerically obtained. Acoustic modes could be easily excited by vortex shedding. The vortex shedding frequency and second acoustic frequency dominate the pressure oscillation characteristics in the chamber. Thermal inhibitor position and gas temperature have little effect on vortex shedding frequency, but have great impact on pressure oscillation amplitude. Pressure amplitude is much higher when the thermal inhibitor locates at the acoustic velocity anti-nodes. The farther the thermal inhibitor is to the nozzle head, the more vortex energy would be dissipated by the turbulence. Therefore, the vortex shedding amplitude at the second acoustic velocity antinode near 3/4L (L is chamber length) is larger than those of others. Besides, the natural acoustic frequencies increase with the gas temperature. As the vortex shedding frequency departs from the natural acoustic frequency, the vortex-acoustic feedback loop is decoupled. Consequently, both the vortex shedding and acoustic amplitudes decrease rapidly.     

低雷诺数下压气机的二次流旋涡结构

为了探究高空低雷诺数条件下跨声速压气机的流动规律,对NASA Rotor37进行单通道数值模拟,探索其在低雷诺数进气条件下二次流的旋涡结构.研究发现:马蹄涡压力面分支诱发压力面角区诱导涡,壁角涡形成了顺流和逆流的两段式结构,脱落涡由叶根角区发展起来后不断从尾缘脱落,泄漏涡近失速点仅局部破裂不是失稳触发的主要原因.通道中的激波系诱发了吸力面和压力面的两个径向涡,压力面径向涡构成闭合的气泡式分离,吸力面径向涡在叶顶的破碎诱导产生分离涡,触发了低雷诺数下压气机的失稳.流场旋涡结构由马蹄涡、壁角涡、径向涡、泄漏涡、分离涡、脱落涡6个大尺度旋涡以及其他小尺度旋涡组成.      

Π形板梁分离流扭转颤振机理数值研究

采用离散涡方法及流场可视化技术识别桥梁Π形板梁断面的旋涡脱落机制.流场可视化显示,来流绕过振动的Π形板梁断面时,在主梁断面前缘下部分离泡形成并发展为主涡结构,主涡结构沿主梁断面下表面漂移并经断面后缘下端进入尾流,主涡结构的形成时间及沿主梁断面的漂移过程对气动力起直接的主导作用;计算结果也表明恰当地布置稳定板能将颤振形态从单自由度扭转颤振转化为两自由度弯扭耦合颤振,颤振临界风速可明显提高.     

大负荷低压涡轮叶型分离转捩流动的大涡模拟

应用动力模式大涡模拟数值方法,对来流无扰动、不可压、雷诺数为5×104(基于进口速度和轴向弦长),定常来流条件下大负荷低压涡轮叶型(Pak B)叶型吸力面非定常分离转捩流动进行了三维数值模拟。在与相关实验数据的对比基础之上,对非定常流动物理信息进行了详细的分析讨论,揭示了计算来流状态下的Pak B叶型吸力面非定常分离转捩流动机理。结果表明,由无粘Kelvin-Helmholtz机制产生的、空间线性增长的初始二维不稳定性在分离剪切内诱导展向旋涡形成并脱落,脱落过程中的展向涡在非线性增长的三维不稳定性作用下发生变形并最终破碎成湍流。计算得到的Kelvin-Helmholtz不稳定性特征频率处于相关实验测量范围内。     

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

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

轴流压气机转子三维非定常流场旋涡流动分析

通过数值模拟,对一个低压压气机转子的非定常流场进行结构分析和频谱分析。通过对于流场时均图和等熵动画图的分析,捕捉到了流场的流动特性和流场中旋涡的存在范围,探索了叶片前缘和尾缘的脱落旋涡以及角涡之间的相互作用,相互影响以及这几种旋涡随时间空间变化的运动规律。对于有脱落现象的旋涡,利用动态演示图估测和频谱分析的方法判定其涡脱落频率,得到了尾缘脱落旋涡和角涡的两种脱落频率。本文的工作为今后采取主动、被动控制,改进压气机的性能提供参考依据。      

机翼尾迹脱落涡噪声的理论研究

Ｍｉｃｈｅｌ等人１９９８年应用平面传声器阵列对飞机过顶噪声进行的测量研究首次发现，机翼尾迹脱落涡噪声是某些类型飞机重要的噪声源。为发展一种预测这种噪声源的理论预测模型，应用ｖｏｎ Ｋａｒｍａｎ涡街模型模拟二维机翼下游尾迹脱落涡，尾迹涡的强度和脱落频率应用这个模型进行计算。基于Ｈｏｗｅ后缘噪声理论，并结合尾迹模型，本文发展了一种预测脱落涡噪声声压级和指向特征的气动声学模型。对６架现代商用飞机的机翼尾迹脱落涡噪声的计算表明，本文理论模型预测的涡脱落频率、声压级以及噪声的指向性等与实验测量结果有较好的一致性。     

用SPIV测量结果分析叶轮机内流动损失的方法初探

基于大涡模拟和PIV(particle image velocimetry)测量的相似性,利用SPIV(stereoscopic PIV)测量结果,可以使用亚格子应力模型求解流场中的湍流耗散率.对比不同亚格子应力模型求出的结果,混合模型求出的结果较为准确,但需要合理选择模型中的系数C_Sm和CAm.利用该方法分析低速大尺寸压气机试验台转子近叶尖区域的SPIV测量结果,发现在设计状态流场中的损失主要源于叶尖泄漏涡,而在近失速状态则主要源于叶尖泄漏涡和角区旋涡.