高强驻波声场中的湍流

反射端声压级高于140dB以后驻波管中的声流结构会发生明显变化。利用流动显示和LDV（Laser Doppler Velocimetry）技术对该问题进行了实验研究。流动显示结果表明，在160dB左右的高强驻波声场中随着旋涡结构的合并、破裂等现象的出现，声致流动有明显的湍流特性。LDV测量结果显示，随着驻波管中声压级的提高，一些测点的速度频谱明显加宽，出现高频速度脉动，而另一些点的速度频谱则没有显著变化。径向湍流强度测量结果显示，此时管中流动已是复杂的三维结构,而非轴对称的二维结构。     

倒置喷管湍流流场计算

本文对倒置喷管的湍流跨音速流场进行了数字模拟，其算法基于ＹａｎＬｅｅｒ的矢通量分裂格式，求解时采用了线性Ｇａｕｓｓ—Ｓｅｉｄｅｌ迭代法，并结合局部时间步长技术对收敛过程进行加速，对雷诺数为１．６６×１０６的倒置喷管的流场进行了数值计算，结果与实验相符合。     

湍流燃烧模拟中化学反应的加速算法研究进展

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法，讨论了动态自适应化学（Dynamic Adaptive Chemistry，DAC）方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中，使用DAC简化可以加速化学反应计算。然而，在并行燃烧数值模拟中，处理器核心的负载极度不平衡，加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心，更有利于整体计算效率的提高。在同等精度下，相比于隐式格式耦合DAC和MTS加速方法，Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。     

数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。     

高雷诺数湍流非预混火焰及NO生成的大涡模拟

采用两种亚网格湍流燃烧模型,即化学建表方法结合假定概率密度模型和稀疏拉格朗日过滤密度函数方法,对高雷诺数湍流非预混火焰Flame D进行数值研究,定量比较不同亚网格模型的差异,并对火焰特征、污染生物生成特性进行分析。结果表明,两类亚网格燃烧模型预测的温度及大组分分布相近,稀疏拉格朗日过滤密度函数方法可以更好地模拟CO质量分数分布。不同的假定概率分布均可合理描述湍流与火焰的相互作用,之间的差别主要体现在NO分布,Dirac函数远高估了NO生成,而Top-hat函数则略低估了NO生成,Beta函数表现最优。Flame D的高温区及NO质量分数均主要分布在当量混合线及富燃侧附近。受高温伴流的影响,NO质量分数与温度一直保持高度正相关,峰值主要集中在标量耗散率很小的区域。不同截面上,反应物中的O_2和生成物中的H_2O均与NO高度相关。     

航空弧齿锥齿轮副风阻功率损失分析与优化

基于CFD理论，利用Fluent求解软件，借助超级计算机强大的并行运算能力对航空弧齿锥齿轮副风阻功率损失进行仿真计算。采用局部综合法建立弧齿锥齿轮副三维模型，选用RNG k-ε湍流模型，考虑平均流动中的旋流流动情况，与标准k-ε模型相比，RNG通过修正湍流黏度并很好地处理了高应变率以及流线弯曲程度较大的流动。齿轮边界运动通过UDF(user-defined functions)函数驱动，同时采用动网格模拟流场形状由于边界运动而随时间改变问题。最后得出无挡风罩和不同挡风罩配置下的齿轮副风阻功率损失，证实了合理安装挡风罩能够有效降低齿轮风阻损失，并分析多组仿真实验间的减速器内流场压力、速度、湍流动能云图变化，得出了最优化的挡风罩配置，以求最小化风阻功率损失，文中减阻效果最好的挡风罩能降低55.3%的齿轮风阻损失，此时挡风罩间隙为1 mm，为工程实际应用挡风罩的设计提供了参考。     

高主流湍流度下密度比对涡轮导叶全气膜冷却特性的影响

为了研究高主流湍流度下二次流密度比对涡轮导叶全气膜冷却特性的影响,使用热色液晶测量了在主流湍流度为15%,二次流密度比为1.0和1.5下三维涡轮导叶的气膜冷却效率和换热系数。二次流与主流质量流量比为7.0%和12.5%。结果表明:二次流密度比增大可以降低冷气射流的动量,小流量比工况下,在叶片前缘和压力面前半段,动量较低的二次流在高主流湍流度的影响下更易耗散,增大二次流密度比使冷却效率明显降低;大流量比工况下,二次流动量降低使气膜孔后区域冷气贴附性增强,气膜冷却效率和冷气覆盖效果均得到提升。小流量比工况下,二次流密度比增大对叶片表面换热的影响较小;大流量比工况下,二次流密度比增大使吸力面中弦区域和压力面后半段的平均换热系数比分别降低15%和25%。     

可控流场尺度预混湍流燃烧器及其火焰结构分析

为了研究单一湍流场参数对预混湍流火焰结构的影响，以及拓宽湍流场的强度和尺度范围，发展了一套可变结构的预混湍流燃烧器。采用恒温型热线风速仪标定流场，得到了一系列湍流参数。流场标定结果表明:该燃烧器能显著拓宽湍流强度和尺度范围，并能利用不同几何结构产生多种可控流场，实现研究单一湍流参数对湍流燃烧速度和火焰结构影响的目的。选用有代表性的15种湍流孔板组合结构，利用OH-PLIF燃烧激光诊断技术，开展了湍流燃烧实验，结果表明:湍流强度的增大（1 < u'/S_L，0 < 10）使得湍流火焰分区扩展到了薄层反应区，火焰面破碎程度明显增强，孤岛结构明显增多。高宏观雷诺数下，积分尺度的增长对湍流燃烧速度起抑制作用，可能存在临界宏观雷诺数Re_c，能够表现流体惯性力占主导地位的程度，决定积分尺度对湍流燃烧速度的影响效果。积分尺度能量大，扰动能力强，故积分尺度越大，火焰体积越大；但过高的湍流强度会使火焰面褶皱更加剧烈，小尺度叠加在大尺度上的程度增强，最终也使火焰体积显著增大，掩盖了积分尺度对火焰体积的影响，说明积分尺度（表征大尺度）不如湍流强度（表征叠加小尺度的程度）对火焰放热率影响大。     

均匀各向同性大气湍流对声爆传播特性的影响

基于离散Fourier模态有限和生成的随机大气湍流场，采用修正波形参数方法，开展了均匀各向同性大气湍流对典型超声速客机声爆传播特性的影响分析。计算采用的超声速客机模型为自行生成的简化超声速公务机模型。首先，应用航空工业空气动力研究院自主研制的CFD软件ARI_Overset在三维空间求解Navier-Stokes方程，得到作为声爆远场传播初始值的近场空间压力分布；其次，基于Von Karman能量谱，采用离散Fourier模态有限和形式生成随机均匀各向同性大气湍流场；最后，采用修正波形参数方法模拟了声爆信号在随机速度湍流场中的传播过程。数值结果表明：各向同性大气湍流对于地面声爆特征有重要影响。与无湍流状态相比，均匀各向同性大气湍流使得地面声爆特征增强的概率约为55%，使得地面声爆特征减弱的概率约为45%，故总体而言大气湍流效应更倾向于增强地面声爆特征；均匀各向同性大气湍流对于声爆传播路径影响相对较小，但是这一变化仍会导致地面信号接收点的不确定性。