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

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

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

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

涡轮导向叶片综合冷却效率实验研究

为了研究流动参数对涡轮导向叶片综合冷却效率的影响,采用红外热像仪对叶片表面的温度分布进行了测量,得到了叶片的综合冷却效率随流量比、温比、主流进口雷诺数和湍流度的变化规律。实验过程中,次流与主流的流量比分别为0.15,0.18,0.20,0.22和0.24;主次流温比分别为1.4,1.7,1.93和2.2;主流进口雷诺数分别为1.0×10~5,1.1×10~5,1.2×10~5,1.3×10~5和1.4×10~5;主流进口湍流度分别为0.506%,8.156%,14.92%。结果表明,综合冷却效率在前缘处最低,沿流向逐渐升高;增大流量比会显著提高叶片的综合冷却效率,在温比为1.93时,流量比由0.15增大至0.24,综合冷却效率平均增加29.3%;温比和主流进口湍流度的增大均不利于综合冷却效率的提升,流量比为0.20时,温比由1.4增大至2.2,综合冷却效率平均下降46.5%,湍流度由0.506%增大至14.92%,综合冷却效率平均降低15.5%;主流进口雷诺数对叶片综合冷却效率的影响很小。     

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

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

湍流边界层外区超大尺度相干结构相位平均波形

用IFA300恒温式热线风速仪和X形二分量热线探针精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、法向速度分量的时间序列信号。通过对不同速度分量进行子波变换，得到了近壁区和外边界区流向和法向脉动速度分量的分尺度湍动能随尺度的分布，在外边界区发现流向脉动速度和法向脉动速度的能量最大尺度随法向坐标远离壁面增大，相干结构的平均猝发周期也相应增加，说明在外区存在超大尺度的相干结构。采用条件采样和相位平均技术，提取了近壁区和外边界区能量最大尺度相干结构流向及法向速度的喷射和扫掠的条件相位平均波形，发现流向和法向脉动速度及雷诺应力喷射和扫掠事件的条件相位波形较近壁区都发生了改变。雷诺应力负的幅值显著降低，雷诺应力表现为正负交替的波动现象。     

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

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

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

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

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

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

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

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

超声速混合层中PIV粒子的湍流变动作用研究

对二维超声速气固两相混合层进行双向耦合，研究了粒子图像测速技术（PIV）中示踪粒子对超声速混合层的湍流变动作用。超声速气固两相混合层的气相采用大涡模拟，离散相采用拉格朗日颗粒轨道模型求解。结果表明:与无负载示踪粒子时的超声速混合层相比，小Stokes数示踪粒子在超声速混合层中的布撒减弱了流向湍流，而强化了法向湍流，使雷诺应力峰值增大了9.68%；大Stokes数示踪粒子对混合层的湍流脉动起到了一定的削弱作用，最大雷诺应力值只有无负载时的41.74%。大质量载荷时，大量示踪粒子的运动尾迹抹平了部分法向速度脉动，使最大法向速度脉动只有无负载粒子时的38.63%；中等质量载荷时，超声速混合层的法向速度脉动和雷诺应力峰值与无负载粒子时相近；而小质量载荷时，超声速混合层中心线及其附近的法向速度脉动得到较小的增强，而最大流向速度脉动却被削弱了19.29%。小Stokes数和中等质量载荷示踪粒子对原始无负载粒子时的流场影响相对较小，研究结论对高速流动PIV测试有着重要的参考价值。