微通道流动转捩数值模拟

应用CFD计算软件ANSYS CFX对水力直径为0.4mm的矩形截面微尺度通道中的气体流动进行了数值模拟.采用了3种流动模型,即γ-Re_θt转捩模型、层流模型和剪切应力输运（SST）模型,并将模拟结果与实验结果进行了对比.模拟结果表明:层流模型及SST模型不具有预测转捩的能力;γ-Re_θt转捩模型捕捉到了临界雷诺数且摩擦因数相对误差在20%以下.在未达到临界雷诺数之前,通道后部已经出现局部湍流区,且局部湍流区随着雷诺数的增大而增大.壁面摩擦因数的最低点对应的流向位置与通道中前后湍流区的交结位置几乎一致.      

高压燃气涡轮径向内冷叶片气热耦合的数值分析

采用气热耦合方法对采用径向内冷方式的MarkⅡ型跨声速高压燃气涡轮金属导叶进行数值模拟,通过分析叶片通道内的传热和流动过程发现叶片表面附面层内流动非常复杂,包含层流流动、转捩和湍流流动状态,所以只有使用转捩模型计算的叶片附面层内流动与实际情况相符,叶片壁面温度和换热系数分布与实验结果吻合的较好,使用其他湍流模型由于不能准确描述附面层内流动而使得计算结果误差相对较大,但是所有的湍流模型都能很好的模拟附面层以外流动.      

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

应用Transition k-kl-ω转捩模型对内冷叶片气热耦合的数值模拟

选取NASA-Mark Ⅱ跨声速叶片为算例,研究了Transition k-kl-ω转捩模型在内冷叶片气热耦合计算中的应用,探讨了整场耦合与冷却通道内采用对流换热系数准则耦合的差异。结果表明,该转捩模型相比其它全湍流模型能够更准确预测附面层内的层流和转捩状况;由于Transition k-kl-ω转捩模型转捩前期采用层流动能来描述扰动的发展,避免了使用含有来流湍流度的经验公式,引入了"分裂机制"来描述层流与湍流脉动间的相互作用,并且在旁路转捩和自然转捩源项模化中加入了Tollmien-Schlichting波的影响,对强激波后的温度计算相比常用的间歇因子转捩模型与实验值更吻合;换热系数准则耦合用于冷却通道传热计算,避免了冷却通道边界条件带来的误差,计算结果与实验吻合较好,更易于工程应用。     

壁面温度对斜爆震发动机进气道流动特性影响的数值研究

为研究壁面温度条件对层流、转捩、湍流状态下斜爆震发动机进气道流场结构、流场参数的影响,选取Ma10级、具有曲面压缩段的斜爆震发动机进气道为研究对象进行数值模拟,对进气道壁面附近激波诱导分离区、热边界层的变化进行了深入探讨。数值模拟结果表明,进气道肩部圆弧过渡段出现的再层流化现象,壁面冷却对其起抑制作用,绝热壁面条件下再层流化程度最为严重。壁面温度的增加有利于延缓流动转捩,同时也导致了分离区尺寸的增加以及转捩、湍流状态下分离区主体位置逐渐前移,进气道内通道的转捩为分离诱导转捩,转捩位置主要受到分离点位置的影响,整体表现为壁面温度增加转捩位置前移。进气道出口顶板侧热边界层厚度随着壁面温度的增加逐渐变厚,转捩状态下热边界层厚度变化可达5%,温度峰值也随着壁面温度的增加逐渐增加,且峰值位置逐渐靠近壁面。壁面温度条件相同时,层流状态下热流、热边界层厚度均较小。转捩、湍流状态下进气道出口顶板侧热边界层较厚,约为层流状态3倍,同时转捩、湍流状态下热边界层厚度相差可达2%。     

高超声速圆锥边界层转捩数值研究

以类似于处理湍流脉动的方式来模型化非湍流脉动,考虑湍流的间歇性特征及引起超声速流动失稳机制,研究包括转捩起始位置及转捩区范围的高超声速流动转捩问题,发展可反映扰动模态和可压缩性影响的两方程SST转捩/湍流模式.针对尖锥模型,在马赫数6、攻角O°～2°时,采用第二模态(Second-mode)机制,预测并比较完全层流态、完全湍流态,以及自然转捩发展过程中的表面热流、表面摩阻等气动特性.     

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

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

飞行试验热流辨识和边界层转捩滞后现象

对中国空气动力研究与发展中心马赫数为5左右的球锥模型在首次航天模型飞行试验中的温度测量数据进行了分析,通过辨识获得热流分布,发现飞行试验的测热数据后处理方法与地面风洞试验有很大差别,必须考虑温度变化历史,并考虑测温单元与周围飞行器壳体的三维传热才能得到正确的热流结果。采用工程计算方法对模型表面热流分布进行了计算,通过与飞行试验测量结果对比分析,发现测温点在发射上升段由湍流完全变为层流和在再入下降段由层流向湍流转捩具有不同的转捩准则数,边界层转捩存在滞后现象;根据地面风洞试验拟合出的转捩准则受到风洞噪声等因素的影响,预测的转捩位置比实际情况靠前;对于球钝锥外形,当x/R50时,流场和热流趋于锥形流结果。本次模型飞行试验还首次验证了气动热工程方法对于马赫数小于5情况的适应性。     

真实气体效应对双楔绕流流动特性的影响

为研究高马赫数下真实气体效应对双楔绕流流动特性的影响，采用量热完全气体模型和化学非平衡气体模型两种气体模型对层流流态、转捩流态、全湍流流态的双楔绕流流场进行了数值模拟，其中转捩流态计算采用基于间歇因子的γ-Re_θt转捩模型，全湍流计算中采用k-ωSST湍流模型。结果表明：双楔拐角附近发生流动分离，与量热完全气体模型相比，相同流态下，化学非平衡气体模型计算得到的分离区尺寸更小。在边界层转捩位置的分析中发现，双楔拐角附近的分离区对边界层转捩有重要影响，分离激波产生的位置与边界层开始转捩的位置高度吻合，与量热完全气体模型相比，化学非平衡气体模型时的边界层转捩起始位置更靠后。针对真实气体效应下单位雷诺数、拐角角度对边界层转捩的影响规律研究发现，真实气体效应在一定程度上抑制了分离区内的转捩流动，随着单位雷诺数减小、拐角角度增大，真实气体效应对边界层转捩的影响随之增大。     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

Experimental investigation on aero-heating of rudder shaft within laminar/turbulent hypersonic boundary layers

The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m~(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test results is also validated by numerical simulations. Our analysis suggests that this phenomenon is due to the difference of boundary layer velocities caused by different thickness of boundary layer between laminar and turbulent flows, as well as the restricted flow within the rudder gap. When the turbulent boundary layer is more than three times thicker than that of the laminar boundary layer, the heat flux of the rudder shaft under the laminar flow condition is higher than that under the turbulent flow condition. Discovery of this phenomenon has great importance for guiding the design of the thermal protection system for the rudder shaft of hypersonic vehicles.     

考虑转捩的跨声速气冷涡轮叶片气热耦合计算

为了研究转捩对气热耦合计算的影响,在B-L代数模型与SST(shear stress transport)k-ω二方程模型的基础上,增加了两类基于间歇因子的转捩模型:代数AGS(Abu-Gharmam Shaw)模型与一方程间歇因子输运方程.选取NASA-MARKⅡ跨声速叶片为算例,分别采用全湍流模型与加入转捩的模型进行气热耦合计算.数值计算结果与试验对比表明由于能够预测附面层中的转捩过程,采用转捩模型的耦合计算得到的结果与试验吻合最好,由于在叶片壁面附近的网格较粗,采用间歇因子输运方程的转捩模型计算的结果要逊于采用代数转捩模型的结果.      

带有前缘气膜冷却燃气轮机叶片的气热耦合研究(英文)

使用气热耦合数值计算方法对带有"淋浴喷头"式前缘气膜冷却的C3X燃气轮机叶片进行数值模拟。通过对使用不同湍流模型的计算结果与实验测量值进行对比,使用转捩模型能够比较准确的模拟带有气膜冷却射流的边界层内的流动和传热过程。气膜冷却射流在叶片吸力面上游层流边界层区域内表现出比较明显的三维湍流特性,但在吸力面下游边界层转捩为湍流以后,湍流使得气膜冷却射流与主流燃气掺混过程加速而使其三维流动特性迅速消失;在吸力面由于边界层流动始终为层流状态,冷却射流与主流掺混较缓慢而使其仍能保持较强的三维湍流特性。边界层内由高温燃气与被冷却叶片之间温度差形成的温度梯度有益于增强边界层流动的稳定性。     

环形燃烧室冷态流场数值模拟中的数学方法

针对带双级径向旋流器环形燃烧室火焰筒,研究不同紊流模型(standardk-εmodel,RNG(renormalization group)k-εmodel和代数应力模型(algebraic stress model))及算法(SI MPLE(semi-i mplicitmethod for pressure-linked equations),PISO(pressure i mplicit split-operator))对复杂流场求解准确性及效率的影响.火焰筒三维贴体网格采用型线定点与偏微分方程相结合的方式生成,用壁面函数法处理近壁区域流动.计算结果表明,PISO算法提高求解效率;不同紊流模型预测结果均与试验数据符合较好,其中代数应力模型最为准确,更适用于复杂流场的数值预测.      

数值模拟纵向隔热屏加力室热态流场

数值研究纵向隔热屏加力燃烧室的在三维贴体坐标系下燃烧整体流场,所用的数学模型有:k-ε模型、EBU-Arrehenius和二阶矩紊流燃烧模型及六通量模型等.数值分析两种不同几何形状、进口气流参数分布及紊流燃烧模型等对加力室内各气流参数、隔热屏和加力室筒体壁面温度分布的影响,计算结果与试验数据比较表明:不同几何形状对加力室紊流燃烧流场的影响要比进口气流参数分布大些,二阶矩模型更适用模拟紊流燃烧流动.      

The laminar-turbulent transition zone in the boundary layer

The flow during transition from the laminar to a turbulent state in a boundary layer is best described through the distribution of the intermittency. In constant-pressure, two-dimensional flow, turbulent spots appear to propagate linearly; the hypothesis of concentrated breakdown, together with Emmons's theory, leads to an adequate model for the intermittency distribution over flow regimes ranging all the way from low subsonic to hypersonic speeds. However, when the pressure gradient is not zero, or when the flow is not two-dimensional, spot propagation characteristics are more complicated. The resulting intermittency distributions often show peculiarities that may be best viewed as ‘subtransitions’. Previous experimental results in such situations are reviewed and recent results and models are discussed. The problem of predicting the onset of transition remains difficult, but is outside the scope of the present article.Although this paper is intended to be chiefly a survey, several new results in various stages of publication are also included.     

Prediction of Boundary Layer Transition Based on Modeling of Laminar Fluctuations Using RANS Approach

This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...     

旋流燃烧室内多组分气体湍流混合的数值模拟

为合理地描述旋流条件下组分质量的湍流输运,将考虑湍流-旋流相互作用的新代数质量流通量模型与新代数Reynolds应力模型结合,构成了新代数应力/质量流通量模型。应用该模型对旋流燃烧室内两股不同组分射流的相互作用与湍流混合进行了数值模拟,得到了与实验相符合的气体轴向与切向速度、轴向与切向脉动速度均方根值和氦气浓度分布。模拟结果与k-ε模型的结果进行了对比。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

加力室液雾燃烧的模型与应用

提出了描写航空发动机加力燃烧室内气液两相湍流燃烧的控制方程,采用了双流体模型,高温辐射通量模型,双时间尺度湍流K-ε模型,应用SIMPLE方法解这组方程,对二维轴对称的涡喷加力燃烧室进行计算,证明应用本文的液雾燃烧模型可以得出合理的计算结果。