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 vortex axis and vortex core border grid adaptation algorithm

In the design process of hydrodynamical and aerodynamical technical applications, the numerical simulation of massively separated vortical flow is crucial for predicting, for example, lift or drag. To obtain reliable numerical results, it is mandatory to accurately predict the physical behavior of vortices. Thus, the dominant vortical flow structures have to be resolved in detail, which requires a local grid refinement and certain adaptation techniques. In this paper, a vortex flow structure adaptation algorithm is presented, which is particularly designed for local grid refinement at vortex axes positions and associated vortex core border locations. To this end, a fast and efficient vortex axis detection scheme is introduced and the algorithm for the vortex core border determination is explained. As the interaction between vortices makes the assignment of grid points to a certain vortex axis difficult, a helicity‐based vortex distinction approach in combination with a geometrical rotational sensor is developed. After describing the combined different techniques in detail, the vortex feature adaptation algorithm is applied to analytical and more realistic examples, which show that the described grid adaptation algorithm is able to enhance the grid cell resolution locally such that all significant vortical flow phenomena are resolved. Copyright © 2008 John Wiley & Sons, Ltd.     

Population,characteristics and kinematics of vortices in a confined rectangular jet with a co-flow

Vortex behavior and characteristics in a confined rectangular jet with a co-flow were examined using vortex swirling strength as a defining characteristic. On the left side of the jet, the positively (counterclockwise) rotating vortices are dominant, while negatively rotating vortices are dominant on the right side of the jet. The characteristics of vortices, such as population density, average size and strength, and deviation velocity, were calculated and analyzed in both the cross-stream direction and the streamwise direction. In the near-field of the jet, the population density, average size and strength of the dominant direction vortices show high values on both sides of the center stream with a small number of counter-rotating vortices produced in the small wake regions close to jet outlet. As the flow develops, the wake regions disappear, these count-rotating vortices also disappear, and the population of the dominant direction vortices increase and spread in the jet. The mean size and strength of the vortices decrease monotonically with streamwise coordinate. The signs of vortex deviation velocity indicate the vortices transfer low momentum to high-velocity region and high momentum to the low velocity region. The developing trends of these characteristics were also identified by tracing vortices using time-resolved particle image velocimetry data. Both the mean tracked vortex strength and size decrease with increasing downstream distance overall. At the locations of the left peak of turbulent kinetic energy, the two-point spatial cross-correlation of swirling strength with velocity fluctuation and concentration fluctuation were calculated. All the correlation fields contain one positively correlated region and one negatively correlated region although the orientations of the correlation fields varied, due to the flow transitioning from wake, to jet, to channel flow. Finally, linear stochastic estimation was used to calculate conditional structures. The large-scale structures in the velocity field revealed by linear stochastic estimation are spindle-shaped with a titling stream-wise major axis.     

三角翼大迎角流场结构和气动力特性的数值分析研究

采用数值计算方法对亚音速三角翼纵向及带有小侧滑情况下的流场结构和气动力特性进行了计算。文中给出了三角翼大迎角纵向情况下气动力、机翼前缘分离涡轴线位置和旋涡破裂位置随迎角的变化规律，以及带有横侧小扰动和小侧滑情况下流场结构的非对称性对气动力的影响。计算结果表明与实验结果符合较好。     

Volumetric three-component velocimetry measurements of the turbulent flow around a Rushton turbine

Volumetric three-component velocimetry measurements have been taken of the flow field near a Rushton turbine in a stirred tank reactor. This particular flow field is highly unsteady and three-dimensional, and is characterized by a strong radial jet, large tank-scale ring vortices, and small-scale blade tip vortices. The experimental technique uses a single camera head with three apertures to obtain approximately 15,000 three-dimensional vectors in a cubic volume. These velocity data offer the most comprehensive view to date of this flow field, especially since they are acquired at three Reynolds numbers (15,000, 107,000, and 137,000). Mean velocity fields and turbulent kinetic energy quantities are calculated. The volumetric nature of the data enables tip vortex identification, vortex trajectory analysis, and calculation of vortex strength. Three identification methods for the vortices are compared based on: the calculation of circumferential vorticity; the calculation of local pressure minima via an eigenvalue approach; and the calculation of swirling strength again via an eigenvalue approach. The use of two-dimensional data and three-dimensional data is compared for vortex identification; a ‘swirl strength’ criterion is less sensitive to completeness of the velocity gradient tensor and overall provides clearer identification of the tip vortices. The principal components of the strain rate tensor are also calculated for one Reynolds number case as these measures of stretching and compression have recently been associated with tip vortex characterization. Vortex trajectories and strength compare favorably with those in the literature. No clear dependence of trajectory on Reynolds number is deduced. The visualization of tip vortices up to 140° past blade passage in the highest Reynolds number case is notable and has not previously been shown.     

Instantaneous and time-averaged flow fields of multiple vortices in the tip region of a ducted propulsor

The instantaneous and time-averaged flow fields in the tip region of a ducted marine propulsor are examined. In this flow, a primary tip-leakage vortex interacts with a secondary, co-rotating trailing edge vortex and other co- and counter-rotating vorticity found in the blade wake. Planar particle imaging velocimetry (PIV) is used to examine the flow in a plane approximately perpendicular to the mean axis of the primary vortex. An identification procedure is used to characterize multiple regions of compact vorticity in the flow fields as series of Gaussian vortices. Significant differences are found between the vortex properties from the time-averaged flow fields and the average vortex properties identified in the instantaneous flow fields. Variability in the vortical flow field results from spatial wandering of the vortices, correlated fluctuations of the vortex strength and core size, and both correlated and uncorrelated fluctuations in the relative positions of the vortices. This variability leads to pseudo-turbulent velocity fluctuations. Corrections for some of this variability are performed on the instantaneous flow fields. The resulting processed flow fields reveal a significant increase in flow variability in a region relatively far downstream of the blade trailing edge, a phenomenon that is masked through the process of simple averaging. This increased flow variability is also accompanied by the inception of discrete vortex cavitation bubbles, which is an unexpected result, since the mean flow pressures in the region of inception are much higher than the vapor pressure of the liquid. This suggests that unresolved fine-scale vortex interactions and stretching may be occurring in the region of increased flow variability.     

Unsteady characteristics of inlet vortices

An experimental study of the unsteady characteristics of inlet vortices has been conducted using a high-frame rate digital particle image velocimetry system. The results revealed the formation of a pair of counter-rotating inlet vortices for the no-wind configuration and one single inlet vortex when there was crosswind. In all measurement planes, from near the ground to the inlet, evidence of vortex meandering with quasi-periodicity was found. The vortex meander is dominant in the direction of the crosswind, and its amplitude increases with crosswind velocity. The proper orthogonal decomposition analysis of the instantaneous velocity field suggested that the most energetic mode was a helical displacement wave, corresponding to the first helical mode. Similarities with the meandering of the trailing vortices from wings were noted. The present results also suggest that the unsteady characteristics of the focus of separation formed on the ground might be responsible for the unsteady nature of the inlet vortex.     

Wind Impact on Single Vortices and Counterrotating Vortex Pairs in Ground Proximity

Wall-resolved large eddy simulations are employed to investigate the behaviour of wake vortices and single vortices in ground proximity at a variety of wind conditions. The six considered strengths of wind, ranging between 0.5 and 4 times the initial wake vortex descent speed, w₀, include practically and theoretically significant wind speeds. A crosswind of 0.5 w₀ may lead to windward stall posing a potential hazard to subsequently landing aircraft, whereas theoretical considerations predict that at 4 w₀ the rebound of the luff vortex is completely suppressed. The same range of wind speeds is also used to investigate the effects of headwind and diagonal wind in order to discriminate between effects of environmental turbulence increasing with wind speed and the direction of the wind shear. The study has been complemented by a number of single vortex computations in order to differentiate between effects related to the mutual interaction of the vortex pair and the individual vortices with the turbulent boundary layer flow. It is shown that vortex ascent, descent, rebound and decay characteristics are controlled by (i) the interaction of the vortices with secondary vorticity detaching from the ground, (ii) the redistribution of vorticity of the boundary layer which is altering the path of the primary vortices by mutual velocity induction, and (iii) the interaction of the vortices with the environmental turbulence.     

Computation of field structure and aerodynamic characteristics of delta wings at high angles of attack

A numerical investigation of the structure of the vortical flowfield over delta wings at high angles of attack in longitudinal and with small sideslip angle is presented. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics that are dominated by the effect of the breakdown of the leading-edge vortices. The methods that analyze the flowfield structure quantitatively were given by using flowfield data from the computational results. In the region before the vortex breakdown, the vortex axes are approximated as being straight line. As the angle of attack increases, the vortex axes are closer to the root chord, and farther away from the wing surface. Along the vortex axes, as the adverse pressure gradients occur, the axial velocity decreases, that is, A is negativee, so the vortex is unstable, and it is possible to breakdown. The occurrence of the breakdown results in the instability of lateral motion for a delta wing, and the lateral moment diverges after a small perturbation occurs at high angles of attack. However, after a critical angle of attack is reached the vortices breakdown completely at the wing apex, and the instability resulting from the vortex breakdown disappears.     

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

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

Bubble clustering and trapping in large vortices. Part 2: Time-dependent trapping conditions

In turbulent, periodically excited jets, interactions between bubbles and large coherent vortices are quantitatively studied. Simultaneous, two-phase PIV (particle image velocimetry) and photographic recordings were applied for tracking the large vortices and bubble structures and for investigating trapping phenomena. In order to quantify the interaction between bubbles and the large vortices that are formed in the shear layer, characteristic phase-averaged quantities were determined by PIV. The time-dependent vortex radius, the vorticity at the vortex centre and the time-dependent trapping conditions, obtained from the simulation of the vortex development, were tested against the experimental data.     

Granular vortices: Identification,characterization and conditions for the localization of deformation

We relate the micromechanics of vortex evolution to that of force chain buckling and, on this basis, formulate the conditions for strain localization in a continuum model of dense granular media. Using the traditional bifurcation analysis of shear bands, we show that kinematic vortex fields are in fact solutions to the boundary value problem satisfying null boundary conditions. To establish an empirical basis for our study, we first develop a method to identify the location of the core and boundary of each vortex from a given displacement field in two dimensions. We then employ this method to characterize the residual deformation field (i.e., the deviation of particle motions from the continuum deformation) in a physical experiment and a discrete element simulation of dense granular samples submitted to biaxial compression. Vortices in the failure regime are essentially confined to the shear band. Primary vortices, the clear majority, rotate in the same direction as the shear band; secondary vortices, the so-called wakes, rotate in the opposite direction. Primary vortices align in spatial succession along the central axis of the band; wakes form next to the band boundaries, in between and beside two adjacent primary vortices. Force chain buckling, the governing mechanism for shear bands, is responsible for vortex formation in the failure regime. Vortex dynamics are consistent with stick-slip dynamics. From quiescent conditions of jamming or stick, vortical motions arise from force chain buckling and associated relative particle rotations and sliding; these in turn precipitate intermittent periods of unjamming or slip, evident in the attendant drops in stress ratio and bursts in both kinetic energy and local nonaffine deformation. A kinematic vortex field inside shear bands is proposed that is consistent with the equations of continuum mechanics and the underlying instability of force chain buckling: such a field is periodic with a repeating unit cell comprising a primary vortex at the center of the band, with two trailing wakes close next to the band boundaries.     

The flow in a planar-radial vortex chamber. 2. Vortex structure of the flow

A vortex structure of an air flow with a characteristic size of vortices comparable with the primary vortex size was observed in a vortex chamber of planar-radial geometry for the first time. The vortex component of the flow velocity along the chamber radius and its axis was calculated. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 1, pp. 41–49, January–February, 2000.     

轴对称射流场涡结构的离散涡段方法研究

本文用三维离散涡方法,模拟了轴对称圆射流涡结构的发展。     

Numerical studies on evolution of secondary streamwise vortices in transitional boundary layers

Formation and evolution of secondary streamwise vortices in the compressible transitional boundary layers over a flat plate are studied using a direct numerical simulation method with high-order accuracy and highly effective non-reflecting characteristic boundary conditions. Generation and development processes of the secondary streamwise vortices in the complicated transitional boundary flow are clearly analyzed based on the of numerical results, and the effects on the formation of the ring-like vortex that is vital to the boundary layer transition are explored. A new mechanism forming the ring-like vortex through the mutual effect of the primary and secondary streamwise vortices is expressed.     

Numerical simulation of jet-forced flow in a circular reservoir using discrete and random vortex methods

This paper describes a Biot–Savart discrete vortex model for simulating the flow patterns which occur when a single high-velocity inflow jet is used to stir the fluid within a circular container. The first stage of the model consists of conformally mapping the circular perimeter of the container onto a rectangle by means of a Schwarz–Christoffel transformation. A potential flow solution is then obtained for the flow inside the rectangle and this is transformed to give the potential flow inside the circle. In the second stage of the simulation, discrete vortices are added at the inlet of the physical system in order to model the inflow shear layers. Velocity components resulting from the discrete vortices and their images in the walls of the cylinder are superimposed on the uniform potential flow solution. The positions of the vortices are updated using a Lagrangian tracking procedure. Viscous effects are incorporated through the use of random walks. From the results it is shown that the discrete vortex method does predict qualitatively the important features of jet-forced reservoir flow.     

Analytical model of motion of turbulent vortex rings in an incompressible fluid

An analytical model describing the motion of vortex rings in an incompressible fluid is constructed. The model is valid both for homogeneous and inhomogeneous vortices buoyant in the gravity field, as well as for combined vortices. The expansion angle of a buoyant vortex is found from the characteristic parameters that define the flow rather than specified on the basis of experiments. Significant differences in the expansion angles of homogeneous and buoyant vortex rings are explained. The calculation results for the proposed model are compared with the results of laboratory experiments and data on the rise of the cloud produced by an atomic explosion.     

Diffusion via splitting and remeshing via merging in vortex methods

The technique of splitting a fat vortex element (with a core width larger than some threshold) into some thin ones in order to fix the convergence problem of the core‐spreading vortex methods is convenient and efficient. In particular, it keeps the method purely Lagrangian. In the present investigation, the splitting process is further viewed as part of the physical diffusion process. A new splitting method in which several weaker child vortices surround a thinned but still strong parent vortex is proposed. It is found that because of the survival of the parent vortex, the error arising from the splitting events can be largely reduced. The computational amount on the other hand is kept reasonably large by merging similar and close‐by vortices. The merging scheme designed herein not only involves fewer restrictions but also allows merging vortices of opposite rotations through the viewpoint of remeshing. The validity and accuracy of these techniques, proposed particularly for simulations undergoing lots of splitting and merging events, are verified by successfully simulating the interactions between two Burgers vortices under an external straining field. Copyright © 2005 John Wiley & Sons, Ltd.