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

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法，讨论了动态自适应化学（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高度相关。     

煤油两相连续旋转爆震燃烧室工作特性试验研究

为了深入分析煤油燃料两相连续旋转爆震燃烧室的工作特性,采用富氧空气或氧气为氧化剂,通过试验得到爆震波的时域、频域特征,对两相连续旋转爆震燃烧室中爆震波的起爆过程和稳定后的传播过程进行研究。利用基于激光散射相位多普勒分析(PDA)技术对雾化流场进行了测量,得到喷注器出口不同平面处煤油液滴速度与直径的统计分布。试验结果表明当煤油流量为78g/s,氧气流量为224.0g/s,空气流量为72.5g/s,当量比为1.083时,燃烧室在单波模态下工作,爆震波传播频率为0.904k Hz,平均转速为649m/s。使用氧气作为氧化剂,当煤油流量为81.8g/s,氧气流量为231.8g/s,当量比为1.222时,燃烧室在双波模态下工作,爆震波传播频率为5.882k Hz,平均转速为1848m/s,传播过程中表现出很强的非定常性。在当量比为0.805~0.908的富氧工况下,随着氧化剂中含氧量的增加,爆震波的速度逐步增大,最终达到2440m/s;在当量比为1.057~1.220的富燃工况下,随含氧量的增加爆震波速度呈现线性增长的特征。     

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

基于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，能够表现流体惯性力占主导地位的程度，决定积分尺度对湍流燃烧速度的影响效果。积分尺度能量大，扰动能力强，故积分尺度越大，火焰体积越大；但过高的湍流强度会使火焰面褶皱更加剧烈，小尺度叠加在大尺度上的程度增强，最终也使火焰体积显著增大，掩盖了积分尺度对火焰体积的影响，说明积分尺度（表征大尺度）不如湍流强度（表征叠加小尺度的程度）对火焰放热率影响大。