On variational features of vortex flows

Ideal incompressible fluid is a Hamiltonian system which possesses an infinite number of integrals, the circulations of velocity over closed fluid contours. This allows one to split all the degrees of freedom into the driving ones and the “slave” ones, the latter to be determined by the integrals of motions. The “slave” degrees of freedom correspond to “potential part” of motion, which is driven by vorticity. Elimination of the “slave” degrees of freedom from equations of ideal incompressible fluid yields a closed system of equations for dynamics of vortex lines. This system is also Hamiltonian. The variational principle for this system was found recently (Berdichevsky in Thermodynamics of chaos and order, Addison-Wesly-Longman, Reading, 1997; Kuznetsov and Ruban in JETP Lett 67, 1076–1081, 1998). It looks striking, however. In particular, the fluid motion is set to be compressible, while in the least action principle of fluid mechanics the incompressibility of motion is a built-in property. This striking feature is explained in the paper, and a link between the variational principle of vortex line dynamics and the least action principle is established. Other points made in this paper are concerned with steady motions. Two new variational principles are proposed for steady vortex flows. Their relation to Arnold’s variational principle of steady vortex motion is discussed.      

Topological structures of vortex and helicity analysis

In this paper, we discuss the topological structures of the vortex filaments and vortex tubes with an exact solution of a straight spiral vortex tube. We find that there are some confusions about the calculation of the helicities of a knotted vortex filament and some linked vortex filaments by using different methods. We explain how to unify these methods and give the right results. This work is supported by the National Basic Research Important Project “Nonlinear Science”     

A model of stratified entrainment using vortex persistence

A new model is proposed for the entrainment rate by vortices across stratified interfaces. In the model, different entrainment regimes are distinguished by the conventional parameters Richardson, Reynolds, and Schmidt number as well as a new parameter, the “vortex persistence”. Vortex persistence is defined as the number of rotations a vortex makes during the time it moves its own diameter with respect to the interface. It is further proposed that the concept of vortex persistence is important whenever a vortex is near any kind of surface, either stratified or solid. The model is in accord with most field and laboratory observations in a variety of stratified and bounded flows, including measurements of wall heat transfer and vortex formation in starting jets.     

The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

Dual leading-edge vortex structure for flow over a simplified butterfly model

The dye visualization experiments show that a dual leading-edge vortex (LEV) structure exists on the suction side of a simplified butterfly model of Papilio ulysses at α = 8°−12°. Furthermore, the results of particle image velocimetry (PIV) measurement indicate that the axial velocity of the primary (outer) vortex core reaches the lower extreme value while a transition from a “wake-like” to a “jet-like” axial velocity profile occurs. The work reveals for the first time the existence of dual LEV structure on the butterfly-like forward-sweep wing configuration.     

The nonlinear interaction of vortex rings with a free surface

Nonlinear interactions of vortex rings with a free surface are considered in an incompressible, ideal fluid using the vortex contour dynamics technique and the boundary integral equation method. The flow is axisymmetric and the vorticity is linearly distributed in the vortex. Effects of the gravity and the surface tension as well as the initial geometric parameter of the vortex on the interaction process are investigated in considerable detail. The interaction process may be divided into three major stages: the vortex free-traveling stage, the collision stage, and the vortex stretching and rebounding stage. Time evolutions of both the vortex and free surface under various conditions are provided and analyzed. Two kinds of waves exist on the free surface during interaction. In a special case where the gravity and surface tension are very weak or the vortex is very strong, an electric-bulb-like ‘cavity’ is formed on the free surface and the vortex is trapped in the ‘cavity’ for quite a long time, resulting in a large amount of fluid above the mean fluid surface. The project supported by the National Education Commission of China and NASA under cooperative grant agreement # NCC5-34     

On the formation of vortex rings and pairs

The axisymmetric vortex sheet model developed by Nitsche & Krasny (1994) has been extended to study the formation of vortex rings (pairs) at the edge of circular (2D) tube and opening. Computations based on this model are in good agreement with the experiments (Didden (1979) for circular tube and Auerbach (1987) for 2D tube and opening). Using this new model, evidences are provided to show that the main failure of the similarity theory (the false prediction of axial trajectory of vortex ring) is due to its ignorance of the self-induced ring velocity (mutual induction for vortex pair). We further reason why the similarity theory succeeds in its prediction of radial movement of vortex ring. The effects of various parameters such as turning angle α and piston speedU _p (t) on the formation of vortex ring are investigated. Numerical result shows that turning angle α has no effect on circulation shed τ. We also discuss Glezer (1988)'s summary on the influence ofU _p upon the shedding circulation, and finally give the variation of core distribution of vortex ring with α andU _p (t). The project is supported by National Natural Science Foundation of China and Doctoral Program of Institution of Higher Education     

Flow structures and force characteristics for flat plate in oscillatory flows withKc number from 2 to 40 and in combined flows

The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming flow cases, wyenKc number varies from 2 to 40, the vortex pattern changes from a “harmonic wave” shaped (in a range of smallKc number) to a slight inclined “harmonic wave” shaped (in a range of moderateKc numbers), then to inclined vortex clusters with an angle of 50° to the oncoming flow direction (atKc=20), at last, asKc number becomes large, the vortex pattern is like a normal Karman vortex street. The well predicted drag and inertia force coefficients are obtained, which are more close to the results of Keulegan & Carpenter's experiment as compared with previous vortex simulation by other authors. The existence of minimum point of inertia force coefficientC _m nearKc=20 is also well predicted and this phenomenon can be interpreted according to the vortex structure. For steady-oscillatory in-line combined flow cases, the vortex modes behave like a vortex street, exhibit a “longitudinal wave” structure, and a vortex cluster shape corresponding to the ratios ofU _m toU ₀ which are ofO (10⁻¹)O(1) andO(10), respectively. The effect on the prediction of forces on the flat plate from the disturbance component in a combined flow has been demonstrated qualitatively. In addition to this, the lock in phenomenon of vortex shedding has been checked. The project supported by National Natural Science Foundation of China & LNM, Institute of Mechanics, CAS     

Using the perturbation method to solve the problem of separated incompressible flow past thin airfoils

A model for separated incompressible flow past thin airfoils in the neighborhood of the “shockless entrance” condition is constructed based on the averaging of the vortex shedding flow past the airfoil edges. By approximation of the vortex shedding by two vortex curves, determination of the average hydrodynamic parameters is reduced to a twofold solution of an integral singular equation equivalent to the equation describing steady-state nonseparated airfoil flow. In this case, the calculation time is two orders of magnitude smaller than the time required for the solution of the corresponding evolution problem. The results of a test calculation using the proposed method are in fair agreement with available results of calculations and experiments. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 49–63, May–June, 2006.     

Model of vortex ring formation

A semi-empirical model of vortex ring formation during exhaustion of a pulsed submerged jet from a circular orifice is presented. Formulas for determining the parameters of the vortex ring, depending on the conditions of formation of the latter, are derived. The calculated characteristics of the vortex ring as functions of criteria determining the vortex formation process are demonstrated to be in good agreement with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 6, pp. 25–36, November–December, 2008.     

Primary and secondary vortical structures contribution in the entrainment of low Reynolds number jet flows

Particle image velocimetry measurements and time-resolved visualization are used for the reconstruction of the Kelvin–Helmholtz vortex passing in the near field of a round jet and of a lobed jet. For the round jet, the entrainment is produced in the braid region, where streamwise structures develop. In the Kelvin–Helmholtz ring, entrainment is dramatically affected by the attenuation of the streamwise structures. As for the lobed jet, the special geometry introduces a transverse shear leading to a breakdown of the Kelvin–Helmholtz structures into “ring segments.” Streamwise structures continuously develop at the resulting discontinuity regions and control the lobed jet self-induction. In this case, the entrainment rate is less affected by the primary structures dynamics.     

Modeling of jet flows of a viscous fluid by the discrete vortex method

Results obtained previously by the discrete vortex method with a “viscous” correction are generalized. The boundaries of applicability of this method are determined. Previous results obtained for a flow past a flat plate are supplemented with solution convergence estimates. Exhaustion of a plane jet of a viscous incompressible fluid into the ambient space is modeled. The geometric parameters of the jet (its half-width, shapes of the streamwise velocity profiles, and intensity of oscillations) are analyzed. The calculated results are found to agree well with experimental data and with results calculated by other methods.     

Loss of energy during the motion of a vortex ring

Experimental estimates of energy and energy dissipation of a vortex ring are presented. The energy losses during the motion of a vortex ring and a streamlined solid are compared. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 1, pp. 24–30, January–February, 2008.     

Measurements of Tip Vortex Characteristics and the Effect of an Anti-Cavitation Lip on a Model Kaplan Turbine Blade

Motivated by the problem of cavitation erosion, the position and strength of the roll-up vortex on the suction side of a Kaplan-type turbine blade with various endings is investigated. Measurements are made on different models with simplified two-dimensional geometries using mainly particle image velocimetry. It is found that the greatest danger of cavitation erosion exists when the casing of the turbine is of the “semi-spherical” type, so that there is a large and closing clearance gap at the leading edge of the blade. Furthermore, the anti-cavitation lip used in this study is shown to be ineffective at increasing the distance of the vortex from the blade, although it does reduce the circulation of the vortex and presumably the danger of cavitation. Our measurements are found to be in good agreement with existing models for the position of the vortex when appropriately interpreted. This revised version was published online in July 2006 with corrections to the Cover Date.     

Rising vortex ring in a viscous fluid

Most theoretical results for thermals, whose motion is determined by the complex interaction between dynamics and buoyancy, have been obtained numerically [1–4]. The analytic solutions for a convection element have been limited to consideration of the self-similar regime [5]. At the same time, the preself-similar stage of development of a vortex ring of dynamic origin has been described analytically [6]. This approach is now extended to a rising vortex ring. In this case a modification of the traditional formulation of the problem makes it possible to obtain an analytic solution of the problem of a weak thermal in the form of unsteady temperature, vorticity and stream function fields that tend in the limit to the self-similar regime. The rate of ascent of the convective vortex ring is found. A solution is obtained for the two-dimensional analog of the problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–48, May–June, 1989.     

Impulse generated during unsteady maneuvering of swimming fish

The relationship between the maneuvering kinematics of a Giant Danio (Danio aequipinnatus) and the resulting vortical wake is investigated for a rapid, ‘C’-start maneuver using fully time-resolved (500 Hz) particle image velocimetry (PIV). PIV illuminates the two distinct vortices formed during the turn. The fish body rotation is facilitated by the initial, or “maneuvering” vortex formation, and the final fish velocity is augmented by the strength of the second, “propulsive” vortex. Results confirm that the axisymmetric vortex ring model is reasonable to use in calculating the hydrodynamic impulse acting on the fish. The total linear momentum change of the fish from its initial swimming trajectory to its final swimming trajectory is balanced by the vector sum of the impulses of both vortex rings. The timing of vortex formation is uniquely synchronized with the fish motion, and the choreography of the maneuver is addressed in the context of the resulting hydrodynamic forces.     

Thin‐tube vortex simulations for sinusoidal instability in a counter‐rotating vortex pair

A thin‐tube vortex method is developed to investigate the intrinsic instability within a counter‐rotating vortex pair system and the effects from the core size and the wavenumbers (or wavelengths). The numerical accuracy and the advantages of the scheme are theoretically estimated. A nearest‐neighbour‐image method is employed in this three‐dimensional vortex simulation. Agreement with Crow's instability analysis has been achieved numerically for the long‐wave cases. A short‐wave instability for the zeroth radial mode of bending instability has also been found using the thin‐tube vortex simulations. Then, the combinations of long‐ and short‐wave instability are investigated to elucidate the non‐linear effects due to the interactions of two different modes. It is shown that instability is enhanced if both long‐ and short‐wave instabilities occur simultaneously. Although the method used in the paper is not capable of including effects such as axial flow, vortex core deformation and other complicated viscous effects, it effectively predicts and clarifies the first‐order factor that dominates the sinusoidal instability behaviour in a vortex pair. Copyright © 2002 John Wiley & Sons, Ltd.     

A study on the mechanism of high-lift generation by an airfoil in unsteady motion at low reynolds number

The aerodynamic force and flow structure of NACA 0012 airfoil performing an unsteady motion at low Reynolds number (Re=100) are calculated by solving Navier-Stokes equations. The motion consists of three parts: the first translation, rotation and the second translation in the direction opposite to the first. The rotation and the second translation in this motion are expected to represent the rotation and translation of the wing-section of a hovering insect. The flow structure is used in combination with the theory of vorticity dynamics to explain the generation of unsteady aerodynamic force in the motion. During the rotation, due to the creation of strong vortices in short time, large aerodynamic force is produced and the force is almost normal to the airfoil chord. During the second translation, large lift coefficient can be maintained for certain time period and , the lift coefficient averaged over four chord lengths of travel, is larger than 2 (the corresponding steady-state lift coefficient is only 0.9). The large lift coefficient is due to two effects. The first is the delayed shedding of the stall vortex. The second is that the vortices created during the airfoil rotation and in the near wake left by previous translation form a short “vortex street” in front of the airfoil and the “vortex street” induces a “wind”; against this “wind” the airfoil translates, increasing its relative speed. The above results provide insights to the understanding of the mechanism of high-lift generation by a hovering insect. The project supported by the National Natural Science Foundation of China (19725210)     

Compressible and confined vortex flow

The internal compressible flow of a thin vortex chamber was investigated experimentally by measuring the radial distribution of temperature and pressure, from which the velocity field was calculated. The bulk of the internal vortex was found to be described by u_θr^0.69 = constant. The total resistance of the vortex chamber to the flow was also investigated in the context of fluidic vortex diode behavior under conditions of compressible and choked flow. It was found that the vortex chamber choked at an upstream-to-downstream pressure ratio of about 6 and in doing so passed a mass flow rate of 28% of the equivalent one-dimensional ideal nozzle. The resistance of vortex chambers is known to be strongly influenced by the presence of reversed flow in the exit due to vortex breakdown. Schlieren photography of the swirling exhaust flow was used to show that, while vortex breakdown does occur, it can only do so after the flow has become subsonic downstream of the exit and cannot therefore influence the vortex chamber resistance.     

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.