飞机增升装置的数值模拟研究

增升装置的设计是飞机设计的一项极其重要的内容.本文通过多块网格技术生成计算网格,运用全湍流SST(Shear stress transport)模型和带有转捩修正的SST湍流模型对飞机三维增升装置的气动力特性和流场特征进行了数值模拟研究,通过对数值计算结果与实验结果的对比分析,研究了两种湍流模型的特点及其在飞机增升装置...     

基于S-A与SSG/LRR-ω两种湍流模型的CHN-T1标模计算与分析

通过求解雷诺平均Navier-Stokes方程,对第一届AeCW-1研讨会上要求的所有工况进行了计算,着重研究S-A湍流模型和雷诺应力模型在CHN-T1标模上的表现,评估两种湍流模型的网格收敛性,对抖振特性的模拟能力以及对雷诺数的敏感性。根据与实验数据的对比分析,得出主要结论如下:两种湍流模型都能趋向于网格收敛,雷诺应力模型的网格收敛更快。在较大迎角的工况,两种湍流模型对机翼上激波诱导的分离流有较大的差别。CHN-T1标模的尾撑,改变了飞机平尾的压力分布,并且在飞机机尾上产生复杂的三维分离流动,对俯仰力矩系数的准确计算产生了较大的影响。考虑静气弹效应,两种湍流模型计算的力矩系数比不考虑静气弹效应更大,雷诺应力模型的偏差更大。在高雷诺数工况,两种湍流模型的计算结果趋于一致。     

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

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

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

航空发动机燃烧室流动数值计算中湍流模型的比较

分别采用标准k-ε湍流模型、RNGk-ε湍流模型、Realizablek-ε湍流模型以及雷诺应力模型,对某型航空发动机燃烧室流动进行了数值计算.近壁处理采用标准壁面函数法,计算得到速度矢量分布以及质量流量、湍流粘度比和湍流强度等参数.四种湍流模型计算的总体流动差别较小,但流场的细节有较明显的不同.标准k-ε模型、Realizable k-ε模型和雷诺应力模型的湍流粘性计算结果较为接近,而RNGk-ε模型计算的湍流粘性较小.     

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

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

SSG模型在压缩拐角激波/湍流边界层干扰模拟中的应用

采用非线性方案模化雷诺应力再分配项的二阶矩雷诺应力输运模型——SSG模型和两方程线性涡粘模型SST模型,对24°高雷诺数二维压缩拐角激波/湍流边界层相互干扰流动进行数值模拟。通过在壁面压力分布、壁面摩擦阻力系数分布、截面速度分布及分离区大小等方面与实验数据比较,综合评估了上述两种模型模拟该问题的能力。结果表明:两种湍流模型在上述各方面都取得了与实验较为一致的结果,但相比较而言,除SST模型在流动再附后对壁面摩擦阻力系数分布的模拟更接近实验值外,SSG模型结果在上述其他方面均更接近实验值。本文工作为压缩拐角激波/湍流边界层干扰这类流动数值模拟的湍流模型选择提供了参考。     

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

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

应用二阶雷诺应力封闭的湍流旋流计算

本文叙述了三种不同的湍流旋流流动有限差分的计算结果,即:涡流管中的流动、强自由旋转射流、轴对称,旋转突然扩张的湍流流动。湍流的模化是求解微分形式的雷诺应力分量输运方程的二阶封闭。计算的结果均与试验数据作了比较。比较的结果表明雷诺应力模型明显地优于常用的K—ε模型。而预示与试验数据的吻合程度也是令人满意的。     

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

NASA TrapWing高升力标模数值模拟研究

针对NASA TrapWing高升力全展襟翼构型，采用计算流体力学（ CFD）方法进行三维复杂流场仿真模拟，考察网格尺度和湍流模型对高升力模型气动特性的影响。采用“超立方体”概念，生成绕TrapWing模型的不同网格密度高质量多块结构网格，通过求解雷诺平均Navier-Stokes方程，研究网格尺度对高升力模型气动特性的影响。在此基础上选取中等规模计算网格，考察Spalart Allmaras和Menter k-ωSST湍流模型对高升力全展构型流场模拟能力，分析湍流模型对气动特性的影响。研究结果表明：网格策略具有较好的网格收敛性；湍流模型对翼稍附近位置上翼面压力系数的预测稍有影响，SA湍流模型预测的压力系数较SST更接近实验值。确认研究工作为大型飞机增升装置数值模拟提供了一定的参考。     

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

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

扩压器内跨音速湍流的数值模拟

采用Johnson-King非平衡代数雷诺应力湍流模型(J-K模型)和Baldwin-Lomax零方程湍流模型(B-L模型),数值模拟较强激波/边界层相互作用时扩压器内的分离流动。计算结果与实验值进行了比较,表明J-K模型比B-L代数湍流模型可较好地计算出分离流动的再附点位置,并且可更好地计算出激波强度和沿流程的压力分布,仅增加很少的计算量,并更易推广应用于三维湍流问题的数值模拟。      

包含旋转矢量的三维湍流模型

对雷诺应力方程中的旋转项进行分析后给出比较合理的模化关系。通过湍流的生成项、耗散项以及旋转项的局部平衡,建立了一个包含旋转矢量在内的各向异性的三维湍流模型。该模型对三维性较弱的剪切流是适切的。对旋转螺旋面及压气机转子叶片上的三维旋转湍流边界层进行了计算。计算结果与实验和C-S模型计算结果比较表明,该模型对近壁参数的预测能力有所改善。另外,该模型形式简单,从表达式本身来看,它相当于Bradshaw模型在三维问题上的推广。     

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

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

飞机增升装置气动力特性计算方法研究

本文给出了飞机增升装置气动力特笥的数值计算方法和计算结果，为了计算增升装置的缝隙效应和剪刀口效应，求解了粘性三维Navier-Stokes方程。为了正确 模拟大舵偏和阻力板后的大面积分离，采用了反映涡流特笥的新类型的两方程湍流模型。     

基于CFD技术的脊状表面湍流流场特性研究

在湍流状态下,利用计算流体力学(CFD)技术,对一系列布置不同尺寸顺流向脊状结构的表面进行了数值模拟。模拟过程采用Reynolds应力方程湍流模型,并且编制程序设定了充分发展的湍流入口边界。通过与光滑表面比较发现,脊状结构附近产生的流向涡是脊状表面减阻或者增阻的重要原因。流向涡改变了近壁面流场结构,引起了边界层近壁区雷诺应力和壁面剪应力的变化,脊状表面也因此出现了减阻或者增阻的现象。     

运用低雷诺数非线性k-ε模型研究湍流分离流动

给出了一个低雷诺数非线性ｋ－ε湍流模型，其中采用了一个非线性雷诺应力表达式并引入了两个考虑湍流近壁效应的修正函数。应用该模型和标准ｋ－ε模型计算了两种不同膨胀比的后台阶分离流。与实测结果比较表明，低雷诺数非线性模型结果较标准ｋ－ε模型结果有明显改进，它对流场的各特征量包括分离—再附长度、时均速度及各雷诺应力分量均给出了较好的预测。     

S-A模型在分离流动模拟中的改进

基于湍流不平衡性的分析,通过修正雷诺应力以及模型系数的计算,提高了S-A模型在二维边界层以及三维叶栅的分离流动上的预估能力。修正后的S-A模型给出了与实验吻合较好的结果,在对于分离流动的预估能力上优于原始模型。     

圆柱尾流中眼镜蛇探针与热线的性能对比

在雷诺数为7200的圆柱尾流中对眼镜蛇探针与热线的测量结果进行了详细对比,主要包括：时均速度、雷诺应力、概率密度函数、自相关函数和能谱等。从对比结果可以发现眼镜蛇探针数据的拒绝率随湍流度增大而增大,但在本实验中其对时均测量结果无明显影响。尽管眼镜蛇探针数据有效率在95%以上,但相对于热线其明显高估了雷诺应力。尤其是主雷诺剪切应力u＊v＊,眼镜蛇探针与热线测量结果的差异要显著大于两者在雷诺正应力上的测量差异。由于无法反映高频随机脉动的影响,眼镜蛇探针所测信号自相关系数的周期性明显强于热线。眼镜蛇探针对速度u和v的响应存在区别,使得两者能谱在惯性子区内有显著差异,不能像热线一样反映出惯性子区内小尺度湍流的各向同性。