弯曲管道内湍流流动的数值模拟

为了对弯曲管道内的湍流流动特性进行有效的预测，用显示代数应力模式对－典型的弯曲管道内流动进行了数值模拟和分析，计算结果表明：（1）显式代数应力模式能对弯曲管道内流动的平均量（如速度、压力、摩阻等）进行较为准确的预测；（2）对模式中的ε方程做适当修正后，该模式能较好地模拟出流动曲率效尖对湍流特性的影响，是对雷诺应力模式的一个比较好的近似。此外，该模式在形式上与标准k－ε两方程模式类似，使用方便，具有较强的工程实用价值。     

用非线性涡粘性模式计算三维湍流边界层

针对非线性涡粘性模式在求解三维湍流边界层流动时的不足，本文从压力－变形率关联项中的快速项出发，考虑流动非均匀影响，在显式代数应力模式的雷诺应力表达式中引入了反映雷诺应力“松驰”效应的速度二阶导数项，构造了一个新的非线性涡粘性模式。通过对典型算例－翼体角隔流动的计算结果表明，新模式能较好地再现出三维湍流边界层内雷诺切应力方向的发展滞后于速度梯度方向发展这一流动特性。     

修正曲率及雷诺正应力对分离流的影响

应用Ｓｔｒａｗｎ＆．Ｋｌｉｎｅ的准同时粘性／无粘流匹配模型，于动量积分布方程中加入了曲率和雷诺正应力两个二阶修正项，计算了二维不可压湍流附面层分离流地。文中首先分析推导了新的曲率修正公式，给出了两个算例，比较了修正项对计算结果的影响，并与实验做了对比。     

两种雷诺应力湍流模式的比较及其在涡流计算中的应用

本文将要讨论的两种雷诺应力湍流模式为:(1)雷诺应力代数模式;(2)雷诺应力输运模式。这两种模式分别应用于弱强涡射流的计算中,其结果表明,代数模式并不能很好地近似输运模式,特别是在一些重要的流场区域,如对称轴附近,前者给出非真实的物理结果。本文将探讨代数模式的这一缺陷,同时提出一新的输运模式。     

湍流应力模式对计算有分离三维边界层的效果研究

为计算有分离三维边界层的特性，探索适用的应力模式，对无限展长的后掠翼（Ｘ＝３５°），选择了几种湍流应力模式进行计算，并对结果作了分析，所采用应力模式有代数模式，它包括常用的Ｃｅｂｅｃｉ－Ｓｍｉｔｈ（Ｃ－Ｓ）模式，Ｂａｌｄｗｉｎ－Ｌｏｍａｘ（Ｂ－Ｌ）模式、修正外区粘性系统的Ｃ－Ｓ模式修正形式以及北京航空航天大学近年研究的修正由压强梯度和湍能非平衡引起影响的模式，所采用的还有考虑壁面效应的低Ｒｅ数Ｋ－     

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

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

应用于翼型绕流的线性/非线性湍流模式的研究

本文选取了四个线性湍流模式、四个非线性涡粘性湍流模式和一个显式代数应力模式对绕翼型的不可压缩分离流动进行了数值模拟.因计算鲁棒性的需要,其中部分模式在壁面附近耦合了一方程模式.通过与实验结果的比较,对翼型在大攻角情况下流动产生分离的气动特性进行了评估.计算结果表明,非线性模式能够较好地反映湍流的各向异性和曲率影响.     

新型单方程湍流模型构造及其应用

为了提高雷诺应力本构关系式对于非平衡湍流的预测精度并且兼顾求解效率,发展了一种基于湍动能k的单方程(KDO)湍流模型。其主要思路为:采用平板直接数值模拟(DNS)数据对原始Bradshaw假设进行重新标定,使得当地湍动能和雷诺主应力之比能够根据当地流动条件进行自适应调节;同时,对标准k-ε模型中的湍流耗散率输运方程采用代数形式进行模化,进而形成一种一方程湍流模型。算例结果表明:KDO湍流模型对于对数率能够准确反馈,而在带有激波或部件干扰等流动现象的RAE-2822、ONERA-M6和DLR-F6算例中,KDO湍流模型能够准确控制湍动能的增长和衰减,相比于Spalart-Allmaras和Menter k-ω剪切应力输运(SST)模型,KDO湍流模型的计算结果有了较为明显的改善。     

适用于旋转通道各向异性湍流模型

基于代数雷诺应力方程的简化模型,并结合充分发展旋转通道的直接数值模拟（DNS）湍流脉动数据,发展了一种适用于旋转通道的各向异性k ω模型。采用该模型对进口雷诺数为6000,旋转数为0~0.26的旋转直通道进行模拟,结果表明:将旋转修正系数乘以传统的湍流黏度发展的各向异性k ω模型,能够准确地描述旋转状态下前缘面和后缘面的换热情况,是一种有效的各向异性湍流模型修正方法;旋转修正系数是否合理的关键是对雷诺应力比值进行准确地描述;通道的换热效果与旋转数和流向沿程无量纲位置有关,前缘面的换热随旋转数和无量纲流向距离的增加而减小,而后缘面的换热随旋转数和无量纲流向距离的增加而增大。      

不可压缩平板边界层转捩中突变过程的数值模拟

采用时间模式对不可压平板边界层从层流到湍流的转捩过程进行了直接数值模拟.分析结果表明:在层流-湍流转捩的突变过程开始之前,已经产生非线性作用;二维波与三维波的幅值实际增长率也偏离线性状态,而且在二者之间存在能量转移的迹象.在突变过程中,速度变形率对总应力的影响较大,在转捩结束后雷诺应力对总应力的影响比较显著.研究表明:正是由于大量谐波被快速激发并增大,导致了雷诺应力及其法向梯度也快速增大.这使得平均流剖面快速被修正,流动快速转为湍流.      

改进S-A湍流模型对横向射流的CFD模拟

通过在Spalart-Allmaras(S-A)湍流模型中加入应变率张量,来改善模型对横向射流(JICF)的计算精度。模拟了速度比分别为0.5和1.5工况下的横向射流。经实验值验证,结果表明:改进后的S-A模型预测精度有所提升;并对湍流模型改进部分对横向射流的作用机制做出了分析。将改进后S-A模型的计算结果与k-ε模型,雷诺应力模型(RSM)和大涡模拟(LES)的结果进行比较,发现S-A模型在速度比为0.5时精度最高,而在速度比为1.5工况下没有明显优势;并对这一现象做出了合理分析。      

应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

湍流模型对三维翼梢涡流场数值模拟的影响

 梢涡流场数值模拟对减少梢涡的不良影响意义深远。应用雷诺平均Navier-Stokes(RANS)方法建立了三维翼梢涡流场数值模型,选用的湍流模型包括两种经典两方程模型和一种代数雷诺应力模型。计算内容包括翼梢涡流场速度分布、压力分布及涡核位置,为了准确模拟梢涡涡核内系统旋转和流线曲率影响,将旋转和曲率修正方法应用于选用的湍流模型中。计算时对翼梢涡涡核处网格进行了加密处理。将计算结果与实验结果对比后发现,修正后的湍流模型较标准湍流模型明显提高了预报精度,与实验结果吻合较好,具有工程应用价值。     

矩形横截面螺旋管湍流流动与传热特性的数值研究

 采用了剪切应力输运(SST)k-ω两方程湍流模型并考虑近壁低雷诺数的修正对矩形横截面螺旋管内冷却水流动和传热特性进行了数值研究。数值分析了在不同入口雷诺数、曲率半径以及扭距条件下,螺旋管内的温度、速度场以及流线的变化,讨论了螺旋管内、外壁面对流传热系数的差异及产生机理,同时与直通管道传热性能进行了比较。研究发现由于离心力的作用,螺旋管内存在显著的二次流动,管内、外侧壁面对流传热存在差异。旋转一周后,螺旋管即进入了流动稳定状态,入口雷诺数可以显著提升螺旋管整体的对流换热效率,扭矩和曲率对内外壁面传热效果的影响不大,而窄高型的横截面构型可以显著改善螺旋管的传热效果。研究结果对应用矩形横截面螺旋管的冷却设计提供参考。     

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

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

曲壁气膜冷却数值计算的研究

 本文在一般二维边界层流动和换热计算方法的基础上,引入湍流混合长度的修正模型,用以考虑壁面曲率和注入冷气对换热的影响,研究计算了壁面的弯曲部分和曲率恢复部分换热的变化,并对曲面上的气膜冷却问题作了计算,取得较满意的结果。     

不可压缩湍流的色散模型

分析了广泛存在于湍流运动中的能量逆转现象,揭示了其产生的原因,提出了湍流的色散性质.在此基础上,修改了Boussinesq假设,建立了包含色散效应的新的雷诺应力封闭式和湍流色散系数,给出与不同模型相结合时,湍流色散系数所具有的不同形式,并阐述了湍流运动中能量的传递方向及条件.采用不可压缩平板边界层流动和平面后台阶流动验证了其可信性和优越性.平板摩擦阻力系数及边界层速度型与实验结果吻合良好,平面后台阶流动的流向再附长度、台阶边压力系数及湍流强度等参数均比标准k-ε模型更接近实验结果.结果表明:色散项的加入可以在不显著增加计算量的同时显著改善预测精度,模型具有一定的工程应用价值.      

S弯扩压器中四种湍流模型的比较

采用Realizableκ-ε,S-A(Spalart-Allmaras),SST(Shear stress transport剪切应力输运),RSM(Reynolds stress equation model雷诺应力模型)四种湍流模型对S弯扩压器内的分离流动进行了模拟,分析了不同湍流模型对计算结果的影响.将不同轴向位置的表面压力系数、分离区的大小、马赫数等三个量的实验值与计算值进行了比较,结果表明:分离区湍流模型对表面压力分布影响较大,未分离区不同湍流模型模拟精度相当;S-A模型对分离区的模拟较保守,RSM模型对分离的模拟精度较高,SST与RSM模拟结果很接近;各湍流模型计算的高速区域均比实验值大.扩压器的平均总压恢复系数受湍流模型影响不大,和实验结果吻合较好.      

标准的与非线性的k-ε模型在二维剪切流中的应用

 采用标准 k-ε和非线性 k-ε模型对二维直通道及小曲率弯道中的剪切流进行了数值计算。与实验结果对比后得到 :由于非线性 k-ε模型考虑了流体变形与应力关系中的非线性项,因而能正确地反映湍流各向异性,能正确预测二维剪切流中湍流强度及湍流动能等参数的分布