超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

数值模拟方法对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高度相关。     

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

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

圆形和矩形燃气射流在受限液体中扩展特性的对比研究

为研究不同喷孔结构的燃气射流在受限液体工质空间的扩展特性,采用VOF(Volume of fluids)模型分别对圆柱形充液室中圆形射流和矩形射流扩展过程进行气水耦合数值求解。考虑了燃气可压缩性、流体粘性和热量交换因素对射流扩展特性的影响。通过数值模拟,获得圆形射流扩展形态及轴向速度计算值,与实验结果吻合较好。在此基础上,进一步对比分析了矩形射流和圆形射流的扩展特性,获得射流场中密度、压力、温度和速度的分布图以及涡的演化过程。计算结果表明:燃气自喷孔喷出,由于出口处压力较大,高温燃气继续膨胀产生膨胀波,膨胀波在气液界面反射形成压缩波,圆形射流的膨胀波和压缩波均比矩形射流强。燃气射流扩展过程中,由于周围液体惯性效应,气液卷吸掺混效应以及膨胀压缩波作用,圆形射流的压力场和温度场分布较矩形射流更加复杂,同时刻的圆形射流轴向扩展速度也比矩形射流小。从速度云图发现矩形射流的速度核心区比圆形射流速度核心区要短;同时从圆形射流和矩形射流不同截面流向速度分布图发现流向速度峰值均发生了偏移,但矩形射流发生偏移的区间更长,偏离距离更远。     

磁流体湍流中的结构和能量传输

磁流体湍流中存在多种相干结构,包括电流和涡旋结构。文章主要对近几年有关磁流体湍流相干结构及其相关的能量传输的一些工作进行了介绍。不同于中性流体湍流中的管状结构,磁流体中的相干结构多呈现为片状,强电流片附近也存在强涡片结构。磁流体湍流中的能量级串使能量跨尺度从大尺度传递至小尺度;与中性流体湍流能量传输主要集中在相近尺度上不同,磁流体湍流除了动能传输,还有磁能传输,以及动能与磁能之间的跨尺度互相转换。研究表明,磁流体湍流中的能量传输及耗散具有强间歇性,集中在较小区域范围,且与相干结构存在关联。在存在背景磁场的情况下,磁流体湍流中的相干结构在背景磁场方向上被拉长,而同样方向上的能量传输受到抑制。在高马赫数的情况下,激波的产生会使能量传输明显增强,同时带来很强的间歇性,结构函数标度律远远偏离线性标度关系,并且随着阶数的增加,标度指数会达到饱和值。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

空心气冷低压涡轮动叶气热耦合数值模拟

以空心气冷低压涡轮动叶为研究对象,采用高质量的流体域和固体域网格控制技术,带转捩模型的双方程SST湍流模型,开展了基于CFD方法的叶片气热耦合问题研究。获得了不同冷气流量比(分别为1.0%,1.38%,1.8%和2.2%)、温比(分别为2.1,2.25,2.3,2.4和2.5)和压比(分别为1.4,1.6和1.8)对叶片换热特性的影响规律,设计状态中截面按弧长平均的叶片壁面金属温度计算值较试验值偏小0.3%,气热耦合计算的叶片壁面温度分布与试验结果吻合良好,验证了气热强耦合计算方法的精度,为涡轮叶片温度场分析提供了一种有效的方法。