基于改进SST模型的分离流动数值模拟

对SST湍流模型进行了改进,并通过对典型分离流动进行数值模拟,来检验和提高模型预测分离流动的能力.基于亚声速分离流动,提出了提高雷诺应力的模拟精度和分离流修正的改进方法,并进行了对比研究得出结论;在亚声速分离流动分析结论基础上,采用了可压缩性修正方法,开展了跨声速、超声速和高超声速激波/边界层干扰分离流动的数值模拟研究.结果表明:提高雷诺应力的模拟精度和采用分离流修正明显地提高了SST湍流模型对分离流动的模拟能力;适当地可压缩性修正对超声速和高超声速分离流动的计算精度有改善作用.     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

跨声速和超声速流中激波/边界层干扰数值模拟

对SST湍流模型中的Bradshaw常数a1进行了修正,并对跨声速和超声速流中激波/边界层干扰进行了数值模拟研究,空间离散采用二阶精度差值的低耗散通量分裂格式(LDFSS),时间离散采用对称高斯-赛德尔(SGS)算法。结果表明:在跨声速流动中,计算得到的壁面压力分布、分离区长度和速度剖面都与实验值吻合较好,而且很好地模拟了典型的λ激波结构;在超声速流动中,修正后模型的计算精度较原始模型有了较大改善,计算得到壁面压力分布和分离点的位置都和实验值吻合较好。     

SST湍流模型在高超声速绕流中的改进

为模拟高超声速湍流问题,对剪切应力输运(SST)湍流模型系数进行了修正。数值格式采用改进的总变差递减(TVD)格式,并对湍流模型的负值强制项进行了隐式处理。在此基础上计算了绕平板以及具有分离、再附、激波/边界层干扰等复杂流动结构的压缩拐角的高超声速流动。计算结果与试验数据及半经验公式的对比表明:SST湍流模型引入的雷诺剪切应力与湍动能之比为常数(Bradshaw数)在高超声速绕流中并不成立。Bradshaw数修正后的SST湍流模型与原模型相比,所计算的壁面压力、摩擦阻力和壁面热流分布更接近试验结果。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

NPLS技术及其在高速飞行器气动研究中的应用

近年来,与高速飞行器相关的超声速/高超声速流动受到了极大关注。这类流动所具有的非定常性、强梯度和可压缩性对试验研究提出了挑战。纳米示踪的平面激光散射技术（NPLS）是2005年由作者所在的研究团队研发的非接触光学测试技术。它能够获得超声速三维流场的某个剖面的瞬态流动结构,并且具有较高的时空分辨率。目前,许多研究结果表明NPLS是研究超声速湍流的一项非常有效的技术。近年来,作者应用 NPLS 技术在超声速湍流研究中取得了较大的进展,并且基于NPLS开发了其它几种技术,比如基于 NPLS 的密度场测量技术（NPLS-DT）,能够获得超声速流动的密度场信息并还能进一步得到雷诺应力分布。本文介绍了NPLS技术并回顾了其在超声速边界层、激波/边界层相互作用等流动中的应用。由于能够获得雷诺压力和湍动能等统计量, NPLS技术有望在发展可压缩湍流模型的研究中发挥作用。     

来流边界层拟序结构对超声速混合层流场结构的影响研究

超声速混合层涉及可压缩湍流的根本问题,具有重要的应用背景。通过设计超声速混合层实验装置、应用新近提出的高分辨率NPLS测试技术,拍摄了来流边界层分别为层流和湍流流态下混合层的流向和展向流动图像。根据流动图像的特征,分析了混合层的流向与展向流场切面中拟序结构的成因;深入讨论了来流边界层中拟序涡结构与混合层涡结构的相互作用问题;比较了层流和湍流来流条件下混合层拟序结构的异同及其对混合效率的影响。结果表明：当来流边界层为湍流时,对应的混合层具有较高的混合效率。     

基于SST k-ω的超声速气膜冷却湍流模型温度修正

针对现有湍流模型无法准确预测非等温超声速气膜冷却行为的问题,在现有SST k-ω可压缩修正模型基础上,以总温梯度为变量,完成了湍流模型的温度修正,并首先通过非等温可压缩自由剪切流动实验数据初步验证了其修正效果,在此基础上对温度修正模型预测超声速气膜冷却传热的准确性进行了验证。结果表明,基于剪切层总温变化的湍流模型修正效果显著,可准确预测大温度梯度下的自由剪切流动轴向速度分布;修正模型计算得到超声速气膜冷却壁面热流分布与对应的实验结果吻合;当用于剪切层温度大梯度变化的超声速气膜冷却数值模拟时,温度修正后的SST k-ω模型与可压缩修正的k-ω模型、SST k-ω模型相比,具有显著优越性。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用。与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要。因此,基于两步四阶时间离散和高精度WENO重构,发展了具有四阶时间精度的气体动理学格式。该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力。目前,两步四阶格式已经成功地应用到低雷诺数湍流直接数值模拟和高雷诺数工程湍流RANS模拟中,包括低速槽道湍流、超声速均匀各向同性衰减湍流、二维亚声速翼型湍流和三维跨声速翼身湍流等。数值结果表明该格式对湍流直接数值模拟和湍流RANS模拟具有高数值精度和高数值稳定性。下一步将利用高阶气体动理学格式研究更具有挑战性的可压缩湍流问题,例如超声速湍流边界层和激波边界相互作用等。     

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

高超声速边界层转捩的若干问题及工程应用研究进展综述

高超声速飞行器绕流存在着激波、边界层、流动分离、稀薄气体效应和高温气体效应等多种复杂流动现象的空气动力学问题,其中高超声速边界层转捩既是空气动力学的基础问题,也是高超声速流动研究的热点和难点。若能对边界层转捩进行准确预示及有效控制,则可以实现对飞行器气动力热特性的精细设计,改进飞行器性能,提高任务执行能力。文章针对工程中具有复杂外形飞行器存在的典型失稳特征进行了研究进展回顾,提出了工程实际中亟需解决的复杂边界层转捩问题,明确了高超声速边界层转捩研究的工程应用方向。文章最后还对高超声速边界层的流动控制进行了回顾,以期在今后高超声速飞行器设计中实现对边界层的流动控制,提高飞行器的飞行性能。     

几种常用格式和湍流模型在高速流动中的适用性研究

针对高超声速飞行器中存在的强激波、激波/边界层干扰、分离、湍流等复杂流动现象,对比分析了当前计算流体力学中的主要空间离散格式及湍流模型,发现不同格式对强激波的分辨率基本相同,Roe和LDE(low diffusion E-CUSP(convective upwind split pressure))格式对摩擦因数和传热系数的模拟优于其他格式;S-A(Spalart-Allmaras)一方程湍流模型计算的摩擦因数比k-ωSST(shear stress transport)两方程湍流模型高10%左右,而后者预测的分离区约为前者的2倍,且分离点靠前.     

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟，评估了Spalart-Allmaras， 标准k-ε， RNG (re-normalization group) k-ε，realizable k-ε，标准k-ω，SST(shear stress transport) k-ω，v2-f，Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确，其中SST k-ω模型比v2-f更能准确地预测壁面压力，然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

激波/平板湍流边界层干扰的数值模拟

采用理性GAO-YONG可压缩湍流模型,模拟了激波/平板湍流边界层干扰现象,结果表明:计算所得到的壁面压力分布、摩阻系数分布和速度型分布均与实验值吻合很好,并且比较准确地预报出了入射斜激波与平板湍流边界层干扰所引起的边界层分离点和再附点等流动特性.由此表明GAO-YONG可压缩湍流模型能够准确地用来模拟激波/平板湍流边界层的干扰流动.      

高超声速边界层感受性问题高精度数值模拟

通过五阶WENO格式和六阶对称紧致格式以及三阶TVD R-K结合的方法,对存在强激波和小扰动相互干扰的高超声速边界层感受性问题进行研究.结果表明:此方法能够模拟边界层内不稳定波的产生和发展以及小扰动和激波相互干扰等现象,因此能广泛应用于含激波的感受性等问题的可压缩湍流直接数值模拟以及具有激波和边界层干扰等复杂流场的计算.     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

超声速压缩拐角激波/边界层干扰动力学模态分解

压缩拐角激波与边界层干扰问题广泛存在于高速飞行器的外部和内部流动中,其非定常复杂流场结构对飞行器气动性能影响显著。动力学模态分析将有助于进一步加深理解激波与边界层干扰流场不同特征频率对应的流动结构及动力学特性,为揭示其复杂流动机理提供参考。本文采用动态模态分解(DMD)方法对来流马赫数为2.9、24°压缩拐角内激波与超声速边界层干扰下的非定常流动进行了模态分析。评估了稀疏改进动态模态分解方法在压缩拐角流动中的适用性,研究了湍流干扰和转捩干扰下典型特征频率对应的动力学模态空间结构差异及其原因,分析了转捩边界层展向非均匀性对低频/高频模态动力学机制的影响规律。研究发现,湍流干扰与转捩干扰下拐角干扰区内均存在两类截然不同的动力学模态:低频模态和高频模态。低频模态结构集中在分离激波及分离泡剪切层的根部,表征为分离泡的大尺度膨胀和收缩运动;高频模态空间分布则以平均声速线附近正负交替结构为主,对应为边界层内不稳定波沿剪切层往下游的传播。转捩边界层的展向结构对低频模态运动特性影响明显,而对高频模态的影响则相对较小。     

涡轮叶栅超声速流场流动特征与气膜冷却特性

应用shear strain transport(SST) k-ω 两方程湍流模型,对超声速涡轮叶栅通道内气膜冷却特性进行数值研究,得到不同气膜孔倾角和吹风比下叶栅通道内流场流动特征以及气膜冷却效率的变化规律.在激波入射点附近的气膜射流能够向分离区边界层中补充动量,克服逆压力梯度,有效改善由于激波引起的局部过热.亚声速流动状态下的气膜入射角度对冷却效率的影响能够在较大吹风比下得以体现,而超声速主流状态下,气膜冷却效率与入射角度基本无关,说明亚声速的气膜冷却射流对超声速主流的穿透力要弱于对亚声速主流的穿透力;超声速主流条件下,在激波入射位置的气膜冷却效率要高于激波入射位置下游的气膜冷却效率,这与气膜孔出流在当地的湍流度有关.