应用GAO-YONG可压缩湍流模式数值模拟RAE2822翼型绕流

应用Gao-Yong可压缩湍流模式,数值模拟RAE2822二维翼型在两种不同来流情况下的跨音速粘性绕流问题.湍流模式的对流项用ROE格式离散,扩散项用中心差分格式离散,空间离散后的控制方程用多步Runge-Kutta显式时间推进格式求解.计算结果预测了翼型表面的压力系数的分布、平均速度剖面、激波的位置、马赫数等值线等情况.同时,对翼型表面激波与边界层相互干扰以及转捩问题进行分析计算,结果表明,Gao-Yong可压缩湍流模式结合适当的数值方法能够成功地模拟翼型跨音速粘性流动.最后,基于Gao-Yong可压缩湍流模式各项异性湍流粘性的机理,初步提出一种预测转捩起始位置的方法.     

转捩与湍流对激波边界层干扰及底部流动结构的影响

为研究转捩与湍流对激波边界层干扰及底部流动结构的影响,文章选取了二维与三维高超声速双斜面进气道模型与大钝头着陆器模型,并使用γ-Re_θ转捩模型开展数值模拟研究.研究表明,对于二维进气道模型,随着前缘钝度的增加,激波边界层干扰位置前移,分离区变大,与层流流动情况相比,有转捩流动发生时,激波边界层干扰位置后移,同时分离流动强度变弱,分离区缩小;对于三维进气道模型,其拐角附近的分离流动呈现明显的三维特征,转捩流动也存在三维流动结构,与静风洞状态相比,噪音风洞状态下,有转捩流动发生,对壁面热流影响较大,对激波系影响很小.对于着陆器模型,底部流动发生转捩,使得底部流动由不稳定非定常的流动结构变为稳定定常的流动结构,这有益于姿态控制设计.      

一种数据驱动的翼型流动转捩预测方法

研究翼型绕流的转捩预测方法,对于翼型流动细节的精确模拟和气动力的准确计算以及精细化设计均具有十分重要的意义.采用动模态分解（dynamic mode decomposition,DMD）代替线性稳定性理论（linear stability theory,LST）与eN方法结合,不需要求解稳定性方程,成为一种数据驱动的翼型边界层转捩预测新方法,称为DMD/eN方法.在原有方法的基础上,改进了DMD网格线生成方法和扰动放大N因子的积分策略,并将RANS求解器与改进的DMD/eN方法进行耦合,实现了翼型定常绕流转捩预测自动化.采用该方法对LSC72613跨声速自然层流翼型以及NLF0416低速自然层流翼型在不同攻角下的绕流进行转捩预测,转捩点计算结果均与实验值和LST/eN方法吻合良好.该方法计算得到的N值增长曲线与LST/eN方法的包络线也较为吻合,进一步验证了积分策略的正确性.改进的DMD/eN方法可作为自然层流翼型设计的新的有力工具.      

压缩拐角激波/湍流边界层干扰特性分析

用大涡模拟方法对Mach数3.0下的压缩拐角激波/湍流边界层干扰问题进行数值研究.对拐角上游平板区域边界层转捩及湍流进行模拟,设定平板区域长度,使得转捩过程于平板区域发生并充分完成,从而在拐角处产生激波/湍流边界层相互干扰,研究激波/湍流边界层的作用机理.研究表明:流场能够在非定常扰动激励下迅速转捩,并于平板区发展为完全湍流;湍流边界层与激波相互作用过程中,拐角附近分离区较层流情况明显减小;展向不同区域分离区大小差异较大,局部区域分离现象消失.     

基于层流动能湍流模型的数值模拟方法

从层流动能出发提出数值模拟原则;综合考虑自然、旁路和分离流转捩的因素构建实用的层流动能湍流模型,结合预处理和基本求解技术发展出适于转捩流动的数值模拟方法和程序.针对预处理技术,以Weiss-Smith矩阵为基础,考虑湍流粘性的影响;针对基本离散格式和边界条件,结合模型方程进行对角占优强化等特殊处理.最后通过平板边界层和典型翼型,特别是低雷诺数翼型的数值模拟,验证数值方法的有效性和鲁棒性.算例表明本文的方法能够为求解更复杂的流动提供参考.     

湍流模式对转捩的初步研究

采用二阶时间精度的LU-SGS-τTS目的进行时间推进求解Navier-Stokes方程组,用基于Roe的六阶对称TVD格式离散对流项,并采用Baldwin-Lomax和Spalart-Allmaras模式对平板边界层、Aerospatial A-profile翼型低速大迎角及RAE-2822翼型跨声速激波/边界层干扰流动进行数值分析,自动捕捉流动的转捩过程,获取了有意义的结果.     

超声速层流/湍流压缩拐角流动结构的实验研究

在Ma=3.0的超声速风洞中, 分别对上游边界层为超声速层流和湍流, 压缩角度为25°和28°的压缩拐角流动进行了实验研究. 采用纳米粒子示踪平面激光散射(NPLS)技术获得了流场整体和局部区域的精细结构, 边界层、剪切层、分离激波、回流区和再附激波等典型结构清晰可见, 测量了超声速层流压缩拐角壁面的压力系数. 从时间平均的流场结构中测量出分离激波、再附激波的角度和再附后重新发展的边界层的增长情况, 通过分析时间相关的流场NPLS图像, 可以发现流场结构随时间的演化特性. 实验结果表明: 在25°的压缩角度下, 超声速层流压缩拐角流动发生了典型的分离, 边界层迅速增长失稳转捩, 并引起一道诱导激波, 流场中出现了K-H涡、剪切层和微弱压缩波结构, 而超声速湍流压缩拐角流动没有出现分离, 湍流边界层始终表现为附着状态; 在28° 的压缩角度下, 超声速层流压缩拐角流动进一步分离, 回流区范围明显扩大, 诱导激波、分离激波向上游移动, 再附激波向下游移动, 分离区流动结构复杂, 相比之下, 超声速湍流压缩拐角流动的回流区范围明显较小, 边界层增长缓慢, 流场中没有出现诱导激波、K-H涡和压缩波, 流动分离区域的结构也相对简单, 但分离激波的强度则明显更强. 关键词： 压缩拐角 层流 湍流 流动结构     

后掠机翼的横流不稳定性分析及转捩预测

后掠机翼边界层的流动稳定性及转捩对翼型的设计及优化有着重要的参考价值,而横流失稳是引起后掠机翼边界层转捩的关键因素之一.以NLF（2）-0415翼型为研究对象,采用三维可压缩Navier-Stokes方程并结合γ-Re_θt转捩模式计算了展向无限长后掠机翼的基本流场.由于原始γ-Re_θt模式只能预测流向边界层转捩,因此在原始转捩模式中添加横流间歇因子项,进而对复杂构型进行横流不稳定性转捩预测.计算结果显示,利用改进后γ-Re_θt转捩模式预测得到的后掠翼型的转捩位置与实验数据吻合较好,证明了修正的转换模式的合理性和实用性.      

激波与层流/湍流边界层相互作用实验研究

在超声速风洞中,分别对层流和湍流来流条件下的边界层和斜激波(激波强度足以引起流动分离)相互干扰进行了实验研究,利用纳米粒子示踪平面激光散射(NPLS)技术获得了两种条件下流场的精细结构图像;利用粒子图像测速(PIV)技术获得了两种条件下流场的速度场和涡量场;综合运用NPLS结果和PIV结果对比分析了两种流动的瞬时流动结构和时间相关性,实验结果表明:层流边界层内的分离区呈现出狭长的条状,而湍流边界层内分离区呈现出较规则的椭圆;在入射激波上游距入射点较远的位置,层流边界层外围拟序结构会诱导出一系列压缩波系,进而汇聚成空间位置不稳定的诱导激波,而湍流边界层则是在入射激波上游较近的地方直接形成较强且稳定的诱导激波;在入射激波下游,层流边界层内的膨胀区域较小且急促,膨胀后产生的再附激波很弱,而湍流边界层内的膨胀区域较大,膨胀后产生的激波较强。     

湍流强度对风力机翼型气动性能的影响

本文通过采用Transition-SST湍流模型对UMY02-T01-26风电机组专用翼型绕流流场的数值计算,探究了湍流强度对风力机翼型气动性能的影响。结果表明,随着湍流强度的提高,翼型升力系数由前缘失速转变为混合失速。在一定的攻角范围内,升力系数略有增大。对于攻角处于升力系数非线性增长区域范围内,湍流强度的增大导致翼型壁面最大负压值增大。当湍流强度变化时,其壁面上出现层流分离泡的位置大小随之发生变化。此外,本文通过流场分析进一步确定了层流分离泡的产生与变化。     

Self-sustained oscillations of a shock wave interacting with a boundary layer on a supercritical airfoil

A theory is proposed of the self-sustaining oscillations of a weak shock on an airfoil in steady, transonic flow. The interaction of the shock with the boundary layer on the airfoil produces displacement thickness fluctuations which convect downstream and generate sound by interaction with the trailing edge. A feedback loop is established when this sound impinges on the shock wave, resulting in the production of further fluctuations in the displacement thickness. The details are worked out for an idealized mean boundary layer velocity profile, but strong support for the basic hypotheses of the theory is provided by a comparison with recent experiments involving the generation of acoustic “tone bursts” by a supercritical airfoil section.     

新型可压缩尺度在跨声速流动中的应用

基于von Karman长度尺度和新型Reynolds应力本构关系对k-ε湍流模型重构,将k方程封闭,采用代数形式对湍流耗散项进行模化.在KDO（kinetic dependent only）模型的基础上,引入可压缩von Karman长度尺度,得到一种适用于复杂可压缩流动的新型湍流模型CKDO（compressible kinetic dependent only）,在CKDO模型中没有任何经验系数,仅有两个来自边界层精细化标定的可调参数.对RAE2822翼型、轴对称圆筒管道凸起流动、ONERA-M6机翼跨声速流动等算例进行数值计算,结果显示CKDO湍流模型对上述算例流场的压力系数模拟结果与实验值吻合较好,表明CKDO模型能够对跨声速流场进行较为准确的模拟.      

Features of the generation of a collisionless electrostatic shock wave in a laser-ablation plasma

The appearance of a density bump is experimentally revealed in an electrostatic shock wave during the ablation of an aluminum foil by a femtosecond laser pulse. The numerical simulation shows that this phenomenon can be explained by the generation of a packet of ion acoustic waves under the action of high-energy electron flows in a collisionless plasma. It is found that, for the formation and maintenance of the dense plasma layer in the shock wave, the contributions of accelerated ions overtaking it and wave-captured ions of the background plasma formed by a nanosecond laser prepulse in the process of ablation are significant.     

Numerical Study of Air Flow Induced by Shock Impact on an Array of Perforated Plates

This study is focused on the propagation behavior and attenuation characteristics of a planar incident shock wave when propagating through an array of perforated plates. Based on a density-based coupled explicit algorithm, combined with a third-order MUSCL scheme and the Roe averaged flux difference splitting method, the Navier–Stokes equations and the realizable k-ε turbulence model equations describing the air flow are numerically solved. The evolution of the dynamic wave and ring vortex systems is effectively captured and analyzed. The influence of incident shock Mach number, perforated-plate porosity, and plate number on the propagation and attenuation of the shock wave was studied by using pressure- and entropy-based attenuation rates. The results indicate that the reflection, diffraction, transmission, and interference behaviors of the leading shock wave and the superimposed effects due to the trailing secondary shock wave are the main reasons that cause the intensity of the leading shock wave to experience a complex process consisting of attenuation, local enhancement, attenuation, enhancement, and attenuation. The reflected shock interactions with transmitted shock induced ring vortices and jets lead to the deformation and local intensification of the shock wave. The formation of nearly steady jets following the array of perforated plates is attributed to the generation of an oscillation chamber for the inside dynamic wave system between two perforated plates. The vorticity diffusion, merging and splitting of vortex cores dissipate the wave energy. Furthermore, the leading transmitted shock wave attenuates more significantly whereas the reflected shock wave from the first plate of the array attenuates less significantly as the shock Mach number increases. The increase in the porosity weakens the suppression effects on the leading shock wave while increases the attenuation rate of the reflected shock wave. The first perforated plate in the array plays a major role in the attenuation of the shock wave.     

Gas-dynamic mechanism of the effect of pulsed-periodic energy supply on shock-wave structure of transonic flow around airfoils

Gas-dynamic mechanisms of the shock-wave structure restructuring is investigated in a transonic flow around an airfoil at a pulsed-periodic energy supply. It is shown that it does not reduce only to the effects of the low-density wake forming behind the energy source on the shock wave but represents a combination of several factors, of which the “explosion” character of the disturbance development is determining.     

Experimental Research on Wet Steam Flow With Shock Wave

The article presents research on a wet steam transonic flow in a Laval half-nozzle with a shock wave. The motivation for this research was to investigate the shock wave/liquid phase interaction in the transonic wet steam flow. This phenomenon is responsible for the lack of good correspondence between experimental data and computational fluid dynamics results. For the tests, the geometry of the half-arc nozzle was used. The shock wave formation at the divergent section of the nozzle is caused by too high back-pressure. The observed instabilities in the flow are mainly initiated by the shock wave/liquid film (boundary layer) interaction. The numerical calculations were compared with experimental results with respect to the static pressure distribution along the nozzle.     

激波冲击火焰的流动拓扑研究

为深入研究激波冲击火焰现象的内在机制,采用二维带化学反应的Navier-Stokes方程对现象进行数值研究,通过对速度梯度张量特征方程的分析证明Okubo-Weiss函数适用于可压缩流动,并重点分析火焰区的流动拓扑特性.结果表明,波后火焰区内Okubo-Weiss函数积分量基本守恒,但在火焰区内部和表面具有截然不同的流动状态,且火焰发展基本不受流场可压缩性的影响;波后火焰区的流动拓扑分类主要以焦点和鞍点为主,意味着流场中变形占主导.     

Nonlinear effects of pulsed periodic energy supply on shock wave structure of transonic flow around airfoils

The work deals with an investigation of possibilities of controlling the aerodynamic characteristics of airfoils with the aid of external local pulsed-periodic energy supply at transonic flight regimes. The alteration of flow structure near a symmetric airfoil and its wave drag has been studied on the basis of a numerical solution of two-dimensional nonstationary gas dynamics equations versus the energy supply period, the localization and shape of the energy supply zone. The energy supply upstream of the closing shock in the immediate proximity of the contour in the zones extended along it is found to result in a considerable reduction of the profile wave drag. The nature of such a drag reduction is elucidated. The existence of a limiting frequency of energy supply is established.