详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

SAS与PANS模型在圆柱绕流中的应用比较

分别运用SST(Shear Stress Transport)、SST-SAS(Scale-Adaptive Simulation)、两种变fk(模化湍动能的比例)函数的SST-PANS(Partially Averaged Navier-Stokes)湍流模型对Re=3900的圆柱绕流进行了数值研究,重点从湍流结构捕捉、气动力计算、涡黏性控制等方面,比较了SAS与PANS两类RANS/LES混合模型的计算能力,并通过不同网格计算分析了模型的网格敏感性。数值结果表明:SAS及两种变fk方法的PANS模型均具有求解小尺度涡运动的能力,并能较好地反映出绕流尾迹的三维非定常特性,同时PANS模型能捕捉到更多的非定常结构;SAS模型中自适应尺度Lvk立足于当地流动,对网格依赖较小,计算的湍动粘度分布更合理,能够更好地计算剪切层及回流区;两种PANS模型网格独立性较差,出现了雷诺应力不足的现象;类DES可变fk函数构造相对简单,所得fk分布更准确,使用tanh函数计算的尾迹区fk值偏低,对流场调控能力稍差。     

基于等离子体激励的湍流边界层减阻控制

为了探究介质阻挡放电（dielectric barrier discharge, DBD）等离子体气动激励对平板湍流边界层的减阻情况，在控制来流速度为10.7 m/s的低速风洞中进行等离子体平板湍流边界层减阻控制实验。本实验重点研究了不同激励频率对湍流边界层的减阻控制效果，使用热线风速仪系统采集流向速度信号，获得边界层平均速度分布和脉动速度分布。对实验结果进行对比分析发现，在施加不同频率的等离子体激励之后，边界层内对数区速度明显减小；随着激励频率的增加，局部减阻率呈现出先增大后减小的趋势，在激励频率为200 Hz时，减阻率达到最大为7.4%。     

场协同原理在高超声速化学非平衡流动中的推广

应用理论分析与数值模拟方法，将对流传热的场协同原理从不可压缩流动推广至高超声速化学非平衡流动中。结果表明，高超声速化学非平衡层流与湍流的热流密度取决于流动的当地单位体积的动量与单位质量总焓梯度的协同。用当地单位体积的动量与单位质量总焓梯度的协同研究高超声速化学非平衡流动的壁面传热问题，对层流流动下的对流传热，不但计及了高超声速化学非平衡流的密度变化对热流密度的影响，而且包括了静焓梯度、压力梯度、边界层内的分子黏性剪切效应对热流密度的作用；对湍流问题，除了上述层流流动各项对热流密度的影响外，还计及了雷诺剪切应力对热流密度的作用。考虑到高超声速化学非平衡流静焓的定义，高超声速化学非平衡层流及湍流的场协同同时计及所有组分的平动能、转动能、振动能及电子能等梯度的贡献。     

基于气动外形优化和主动流动控制的减阻技术

增升和减阻是民用飞机设计的永恒目标。本文对民用飞机减阻技术最新进展及未来发展进行综述,在此基础上,报道了计算所团队通过数值优化设计及主动流动控制实现减阻的研究进展。为此,发展了基于壁面模化大涡模拟(WMLES)准确模拟湍流边界层的先进数值方法,以及利用伴随方程方法实现气动外形快速数值优化的技术。在具体实践中,通过组合基于自由变形(FFD)的几何参数化模块、基于线弹性体法的网格变形模块,以及基于伴随方程方法和序列二次规划(SQP)的优化算法模块,搭建了面向工程、适于超大规模变量优化设计的整机气动优化设计平台,在NASACRM模型上实现了设计变量超过600的大变量、跨声速气动数值优化设计,在合理的工程约束条件下,有效削弱了机翼表面激波,总体减阻率达到2%以上;通过对湍流平板边界层外层结构进行射流非定常控制,在中等雷诺数Reτ=4700条件下,实现了5%～6%的当地减阻率。研究证实,数值优化设计和基于外层射流控制的主动流动控制技术将是大型民用飞机减阻中非常有前景和值得重视的两种方法。     

预混湍流火焰面褶皱结构网络拓扑研究

湍流火焰结构是表征湍流与火焰相互作用的组分、速度、温度等标量场信息，理解湍流与火焰相互作用规律，验证和发展湍流燃烧模型的实验基础。针对传统曲率PDF分布反映湍流火焰面褶皱结构失准问题，利用网络拓扑结构方法可以标记系统关键节点和特征结构，构建湍流火焰面的拓扑结构。本文标记了湍流火焰面上的关键褶皱结构，分析了湍流与火焰的作用规律，结果表明:低湍流强度下，湍流火焰面的关键褶皱结构由火焰自身不稳定性引起；当湍流强度增大，湍流火焰面的关键褶皱结构由湍流尺度决定。在本生灯湍流火焰中，火焰自身不稳定性引起的火焰褶皱与火焰发展距离有关。在本生灯火焰底部，火焰自身不稳定性不引起火焰面褶皱，随着火焰向下游发展，其对火焰面影响逐渐增大，火焰褶皱程度增加。     

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

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

基于湍动能输运的一方程转捩模型

传统转捩模型的构造十分依赖转捩机理的发现和理解,对机理的模化也存在误差。本文基于尺度自适应的KDO湍流模型,以输运变量r=μt/μ重新标定模型参数,使模型整体达到了流动结构自适应,可实现湍流/转捩一体化计算。根据长度尺度的不同,可细分为KDO-tran和CKDO-tran模型。CFD计算评估了经典的T3A,T3B平板边界层旁路转捩,T3A-平板边界层自然转捩,Aero-A翼型分离泡转捩,DLR-F5,6:1椭球横流转捩及可压缩平板边界层、尖锥转捩。该模型无法精确捕捉层流-湍流的过渡,但能准确捕捉转捩的起始位置(Transition onset),在高Re数、复杂流动的表现尤其突出。该模型捕捉转捩的机制在于流动结构的自适应和湍流输运特性的保存,不引入任何转捩机理却能捕捉多种类型的转捩现象,具有较好的适用性。     

小型回流燃烧室多工况碳黑分布预测

为了研究小型回流燃烧室不同条件下的工作状态,并分析其产生积碳现象的原因,采用k-ε湍流模型、EDC湍流燃烧等模型对该燃烧室多个工况进行数值计算。结果表明,燃烧室性能参数计算值与台架数据较吻合,各个状态下参数相对误差均小于1.37%。根据各工况的温度分布的变化,发现燃烧室火焰随工况降低逐渐向火焰筒头部收缩。计算了各工况蒸发管内的燃油蒸发率,发现蒸发率在低工况时明显降低。根据出口温度分布情况,计算出口温度分布系数OTDF为0.35。分析碳黑粒子的分布,发现积碳现象主要在火焰筒头部,其原因主要是燃烧不充分和冷却气膜的缺失。     

关于约束大涡模拟方法的一些思考

高雷诺数壁湍流是工程设计和应用中非常重要的问题之一,其高效高精度的模拟方法一直是湍流研究的重要研究方向。约束大涡模拟方法(CLES)是近些年提出的新模拟方法之一,和传统的RANS/LES混合方法不同,CLES在全场做LES计算:在靠近壁面的内区,它采用带约束的亚网格模型,而在外区,它使用传统的亚网格模型。经过近10年的努力,CLES方法已经在不可压/可压缩附着流、不可压缩/可压缩分离流等经典算例中得到了验证,并成功应用于航空气动中复杂流动的模拟。本文在介绍CLES方法基本原理的基础上,对CLES方法应用中的一些问题进行了讨论,最后对CLES方法的未来研究方向也做了一些概括。