圆角化对方柱气动性能影响的流场机理

方形截面柱体的圆角化处理是常用的流动控制方法,但其流场作用机理尚未被澄清.采用大涡模拟方法,在雷诺数为2.2×104时,考虑风攻角的影响,对均匀流作用下的标准方柱和圆角方柱的气动性能和流场特性进行了研究,定量分析了圆角化气动措施和风攻角变化对分离泡特性的影响规律,从流场角度澄清了圆角化气动措施对方柱气动性能的影响机理.研究表明:与标准方柱相比,圆角方柱的表面风压、气动力和涡脱强度呈整体下降的趋势,但圆角方柱的斯特劳哈尔数更高;圆角方柱的"分离泡流态"发生在更小的风攻角范围内,分离泡的出现会进一步造成方柱的尾流变窄,涡脱强度减弱;随着风攻角的增大,分离泡的长度会逐渐减小直至消失,分离泡的中心会逐渐向方柱前角(迎风向)和方柱壁面移动;与标准方柱相比,圆角方柱的气流发生初次分离的位置向下游移动,分离后的剪切层更贴近方柱,因而更易发生再附现象;方柱尾流宽度的减小和涡脱强度的减弱是导致圆角方柱气动力减小和斯特劳哈尔数增大的主要原因.     

带有柔性薄板三维方柱流固耦合的数值模拟研究

采用双向流固耦合方法,对带有柔性薄板三维方柱的流场变化特性进行了研究。通过对比单方柱,分析了带有柔性薄板三维方柱阻力系数、升力系数以及斯特劳哈尔数的变化规律。研究表明,在方柱尾流区域附加一柔性薄板可以使其阻力系数降低34.6%,同时其变化幅值大大减小;其升力系数的均方根减小84.8%,流场脉动大幅度减小;斯特劳哈尔数降低79.5%。研究结果表明,在三维方柱后设置柔性薄板可有效抑制涡脱落,从而改善三维方柱的尾流特性。     

圆角化对受迫振动方柱绕流特性影响机理分析

为分析圆角化对低雷诺数下受迫振动方柱绕流特性的影响机理, 对Ansys Fluent软件进行二次开发, 即通过用户自定义函数中的DEFINE_ CG_MOTION宏对柱体周期性受迫振动的函数进行编程, 并对流场计算域进行区域划分以便利用动网格技术中动态层法实现柱体受迫振动, 从而实现对受迫振动柱体绕流流场的流固耦合模拟.在雷诺数Re = 200时, 考虑方柱截面不同圆角的影响, 对均匀流作用下5种圆角化r/D = 1/2, 1/4, 1/5, 1/8和0受迫振动方柱的绕流进行数值模拟, 分析了这5种参数下受迫振动方柱的升阻力系数、尾流涡量和锁定区间的变化规律, 澄清了圆角化对受迫振动方柱稳定性的影响机理.研究表明: 与尖角方柱相比, 圆角化方柱升阻力系数有了明显的减小, 且升力、阻力系数随圆角增大而减小; 低振幅比下圆角方柱的涡旋脱落模式均为2S模态, 涡旋尾迹变窄; 锁定区间范围基本关于F = 1对称, 锁定区间的变化趋势与圆柱类似.      

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

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

低雷诺数串列双方柱流致振动质量比效应的数值研究

为澄清串列双方柱流致振动的质量比效应, 采用数值模拟方法, 在雷诺数为150时, 研究了质量比($m^{\ast }=3$, 10, 20)对下游方柱振动响应特性的影响规律, 分析了下游方柱尾流模态的演变过程, 探讨了导致下游方柱振动的流固耦合机制. 结果表明: 质量比对下游方柱的流致振动有重要影响, 低质量比($m^{\ast }=3$)时下游方柱的振动响应更为复杂, 随着折减速度的增大, 下游方柱并未出现传统“锁定”现象(即振动频率比$f_{y}$/$f_{\rm n} \approx1$的锁定), 而发生了“弱锁定”现象(即$f_{y}/f_{\rm n}<1$的锁定); 随着质量比的增加($m^{\ast }=10$和20), “弱锁定”现象消失, 而出现传统“锁定”现象, 且下游方柱横流向最大振幅减小. 质量比对串列双方柱的柱心间距有明显影响, 低质量比($m^{\ast }=3$)时的柱间距在振动锁定区内会急剧减小, 而较高质量比($m^{\ast }=10$和20)下的柱间距则变化不大. 此外, 质量比对串列双方柱的尾流模态和流固耦合机制也有显著影响, 其中低质量比($m^{\ast }=3$)下的情况更为多样.      

可渗透壁对翼型声涡相互作用下流场及声场的影响

由仿生学原理构建的可渗透翼型对湍流气动噪声抑制作用已展现良好的应用前景。对NACA 0012可渗透翼型和实体翼型进行了数值计算,得到了声涡相互作用下气动噪声声场和流场,分析了可渗透壁对翼型流场和声场的影响。研究表明,相对实体翼型,可渗透壁通过减小声源强度降低了主纯音噪声声压级幅值和远场总声压级,消除了高阶离散纯音,但对噪声的指向性没有较大改变。进一步的流场分析表明,可渗透壁对翼型气动性能影响不大的情况下能够降低边界层扰动和翼型后缘大尺度涡旋强度,并推迟分离泡转捩和再附位置。     

基于浸入边界算法的振动钝体绕流模拟研究

传统CFD方法在振动钝体绕流计算中常借助动网格技术,网格再生任务繁重。针对于此,本文利用可在静止网格中计算动边界绕流问题的浸入边界算法(IBM),编写数值模拟程序,分别对竖向强迫正弦振动方柱(Re=UD/v=103、振幅恒定、振动频率变化)以及桥梁断面(Re=UB/v=7.5×103、振幅、振动频率均变化)展开气动特性和流场特征结构分析。初步研究结果表明,振幅恒定为方柱高度的14%时,其涡脱锁定区长度为0.06~0.2,锁定区后端(Stc0.2)振动方柱涡脱频率回归静止涡脱频率;不同振幅下的桥梁断面阻力系数均在静止涡脱频率处产生峰值,桥梁断面升力系数则在此处均出现归零效应,且振幅越大,归零效应愈明显。     

钝体尾流控制机理及方法研究进展

首先从涡脱落生成理论出发对钝体尾流控制方法进行了分类,并简单介绍了国内尾流控制研究情况. 之后介绍了我们用窄条或小方柱取代小圆柱后,对Strykowsky和Sreenivasan 控制方法的改进及其在高雷诺数下对圆柱和方柱尾流涡脱落的有效抑制情况, 并探讨了控制件钝度对抑制效果的影响.第3部分用实验数据对各个涡脱落生成模型做了分析与检验, 指出控制件方法的机理与改变钝体分离位置、减小钝体背压吸力、改变流动的展向相关性、 防止钝体两侧剪切层相互作用等无关,而与钝体近尾流速度剖面的局部修正及其稳定性的改变有关. 最后简单介绍了控制件方法今后研究工作展望及其工程应用前景.     

风力机翼型动态失速等离子体流动控制数值研究

针对动态失速引起的风力机翼型气动性能恶化的问题,本文基于动网格和滑移网格技术, 开展了大涡模拟数值计算研究,探索了非定常脉冲等离子体的动态流动控制机理. 结果表明,等离子体气动激励能够有效控制翼型动态失速, 改善平均和瞬态气动力,减小力矩负峰值和迟滞环面积. 压力分布在等离子体施加范围内出现了负压"凸起",上翼面吸力峰值明显增大.脉冲频率和占空比这两个非定常控制参数对流动控制影响显著,无因次脉冲频率为1.5时等离子体控制效果较好,占空比为0.8时即可接近连续工作模式下的气动收益. 翼型深失速状态,等离子体促使流动分离位置明显向后缘移动, 抵抗了大尺度动态失速涡的发生,分离涡结构破碎耗散、重新附着, 涡流影响范围减小; 浅失速状态,等离子体激励具有较强的剪切层操纵能力, 诱导了翼型边界层提前转捩,促进了与主流的动量掺混. 等离子体气动激励诱导出前缘附近贴体翼面"涡簇",起到了虚拟气动外形的作用.不同尺度、频域的动态涡结构与等离子体气动激励的非线性、强耦合作用导致了气动力/力矩的谐波振荡.      

翼型跨声速气动特性的不确定性及全局灵敏度分析

针对马赫数和仰角的随机不确定性会导致气动性能波动的现象, 采用非嵌入式的混沌多项式方法对绕NACA0012 翼型跨声速随机气动特性进行不确定性及全局灵敏度分析. 具体分析了飞行状态的不确定性对气动载荷分布、流场及气动力系数的影响并通过全局灵敏度分析找出重要因素. 不确定性分析结果表明翼型上表面的激波以及激波后分离泡是造成气动性能剧烈波动的主要原因. 灵敏度分析结果表明在跨声速区域马赫数对激波处气动性能影响最大, 此外, 虽然马赫数和仰角相互耦合作用对气动力系数贡献比较小, 但对于激波位置处的流场, 这种互耦合作用不可忽略.      

Parametric study of separation and transition characteristics over an airfoil at low Reynolds numbers

Time-resolved surface pressure measurements are used to experimentally investigate characteristics of separation and transition over a NACA 0018 airfoil for the relatively wide range of chord Reynolds numbers from 50,000 to 250,000 and angles of attack from 0° to 21°. The results provide a comprehensive data set of characteristic parameters for separated shear layer development and reveal important dependencies of these quantities on flow conditions. Mean surface pressure measurements are used to explore the variation in separation bubble position, edge velocity in the separated shear layer, and lift coefficients with angle of attack and Reynolds number. Consistent with previous studies, the separation bubble is found to move upstream and decrease in length as the Reynolds number and angle of attack increase. Above a certain angle of attack, the proximity of the separation bubble to the location of the suction peak results in a reduced lift slope compared to that observed at lower angles. Simultaneous measurements of the time-varying component of surface pressure at various spatial locations on the model are used to estimate the frequency of shear layer instability, maximum root-mean-square (RMS) surface pressure, spatial amplification rates of RMS surface pressure, and convection speeds of the pressure fluctuations in the separation bubble. A power-law correlation between the shear layer instability frequency and Reynolds number is shown to provide an order of magnitude estimate of the central frequency of disturbance amplification for various airfoil geometries at low Reynolds numbers. Maximum RMS surface pressures are found to agree with values measured in separation bubbles over geometries other than airfoils, when normalized by the dynamic pressure based on edge velocity. Spatial amplification rates in the separation bubble increase with both Reynolds number and angle of attack, causing the accompanying decrease in separation bubble length. Values of the convection speed of pressure fluctuations in the separated shear layer are measured to be between 35 and 50% of the edge velocity, consistent with predictions of linear stability theory for separated shear layers.     

翼型绕流的电磁力控制

将表面包覆电磁激活板的翼型,按一定的攻角,置于流动的弱电介质溶液中,电磁激活板可产生作用于流体的切向电磁力(Lorentz力),从而改变流体边界层的结构. 在转动水槽中,对翼型绕流及电磁力控制下的绕流形态进行了实验研究. 结果表明,未加电磁力时,前缘涡的脱落点是不确定的,与流场具体条件有关,而后缘涡仅在尖角处脱落. 前缘涡与后缘涡相互影响,并周期性的脱体,在尾部形成涡街. 施加电磁力后,当力的方向与流动方向相同时,可以在一定程度上抑制分离,消除涡街,其效果与减小攻角类似. 加反向电磁力时,则相当于加大攻角,在翼型体的背风面形成涡街.      

Effect of buoyancy on the wakes of circular and square cylinders: a schlieren-interferometric study

Wakes behind heated cylinders, circular, and square have been experimentally investigated at low-Reynolds numbers. The electrically heated cylinder is mounted in a vertical airflow facility such that buoyancy aids the inertia of main flow. The operating parameters, i.e., Reynolds number and Richardson number are varied to examine flow behavior over a range of experimental conditions from forced to mixed convection regime. Laser schlieren-interferometry has been used for visualization and analysis of flow structures. Complete vortex shedding sequence has been recorded using a high-speed camera. The results on detailed dynamical characteristics of vortical structures, i.e., their size, shape and phase, Strouhal number, power spectra, convection velocity, phase shift, vortex inception length, and fluctuations are reported. On heating, alteration of organized (coherent) structures with respect to shape, size and their movement is readily perceived from instantaneous Schlieren images before they reduce to a steady plume. For both cylinders, Strouhal number shows a slow increase with an increase in Richardson number. At a critical value, there is complete disappearance of vortex shedding and a drop in Strouhal number to zero. The corresponding spectra evolve from being highly peaked at the vortex shedding frequency to a broadband appearance when vortex shedding is suppressed. The geometry of vortex structures transforms to a slender shape before shedding is suppressed. At this heating level, absence of multiple peaks in power spectra at cylinder centerline indicates absence of interaction between opposite shear layers. The convection velocity of vortices increases in stream wise direction to an asymptotic value and its variation is a function of Richardson number. The convection speed abruptly falls to zero at critical Richardson number. The phase difference of shed vortices between upstream and downstream location increases with an increase in Richardson number. Velocity profiles show an increase in fluid speed and beyond the critical point, buoyancy forces add enough momentum to cancel momentum deficit due to the cylinder. Overall, the combined effect of temperature gradient on the separating shear layer velocity profile in near field and vortical structures interaction in far field influences wake instability of a heated cylinder. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Flow development upstream of a fence

The effect of Reynolds number on the flow development upstream of a rigid, non-porous, static fence is investigated experimentally. The flow field is measured using time-resolved, two-component particle image velocimetry at Reynolds numbers based on fence height of 18000, 36000, and 54000. The results show that a laminar separation bubble forms upstream of the junction vortex at the base of the fence. The mean extent of the bubble decreases with increasing Reynolds number, with mean separation moving downstream and mean reattachment moving upstream. In the aft portion of the bubble, shear layer vortices form and are shed at scaled frequencies and wavelengths that are comparable to laminar separation bubble shedding in low Reynolds number airfoils and flat plates with an imposed adverse pressure gradient. The strong periodicity of the associated coherent structures and the proximity of shear layer roll-up relative to the fence should be taken into consideration in the relevant designs due to potential implications to structural loading. A simple flow separation prediction model combining inviscid fence flow solution with Thwaites’ method is introduced and shows good agreement with the experimental results for the Reynolds number range considered.     

Identification of flow regimes around two staggered square cylinders by a numerical study

The flow over two square cylinders in staggered arrangement is simulated numerically at a fixed Reynolds number (\(Re =150\)) for different gap spacing between cylinders from 0.1 to 6 times a cylinder side to understand the flow structures. The non-inclined square cylinders are located on a line with a staggered angle of \(45^{\circ }\) to the oncoming velocity vector. All numerical simulations are carried out with a finite-volume code based on a collocated grid arrangement. The effects of vortex shedding on the various features of the flow field are numerically visualized using different flow contours such as \(\lambda _{2}\) criterion, vorticity, pressure and magnitudes of velocity to distinguish the distinctive flow patterns. By changing the gap spacing between cylinders, five different flow regimes are identified and classified as single body, periodic gap flow, aperiodic, modulated periodic and synchronized vortex shedding regimes. This study revealed that the observed multiple frequencies in global forces of the downstream cylinder in the modulated periodic regime are more properly associated with differences in vortex shedding frequencies of individual cylinders than individual shear layers reported in some previous works; particularly, both shear layers from the downstream cylinder often shed vortices at the same multiple frequencies. The maximum Strouhal number for the upstream cylinder is also identified at \({G}^{*}=1\) for aperiodic flow pattern. Furthermore, for most cases studied, the downstream cylinder experiences larger drag force than the upstream cylinder.     

Flow around four cylinders arranged in a square configuration

This paper presents an experimental study of the flow around four circular cylinders arranged in a square configuration. The Reynolds number was fixed at Re=8000, the pitch-to-diameter ratio between adjacent cylinders was varied from P/D=2 to 5 and the incidence angle was changed from α=0° (in-line square configuration) to 45° (diamond configuration) at an interval of 7.5°. The flow field was measured using digital Particle Image Velocimetry (PIV) to examine the vortex shedding characteristics of the cylinders, together with direct measurement of fluid dynamic forces (lift and drag) on each cylinder using a piezoelectric load cell. Depending on the pitch ratio, the flow could be broadly classified as shielding regime (P/D≤2), shear layer reattachment regime (2.5≤P/D≤3.5) and vortex impinging regime (P/D≥4). However, this classification is valid only in the case that the cylinder array is arranged nearly in-line with the free stream (α≈0°), because the flow is also sensitive to α. As α increases from 0° to 45°, each cylinder experiences a transition of vortex shedding pattern from a one-frequency mode to a two-frequency mode. The flow interference among the cylinders is complicated, which could be non-synchronous, quasi-periodic or synchronized with a definite phase relationship with other cylinders depending on the combined value of α and P/D. The change in vortex pattern is also reflected by some integral parameters of the flow such as force coefficients, power spectra and Strouhal numbers.     

On vortex structures and processes over bluff bodies in impulsive flow

An experimental system that generates and visualizes unsteady motions in air was employed to visualize vortex developments over bluff bodies in impulsive started flow. The study presents photographic sequences of impulsive flow over circular and square cylinders for flow Reynolds number of 200. A comparison of physical and numerical visualization for impulsive flow over a flat plate at 50° angle of attack and for flow Reynolds number of 1,000 is also presented. The visualization examples reveal the details of vortex separation and subsequent developments and may be utilized as a reference to guide and economize computer visualization efforts on unsteady separated flows.Assistant Professor of Aerospace Engineering     

对称翼型低雷诺数小攻角升力系数非线性现象研究

采用Rogers发展的三阶Roe格式,求解非定常不可压N-S方程,时间方向为二阶精度双时间步方法, 数值模拟了对称翼型SD8020低雷诺数(Re=40000,100000)条件下,流场层流分离涡结构和升力系数随攻角的变化．同试验比较证明了数值模拟的正确性．通过对数值模拟时均化流场结果的详细分析,发现对称翼型在小雷诺数0°攻角附近出现的层流分离泡,其内部结构和演化规律都不同于经典层流分离泡模型,从而提出了一种后缘层流分离泡模型．并应用该模型对对称翼型小攻角低雷诺数流场特性以及升力系数非线性效应的形成机理进行了研究和解释．     

Vortex shedding from a square cylinder in presence of a moving wall

A numerical study on the flow past a square cylinder placed parallel to a wall, which is moving at the speed of the far field has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the cylinder to wall separation length. The governing unsteady Navier–Stokes equations are discretized through the finite volume method on a staggered grid system. A SIMPLE type of algorithm has been used to compute the discretized equations iteratively. A shear layer of negative vortex generates along the surface of the wall, which influences the vortex shedding behind the cylinder. The flow‐field is distinct from the flow in presence of a stationary wall. An alternate vortex shedding occurs for all values of gap height in the unsteady regime of the flow. The strong positive vortex pushes the negative vortex upwards in the wake. The gap flow in the undersurface of the cylinder is strong and the velocity profile overshoots. The cylinder experiences a downward force for certain values of the Reynolds number and gap height. The drag and lift are higher at lower values of the Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigation of the vortex induced unsteadiness of a separation bubble via time-resolved and scanning PIV measurements

A transitional separation bubble on the suction side of an SD7003 airfoil is considered. The transition process that forces the separated shear layer to reattach seems to be governed by Kelvin–Helmholtz instabilities. Large scale vortices are formed due to this mechanism at the downstream end of the bubble. These vortices possess a three-dimensional structure and detach from the recirculation region, while other vortices are formed within the bubble. This separation of the vortex is a highly unsteady process, which leads to a bubble flapping. The structure of these vortices and the flapping of the separation bubble due to these vortices are temporally and spatially analyzed at angles of attack from 4° to 8° and chord-length based Reynolds numbers Re _c = 20,000–60,000 using time-resolved PIV measurements in a 2D and a 3D set-up, i.e., stereo-scanning PIV measurements are done in the latter case. These measurements complete former studies at a Reynolds number of Re _c = 20,000. The results of the time-resolved PIV measurements in a single light-sheet show the influence of the angle of attack and the Reynolds number. The characteristic parameters of the separation bubble are analyzed focusing on the unsteadiness of the separation bubble, e.g., the varying size of the main recirculation region, which characterizes the bubble flapping, and the corresponding Strouhal number are investigated. Furthermore, the impact of the freestream turbulence is investigated by juxtaposing the current and former results. The stereo-scanning PIV measurements at Reynolds numbers up to 60,000 elucidate the three-dimensional character of the vortical structures, which evolve at the downstream end of the separation bubble. It is shown that the same typical structures are formed, e.g., the c-shape vortex and the screwdriver vortex at each Reynolds number and angle of attack investigated and the occurrence of these patterns in relation to Λ-structures is discussed. To evidence the impact of the freestream turbulence, these results are compared with findings of former measurements.