应用PIV技术测试涡旋波流场

涡旋波流动作为一种特殊的流动现象,可以使流体在相对较宽的槽道中产生较强的波动和对流混合,从而在小Re数条件下起到强化传质的效果。本文利用PIV流场显示技术,对振荡流在非对称槽道中所形成的涡旋波的产生机理和发展规律进行了实验研究和定量分析,测得了涡旋波流场的速度矢量图,阐明了涡旋波流场周期性变化的特点。分析了Re数和St数对涡旋波流动的影响,并得出了旋涡涡心位置以及涡心处涡量的动态变化规律。     

非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

用PIV测试涡旋波流场的速度和剪应力分布

采用PIV瞬态流场测试技术,对二维槽道中的涡旋波流场不同相位上的速度分布和应力分布进行了测试和计算,本文定量地描述了槽道中涡旋波的形成过程及发展规律。通过调整振荡流的振幅和活塞行程,分析了Re数和Sr数对涡旋波流动的影响,得出了槽道内剪应力的分布状况以及平均剪切应力的周期性变化规律。深入分析了涡旋波流场强化传递现象的本质。     

槽道内涡波流场的POD分析

对槽道内涡波流场的瞬态速度矢量场进行了2DPIV测量实验,将2DPIV测量的矢量场数据进行POD分析,根据POD分解的各阶模态的能量比确定了表征涡波流场主导结构的前15阶模态。结果表明,POD分解的前15阶模态发现槽道内涡波流场是由槽道壁面剪切层诱导的涡列以及伴随的波状主流组成;流场中大尺度的涡旋发展为涡对,对波状主流的脉动频率产生影响;根据涡波流场中的驻点和鞍点,获取了流场的大尺度涡对、平均流场以及Helmholtz涡环等明显特征;最后根据POD分解的前15阶模态对槽道内涡波流场进行重组,重组流场表征了槽道内层流状态下波状主流的形态和涡旋共存的涡波结构以及驻点和鞍点的位置处涡旋的变化等主要特征,有效地剔除了PIV测量流场中的随机信息,保留了PIV测量流场的主导特征。     

槽道内涡波流场的POD分析

对槽道内涡波流场的瞬态速度矢量场进行了2DPIV测量实验，将2DPIV测量的矢量场数据进行POD分析，根据POD分解的各阶模态的能量比确定了表征涡波流场主导结构的前15阶模态。结果表明，POD分解的前15阶模态发现槽道内涡波流场是由槽道壁面剪切层诱导的涡列以及伴随的波状主流组成；流场中大尺度的涡旋发展为涡对，对波状主流的脉动频率产生影响；根据涡波流场中的驻点和鞍点，获取了流场的大尺度涡对、平均流场以及Helmholtz涡环等明显特征；最后根据POD分解的前15阶模态对槽道内涡波流场进行重组，重组流场表征了槽道内层流状态下波状主流的形态和涡旋共存的涡波结构以及驻点和鞍点的位置处涡旋的变化等主要特征，有效地剔除了PIV测量流场中的随机信息，保留了PIV测量流场的主导特征。     

基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

均匀旋转对圆柱水动力及流动结构的影响

基于格子Boltzmann方法 (LBM)对均匀旋转控制下的低雷诺数(Re=100)圆柱绕流问题进行了数值模拟,得到了转速比从0～10变化下,旋转控制对圆柱水动力及流动结构的影响规律.使用动态模态分解(DMD)对流场特征进行提取,并分析了施加旋转控制之后转速比对流场不同模态和增长率的影响.结果表明,随着转速比增大,圆柱下游流动结构依次呈现出卡门涡街、剪切层、反向剪切层、单侧涡和附着涡5种结构;阻力系数时均值先减小,随后在转速进入单侧涡区间后增大,升力系数与力矩系数的时均值均单调增加,同时,在出现涡脱落的两个转速区间内,水动力出现了明显的波动,且二次失稳时波动幅度更大. DMD的结果表明,圆柱下游的流动结构主要受圆柱壁面的旋转影响而发生改变并产生全新流动模态;旋转会对流动稳定性产生影响:在未充分发展阶段,旋转对流动稳定性的影响不显著,而在充分发展后,各转速下的流场不稳定模态数均远少于未充分发展阶段,随着转速比的增大,流动稳定性会产生不同程度的增强或减弱,且无涡脱落时的稳定性高于有涡脱落时,因此,通过旋转控制抑制尾涡脱落可以有效增强流动的稳定性.     

基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

分离涡方法对梢涡中脉动量的数值模拟

分别用RANS-SA方法和DES方法对NACA0012翼端梢涡进行模拟计算,分析了梢涡区域网格局部加密对梢涡计算结果的影响,并与实验结果进行了对比.相比于RANS-SA方法,DES方法在梢涡流场计算中具有更好的适用性,能够得到更准确的流动信息和更精细的涡结构;另外,网格局部加密对脉动量的计算影响很大.通过对脉动量的分析发现,在近尾缘处,几股涡的融合产生了比较强烈的脉动,随着梢涡的逐渐稳定,脉动量也逐渐减小;现有的实验结果显示在偏下游处会产生梢涡的振荡现象,使统计脉动量增大,而本文计算中未发现该现象.     

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

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

Discretization of the vorticity field of a planar jet

In data assimilation, information from sensors is used to correct the state variables of a numerical model. This has been used to great advantage by the weather prediction community in the context of direct numerical simulation (DNS) models, but has seen comparatively little use in point-vortex models. This is due in large part to data-processing issues. In order to keep up with the speeds necessary for effective data assimilation, one must extract and discretize the vortex structures from velocity field data in a computationally efficient fashion—i.e., using as few discrete vortices as possible to model the measured flow. This paper describes a new strategy for accomplishing this and evaluates the results using data from a laboratory-scale vortex-dominated planar jet. Large-scale vortex structures are found using a family of variants on traditional vortex extraction methods. By augmenting these methods with simple computational topology techniques, one obtains a new method that finds the boundaries of the coherent structures in a manner that naturally follows the geometry of the flow. This strategy was evaluated in the context of two standard vortex extraction methods, vorticity thresholding and Okubo–Weiss, and tested upon velocity field data from the experimental fluid flow. The large-scale structures found in this manner were then modeled with collections of discrete vortices, and the effects of the grain size of the discretization and the parameters of the discrete vortex model were studied. The results were evaluated by comparing the instantaneous velocity field induced by the discrete vortices to that measured in the jet. These comparisons showed that the two extraction techniques were comparable in terms of sensitivity and error, suggesting that the computationally simpler vorticity thresholding method is more appropriate for applications where speed is an issue, like data assimilation. Comparisons of different discretization strategies showed that modeling each large-scale vortex structure with a single discrete vortex provided the best compromise between mean-squared error and computational effort. These results are of potential interest in any situation where one must balance accuracy and expense while extracting vortices from a snapshot of a flow field; data assimilation is only one example.     

A numerical exploration of secondary separation in starting flow around a flat plate by the discrete vortex model

In the present work, a further numerical simulation of the starting flow around a flat plate normal to the direction of motion in a uniform fluid has been made by means of the discrete vortex method. The secondary separation occurring at rear surface of the plate is explored, and predicted approximately using Thwait's method. The calculated results show that in the early stages of the flow secondary separation does occur. The evolution of flow field, the vortex growing process and the characteristics of secondary vortices have been described. The time dependent drag coefficients, the vorticity shed from the edges and rear surface, and the separation positions are calculated as well as the distributions of velocity and pressure on the plate. In the case of flow normal to the plate, the calculated secondary vortices are weak. Their existence will change the local velocity distributions and affect pressure distributions. However, the effect on drag coefficient is negligible.     

The effect of fine structure on the stability of planar vortices

This study considers the linear, inviscid response to an external strain field of classes of planar vortices. The case of a Gaussian vortex has been considered elsewhere, and an enstrophy rebound phenomenon was noted: after the vortex is disturbed enstrophy feeds from the non-axisymmetric to mean flow. At the same time an irreversible spiral wind-up of vorticity fluctuations takes place. A top-hat or Rankine vortex, on the other hand, can support a non-decaying normal mode.In vortex dynamics processes such as stripping and collisions generate vortices with sharp edges and often with bands or rings of fine scale vorticity at their periphery, rather than smooth profiles. This paper considers the stability and response of a family of vortices that vary from a broad profile to a top-hat vortex. As the edge of the vortex becomes sharper, a quasi-mode emerges and vorticity winds up in a critical layer, at the radius where the angular velocity of the fluid matches that of a normal mode on a top-hat vortex. The decay rate of these quasi-modes is proportional to the vorticity gradient at the critical layer, in agreement with theory. As the vortex edge becomes sharper it is found that the rebound of enstrophy becomes stronger but slower.The stability and linear behaviour of coherent vortices is then studied for distributions which exhibit additional fine structure within the critical layer. In particular we consider vorticity profiles with ‘bumps’, ‘troughs’ or ‘steps’ as this fine structure. The modified evolution equation that governs the critical layer is studied using numerical simulations and asymptotic analysis. It is shown that depending on the form of the short-scale vorticity distribution, this can stabilise or destabilise quasi-modes, and it may also lead to oscillatory behaviour.     

Vortex shedding flow behind a slowly rotating circular cylinder

This paper investigates flow past a rotating circular cylinder at 3600?Re?5000 and α?2.5. The flow parameter α is the circumferential speed at the cylinder surface normalized by the free-stream velocity of the uniform cross-flow. With particle image velocimetry (PIV), vortex shedding from the cylinder is clearly observed at α<1.9. The vortex pattern is very similar to the vortex street behind a stationary circular cylinder; but with increasing cylinder rotation speed, the wake is observed to become increasing narrower and deflected sideways. Properties of large-scale vortices developed from the shear layers and shed into the wake are investigated with the vorticity field derived from the PIV data. The vortex formation length is found to decrease with increasing α. This leads to a slow increase in vortex shedding frequency with α. At α=0.65, vortex shedding is found to synchronize with cylinder rotation, with one vortex being shed every rotation cycle of the cylinder. Vortex dynamics are studied at this value of α with the phase-locked eduction technique. It is found that although the shear layers at two different sides of the cylinder possess unequal vorticity levels, alternating vortices subsequently shed from the cylinder to join the two trains of vortices in the vortex street pattern exhibit very little difference in vortex strength.     

Aerodynamic forces and three-dimensional flow structures in the mean wake of a surface-mounted finite-height square prism

Different flow models have been proposed for the flow around surface-mounted finite-height square prisms, but there is still a lack of consensus about the origin and connection of the streamwise tip vortices with the other elements of the wake. This numerical study was performed to address this gap, in addition to clarifying the relationship of the near-wake structures with the far wake and the near-wall flow, which is associated with the fluid forces. A large-eddy simulation approach was adopted to solve the flow around a surface-mounted finite-height square prism with an aspect ratio of AR = 3 and a Reynolds number Re = 500. The mean drag and normal forces and the bending moment for the prism were quantitatively compared in terms of skin-friction and pressure contributions, and related to the near-wall flow. Both three-dimensional visualizations and planar projections of the time-averaged flow field were used to identify, qualitatively, the main structures of the wake, including the horseshoe vortex, corner vortices and regions of high streamwise vorticity in the upper part of the wake. These features showed the same qualitative behavior as reported in high Reynolds number studies. It was found that some regions of high streamwise vorticity magnitude, like the tip vortices, are associated with the three-dimensional bending of the flow, and the tip vortices did not continuously extend to the free end of the prism. The three-dimensional flow analysis, which integrated different observations of the flow field around surface-mounted finite-height square prisms, also revealed that the mean near-wake structure is composed of two sections of different origin and location of dominance.     

Two- and three-dimensional numerical simulations of the clap–fling–sweep of hovering insects

The importance of three-dimensional effects for flapping wings is addressed by means of numerical simulation. In particular, the clap–fling–sweep mechanism is examined. The flow at the beginning of the downstroke is shown to be in reasonable agreement with the two-dimensional approximation. After the wings move farther than one chord length apart, three-dimensional effects become essential. Two values of the Reynolds number are considered. At Re=128, the spanwise flow from the wing roots to the wing tips is driven by the centrifugal forces acting on the mass of the fluid trapped in the recirculation bubble behind the wings. It removes the excess of vorticity and delays the periodic vortex shedding. At Re=1400, vortex breakdown occurs past the outer portion of the wings, and multiple vortex filaments are shed into the wake.