不同进口条件对OGV/ 前置扩压器流场大涡 模拟结果的影响

为了研究不同进口条件对压气机出口导向叶片(Outlet Guide Vane,OGV)/前置扩压器流场性能的影响,分别采用合成涡方法与白噪声方法生成进口条件,研究不同进口条件对流场数值模拟结果的影响,并与文献给出的试验结果进行对比。结果表明:采用合成涡方法生成的进口速度场满足试验测得的时均速度分布与脉动均方根分布,且其生成的进口湍流结构可传播到下游较远处;合成涡方法能够较好地预测OGV内的流场,与白噪声方法相比,所预测的OGV吸力面分离区较小;在OGV出口截面、扩压器内和扩压器出口截面,采用2种方法得到的速度分布相差不大,与试验相比分布趋势相同。     

均匀湍流场下矩形断面流场及气动参数数值模拟

采用大涡模拟方法数值模拟了两种矩形断面(B/H=1,B/H=2)在不同湍流来流下的流场及气动参数,采用谐波合成方法生成与目标谱一致的脉动入口风速.数值计算结果表明来流湍流对断面的流场及气动参数有较明显的影响,通过与均匀来流的计算结果作对比可以看出,湍流导致断面的阻力系数减小明显,其主要原因为背压区负压系数的增加,与实验的结果较一致;从流场分布看湍流主要影响断面尾流涡的形成,这种影响会影响尾流涡的形成.     

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

钝头双锥喷流致冷流场结构及密度脉动特性

采用大涡模拟方法对钝头双锥喷流致冷流场开展了数值模拟,研究了超声速喷流混合流场结构特征及密度脉动特性。大涡模拟方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runger-Kutta方法。数值模拟清晰地捕捉到了流场波系结构,精细地预测了流动发生失稳、转捩以及发展为充分发展湍流的物理过程,直接获得了流场密度脉动特性。通过有、无喷流状态对称面流场的对比,发现超声速喷流能够有效冷却光学窗口;喷流与主流形成的混合层不稳定,很快发生失稳和转捩,形成大尺度湍流结构,进而引起强烈的密度脉动。此外,获得了钝头双锥整体模型喷流致冷流场的空间发展形态特征。      

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

用于可压缩壁面湍流DNS及LES入口湍流边界条件的回收-调节方法综述

很多湍流流动问题的大涡模拟和直接数值模拟的入口湍流脉动边界条件对于计算结果的影响很大,本文叙述了用于可压缩壁面湍流入口脉动生成的一种方法——“回收－调节”方法（recycling-rescaling method）,叙述了这种方法的理论基础及简要推导过程,介绍了几种常用的“回收－调节”方法并指出了实现中的几个关键问题,分析了这种方法中存在的不足并对其应用进行了总结评述。本文可为超声速流动、燃烧问题的大涡模拟和直接数值模拟方法的发展提供参考。     

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

微吹气对湍流平板边界层流动特性的影响及其减阻机理

微吹气技术能够改变平板湍流流场结构,减小平板表面的摩擦阻力。采用直接数值模拟方法,计算了来流马赫数0.7条件下,流场流过光滑平板和NASA-PN2多孔平板表面两种情况,通过对比这两个算例的相关流场特征,验证了微吹气控制减阻的有效性,局部最大减阻率达到了45%,并且由于微吹气控制的"记忆"功能,减阻效果在微吹气流域下游仍会持续一段距离,增加了减阻区域的流向面积。壁湍流摩擦减阻的原因在于近壁区域出现了一个低速的"湍流斑",黏性底层厚度增加,速度型曲线被抬升。但与此同时,边界层内湍流速度脉动也得到了增强。进一步对流向脉动涡演化规律分析,发现微吹气对流向脉动涡发挥着多重作用。在增加流向脉动涡强度的同时,还使得流向涡团向远离壁面抬升,这样减小了流向涡与壁面之间直接作用。此外,微吹射流产生的冲击作用会在流向涡表面留下凹痕,使得流向涡分散成相对小的涡团结构。     

分区LES/DES混合方法及缝翼三维流场模拟

提出一种分区大涡/脱体涡模拟混合方法,采用高阶空间有限差分方法和二阶隐式LU-SGS时间推进方法,对有限展长缝翼的三维流场进行了数值计算。该混合方法吸收了脱体涡模拟方法和大涡模拟方法各自的优点,在非核心区域使用脱体涡方法进行计算,相对于大涡模拟方法具有较高的计算效率;相对于分区脱体涡模拟中采用的雷诺平均方法,在计算量没有明显增加的条件下提高了对流场各尺度流动的模拟能力。计算结果与风洞实验数据吻合良好,同时缝翼流场主要观测点的压力脉动数据为后续缝翼噪声分析和低噪声优化提供了基础。     

大型双曲冷却塔表面脉动风压随机特性——风压极值探讨

针对风压信号呈现的非高斯特性和传统极限估计方法的局限,首次提出基于保证率和相关性的极值估计方法——全概率迭代法进行冷却塔表面脉动风荷载极值分析,并和传统的峰值因子法及改进的Sadek-Simiu法计算结果进行对比验证。结果表明：全概率迭代法避开了对随机过程的高斯分布假定,相比传统的极值估计方法其结果更加真实可靠;表达风压极值中脉动分量的峰值因子数值沿着环向和子午向变化显著,如取为同一数值则偏于危险或过于保守;采用全概率迭代法得到的表面风压系数极值分布曲线与规范取值相比,迎风面和负压峰值区域极值偏小,背风区域极值偏大,且最小值对应角度相差约10°。     

旋流预混燃烧室流动混合的大涡模拟

采用大涡模拟的方法研究旋流预混燃烧室中的流动及燃料空气混合过程,选取亚格子湍动能模式为基准的涡黏模型,通过模拟涡核的运动过程定性的表征了拟序结构的发展;引入混合不均匀度的概率密度函数,研究入口边界条件尤其是入口速度对于涡团破碎及混合程度的影响规律;通过研究速度分量与燃料质量分数脉动的相关性,分析了不同方向的速度脉动对于混合过程的影响因子.这些均为进一步研究旋流预混燃烧室的燃烧不稳定性及污染物排放提供了理论上的指导.      

航空发动机进气总压畸变地面试验数据处理方法综述

在航空发动机进气总压畸变试验研究中,带有畸变发生器的地面模拟试验在发动机稳定性评定中具有重要作用。试验中通常采用多孔探针/探针梳/探针耙/PIV等测试方法对AIP及其他截面的流场参数进行测试。为了准确描述进气总压畸变的稳/动态特性,需要采用合适的数据处理方法对试验数据进行分析。基于国内外关于进气总压畸变地面模拟试验数据处理方法的研究现状,总结了试验中涉及的相关数据处理方法,分析了不同畸变指标、稳定性评定参数的特点,总结了滤波截止频率选取的方法及原则、相似原理在缩比模型试验中的应用以及总压畸变图谱的类型,对比分析了功率谱密度函数、相关性分析、幅值概率密度函数、小波分析、神经网络、POD分析及DMD分析等统计分析方法的应用特点。最后对航空发动机进气总压畸变地面试验数据处理方法的发展趋势进行了展望。     

三维进气道湍流流场数值模拟

结合实验数据,采用有限容积法,对超音速S形进气道三维湍流流场进行了数值模拟,分析了0°攻角、10°攻角、10°侧滑角三种工况下机头激波对入口气流的影响、进气道入口激波分布和进气道出口压力分布情况.计算结果表明:0°攻角和10°攻角时,机头激波对入口气流影响不大,出口气流总压分布相对均匀,且高压区面积较大;而10°侧滑角时,受机头激波的影响,进气道前方出现小范围的低压区,且出口气流低压区较大.     

直升机附件舱温度场的数值模拟

为考查某型直升机附件舱的热控制能力,采用有限元软件分析了其在强制冷却和自然冷却情况下的温度场分布.应用非结构化网格和有限体积法进行计算区域和控制方程的离散,采用标准k-ε模型封闭湍流控制方程.针对3种冷却系统设计方案,分析了不同冷却气流进口尺寸、相对位置以及不同附件散热量对附件舱温度分布的影响,计算了自然冷却情况下附件舱密闭时的温度场.计算和试验结果均表明,所采用的冷却方案可以满足附件舱冷却的需要.      

发动机喘振-爆燃故障的信号分析

本文着重介绍某双轴发动机喘振—爆燃故障用 FFT信号处理机进行诊断和信号分析的方法。文章结合发动机在湍流动态畸变流场下的喘振—爆燃试验研究,介绍了试验设备、测量系统和对发动机“故障”信号分析的结论。湍流型动态畸变流场下的压气机喘振多为“漂移型”喘振,表现有随机性。发动机对动态畸变流场的敏感度远高于稳态畸变。利用 FFT信号分析技术进行故障诊断,有广泛应用价值。它较其它方法精确、迅速、简便。      

二维进气道不启动流场非定常特性的混合LES/RANS模拟

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)方程模拟方法结合5阶WENO(weighted essentially non-oscillatory)格式对马赫数为3的来流中、内收缩比为1.5的不启动状态下的二维进气道进行了计算,再现了不启动进气道中的非定常流场.计算结果表明:所采用的模拟方法对入口处的平均绝热壁温、摩擦速度和雷诺应力的计算精度较好,进气道不启动流场中激波波系和分离区存在大时空尺度的低频运动,其占主导的特征频率和典型的激波/湍流边界层干扰问题中激波和分离区的低频频率接近,且进气道出现了间歇性的启动状态.      

气膜出流对叶片内表面换热系数影响的实验研究和计算

本文用实验和数值计算的方法研究了气膜出流对通道内壁换热系数分布的影响,研究的重点是在三维湍流状况下气膜孔间距比的变化和进气雷诺数等因素的影响规律。研究结果表明:在所研究的范围内,湍流对流换热系数比层流换热系数有数量级上的增加,换热系数随着雷诺数的增大而增大,孔间距比的变化对换热系数的影响不是单调的,内壁的换热特性在孔的间距比变化时并没有明显的变化。      

Effect of drag reducing riblet surface on coherent structure in turbulent boundary layer

The characteristics of turbulent boundary layer over streamwise aligned drag reducing riblet surface under zero-pressure gradient are investigated using particle image velocimetry. The formation and distribution of large-scale coherent structures and their effect on momentum partition are analyzed using two-point correlation and probability density function. Compared with smooth surface, the streamwise riblets reduce the friction velocity and Reynolds stress in the turbulent boundary layer, indicating the drag reduction effect. Strong correlation has been found between the occurrence of hairpin vortices and the momentum distribution. The number and streamwise length scale of hairpin vortices decrease over streamwise riblet surface. The correlation between number of uniform momentum zones and Reynolds number remains the same as smooth surface.