基于S-A湍流模型和间歇因子输运方程的转捩流数值模拟

考虑到γ-Re_θt转捩模型中间隙因子输运方程并不依赖于具体的湍流模型,因而提出了耦合S-A湍流模型和间歇因子输运方程的转捩流模拟方法.其主要思想是构造与间歇因子相关的两个耦合函数,分别作用到S-A模型的生成项和耗散项用以控制湍流的产生与发展,从而实现层流到湍流的数值转捩.结果表明:该方法中模型方程具有完全当地性的特点,易于实施,可直接用以求解三维转捩流动.4个典型算例表明该方法模拟结果与γ-Re_θt模型准确度相当,但少求解了两个输运方程,计算耗时减少了15%.      

真实气体效应对双楔绕流流动特性的影响

为研究高马赫数下真实气体效应对双楔绕流流动特性的影响，采用量热完全气体模型和化学非平衡气体模型两种气体模型对层流流态、转捩流态、全湍流流态的双楔绕流流场进行了数值模拟，其中转捩流态计算采用基于间歇因子的γ-Re_θt转捩模型，全湍流计算中采用k-ωSST湍流模型。结果表明：双楔拐角附近发生流动分离，与量热完全气体模型相比，相同流态下，化学非平衡气体模型计算得到的分离区尺寸更小。在边界层转捩位置的分析中发现，双楔拐角附近的分离区对边界层转捩有重要影响，分离激波产生的位置与边界层开始转捩的位置高度吻合，与量热完全气体模型相比，化学非平衡气体模型时的边界层转捩起始位置更靠后。针对真实气体效应下单位雷诺数、拐角角度对边界层转捩的影响规律研究发现，真实气体效应在一定程度上抑制了分离区内的转捩流动，随着单位雷诺数减小、拐角角度增大，真实气体效应对边界层转捩的影响随之增大。     

高超声速三维边界层转捩数值研究进展及预测软件

高超声速边界层转捩问题已经成为飞行器设计需要考虑的一个关键因素，也是高超声速流动研究的热点和难点。三维边界层中存在横流、流向涡、接触线三种典型的流动结构和失稳特征，他们诱导转捩的机理不同。本文总结了低雷诺数湍流、γ-Re_θt、k-ω-γ、层流动能4种常用转捩模型的构造、适用性及应用情况，介绍了适用于三维边界层的BiGlobal、 PSE3D、 TriGlobal三种全局稳定性分析方法，分析了基于转捩预测e^N方法的e Malik、 LASTRAC、 LILO、HyTEN等转捩软件的优缺点，展望了高超声速三维边界层转捩研究及软件开发的未来。     

高雷诺数下跨声速自然层流短舱优化设计

为解决高雷诺数下大涵道比发动机自然层流短舱高维优化设计问题,提取短舱3个基准面实现轴对称自然层流优化设计。通过获得较大范围层流区域,从而降低短舱表面摩擦阻力。采用类别形状函数（CST）参数化、γ-Re_θt转捩模型和遗传算法建立自然层流（NLF）短舱自动优化流程。表明通过优化短舱基准面建立三维NLF短舱设计方法的可行性。进而通过CATIA二次开发构建三维非轴对称NLF短舱,解决了基准面优化后大量数据点的高效导入和曲面生成问题。针对设计的非轴对称NLF短舱,进行了设计点附近迎角、侧滑角、来流马赫数以及湍流度的转捩敏感性分析。结果表明:跨声速状态下,迎角增大层流范围减小;来流马赫数增大层流范围扩大;侧滑角和湍流度对层流范围影响很小。      

γ-Reθt转捩模型在微通道流动特性研究中的应用

应用ANSYS CFX对矩形微通道中的流动进行了数值研究。通过与实验数据的对比检验了层流模型、γ-Reθt转捩模型和剪切应力输运模型在微通道流动转捩模拟中的适用性,并在此基础上探究了微通道长径比以及截面宽高比对临界雷诺数的影响。结果表明:3种模型中,只有γ-Reθt转捩模型捕捉到了转捩的趋势,并较为准确地捕捉到了临界雷诺数。当微通道的长径比不小于200,进口段对整体流动特性产生的影响可以忽略。当微通道的当量直径和长径比一定,宽高比由1增加到5时,临界雷诺数呈逐渐增长趋势。      

粗糙壁面诱导的流动转捩数值模拟方法

为预测粗糙壁面诱导的流动转捩,对Langtry提出的γ-Re_θ转捩模型增加粗糙壁面的比耗散率和涡黏性的边界条件,并对模型中的经验关联函数——转捩动量厚度雷诺数进行修正,引入等效沙粒表面粗糙度作为变量,使模型具有一定的预测粗糙壁面诱导的流动转捩的能力.对粗糙平板自然转捩算例和变压力梯度平板绕流算例进行数值模拟,计算结果与风洞实验数据符合的较好.主要结论如下:表面粗糙度一般会增加边界层内的传热系数和阻力系数,同时使层流到湍流的转捩位置提前;自然转捩的转捩位置受表面粗糙度影响较大,与光滑壁面相比,平板算例中0.15mm的表面粗糙度使转捩位置提前40%;分离诱导转捩的转捩位置受表面粗糙度影响较弱,随着表面粗糙度逐渐增加,转捩位置和分离泡位置略有后移,且分离泡强度逐渐减弱,分离泡之后的阻力系数增加.      

低雷诺数下翼型层流分离泡及吹吸气控制数值研究

在低雷诺数下Eppler387翼型表面会出现层流分离泡现象,因此本文使用Fluent求解器开展吹/吸气控制翼型表面层流分离泡的数值研究,主要探究了射流位置、射流角度、射流速度比三个控制参数对层流分离泡控制的影响规律.研究结果表明,采用与SSTκ-ω湍流模型耦合的γ～(Reθt)转捩模型可以准确预测层流分离泡的位置;吹/吸气可以有效抑制低雷诺数下层流分离泡的发展,明显提高低雷诺数下翼型升阻比.固定射流位置,较大吸气速度比和较小吹气速度比可分别获得较好的流动控制效果,且吸气控制比吹气控制对层流分离泡的抑制作用更加有效.     

固定转捩和自由转捩流动预测方法验证

采用固定转捩位置和基于SST k-ω湍流模型的γ-Reθt转捩模型自由转捩模拟两种方法研究了RAE2822翼型流动转捩问题。先采用γ-Reθt转捩模型确定转捩位置,再采用所获得的转捩位置实施固定转捩位置计算,并与全湍流计算结果和其他软件计算结果进行了对比。表明两种方法均能有效捕捉超临界翼型流动转捩现象,固定转捩位置模拟结果基本上重现了自由转捩预测结果。     

用转捩模型预测转捩及确定最佳粗糙带高度

采用SST k-ω二方程湍流模型,附加γ-Reθ转捩模型对不同迎角下的NACA0012翼型绕流进行了转捩预测,计算预测的转捩位置与实验结果符合较好,说明γ-Reθ转捩模型对转捩位置具有较好的预测能力。进一步用模型对带有粗糙带的翼型进行转捩预测,获得了不同雷诺数下刚好激发转捩又不引起过大局部压力扰动和附加阻力的最低粗糙带高度,对固定转捩风洞实验有指导意义。     

基于SA和SST γ-Reθt模型对边界层转捩预测的比较

利用作者共同研发的In-house代码TRANS3D平台,在SA和SST两种常用湍流模型框架下构建了γ-(Reθt)转捩模型,并以二维低速S809层流翼型和三维小后掠跨音速F5层流机翼为对象,比较了两种不同湍流模型构架下γ-(Reθt)转捩模型预测的气动力特性和流场分布.结果表明:两种转捩模型均能预测航空转捩计算中常见的分离流转捩和自然转捩类型,明显改善了中低雷诺数流场下的预测精度,但由于选取基准湍流模型的不同,基于SA和SST的γ-(Reθt)转捩模型在流动细节上依然存在着一些差异.     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

微小圆管内煤油流动研究

通过实验研究了RP-3航空煤油在微小圆管内的流动情况.分别研究了内径为1.0,0.8 mm和0.5 mm的微小圆管中的流动压降和摩擦阻力系数.实验结果表明:摩擦阻力系数(层流损失系数)随雷诺数的增加而减小,而且微小圆管内的摩擦阻力系数比常规尺寸圆管中摩擦阻力系数小;微小圆管中从层流转变成湍流的临界雷诺数也远小于常规尺寸圆管中的临界雷诺数2 300.此外,进、出口压差随雷诺数的增大而增大;相同雷诺数、相同管长的情况下,管内径越大,进、出口压差越小;管内径越小,进、出口压差越大.      

基于SA和SSTγ-■模型对边界层转捩预测的比较

利用作者共同研发的In-house代码TRANS3D平台,在SA和SST两种常用湍流模型框架下构建了γ-■转捩模型,并以二维低速S809层流翼型和三维小后掠跨音速F5层流机翼为对象,比较了两种不同湍流模型构架下γ-■转捩模型预测的气动力特性和流场分布。结果表明:两种转捩模型均能预测航空转捩计算中常见的分离流转捩和自然转捩类型,明显改善了中低雷诺数流场下的预测精度,但由于选取基准湍流模型的不同,基于SA和SST的γ-■转捩模型在流动细节上依然存在着一些差异。     

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

槽流的亚临界转捩与局地湍流

在中低雷诺数时,处于线性稳定的槽道流若受到强扰动可发生亚临界转捩变为湍流。近年人们对该类转捩机理的研究取得重要进展。对于平板泊肃叶流,亚临界转捩之初是稀疏湍流态,其特征结构是远间隔的包含小尺度涡和高低速条带的大尺度的湍流带,可倾斜伸长。该阶段的湍流占比有上限但并非雷诺数的单值函数。随着雷诺数的增加转捩进入平衡局地湍流态,即存在统计定常态,其湍流占比是雷诺数单值函数,可由定向逾渗模型描述。进一步增加雷诺数,湍流带的分裂愈发频繁,最终流场会布满湍流带,在更高的雷诺数时变为均匀湍流。论文概述了为比拟亚临界转捩过程所提出的动力学模型,以及为定量表征管流、平板库艾特流和平板泊肃叶流转捩过程的相似性所提出的局地稳定性参数,并在最后对槽道流亚临界转捩研究的发展做了简要展望。     

转捩对平流层飞艇外形优化的影响

采用梯度算法优化了平流层飞艇外形以减少气动阻力,目标函数飞艇总阻力系数通过求解Navier-Stokes方程获得.优化过程中采用了SSTκ-ω转捩湍流模型和Realizable κ-ε全湍流模型求解绕流流场,采用大涡模拟方法详细分析了获得的两种飞艇外形的气动性能和特点,并且比较了两种湍流模拟方法的流场求解对梯度法优化的目标函数以及飞艇形状特点的影响.结果表明:是否考虑转捩对优化的结果有明显的影响,考虑转捩优化方案的气动性能明显优于全湍流流动求解优化方案.采用转捩模型的优化由于层流区域的存在改变了飞艇尾部的形状,从而获得了尾部分离区较小的飞艇外形轮廓.