应用各向异性湍流模型对气膜冷却的数值研究

应用引入各向异性系数的k-w湍流模型和标准k-ε模型,对平板气膜冷却经典实验进行了数值研究,分析了各向异性系数的引入对气膜孔下游流场与温度场的影响机理,探讨了该系数的适用性问题.研究表明:①展向各向异性系数的引入,增加了流体在展向的湍流脉动,强化了气膜冷气沿展向的扩散;②各向异性模型的计算结果在低吹风比(M=0.5)的工况下与实验值符合得较好;③各向异性系数在气膜脱离壁面与主流掺混的区域适用性较差.      

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

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

区域分裂法在湍流流场数值解中的研究

应用区域分裂算法(简称DDM)对湍流流场进行数值研究。对湍流流动控制方程建立一种通用的DDM-1计算格式S_{I,ω,给出重叠区域内变量处理方法。DDM-1的SI,ω计算格式的适用性由突扩圆管分离流流场数值计算所验证,数值计算结果与实验数据比较表明:采用SI,ω计算格式所预测的再附点位置及平均轴向速度、湍动能强度比在整个计算区域内采用单一模型的方法能够更好地与实验值吻合。}      

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

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

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用.与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要.因此,基于两步四阶时间离散和高精度WE-NO重构,发展了具有四阶时间精度的气体动理学格式.该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力.目前...     

GAO-YONG湍流模型对湍流传热的数值模拟

将GAO-YONG湍流模型应用于湍流传热的研究,分别计算了平板剪切湍流和二维平面冲击射流的湍流传热问题.边界层剪切湍流流动与换热的计算表明:与传统的湍流模型不同,GAO-YONG湍流模型不需要对近壁区域做任何特殊处理(比如壁面函数、低Reynolds数修正等)即可模拟出从壁面到主流区的全部流动与传热情况;另外,对于冲击射流Nusselt数的模拟也得到了与实验符合较好的计算结果,准确地捕捉到了2种冲击高度下流场换热的不同特征,表明了GAO-YONG湍流模型能够较高精度地计算湍流换热.      

湍流燃烧模型在燃烧室数值模拟中的对比分析

以某单头部矩形燃烧室为研究对象,采用多种湍流模型和燃烧模型进行组合计算,模拟燃烧室内部的速度场和温度场,并对计算结果进行对比分析。结果表明:湍流模型主要影响火焰筒内部主燃孔横截面上游的速度分布,Standard k-ε和Realizable k-ε模型的速度场计算结果差异相对较小;湍流动能预测受湍流模型的影响较大,并具有一定规律性;不同组合模型对燃烧室内部和出口温度分布的局部细节模拟差异较大,燃烧模型影响最大;PDF模型计算的温度值较合理,另外三种燃烧模型在单步完全反应燃烧机理下的计算值偏高,计算获得的OTDF也相差较大。     

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。     

两种k-ε模型用于扩压器内流数值模拟的比较

同时采用标准k-ε及非线性k-ε模型在附体坐标系下求解定常不可压N-S方程组计算扩压器内湍流流场, 得到了扩压器的压力恢复系数和速度分布, 与实验结果对比发现非线性k-ε模型预测的结果更吻合实验值, 且能测出临界扩压角下的局部分离流。文中着重研究了非对称扩压器, 并给出了用两种模型计算的湍流量。      

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

对 k-ω 湍流模型的改进

为了更适合于较复杂几何外形的数值计算，本文提出仅引入１个参数来改进ｋ－ω湍流模型并采用了一种双网格系统的有限体积格式来求解Ｎ－Ｓ方程与改进的的ｋ－ω模型方程。通过平板附面层的计算结果表明了用本文提出的对ｋ－ω湍流模型的改进对控制转捩是有效的。     

低雷诺数κ-ε封闭模式在三维湍流边界层计算中的应用

本文采用低雷诺数形式的k-ε湍流模型,计算具有大横向流的不可压缩三维湍流边界层。在湍流模型中,应用了涡粘性的各向异性系数。本文所得结果与用混合长模型所得结果以及量测结果的比较表明:对于边界层的不同特征量,湍流模型和涡粘性各向异性对计算结果的影响程度是不同的。总的来说,应用考虑涡粘性各向异性的k-ε模型,较之应用混合长模型,所得计算结果有所改善,但还不是根本性的改善。     

k-ε-g紊流燃烧模型数值研究

本文对k-ε-g紊流燃烧模型中模型系数C_R,C_(g1),C_(g2)的选定和计算浓度脉动均方值g的不同方法进行数值研究。计算结果表明模型系数值的选定直接与化学反应速率有关,而其中C_R影响最大。g的计算虽然可有不同方法,而其中微分方程求解最佳。并用实验数据来验证计算结果。     

基于火焰面模型的超声速湍流燃烧数值模拟研究

在非结构混合网格有限体积框架下,将RANS方程和层流火焰面模型相结合,开展超声速湍流燃烧流动数值模拟研究,时间和空间离散分别采用LU-SGS格式和HLLC格式。用SST k-ω模型模拟湍流,FlameMaster 3.9生成火焰面数据库,并采用β-PDF和δ-PDF分布来考察湍流和化学反应的相互作用。采用该数值方法,对DLR氢燃料超燃燃烧室和北航乙烯燃料气动斜坡超燃燃烧室开展了数值模拟研究,并将数值计算结果与实验数据进行对比,二者符合较好。研究结果表明:在主流中心支板顺流喷射超燃燃烧室内,燃烧主要发生在支板后的燃料/空气混合层内,并且火焰核心落在恰当混合物分数Zst附近区域。      

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

针对叶轮机流场分析的不同Spalart-Allmaras湍流模型评价

Spalart-Allmaras（SA）湍流模型在Reynolds-Average Navier-Stokes（RANS）方程求解中得到了广泛的应用。针对叶轮机数值模拟,国内外学者提出了很多修正的SA模型,但它们之间的对比少有研究。为此,综述了原始的SA（SA-standard）,SA-neg,SA-Helicity,SA-Ning及SA-R五种不同SA湍流模型。针对NASA Rotor 67和NASA Rotor 37,考察这五种湍流模型对数值计算稳定性的影响及对叶轮机内部流动细节捕捉的能力。研究表明：SA-standard和SA-neg湍流模型计算结果几乎相同;相比于SA-standard湍流模型,SA-Ning湍流模型计算的压比、效率和堵塞流量均较大;SA-R湍流模型模拟的激波/边界层干涉更严重,分离区域更大,计算的压比小,模型修正效果不理想;SA-Helicity湍流模型不仅能极大地提高计算的压比和流量范围,而且能提高计算的稳定性和失速工况下的计算准确性,但计算的堵塞流量和效率较小。     

三维进气道湍流流场数值模拟

结合实验数据,采用有限容积法,对超音速S形进气道三维湍流流场进行了数值模拟,分析了0°攻角、10°攻角、10°侧滑角三种工况下机头激波对入口气流的影响、进气道入口激波分布和进气道出口压力分布情况.计算结果表明:0°攻角和10°攻角时,机头激波对入口气流影响不大,出口气流总压分布相对均匀,且高压区面积较大;而10°侧滑角时,受机头激波的影响,进气道前方出现小范围的低压区,且出口气流低压区较大.     

发动机加力燃烧室湍流流场数值计算

采用非正交曲线坐标系下非交错网格的ＳＩＭＰＬＥ方法,对航空发动机加力燃烧室气相燃烧的湍流流场进行了数值计算,湍流模型采用ｋ－ε双方程模型,平均化学反应速率采用涡旋破碎模型（ＥＢＵ）计算,对ＥＢＵ模型的缺陷作了讨论。差分网格采用分区方法生成,计算时对整个流场进行分区迭代直至得到收敛结果,数值计算结果合理。     

超声速H_2/Air湍流扩散燃烧RANS数值模拟

为了研究修正的火焰面反应进度变量燃烧模型在超声速湍流扩散燃烧问题中的适用性,对德国宇航中心(DLR)超声速燃烧室开展RANS数值模拟。基于Open Foam软件平台中密度求解器分别对三维冷态场和燃烧场进行模拟分析。将网格自适应加密技术用于流场的计算;燃烧场计算中,通过分析不同压力下层流火焰面数据库,引入了反应进度变量源项的压力修正系数,压力修正系数α等于2.2。计算结果表明,冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合较好。压力修正方法能够较好地解决超声速湍流扩散燃烧问题。湍流Schmidt数敏感性分析表明,湍流Schmidt数Sct对湍流火焰结构有较大影响,文中Sct等于0.7时能得到与实验值较为一致的分布。     

不同湍流模型对插板式进气畸变的数值模拟

为得出一种较为准确的湍流模型,应用三种湍流模型对插板式进气畸变进行数值模拟,通过与畸变实验实测数据的对比,研究不同湍流模型在插板式进气畸变计算中的预测能力。不同湍流模型计算出的出口流场畸变不同,重整化群k-ε模型得出高压区范围较大,S-A模型过渡区域范围大。结果表明:插板高度为畸变发生器内径的48%时,低压区范围角与实验值较为吻合,插板高度为畸变发生器内径的20%时,低压区范围角偏小。畸变指数的计算值以可实现化k-ε模型最为准确。