LARGE EDDY SIMULATION ON THE EVOLUTION AND THE FAST-TIME PREDICTION OF AIRCRAFT WAKE VORTICES

随着我国人民生活水平的提高,航空运输的重要性与日俱增,航班延误问题也日益严重.尾流间隔（保障后机不受前机尾流影响的最小安全间隔）是制约机场效率的关键因素.针对这一工程应用问题,采用大涡模拟方法研究飞机尾涡在大气中的演变特性.研究工作首先发展了飞机尾涡演变的大涡模拟方法,将自适应网格技术应用于飞机尾涡演变的大涡模拟,大幅减少所需的网格量,提高计算效率.提出了升力面尾涡生成方法,在不增加计算量的情况下实现了尾涡卷起过程和远场衰减的组合模拟.在系列算例分析研究基础上,创建了基于大涡模拟计算结果的尾流间隔快速预测系统.该系统可以根据实时大气风场和进出港的前后飞机机型,快速预测并输出所需的尾流间隔.经过与场地测试数据比较表明,在北京市2014年的平均风速条件下,本系统预测的尾流间隔可在现有标准基础上缩减7%~50%,能够有效提高机场容量.     

Spatial perturbation of a wing-tip vortex using pulsed span-wise jets

The separation distance required between transport aircraft to avoid wake vortices remains a limiting factor on airport capacity. The dissipation of the wake can be accelerated by perturbing co-operative instabilities between multiple pairs of vortices. This paper presents the results of a preliminary experimental investigation into the use of pulsed span-wise air jets in the wing tip to perturb a single tip vortex in the very near field. Velocity measurements were made using PIV and hot-wire anemometry. The results demonstrate that the vortex position can be modulated at frequencies up to 50 Hz and, as such, the method shows promise for forcing instability in multiple vortex wakes.     

水下发射航行体尾涡不稳定性分析

在水下连续发射过程中前一发航行体尾流会对后一发航行体运动姿态稳定性产生流动干扰现象. 因此, 研究尾流中涡旋结构演变机理对解决多弹体水下连续发射流动干扰难题具有重要的意义. 本文采用改进型分离涡模型与能量方程, VOF多相流模型与重叠网格技术相结合方法, 对航行体水下发射尾流演变过程开展精细化模拟研究, 其中模拟结果和实验吻合度较好, 验证了本文数值方法的有效性. 以航行体尾流区域为重点研究对象, 分析了尾流区瞬态流场分布, 讨论了横流强度和雷诺数对尾涡结构演变以及脉动压力分布特性的影响. 结果表明: 由于尾流区高速流体核心区与低速自由流相互作用导致Kelvin-Helmholtz不稳定现象出现, 可以清晰地发现涡旋结构在剪切力的作用下发生脱落. 在横流条件下, 航行体尾端脱落的涡环与涡腿形成发卡涡, 而多个发卡涡沿轴向间隔排列组成发卡涡包存在于尾流中. 随着横流强度增大, 形成多级发卡涡包结构, 而导致脉动压力二次峰值均出现的主要原因是尾流涡旋流场演变引起的. 随着雷诺数的增大, 尾流中由圆柱形涡和U型涡组成的二次涡结构逐渐明显, 不稳定性加强.      

Numerical study on wake characteristics of high-speed trains

Intensive turbulence exists in the wakes of high speed trains, and the aerodynamic performance of the trailing car could deteriorate rapidly due to complicated features of the vortices in the wake zone. As a result, the safety and amenity of high speed trains would face a great challenge. This paper considers mainly the mechanism of vortex formation and evolution in the train flow field. A real CRH2 model is studied, with a leading car, a middle car and a trailing car included. Different running speeds and cross wind conditions are considered, and the approaches of un- steady Reynold-averaged Navier-Stokes （URANS） and de- tached eddy simulation （DES） are utilized, respectively. Re- suits reveal that DES has better capability of capturing small eddies compared to URANS. However, for large eddies, the effects of two approaches are almost the same. In conditions without cross winds, two large vortex streets stretch from the train nose and interact strongly with each other in the wake zone. With the reinforcement of the ground, a complicated wake vortex system generates and becomes strengthened as the running speed increases. However, the locations of flow separations on the train surface and the separation mechanism keep unchanged. In conditions with cross winds, three large vortices develop along the leeward side of the train, among which the weakest one has no obvious influence on the wake flow while the other two stretch to the tail of the train and combine with the helical vortices in the train wake. Thus, optimization of the aerodynamic performance of the trailing car should be aiming at reducing the intensity of the wake vortex system.     

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.      

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

Natural and forced evolution of a counter rotating vortex pair

Naturally occurring instabilities such as the long wave Crow instability are one of the mechanisms that can bring about the break-up of aircraft wake vortex pairs. Encouraging the early onset of such instabilities by active or passive means offers the possibility of accelerating wake break-up and dissipation with favourable consequences in terms of safe aircraft separation and airport capacity. This paper describes an experimental investigation of the evolution of the Crow instability in a counter rotating vortex pair perturbed by a pulsed air sheet. Flow measurements were made using two-dimensional particle image velocimetry and two methods devised to detect wake break-up from the two-dimensional data. Pulsed excitation at a wavelength within the range amplified by the Crow instability resulted in wake break-up and the formation of vortex rings at a downstream distance reduced by up to 37% compared to the baseline case. For excitation at a wavelength outside the excited range early break-up was also observed but the extent of the accelerated break-up was reduced as the amplitude of the excitation was increased.     

The Vortical Structures in the Rear Separation and Wake Produced by a Supersonic Micro-Ramp

The vortical structures in the rear separation and wake region produced by a micro-ramp that immersed in a supersonic turbulent boundary layer are investigated. The small scale separation close to the trailing edge was revealed and this confirms the previous experimental observation. Between the reverse region and surrounding fast moving flow, a three-dimensional shear layer was formed, and vortices are generated. By using vortex line method, the spiral points were illustrated as the cross-sections of the Ω-shaped vortices that follow the shape of the separation. The vortical structure was analogous to that in the wake region, where similar Ω-shaped vortex which follows the deficit region caused by the micro-ramp can be observed. Finally, the revealed flow topology was conceived new and beneficial to the studying of wall bounded turbulence which involves similar vortical structures but in a smaller scale, compared with the vortical pattern in the current micro-ramp wake.     

Effect of homogeneous stratification of the fluid on the dynamics of a kármán vortex street behind a circular cylinder

The experimental results of studying the effect of homogeneous stratification of the fluid on the conditions of generation of a Kárman vortex street [1] developing in the wake of a cylinder in steady horizontal motion are described. In a homogeneous medium at Reynolds numbers Re >5 two symmetrical regions of vorticity of opposite sign are formed behind the cylinder and move together with the latter. As the speed of the cylinder increases, the link between the vortices and the cylinder grows weaker, the vortices are stretched out along the flow and at Re > 40 begin to separate alternately, forming a vortex street in the wake. At first, the frequency of vortex separation increases sharply with increase in Re, but then levels off. It is found that in a uniformly stratified fluid the onset of vortex separation from the moving cylinder is delayed. The dependence of the critical Reynolds number (onset of vortex separation) on the internal Froude number is obtained. The effect of stratification of the fluid on the frequency of separation of the vortices in the Kármán street is investigated. The effect of the Froude number on the dependence of the Strouhal number on the Reynolds number is established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 83–86, January–February, 1986.In conclusion the authors wish to thank A. T. Onufriev for his interest in their work and useful discussions of the results.     

Analysis of tornado and near-ground turbulence using a hybrid meteorological model/engineering LES method

A new hybrid meteorological model/engineering LES method was used to analyze a tornado and near-ground turbulence under realistic conditions. The hybrid method estimates high frequency turbulence regeneration in an actual severe storm that has complex thermal conditions in order to alleviate the problem of high frequency component dissipation in meteorological models and generate appropriate field data to investigate near-ground flow and wind loading on buildings. The high frequency regeneration method was validated in the convective boundary layer and the energy spectrum of the velocity field was successfully extended to higher frequency region in good agreement with −5/3 Kolmogorov law.Multi-scale simulations of a tornado were conducted based on meteorological model outputs with three different terrain conditions, including an actual urban building geometry and terrain, and the changes in the main tornado vortex and the interactions with near-ground turbulent field were examined. On roughened surfaces with uniform blocks and buildings, the main tornado vortex was connected to near-ground vertical vortices arising due to the urban geometry. The structure was significantly deformed, splitting into numerous finer vortices below several times the building height. These near-ground vortices create localized sharp pressure drop patterns and also affect the movement of the tornado. These effects of low pressure on buildings were quantitatively analyzed.     

斑马鱼直线加速游动试验研究

海洋生物低噪音、高速、髙效游动能力是任何人造水下航行器所无法比拟的. 借助时间解析粒子图像测速技术对斑马鱼直线加速游动过程进行精细流场测量, 对其运动学行为特性和动力学机理进行分析. 同时应用双正交分解对涡量场进行模态分解, 获取流场的时间演化和空间分布特征. 从流动机理的角度探究斑马鱼游动过程的流动结构特征及旋涡动态演化特性. 试验结果表明: 流动可视化展现了整体涡流尾迹的结构分布,方便探究运动特性与旋涡尾迹之间的耦合关系. 斑马鱼从运动开始时体干保持着鲹科式的运动规律, 游动时的动能主要由前几次大幅的摆尾过程提供, 后续的摆尾主要调整方向及姿态. 两次不同方向的摆尾动作会形成一对方向相反的旋涡, 并在时序下旋涡逐渐脱落. 同时尾流的涡量变化在一定程度上反映鱼体的游向的变化. 基于双正交分解分解后的时间演化结果验证本次试验在时间上涡量场具有合理的恒定幅度, 空间分布表明低阶空间模态表征斑马鱼游动的主要涡流动结构, 高阶空间模态表征涡流动的细节结构. 研究鱼类游动时的摆尾推进机制与动力学特性能够为高效率的仿鱼类推进装置设计提供一定科学参考.      

Large eddy simulation of aircraft wake vortex with self-adaptive grid method

A self-adaptive-grid method is applied to numerical simulation of the evolution of aircraft wake vortex with the large eddy simulation(LES). The Idaho Falls(IDF)measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V¨ais¨al¨a(BV) frequencies are computed with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, secondary structures such as baroclinic vorticity and helical structures are also simulated.Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.     

横摆角下汽车尾部湍流特征量的PIV试验分析

汽车尾部湍流场是汽车压差阻力的主要来源,在HD-2汽车模型风洞中,首先使用测力天平和测压系统,对横摆角工况下汽车模型的气动六分力和纵对称截面48个测点的表面压力进行了测量,然后利用PIV测量技术对模型在横摆角分别为0°、15°的尾部湍流场进行了测量,获得该模型尾流场的速度场、涡量场和雷诺应力流场信息,通过计算得出尾流场区域空间相关系数和湍流积分尺度。结果表明:在横摆角工况下,汽车模型尾部涡流的结构呈现向上发展的趋势;尾流场拖拽涡的范围和强度的增大导致了模型气动力出现较大的增加;湍流积分尺度的变化表明,尾部涡流区的分离噪声与涡流分离位置有关,在汽车尾部造型设计中,要尽量推迟尾部涡流的分离。     

舵几何特征对桨-舵系统尾流场演化的影响

发生在桨和舵之间的干扰会影响螺旋桨尾流的演化,导致尾流场中的湍流在下游增强,恶化船舶的振动和噪声性能,深入分析舵几何参数对桨-舵系统尾流场演化的影响能够为推进器尾流场的调节和减振降噪提供新思路.因此,从弦长、剖面和梯形舵入手分析不同的舵几何参数对螺旋桨尾流场演化特性的影响,使用大漩涡模拟方法模拟流场中的湍流结构,对不同舵弦长、剖面下的螺旋桨尾涡结构演化进行了分析,在舵弦长、剖面影响螺旋桨尾流场演化的研究的基础上分析了梯形舵对螺旋桨尾涡结构的影响,进一步分析了梯形舵影响下的螺旋桨尾流场中湍动能的分布.结果表明舵的弦长和剖面均会影响螺旋桨尾流场的演化,这种影响表现为更大的弦长和更厚的剖面会促进螺旋桨梢涡在舵压力面上的偏移,更薄的舵剖面会带来更强烈的螺旋桨毂涡偏移;涡管轮廓和舵表面脉动压力的对比均表明梯形舵会促进螺旋桨尾流场沿逆舵梯度方向偏移,从而导致螺旋桨的尾涡结构在舵两侧及下游呈现不对称分布,桨-舵系统下游的湍流结构与螺旋桨尾涡-舵碰撞过程、螺旋桨尾涡-舵随边涡干扰过程、螺旋桨梢涡-螺旋桨毂涡干扰有关,偏移更大的螺旋桨尾涡结构会在尾流场中更早地引起湍动能增强.     

A comprehensive PIV measurement campaign on a fully equipped helicopter model

The flow field around a helicopter is characterised by its inherent complexity including effects of fluid?Cstructure interference, shock?Cboundary layer interaction, and dynamic stall. Since the advancement of computational fluid dynamics and computing capabilities has led to an increasing demand for experimental validation data, a comprehensive wind tunnel test campaign of a fully equipped and motorised generic medium transport helicopter was conducted in the framework of the GOAHEAD project. Different model configurations (with or without main/tail rotor blades) and several flight conditions were investigated. In this paper, the results of the three-component velocity field measurements around the model are surveyed. The effect of the interaction between the main rotor wake and the fuselage for cruise/tail shake flight conditions was analysed based on the flow characteristics downstream from the rotor hub and the rear fuselage hatch. The results indicated a sensible increment of the intensity of the vortex shedding from the lower part of the fuselage and a strong interaction between the blade vortex filaments and the wakes shed by the rotor hub and by the engine exhaust areas. The pitch-up phenomenon was addressed, detecting the blade tip vortices impacting on the horizontal tail plane. For high-speed forward flight, the shock wave formation on the advancing blade was detected, measuring the location on the blade chord and the intensity. Furthermore, dynamic stall on the retreating main rotor blade in high-speed forward flight was observed at r/R?=?0.5 and 0.6. The analysis of the substructures forming the dynamic stall vortex revealed an unexpected spatial concentration suggesting a rotational stabilisation of large-scale structures on the blade.     

圆柱尾迹涡的三维演化及结构特征

应用无粘涡丝运动学理论和局部诱导近似(LIA)方法,以Lagrange观点跟踪涡丝在背景流场中运动,用数值方法研究了中等Re数(≈10~3)下圆柱分离尾迹中Kármán涡和涡辫区涡丝的三维演化的机制和动力学过程,及其结构特征。背景流场考虑为尾迹时间平均速度流场和Kármán涡街流场。初始展向小扰动为指数形式和谐波形式。结果指出:Kármán涡和涡辫区中的涡丝具有展向不稳定性,形成流向涡量。在尾迹的初期输运过程中,表现出有序的大、小尺度涡结构。并进一步分析了其产生的机理。     

Hydrodynamic instability of separated viscous flow around a circular cylinder

A uniform viscous flow around a circular cylinder is studied numerically in the Reynolds number range from 0 to 500. It is shown that the existence and the basic properties of self-oscillating regimes are specified by the evolution of their hydrodynamic instability. It is found that the vortex formation in a near wake is associated with the separation zone dynamics in the main flow. The values of critical Reynolds numbers for the four successive bifurcations of the self-oscillating regimes of flow are obtained. An interpretation of experimental data on the vortices in the near wake is discussed.     

螺旋桨梢涡不稳定性机理与演化模型研究

螺旋桨尾流场的涡流特性是一个基础但又十分复杂的流体力学问题, 它的复杂性源于其蕴含复杂的漩涡系统, 且该漩涡系统会在高速的剪切层流动中不断演化, 其流体动力学行为, 如由稳定态演变为不稳定态的机理以及复杂工况环境中的流动现象, 一直是流体力学领域的难点和备受关注的热点问题. 从工程应用的角度看, 桨后梢涡的演化特性与船舶结构物的宏观特性直接相关, 更好地理解多工况下螺旋桨尾流的动力学特性, 将有助于改善与振动、噪声以及结构问题等相关的推进器性能, 对综合性能优良的下一代螺旋桨的设计和优化有着重要的现实意义. 本文基于延迟分离涡模拟、大涡模拟和无湍流模型模拟方法以及粒子图像测速流场测试分别开展了螺旋桨尾流动力学特性的数值与试验研究, 对螺旋桨尾流不稳定性的触发机理进行了揭示. 基于均匀来流中螺旋桨梢涡的演化机理, 提出了螺旋桨梢涡演化模型. 该模型能够较为准确地模拟螺旋桨梢涡的演化过程, 预测螺旋桨梢涡融合的时间和位置, 对螺旋桨流噪声预报和控制以及性能优良的螺旋桨设计具有重要意义.      

Large-Eddy Simulations of Forced Three-Dimensional Impinging Jets

Large-eddy simulations of the flow field around twin three-dimensional impinging jets were carried out to simulate the near-ground hover configuration of a vertical takeoff and landing (VTOL) aircraft. Both the impinging jet and the upwash caused by the collision of the wall jets are modeled in this study. The evolution of the vortical structures in the impinging jet flow field, due to the introduction of axisymmetric and azimuthal perturbations at the jet exit, has been investigated. The vortical structures formed in the jet shear layer due to azimuthal forcing, show significant three-dimensional vortex stretching effects when compared to the structures formed during axisymmetric forcing. Breakdown of the large-scale structures into smaller vortices also occurs much earlier during azimuthal forcing. When compared to the upwash formed during axisymmetric forcing, the azimuthally perturbed jet forms an upwash that is less coherent and results in a weaker upload or lift-off force on the aircraft undersurface. Comparison with available experimental data indicates good agreement for the centerline velocity decay, the wall pressure variation and the phase speed of the vortical structures.     

Flow control over a NACA 0012 airfoil using dielectric-barrier-discharge plasma actuator with a Gurney flap

Flow control study of a NACA 0012 airfoil with a Gurney flap was carried out in a wind tunnel, where it was demonstrated that a dielectric-barrier-discharge (DBD) plasma actuator attached to the flap could increase the lift further, but with a small drag penalty. Time-resolved PIV measurements of the near-wake region indicated that the plasma forcing shifted the wake downwards, reducing its recirculation length. Analysis of wake vortex dynamics suggested that the plasma actuator initially amplified the lower wake shear layer by adding momentum along the downstream surface of the Gurney flap. This enhanced mutual entrainment between the upper and lower wake vortices, leading to an increase in lift on the airfoil.