两方程湍流模型的可压缩性修正及其应用

本文在不可压两方程湍流模型的基础上进行可压缩性修正,给出了ｋ－ε、ｋ－ε两方程模型的可压缩性修正方程,并对亚声平板湍流流动和跨声喷管流动进行了数值模拟,得到了较好的结果。从而证实了两方程湍流模型经过可压缩性修正,能较好地模拟可压缩流动。     

多尺度湍流模式在压缩拐角绕流中的应用

 应用4 个多尺度涡粘性湍流模式对一个二维压缩拐角流动进行了数值模拟,通过与实验结果以及两个基准模式(一个单尺度k-ε两方程模式和一个零方程模式) 的计算结果的比较,结合压缩性对湍能输运的作用、湍流的非平衡性对多尺度湍流模式在拐角非平衡湍流中的应用进行了评估。     

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

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

一种新的非线性K-ε两方程湍流模型

本文基于理论分析和对标准Ｋ－ε两方程湍流模型数值模拟结果的分析,在标准Ｋ－ε两方程模型的雷诺应力表达式中加入６高阶非线性项;并在模型的涡粘性系数μ１中考虑流动各向异性的影响,将涡粘性系数中参数Ｃμ取为表征流动各向异性的参数Ａ的一个线性函数,从而得到了一个新的非线性Ｋ－ε两方程湍流模型。     

比熵增概念及其在湍流模型中的应用

通过对不同马赫数下(Ma=0.7、2.25、6)平板边界层进行直接数值模拟研究,提出比熵增的概念.比熵增表征单位机械能的耗散,并且在数值上表现出较好的单调性和马赫数无关性,因而能够可靠地表征边界层的范围.应用比熵增概念修正BL零方程湍流模型长度尺度,构造BL-entropy,同时采用原始的BL模型,SA一方程模型对后台阶(Ma=0.128)以及高超声速锥柱裙组合体(Ma=7.05)流动进行数值模拟,并与实验数据进行比较,结果表明BL-entropy能够准确地获取长度尺度进而得到均匀合理的涡粘性分布,极大地提升了原始BL模型模拟复杂流动的能力.     

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

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

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

增程制导航弹绕流场数值模拟

研究某增程制导航弹绕流场的结构,为气动特性分析和外形改进提供依据。以三维N-S方程为主控方程,引入S-A一方程湍流模型,对某增程制导航弹的绕流场进行了数值模拟,给出了流场图像和气动特性计算结果,并将CFD计算结果与试验结果进行了比较分析。CFD计算结果与试验结果吻合良好,验证了所采用CFD方法的准确性和有效性,数值模拟结果已得到应用。     

室内气流分布的大涡模拟

准确地模拟室内气流分布对优化室内通风设计及预测舒适和健康的室内环境有着重要的意义，传统的空内气流分布的数值模拟大都采用к-ε两方程湍流模型，该方法求解出的是流动变量的平均值，无法给出流场结构的详细信息，体现不出湍流流动的瞬时性的特点。本文采用大涡模拟（LES）法求解室内气流分布，使HSMAS法保证计算的稳定性和收敛性。通过对室内等温贴附射流和非等温受限射流两种类型的射流作用下的室内气流分布计算，以及导入к-ε两方程湍流模型对计算结果进行比较与分析表明，大涡模拟求解室内气流分布是一种确定可行和有效的方法。     

绕翼型低速湍流的数值模拟

采用求解低速流动的SIMPLE算法,用标准K-ε两方程模型、非线性K-ε两方程模型和一种改进的K-ε两方程模型来对绕两种翼型（A翼型和GAW-1翼型）在接近失速攻角情况下的低速湍流流动进行了数值模拟。计算结果表明,非线性K-ε两方程模型和改进的K-ε两方程模型较好地模拟出了翼型表面上的分离流动特性,较准确地预测出了分离点位置以及翼型的气动力系数,对翼型与机翼的工程设计具有一定的参考价值。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

高超声复杂流动中湍流模式应用的评估

本文选择几个高超声速基准流动：二维可压缩拐角、锥柱裙组合体绕流、斜激波与湍流边界层干扰，采用几个常见的湍流模式，BL模式，CHκ-ε模式，κ-ε模式，SST模式，CMOTTκ-ε模式，SHIHκ-ε模式，通过将数值计算结果和实验结果进行比较，对有关的湍流模式地 评估，得到一些有意义的结论。     

基于Reynolds平均的高超声速二方程湍流模型

 在Reynolds平均的框架下推导了可压缩湍流Reynolds应力方程和湍动能方程。根据一定的假设和尺度分析简化并封闭了所推导的湍动能方程。在均匀湍流假设下,湍动能耗散率可分解成可压缩性耗散和旋度耗散,并对其中的可压缩性耗散进行了封闭;同时认为旋度耗散不受可压缩性影响,直接引用传统的Reynolds平均不可压缩湍动能耗散率模型方程。由此构造了适用于高马赫数的二方程可压缩湍流模型。应用所发展的模型计算了高超声速平板绕流,并与若干现有模型的计算、实验与半经验公式的计算结果进行了对比,验证了所发展的模型。在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对所发展的湍流模型,以及若干现有模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了湍流模型可压缩修正的必要性,得到了所发展模型的计算结果最接近实验结果的结论。     

关于两方程湍流模型的考虑

通过对两方程湍流模型的分析,提出了在主流区或远壁区用两方程模型,近壁区用代数湍流模型的思路,即所谓的“双层”模型概念。通过对平板湍流附面层的计算表明,这种分层计算措施获得的速度分布、壁面剪切力较单纯应用低雷诺数两方程湍流模型时的结果准确。对非均匀来流条件下90°弯管内高湍流度流场的计算表明,双层模型可获得更加准确的速度分布。从而表明本文提出的分层计算措施是有效的,可满足工程计算要求。      

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。     

基于物理知识约束的数据驱动式湍流模型修正及槽道湍流计算验证

对湍流摩擦阻力的精准预测是学术界和工业界普遍关心的重要问题，而数据驱动式的湍流模型修正方法对此显示出较大的潜力和前景。提出了一种基于物理知识约束的数据驱动式湍流模型修正方法，根据湍流摩擦阻力分解获得先验物理知识，在S-A湍流模型的生成项中引入非均匀分布的修正因子，以修正因子为设计变量，设定包含物理知识约束的目标函数，利用离散伴随方法求解目标函数与设计变量之间的梯度关系，通过高效率的迭代求解获得修正因子的分布。以槽道湍流为例，验证了包含物理知识约束的数据驱动式建模方法的优势，并分析了物理知识约束对湍流摩擦阻力预测精度的影响，结果表明引入物理知识约束可进一步提高湍流摩擦阻力的预测精度。     

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

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

基于遗传算法及转捩模型的层流翼型优化设计研究

讨论了层流翼型对飞机机翼减阻的重要意义,概述了遗传算法的优化原理及其特点。在基准翼型的基础上,以Hicks—Henne型函数的解析函数线性叠加法来描述翼型。采用遗传算法完成了层流翼型优化设计。计算中采用了转捩模型耦合SST两方程湍流模式,通过求解雷诺平均的N-S方程组模拟了翼型的转捩流动。对计算结果进行了分析,分析表明优化翼型在升力特性、阻力特性方面比基准翼型有很大的提升,达到了优化升阻比的目标。     

应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.