Model of two-dimensional vortex motion with an entrainment mechanism

A model for describing the vertically averaged vortex motions of an incompressible viscous fluid with an arbitrary vertical structure of the bottom Ekman boundary layer is proposed. An approach similar to that adopted in [1] is used: the second moments of the deviations from the average velocities required in order to close the vorticity equation are calculated by means of the Ekman solution for gradient flows, which makes it possible to take the integral bottom boundary layer effect into account. As a result, these terms lead to a specific form of nonlinear friction with a coefficient that depends on the vorticity of the average flow. In the particular case of a constant vertical turbulent transfer coefficient the inaccuracies of the model described in [1] can be eliminated. The generalized vorticity equation obtained has solutions of the vorticity spot type with a uniform internal vorticity distribution, which can be effectively investigated by means of appropriate algorithms [2]. The mechanism of entrainment at the vorticity front is illustrated with reference to the example of the evolution of vorticity spots. An exact solution of the problem of the evolution of an elliptic vortex (generalized Kirchhoff vortex), which in the case of fairly strong anticyclonic vorticity degenerates first into a line segment (vortex sheet) and then into a point vortex, is constructed. Equations describing the dynamics of an elliptic vorticity spot in an external field with a linear dependence of the velocity on the horizontal coordinates and generalizing the classical Chaplygin-Kida model [3, 4] are constructed.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 49–56, November–December, 1992.     

用拉格朗日相关结构研究圆盘启动过程的流体输运

利用粒子成像测速(PIV)技术,得到了圆盘启动涡环流场的速度分布和涡量分布.圆盘启动涡环流场的有限时间李雅普诺夫指数场(Finit-time Lyapunov exponents,FTLE)以及拉格朗日相关结构(Lagrangian coherent structures,LCS)被计算出来.基于圆盘启动涡环流场的有限时间李雅普诺夫指数场以及拉格朗日相关结构,通过跟踪流体质点,对圆盘启动涡环流场的输运情况进行了分析.在圆盘启动涡环形成过程中,流体发现被圆盘和相互排斥的拉格朗日相关结构分成三部分.剪切流窗口(vorticity-flux window)被发现,涡量流通过剪切流窗口进入涡核.涡环的非定常边界被确定,它由相互排斥的拉格朗日相关结构背风面、圆盘以及剪切流窗口组成.     

Investigation of vortex development on accelerating spanwise-flexible wings

The effect of spanwise flexibility on the development of leading-edge vortices for impulsively started flat plates at low Reynolds numbers has been investigated. A theoretical model is proposed, based on the Euler–Bernoulli beam theory, coupled with a vortex growth model based on vorticity flux through a leading-edge shear layer. The model was validated for rigid and flexible flat plates undergoing a towing motion at an angle-of-attack of 45°. It is shown that a phase-delay in lift and drag generation occurs between rigid and flexible cases. The model indicates decreasing vorticity along the span as the wing is accelerated and begins to bend. Particle image velocimetry measurements conducted at three different spanwise planes showed a delay in vortex growth along the span, despite a uniform spanwise circulation. This uniform spanwise distribution of circulation is in contrast to the quasi-two-dimensional model, which predicted a reduced circulation near the profile tip where plate motion was delayed. It is therefore concluded that circulation must be dynamically redistributed through vorticity convection during the impulsive motion.     

Numerical study of an inviscid incompressible flow through a channel of finite length

A two‐dimensional inviscid incompressible flow in a rectilinear channel of finite length is studied numerically. Both the normal velocity and the vorticity are given at the inlet, and only the normal velocity is specified at the outlet. The flow is described in terms of the stream function and vorticity. To solve the unsteady problem numerically, we propose a version of the vortex particle method. The vorticity field is approximated using its values at a set of fluid particles. A pseudo‐symplectic integrator is employed to solve the system of ordinary differential equations governing the motion of fluid particles. The stream function is computed using the Galerkin method. Unsteady flows developing from an initial perturbation in the form of an elliptical patch of vorticity are calculated for various values of the volume flux of fluid through the channel. It is shown that if the flux of fluid is large, the initial vortex patch is washed out of the channel, and when the flux is reduced, the initial perturbation evolves to a steady flow with stagnation regions. Copyright © 2008 John Wiley & Sons, Ltd.     

Lateral Spreading Mechanism of a Turbulent Spot and a Turbulent Wedge

We examine the spreading of turbulent spots and wedges into a surrounding laminar Blasius boundary layer. The spreading is not due to the lateral propagation of turbulent eddies but rather to a developing disturbance in the surrounding spanwise vorticity of the laminar boundary layer. We concentrate on the mechanisms for generating streamwise vorticity. In particular, inclined generally streamwise vortex tubes along the spot/wedge boundary tilt mean shear vortex lines either up or down. These lines subsequently tend to either lag back or lead forward. As the leading or lagging vortex lines continue to wrap around and reinforce the causative inclined tube, the lines arch up or down. The outboard portion of the resulting arch must acquire a vertical, ω_y, component of vorticity which induces the rollup of a new inclined tube now outboard of the first. Close to the wall the arching mechanism is inhibited by the no through flow boundary condition while far from the wall the process is inhibited by the lack of sufficient mean spanwise vorticity.     

Parallelization of a vorticity formulation for the analysis of incompressible viscous fluid flows

A parallel computer implementation of a vorticity formulation for the analysis of incompressible viscous fluid flow problems is presented. The vorticity formulation involves a three‐step process, two kinematic steps followed by a kinetic step. The first kinematic step determines vortex sheet strengths along the boundary of the domain from a Galerkin implementation of the generalized Helmholtz decomposition. The vortex sheet strengths are related to the vorticity flux boundary conditions. The second kinematic step determines the interior velocity field from the regular form of the generalized Helmholtz decomposition. The third kinetic step solves the vorticity equation using a Galerkin finite element method with boundary conditions determined in the first step and velocities determined in the second step. The accuracy of the numerical algorithm is demonstrated through the driven‐cavity problem and the 2‐D cylinder in a free‐stream problem, which represent both internal and external flows. Each of the three steps requires a unique parallelization effort, which are evaluated in terms of parallel efficiency. Copyright © 2002 John Wiley & Sons, Ltd.     

A viscous vortex particle method for deforming bodies with application to biolocomotion

Bio‐inspired mechanics of locomotion generally consist of the interaction of flexible structures with the surrounding fluid to generate propulsive forces. In this work, we extend, for the first time, the viscous vortex particle method (VVPM) to continuously deforming two‐dimensional bodies. The VVPM is a high‐fidelity Navier–Stokes computational method that captures the fluid motion through evolution of vorticity‐bearing computational particles. The kinematics of the deforming body surface are accounted for via a surface integral in the Biot–Savart velocity. The spurious slip velocity in each time step is removed by computing an equivalent vortex sheet and allowing it to flux to adjacent particles; hence, no‐slip boundary conditions are enforced. Particles of both uniform and variable size are utilized, and their relative merits are considered. The placement of this method in the larger class of immersed boundary methods is explored. Validation of the method is carried out on the problem of a periodically deforming circular cylinder immersed in a stagnant fluid, for which an analytical solution exists when the deformations are small. We show that the computed vorticity and velocity of this motion are both in excellent agreement with the analytical solution. Finally, we explore the fluid dynamics of a simple fish‐like shape undergoing undulatory motion when immersed in a uniform free stream, to demonstrate the application of the method to investigations of biomorphic locomotion. Copyright © 2008 John Wiley & Sons, Ltd.     

Interaction of Two-Dimensional Trailing Vortex Pair with a Shear Layer

The basic laws governing the interaction of a two-dimensional vortex pair with a shear layer of constant thickness are considered. The main idea of the study is to develop and adapt a simplified representation of a hydrodynamic flow based on a point-vortex model simulating the actual interaction of full-scale vortex patterns over the ground surface. It is shown that vortices with vorticity opposite in sign to the shear layer may stop or even ricochet from this layer, while the other vortex may penetrate through the layer. Numerical results are presented as plots and analyzed     

Control of flow distortion in offset diffusers using trapped vorticity

Controlled concentrations of trapped vorticity within an offset, subsonic (M_AIP ≤ 0.7) diffuser are explored for active suppression of flow distortion in joint experimental and numerical investigations. The coupling between trapped vorticity, used to model boundary-layer separation, and secondary-flow vortices is manipulated using an array of fluidic oscillating jets, which are spanwise distributed just upstream of the trapped vortex. Actuation energizes the separated shear layer, reducing the size of separation and effecting an earlier reattachment of the boundary layer, which favorably effects the flow field downstream of reattachment. It is shown that optimal interactions between actuation and the trapped vortex fully suppress the central vortex pair, and redistributes the residual vorticity around the diffuser's circumference. This results in a 68% reduction in circumferential distortion at the Aerodynamic Interface Plane (AIP), using an actuation mass flow rate that is only 0.25% of the diffuser mass flow rate.     

A visual study of vortex-induced subcritical instability on a flat plate laminar boundary layer

This paper reports results of our experimental investigation on flow instability on a flat plate laminar boundary layer caused by a captive vortex migrating far outside the boundary layer. Results show that the sign of the circulation associated with the vortex is the main determinant for the severity of the boundary layer instability. A captive vortex with an opposite sign to that of the unperturbed shear layer vorticity causes a breakdown ahead of it, while the one with the same sign as the unperturbed shear layer vorticity gives rise to weaker excitation trailing it. Additional parameters that influence the flow instability are the strength and distance of the vortical disturbance from the boundary layer, as well as the translational speed of the vortex. These experimental results compliment the corresponding theoretical analysis of Sengupta et al. (J Fluid Mech 493:277–286, 2003).     

Strong vortex/boundary layer interactions

This is the second of two papers on the interaction between a longitudinal vortex pair, produced by a delta-wing at angle of attack, and a turbulent boundary layer developing on a flat plate. In the first paper only the outer parts of the vortices entered the boundary layer whereas in this paper the vortices merge with it. In the resultant interaction, the boundary layer between the vortices is kept thin by lateral divergence and a three-dimensional separation line is formed outboard of each vortex. Turbulent, momentum-deficient fluid containing longitudinal vorticity is entrained from the boundary layer along these lines and wrapped around the vortices. As a consequence, the turbulent region of the vortices increases in size and the circulation slowly decreases. It is shown that the flow near the separation line and in the vortices is complicated, and this interaction is expected to be more difficult to calculate than the first. Detailed mean flow and turbulence measurements are reported.     

Method of mean-mass parameters for a three-dimensional boundary layer in a flow with vorticity

When a gas flows with hypersonic velocity over a slender blunt body, the bow shock induces large entropy gradients and vorticity near the wall in the disturbed flow region (in the high-entropy layer) [1]. The boundary layer on the body develops in an essentially inhomogeneous inviscid flow, so that it is necessary to take into account the difference between the values of the gas parameters on the outer edge of the boundary layer and their values on the wall in the inviscid flow. This vortex interaction is usually accompanied by a growth in the frictional stress and heat flux at the wall [2, 3]. In three-dimensional flows in which the spreading of the gas on the windward sections of the body causes the high-entropy layer to become narrower, the vortex interaction can be expected to be particularly important. The first investigations in this direction [4–6] studied the attachment lines of a three-dimensional boundary layer. The method proposed in the present paper for calculating the heat transfer generalizes the approach realized in [5] for the attachment lines and makes it possible to take into account this effect on the complete surface of a blunt body for three-dimensional laminar, transition, or turbulent flow regime in the boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 80–87, January–February, 1981.     

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.     

An inclined rectangular jet in a turbulent boundary layer-vortex flow

A model test study was performed on streamwise vortices generated by a rectangular jet in an otherwise flat-plate turbulent boundary layer. The study was conducted in a low speed wind tunnel. The rectangular jet had a cross-section size of 28 mm by 5.5 mm. The oncoming boundary layer had a 99.5 percent thickness of 25 mm. The freestream speed of the oncoming flow was 20 m/s. Measurements were performed with a three-element LDA system. The effects of skew angle and streamwise development of vortex were investigated and the mean flow properties are presented. The study showed that the rectangular jet was able to produce a streamwise vortex of higher strength than that of a round jet, while at the same time keeping the same size and shape as that of a round jet. A 63% increase in the maximum vorticity was found. The 45^∘ skew angle was identified as the optimal skew angle for vortex production. Received: 24 June 1998/ Accepted: 21 May 1999     

Evaluation of unsteady pressure fields and forces in rotating airfoils from time-resolved PIV

The final stages of transitional phenomena in laminar separation bubbles play a key role in their reattachment process, and they condition the boundary layer properties and flow structure after reattachment. In this experimental study, the evolution of the perturbation velocity spectra found in this zone is first presented, showing the nonlinear growth of instabilities in their path to develop fully turbulent spectra. The study of the average flow field allows the scaling of the reattachment region, both in its extension and in the characterization of the integral boundary layer magnitudes. Experimental laws are proposed for the evolution of the momentum thickness and of the shape factor. In addition, a universal, wake-like mean velocity profile is found shortly after the reattachment station. The phase-locked characterization technique allows measurements conditioned to the presence of a fluid event. This technique is used to track the evolution of large-scale structures, whose dynamics is seen to dominate the fluid behavior in the reattachment zone. The simultaneous existence of two vortex blobs is found to characterize this flow region, with the longest lived one being convected toward the wall and stretched. This process results in the fast breakdown of the large-scale vorticity structure and the sudden formation of 3-D, small scales that promote the rapid flow evolution toward a fully developed turbulent state.     

Nonlinear standing waves on the surface of a viscous fluid

Standing surface waves in a viscous infinite-depth fluid are studied. The solution of the problem is obtained in the linear and quadratic approximations. The case of long, as compared with the boundary layer thickness, waves is analyzed in detail. The trajectories of fluid particles are determined and an expression for the vorticity is derived.     

Long-wave instability of a vortex ring

The instability of barrel-shaped vibrations of a vortex ring in an ideal fluid is investigated. These vibrations, stable for a vortex ring with a piecewise-uniform vorticity profile, appear to be unstable for a vortex ring with a smooth vorticity profile. The instability growth rate is found on the basis of the energy balance equation determining the energy transport from perturbations with negative energy in the critical layer to perturbations with positive energy in the rest of the flow. The curvature of the vortex ring, by virtue of which the perturbations with energies of different signs appear to be connected, plays a prominent role in the mechanism under consideration.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 72–78, November–December, 1995.     

Effect of velocity ratio on the streamwise vortex structures in the wake of a stack

The time-averaged velocity and streamwise vorticity fields within the wake of a stack were investigated in a low-speed wind tunnel using a seven-hole pressure probe. The experiments were conducted at a Reynolds number, based on the stack external diameter, of Re_D=2.3×10⁴. The stack, of aspect ratio AR=9, was mounted normal to a ground plane and was partially immersed in a flat-plate turbulent boundary layer, where the ratio of the boundary layer thickness to the stack height was δ/H≈0.5. The jet-to-cross-flow velocity ratio was varied from R=0 to 3, which covered the downwash, crosswind-dominated and jet-dominated flow regimes. In the downwash and crosswind-dominated flow regimes, two pairs of counter-rotating streamwise vortex structures were identified within the stack wake. The tip vortex pair located close to the free end of the stack, and the base vortex pair located close to the ground plane within the flat-plate boundary layer, were similar to those found in the wake of a finite circular cylinder, and were associated with the upwash and downwash flow fields within the stack wake, respectively. In the jet-dominated flow regime, a third pair of streamwise vortex structures was observed, referred to as the jet-wake vortex pair, which occurred within the jet-wake region above the free end of the stack. The jet-wake vortex pair had the same orientation as the base vortex pair and was associated with the jet rise. The peak vorticity and strength of the streamwise vortex structures were functions of the jet-to-cross-flow velocity ratio. For the tip vortex structures, their peak vorticity and strength reduced as the jet-to-cross-flow velocity ratio increased.     

Head on collisions of compressible vortex rings on a smooth solid surface

An experimental study has been conducted on the effects of distance variation on the impingement process of compressible vortex rings on a stationary solid wall. An experimental incident Mach number of 1.61 was used. Qualitative and quantitative studies compared the impingement and interaction flow characteristics of a compressible vortex ring with a stationary, solid, smooth surface at three distances: 1.66, 3.33, and 5.00 inner diameters. The three distances corresponded to an under-developed vortex ring (1.66 inner diameters), a vortex ring at its development threshold (3.33 inner diameters), and a fully developed one (5.00 inner diameters). Qualitative schlieren results showed that the surface distance affected the shock/vortex interaction process along with the impingement process of the vortex ring and the flow structure of its trailing jet. Quantitative data were extrapolated to evaluate the propagation velocity of the vortex ring prior to impingement. The boundary layer thickness was also estimated. Particle image velocimetry studies showed the main and secondary vortices to have opposite vorticity, with the magnitude of the vorticity of the secondary vortices being approximately half of that of the main vortex. Surface pressure results reveal the symmetrical properties of the impinging flow, along with a direct correlation between the maximum pressure measured at the instant the vortex ring impingement and an increase in surface distance.     

A combined vortex and panel method for numerical simulations of viscous flows: a comparative study of a vortex particle method and a finite volume method

This paper describes and compares two vorticity‐based integral approaches for the solution of the incompressible Navier–Stokes equations. Either a Lagrangian vortex particle method or an Eulerian finite volume scheme is implemented to solve the vorticity transport equation with a vorticity boundary condition. The Biot–Savart integral is used to compute the velocity field from a vorticity distribution over a fluid domain. The vorticity boundary condition is improved by the use of an iteration scheme connected with the well‐established panel method. In the early stages of development of flows around an impulsively started circular cylinder, and past an impulsively started foil with varying angles of attack, the computational results obtained by the Lagrangian vortex method are compared with those obtained by the Eulerian finite volume method. The comparison is performed separately for the pressure fields as well. The results obtained by the two methods are in good agreement, and give a better understanding of the vorticity‐based methods. Copyright © 2005 John Wiley & Sons, Ltd.