使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

SA一方程湍流模型参数影响分析与辨识

SA(Spalart-Allmaras)一方程湍流模型是目前工程湍流计算中主要采用的湍流模型之一。首先,针对某典型战斗机的小攻角、中等攻角、大攻角流动工况,利用均匀试验设计方法分析SA一方程模型中8个参数取值对上述工况下飞机升力和阻力系数计算结果的影响规律;然后,建立工程湍流模型参数的辨识方法,并将其用于该战斗机大攻角工况下湍流模型参数cb1值的辨识调整。结果表明:不同工况下,湍流模型参数对计算结果的影响规律不同;在附着流状态下,对升力和阻力影响较大的参数是cv1和cb1;在中等攻角和较大攻角下,对升力和阻力影响较大的参数是cb1;适当减小参数cb1的取值后,升力和阻力系数的计算结果有较明显的改善,这可能与飞机大攻角分离流场中涡粘系数和剪切应力的发展与自由剪切流存在一定差异有关。     

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-ω预测得更准确.      

一种基于湍流模式的转捩预测方法

通过对平板边界层进行计算分析,基于湍流模式来预测流动转捩是可行的,但预测出的转捩位置不准,比实际情况靠前。因此,本文从转捩流动的物理特征出发对湍流模式中的混合尺度作出了修正,提出了一种考虑间歇函数影响的k-ωSST两方程湍流模式。采用该模式来对平板湍流边界层和不同攻角下的NACA0012翼型绕流进行了转捩预测,计算结果表明,改进后的模式对转捩位置具有较好的预测能力。     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

湍流模型对湍流射流CFD模拟的影响

横流中的湍流射流涡旋结构复杂,湍流模型对其模拟结果影响很大。采用S-A,Realizable k-ε,k-ω和SST k-ω四种湍流模型对吹风比为0.5和1.5工况下的横流中射流进行了数值模拟研究。与实验值比较验证了数值模拟结果的可信性。分析了旋转张量和剪切张量对湍流生成的影响。发现四种湍流模型对下游速度分布的预测均较为准确,在小吹风比下S-A模型性能最好;在大吹风比下k-ω模型较为准确。     

两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。     

超声速底部排气弹底部流场与气动特性研究

为了研究超声速底部排气弹气动特性,采用AUSMPW+迎风格式、k-ωSST湍流模型、8组分12反应化学动力学模型和二阶矩湍流燃烧模型耦合求解三维带化学反应的Navier-Stokes方程。在数值方法的有效性和可靠性得到验证的基础上,对超声速底部排气弹底排真实气体流场进行了数值模拟,分析了攻角和船尾角对底部排气弹的底部流场结构和气动特性的影响规律。计算结果表明:攻角对底部流场结构影响较大。随着攻角的增大,迎风面和背风面的初始回流区体积逐渐减小,且迎风面的初始回流区体积始终小于背风区;随着攻角的增大,底部阻力系数、总阻力系数等气动系数均增大;不同工况下存在着相应的最佳船尾角。船尾角的改变会引起底部流场结构的变化,同时影响着富燃气体的二次燃烧区域与强度。有底部排气时对应的最佳船尾角比无底部排气时的最佳船尾角小。     

一种超声速射流计算方法与实验验证

高温射流流场计算是尾焰辐射目标特性计算的前提,然而由于缺少可靠实验数据,针对湍流超声速射流的数值模拟多集中于低温射流,高温射流计算与实验的对比工作还很少见。利用k-ωSST双方程湍流模型,模拟了多个典型超声速射流实验的流场速度与温度分布,通过与实验结果进行对比,建立了一种超声速射流计算方法。首先,通过对比多个低温射流的实验与计算结果,探索了湍流模型中可压缩修正以及来流湍动粘性比对超声速射流计算结果的影响;进而,针对火箭发动机尾焰实验,计算尾焰流场与流场红外辐射,流场辐射计算结果与实验观测结果符合较好,进一步验证了计算方法。最终认为经过可压缩修正的k-ωSST双方程湍流模型结合湍动粘性比取值30可以作为超声射流计算中较为典型的湍流计算方法。     

多尺度湍流模型在圆柱绕流中的应用

使用基于可变时间间隔平均方法的多尺度湍流模型,对高超临界区雷诺数下非定常圆柱绕流进行了数值模拟.结果表明:多尺度湍流模型在升阻力系数预测上与实验值一致,而标准k-ε模型和SST(shear stress transport)k-ω模型却只能兼顾其一;在旋涡脱落特性预测上与实验值相一致,其误差小于SST k-ω模型的5.8%和标准k-ε模型的37.6%;在表面平均压力系数分布预测上其误差仅为3.6%,明显优于SSTk-ω模型的13%和标准k-ε模型的53.7%;在表面摩擦因数分布及分离角度方面与实验结果相吻合,且分离角度误差仅为0.78%,结果优于标准k-ε模型的1.04%和SSTk-ω模型的1.83%,充分验证了该湍流模型应用于复杂湍流模拟预报的潜力.     

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

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

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

Pω增强型k-ω湍流模型在三角翼旋涡流动的应用

在三角翼旋涡绕流数值模拟中,标准 Wilcox k-ω湍流模型生成项未考虑旋度的影响而导致预测的旋涡强度较弱。通过引入探测因子区分剪切层和涡核,在旋涡流动的高旋度区域增加ω方程生成项的方法,基于结构化网格上的 RANS 求解器,加入了 Pω增强型 k-ω湍流模型,对绕尖前缘三角翼亚声速和跨声速旋涡流场进行了数值模拟。计算结果与 NASA 的 NTF 风洞和 DLR 的 DNW-TWG 风洞试验数据进行了对比分析,结果表明：不论在亚声速还是跨声速自由来流条件下,Pω增强型 k-ω湍流模型计算的压力分布、涡破裂位置均与试验数据吻合良好,准确地预测出了三角翼上翼面的主涡、二次涡结构,特别是跨声速条件下激波干扰导致的涡破裂的临界迎角及涡破裂位置,表明 Pω增强型 k-ω湍流模型在绕三角翼旋涡流动数值模拟中具有良好的适用性。     

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

临近空间螺旋桨低雷诺数高效翼型数值分析

针对临近空间翼型边界层变厚更容易发生分离的特点,在修正了湍流模型后,通过分析比较13种低雷诺数翼型的升阻系数得到了比较适合用于临近空间螺旋桨叶素的高效翼型。结果显示,攻角在-6°～8°、8°～14°、14°～22°三种工况下,分别采用SST k-ω、RNG k-ε、Realizable k-ε湍流模型,可以得到与国外实验比较接近的合理结果;工况条件对翼型气动性能好坏的影响很大,翼型S-1223和FX63-137在所研究的工况内都具有较好的气动性能;因此可以选择这两种翼型作为临近空间螺旋桨用高效翼型。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

厚覆冰导线升力突变机理及数值模拟研究

针对新月形厚覆冰导线的升力系数在风攻角15°附近存在突变的问题,分别采用基于k-ωSST湍流模型的雷诺时均法和大涡模拟(LES)的数值方法对新月形厚覆冰导线在风攻角10°～20°范围进行了模拟。通过对比两种数值方法计算得到的覆冰导线气动力系数、流场结构和表面风压,发现LES方法能够更好地捕捉新月形覆冰导线表面的小尺度涡结构,得到的覆冰导线气动力参数计算结果与风洞试验数据高度吻合;而k-ωSST湍流模型难以模拟壁面上小尺度涡,捕捉不到升力系数的突变。根据覆冰导线不同壁面区域的压力分布,发现上侧壁面处的涡结构影响整体流场,并在下侧壁面曲率、来流夹角和壁面切线方向共同作用下导致升力系数突变。LES的气动力参数模拟结果可为覆冰导线防舞提供参考。