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

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

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

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

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

具有无源控制空腔时正激波/湍流附面层干扰的数值模拟

采用雷诺平均N-S方程和B/L代数湍流模型计算了具有无源控制空腔时正激波/湍流附面层干扰流场。计算与实验结果的比较表明,本文方法可较准确地预测激波结构、激波与附面层干扰区流动基本特征及波后流动分离状态、激波位置、波前马赫数等参数。      

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

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

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

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

机器学习辅助湍流建模在分离流预测中的应用

数据驱动湍流建模是近年来发展的提高雷诺平均N-S方程预测精度的有效手段,通过机器学习算法能够从高置信度数据中自动提取特征,建立准确的从平均流动特征到雷诺应力的预测模型。针对高雷诺数下积冰翼型绕流这一类典型的复杂流动分离问题,基于此前研究者提出的机器学习预测框架,从输入输出特征选择和翼型绕流中数据分布特性两个方面出发,对机器学习预测结果的光滑性和准确性进行改善。提出了基于雷诺应力张量分析和流动特征辨识的输入特征选择准则;提出了局部区域建模方法以及基准模型和机器学习预测模型混合的代入计算方法。将改进方法应用于积冰翼型绕流问题之中,结果表明改进的方法能够准确给出训练集和预测集上的雷诺应力结果,并且代入平均流计算可以得到和真实分布更加接近的流动和机翼表面压力分布。     

雷诺应力与涡黏性模型的分离流预测对比分析

现代飞行器设计对流动分离的准确预测需求愈发迫切，但使用广泛的涡黏性模型预测结果却不尽如人意。雷诺应力模型凭借其扎实的理论基础有可能获得可信度更高的结果，但其性能优势仍需进一步评估与发掘。选取了SST模型、Stress-BSL(Baseline)模型分别作为涡黏性和雷诺应力模型的代表，对二维驼峰、二维跨声速凸块、跨声速三维ONERA M6机翼等算例进行了数值仿真。结果表明，相较于实验值，2种模型的预测结果均出现流动分离点提前，再附点滞后的现象，但Stress-BSL模型的预测误差更小，表现出了强逆压梯度下预测分离流动的优势。通过分析发现，2种模型对雷诺应力的低估导致了分离区较大。具体表现为SST模型引入的Bradshaw假设限制了湍动能的生成，使得模型计算的涡黏性系数偏小，强逆压梯度下低估边界层雷诺应力，导致流动分离提前。而分离区上缘处的雷诺应力预测偏小则被认为是流动再附滞后的主要原因。对于雷诺应力模型，误差主要来源于雷诺应力输运方程再分配项的模化不准。最后，针对上述原因，对SST模型关键封闭参数进行了重新标定，并进行了初步验证，结果表明修改后的模型预测表现好于原模型结果。     

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

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

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

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

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

计算流体力学中有限元的强间断处理

本文首次提出了无波动，无自由参数的有限元离散方法，它是基于矢通量分裂后其物理传播特性，利用差分格式构造过程中对激波上下游不同性质的分析，对一维Ｅｕｌｅｒ方程进行了仔细分析，并将其推广到二维，轴对称的无粘，有粘方中，从计算实践看，这种方法简便、实用，对激波，分离，激波与边界层干扰的复杂流场有较好的分辨能力和计算稳定性，比以往的迎风有限元，加入Ｌａｐｄｉｕｓ人工粘性修正的有限元有更好的流场分辨能力，使     

Stagnation temperature effect on the conical shock with application for air

The aim of this work is to realize a new numerical program based on the development of a mathematical model allowing determining the parameters of the supersonic flow through a conical shock under hypothesis at high temperature, in the context of correcting the perfect gas model. In this case, the specific heat at constant pressure does not remain constant and varies with the increase of temperature. The stagnation temperature becomes an important parameter in the calculation. The mathematical model is presented by the numerical resolution of a system of first-order nonlinear differential equations with three coupled unknowns for initial conditions. The numerical resolution is made by adapting the higher order Runge Kutta method. The parameters through the conical shock can be determined by considering a new model of an oblique shock at high temperature. All isentropic parameters of after the shock flow depend on the deviation of the flow from the transverse direction. The comparison of the results is done with the perfect gas model for low stagnation temperatures, upstream Mach number and cone deviation angle. A calculation of the error is made between our high temperature model and the perfect gas model. The application is made for air.     

三维非对称斜劈高速流动计算方法研究

交叉激波是高超声速进气道发展的关键气动问题之一。为实现对此类高速流动的准确数值模拟,基于RANS方程的自研多块结构网格计算软件,对经典三维非对称斜劈构型的斜激波交叉干扰流动特征开展了计算研究。通过对比分析,系统评估了计算网格、湍流模型及高精度格式对空间流场、底板摩擦力流线、压强分布、绝热壁温分布等的影响,获得了适用于三维非对称高速流动问题的网格生成准则和计算策略,同时验证了自研软件的模拟能力。     

两种紊流模型在通流计算中的应用研究

为验证紊流模型在基于周向平均Navier-Stokes方程的通流计算中的适用性，本文给出了周向平均形式的Spalart-Allmaras紊流模型，并选择Baldwin-Lomax模型作为对比。通过两个经典简单算例以及两个典型叶轮机械算例对采用了上述两种紊流模型的通流模型进行验证与对比。两个经典简单算例结果表明，上述两种紊流模型在无叶片的通流模型中均能够准确预测附面层发展，捕获激波，预测分离流动。两个典型叶轮机械算例结果表明两种紊流模型对轴流式叶轮机械的子午流面气动布局以及总体气动性能具备一定的的预测能力。在单转子算例中，S-A模型的计算结果略优于B-L模型的计算结果；在双级高压涡轮算例中两种紊流模型对总体性能参数以及气流角的计算结果相差不大，且S-A模型计算的马赫数比B-L模型的偏小。     

稀薄流到连续流的气体运动论统一数值算法初步研究

从非线性模型Ｂｏｌｔｚｍａｎｎ方程出发,引入简化速度分布函数、使用离散速度坐标法对速度空间进行离散、降维,去掉分布函数对速度分量的连续依赖性;采用时间分裂法,将简化速度分布函数松驰变化方程分解为源项碰撞变化方程、对流运动方程,进行耦合计算,应用ＮＮＤ耗散差分方法直接模拟气体分子速度分布函数;发展离散速度数值积分法,通过宏观取矩获取物理空间各点的流动参数,从而建立一套能有效模拟各流域气动问题的简化的     

唇口对侧压式高超声速进气道及等直隔离段影响的数值分析

用N-S方程和RNGk-ε紊流模型计算了RBCC用侧压式高超音速进气道及等直隔离段三维内流场,重点分析了唇口位置对来流Ma=6流场的影响,并分析了有关的现象。结果表明:唇口位置前移使得出口截面上气流的均匀程度、对气流的压缩程度、唇口激波的强度、进气道的流量捕获率和总压恢复系数等都比唇口位置的情况更大,而在喉部时介于二者之间。      

考虑间隙内流的二元机翼跨声速气动力分析

基于Navior-Stokes (N-S)方程和Spalart-Allmaras湍流模型,建立了含间隙二元机翼的跨声速绕流场计算模型和间隙内流计算模型.通过FLUENT数值仿真,分析了二元机翼的激波压力分布和马赫数分布,并与风洞试验 结果进行了对比验证.数值结果表明,连通波后高压区和波前低压区的间隙内流会提高机翼绕流场中激波前后的速度梯度和压力梯度.局部间隙内流会对机翼结构带来“向外吸”和“向前顶”的气动载荷作用.     

基于时间推进的涡扇发动机整机通流数值模拟

为实现航空发动机整机流场与性能快速分析,基于准三维控制方程模型,开发了整机通流数值模拟程序。程序基于计算流体力学理论,并考虑了引气、冷却和喷油等发动机中的物理现象。为模拟叶片对气流偏转作用,推导了一种鲁棒的无粘叶片力模型,同时采用粘性力模型模拟粘性损失效应。采用该程序对某型双轴分排涡扇发动机一个地面试车工况点进行了整机数值仿真,并与实验数据进行了对比验证。结果表明,所发展程序求解稳定,能够自动捕捉激波,并且可以快速获得完全收敛的航空发动机整机准三维流场,单工况点计算时长不大于30min。其中计算截面平均参数与实验测量偏差不大于8%,符合工程应用精度需求。     

高速粘性绕流的有限体积TVD格式

本文通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,解决二维超声速及高速声速流动问题。通过采用坐标交换技术,可以方便地将有限体积ＴＶＤ格式应用到各种形状的流动问题中,各种典型算例的计算表明该方法数值计算稳定,捕捉激波分辩率高,能很好地模拟超声速和高速超声流动中激波反射及激波与边界层的干扰等各种复杂流动现象。湍流计算中采用ＢａｌｄｗｉｎＬｏｍａｘ湍流。数值结果与实验结果及理论解的对比令人满意。