Flow past a sphere in density-stratified fluid

Stratified flow past a three-dimensional obstacle such as a sphere has been a long-lasting subject of geophysical, environmental and engineering fluid dynamics. In order to investigate the effect of the stratification on the near wake, in particular, the unsteady vortex formation behind a sphere, numerical simulations of stratified flows past a sphere are conducted. The time-dependent Navier–Stokes equations are solved using a three-dimensional finite element method and a modified explicit time integration scheme. Laminar flow regime is considered, and linear stratification of density is assumed under Boussinesq approximation. The effects of stratification is implemented by density transport without diffusion. The computed results include the characteristics of the near wake as well as the effects of stratification on the separation angle. Under increased stratification, the separation on the sphere is suppressed and the wake structure behind the sphere becomes planar, resembling that behind a vertical cylinder. With further increase in stratification, the wake becomes unsteady, and consists of planar vortex shedding similar to von Karman vortex streets.     

Experiments on the formation of a recirculation zone in swirling coaxial jets

 The near flow field of coaxial air jets, with swirl in the outer one, was studied using flow visualization and hot-wire anemometry. The flow is sensitive to both the swirl number and the mass flow ratio between the outer and inner jets. A necessary condition for the formation of an internal recirculation zone (IRZ) is that the swirl number must exceed a minimum value which depends on the mass flow ratio. Spectral analysis of the velocity fluctuations indicates that the formation of an IRZ in the present flow does not appear to be related to the growth of convective flow instabilities. Analysis of the flow visualization and X-wire data indicates that the vorticity dynamics model for vortex breakdown proposed by Brown and Lopez [J Fluid Mech (1990) 222: 553] provides a plausible mechanism for the formation of an IRZ in this flow. Received: 14 June 1999/Accepted: 7 December 1999     

Acoustics and dynamics of coaxial interacting vortex rings

Using a contour dynamics method for inviscid axisymmetric flow we examine the effects of core deformation on the dynamics and acoustic signatures of coaxial interacting vortex rings. Both “passage” and “collision” (head-on) interactions are studied for initially identical vortices. Good correspondence with experiments is obtained. A simple model which retains only the elliptic degree of freedom in the core shape is used to explain some of the calculated features.     

Stability of a viscous subsonic swirl flow

The results of calculating the stability of a three-dimensional swirl flow of a viscous heat-conducting gas are presented. The stability characteristics are determined using the linear time-dependent theory of plane-parallel flow stability. The main undisturbed axisymmetric vortex flow was determined numerically using a quasi-cylindrical approximation for the complete set of Navier-Stokes equations. The circulation of the peripheral velocity in the cocurrent flow surrounding the viscous vortex core was assumed to be constant. In analyzing the stability, nonaxisymmetric perturbations in the shape of waves traveling along the vortex axis with both positive and negative wavenumbers were considered; in these two cases the perturbation rotation is either the same or opposite in sense to the rotation in the vortex core. Neutral stability curves are determined for various values of the swirling parameter and the cocurrent flow Mach number. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 50–59, May–June, 1998.     

Measurements of a stator-induced circumferentially varying flow

The flow physics associated with the generation of both axisymmetric and non-axisymmetric swirl by various deflection patterns of a stator array was investigated experimentally through surface pressure and Stereoscopic Particle Image Velocimetry measurements. A three-dimensional rendering technique was developed to reconstruct the flow field around the model and in its wake. The three-dimensional fluid volume was reconstructed from multiple two-dimensional measurement planes. A cyclic distribution of the stators’ deflections resulted in non-axisymmetric distributions of the surface pressure and the flow field downstream of the stator array. The addition of a shroud had an amplifying effect: accelerating the flow through the stator array while reducing the non-uniform tangential velocity component generated by the stators. In the model near wake the flow field is associated with secondary flow patterns in the form of coherent streamwise vortical structures that can be described by potential flow mechanisms. The collective pitch distribution of the stators produces a flow field that resembles a potential Rankine vortex, whereas the cyclic pitch distribution generates a flow pattern that can be described by a potential vortex pair in a cross-flow.     

Reynolds stress modelling of jet and swirl interaction inside a gas turbine combustor

Influences of the inlet swirl levels on the interaction between the dilution air jets and the swirling cross‐flow to the interior flow field inside a gas turbine combustor were investigated numerically by Reynolds stress transport model (RSTM). Due to the intense swirl and jet interaction, a high level of swirl momentum is transported to the centreline and hence, an intense vortex core is formed. The strength of the centreline vortex core was found to depend on the inlet swirl levels. For the higher swirling inlet, the decay of the swirling motion causes strong streamline variation of pressure; and consequently leads to an elevated level of deceleration of its axial velocity. Predictions contrasted with measurements indicate that the stress model reproduces the flow correctly and is able to reflect the influences of inlet swirl levels on the interior flow structure. Copyright © 1999 John Wiley & Sons, Ltd.     

Near-ground tornado-like vortex structure resolved by particle image velocimetry (PIV)

The near-ground flow structure of tornadoes is of utmost interest because it determines how and to what extent civil structures could get damaged in tornado events. We simulated tornado-like vortex flow at the swirl ratios of S = 0.03–0.3 (vane angle θ_v = 15°–60°), using a laboratory tornado simulator and investigated the near-ground-vortex structure by particle imaging velocimetry. Complicated near-ground flow was measured in two orthogonal views: horizontal planes at various elevations (z = 11, 26 and 53 mm above the ground) and the meridian plane. We observed two distinct vortex structures: a single-celled vortex at the lowest swirl ratio (S = 0.03, θ_v = 15°) and multiple suction vortices rotating around the primary vortex (two-celled vortex) at higher swirl ratios (S = 0.1–0.3, θ_v = 30°–60°). We quantified the effects of vortex wandering on the mean flow and found that vortex wandering was important and should be taken into account in the low swirl ratio case. The tangential velocity, as the dominant velocity component, has the peak value about three times that of the maximum radial velocity regardless of the swirl ratio. The maximum velocity variance is about twice at the high swirl ratio (θ_v = 45°) that at the low swirl ratio (θ_v = 15°), which is contributed significantly by the multiple small-scale secondary vortices. Here, the results show that not only the intensified mean flow but greatly enhanced turbulence occurs near the surface in the tornado-like vortex flow. The intensified mean flow and enhanced turbulence at the ground level, correlated with the ground-vortex interaction, may cause dramatic damage of the civil structures in tornadoes. This work provides detailed characterization of the tornado-like vortex structure, which has not been fully revealed in previous field studies and laboratory simulations. It would be helpful in improving the understanding of the interaction between the tornado-like vortex structure and the ground surface, ultimately leading to better predictions of tornado-induced wind loads on civil structures.     

DNS of flows with helical symmetry

Some flows such as the wakes of rotating devices often display helical symmetry. We present an original DNS code for the dynamics of such helically symmetric systems. We show that, by enforcing helical symmetry, the three-dimensional Navier–Stokes equations can be reduced to a two-dimensional unsteady problem. The numerical method is a generalisation of the vorticity/streamfunction formulation in a circular domain, with finite differences in the radial direction and spectral decomposition along the azimuth. When compared to a standard three-dimensional code, this allows us to reach larger Reynolds numbers and to compute quasi-steady patterns. We illustrate the importance of helical pitch by some physical cases: the dynamics of several helical vortices and a quasi-steady vortex flow. We also study the linear dynamics and nonlinear saturation in the Batchelor vortex basic flow, a paradigmatic example of trailing vortex instability. We retrieve the behaviour of inviscid modes and present new results concerning the saturation of viscous centre modes.     

The nonlinear evolution of swirling jets

We investigate the mechanisms of vorticity concentration, reorientation and stretching in a swirling jet, whose dynamics is dominated by the competition of a Kelvin-Helmholtz-type vortex sheet instability and a centrifugal Rayleigh instability. To this end, we employ an inviscid Lagrangian vortex filament technique. It is found that the axial jet velocity profile breaks the symmetry of the pure swirling flow. Conversely, the swirl is seen to modify the case dominated by a Kelvin-Helmholtz instability in that it results in the formation of counterrotating vortex rings. A pinch-off mechanism is observed which leads to a dramatic decrease in the local jet diameter. Furthermore, the vortex ring circulation is seen to be time dependent.

---

Sommario In questo lavoro si analizza la dinamica della vorticità in un setto rotante in cui siano presenti, ed in competizione reciproca, fenomeni di instabilità di Kelvin-Helmholtz e di Rayleigh. A tale scopo si adotta una metodologia di soluzione non viscosa, Lagrangiana a filamenti vorticosi. Viene mostrato come il profilo di velocità assiale del getto altera la simmetria del moto di pura rotazione. Viceversa, la presenza della rotazione modifica il flusso dominato dall'instabilità di Kelvin-Helmholtz attraverso la formazione di anelli vorticosi controrotanti. L'interazione di questi due campi di velocità porta sia ad una considerevole riduzione del diametro locale del getto, sia ad una variazione temporale della circolazione degli anelli vorticosi.

     

圆环旋转黏性液体射流空间不稳定性研究

利用线性稳定性理论, 进行了液体黏性对不同旋转强度下圆环旋转液体射流 空间不稳定性影响的研究. 在推导出的三维扰动下具有固体涡核型旋转速度分布的圆环旋转 黏性液体射流色散方程的基础上, 针对中低速射流, 进行了类反对称模式与类对称模式下圆 环旋转黏性液体射流的空间不稳定性分析. 研究结果表明, 对于旋转强度较大的圆环旋转液 体射流, 液体黏性的增加, 不利于射流的破碎; 随着液体黏性的增加, 射流的特征频率和最 不稳定波数减小. 然而, 对于旋转强度较小的圆环旋转液体射流, 液体黏性的增加, 有利于 射流的破碎; 随着液体黏性的增加, 类反对称模式下射流特征频率先减小后增大, 类对称模 式下射流特征频率增大; 随着液体黏性的增加, 类反对称模式下射流最不稳定波数先减小后 增大, 类对称模式下射流最不稳定波数增大.     

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.     

An analysis of unsteady highly turbulent swirling flow in a model vortex combustor

This paper reports an experimental investigation of a non-reacting turbulent swirling flow in a practical vortex combustor. The flow was examined for the conditions characteristic of the presence of a breakdown zone and a strong flow instability appearing at swirl numbers S>0.5. Flow visualization techniques, LDA measurements and acoustic probes were employed to study the unsteady flow characteristics. Based on the experimental results a positive first helical mode of instability was identified with a wavelength and frequency depending on swirl. The wavelength was confirmed to grow monotonically with S, while the dominant frequency of the flow pulsations was found to have an unusual parabolic evolution with swirl, with a minimum at S _min=0.88. This finding was interpreted using a proposed kinematic model based on the contribution of two mechanisms: rotation and axial motion of the helical vortex. It was concluded that for S<S _min the instability frequency is essentially dominated by the axial translation of the spiral vortex being inversely proportional to S and therefore giving a decreasing trend. For S>S _min the frequency of the flow precession is more dependent on the angular transportation of the vortex core, which resulted in the expected growing dependence on S.     

Influence of flow field scaling on flashback of swirl flames

In this paper, the effect of geometrical scaling on the onset of flashback into a cylindrical premixing zone of a swirl flame is investigated. We discriminate two types of flashback. In the first type of flashback the flame propagates upstream inside an already present axial recirculation zone. This flashback is caused by turbulent burning along the vortex axis (TBVA¹) and is controlled by flame extinction inside the recirculation zone. The second type of flashback is caused by combustion induced vortex breakdown (CIVB²). This type of flashback is characterised by the aerodynamic influence of the combustion heat release that leads to propagation of the axial recirculation zone and the flame in upstream direction.To study the effects of geometrical scaling on the flow fields and the two types of flashback, the operation of two geometrically scaled burners are compared at equal Reynolds number. By this method it is possible to observe the flashback phenomena in similar swirl flow fields but with different turbulent scales affecting the combustion process. To check flow field similarity and to indentify the flashback type, the non-reacting and reacting flow fields have been examined by planar particle imaging velocimetry and simultaneous recording of the flame luminescence.It is shown that geometrical scaling of the burner shifts the equivalence ratio at which flashback occurs and that this shift is different for the two types of flashback. Consistency and inconsistency with known scaling and stability criterions is discussed. Analysing the fluid dynamics and turbulent combustion gives a first explanation of why CIVB and TBVA are affected differently by geometrical scaling at constant Reynolds number which is in good agreement with the experimental observations.     

Some features of the flow around rapidly rotating bodies made of cellular-porous materials

Two types of gas flows arising near a rapidly rotating cellular-porous disk are studied numerically and experimentally. Steady-state limits for the flow around a disk rotating in free space and the type and scenario of the loss of stability are determined. Transitional flows are characterized by formation of a vortex sheet at the boundary of the exhausting jet. Numerical simulations of the flow around a cellular-porous disk rotating near a flat screen show that it is possible to form a closed swirl flow responsible for redistribution of swirl in the gap between the disk and the flat screen. The computed results offer an explanation for the experimentally observed excess of tangential velocity of the flow in the gap over the velocity of disk rotation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 86–96, January–February, 2007.     

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.     

Obstacle-induced spiral vortex breakdown

An experimental investigation on vortex breakdown dynamics is performed. An adverse pressure gradient is created along the axis of a wing-tip vortex by introducing a sphere downstream of an elliptical hydrofoil. The instrumentation involves high-speed visualizations with air bubbles used as tracers and 2D Laser Doppler Velocimeter (LDV). Two key parameters are identified and varied to control the onset of vortex breakdown: the swirl number, defined as the maximum azimuthal velocity divided by the free-stream velocity, and the adverse pressure gradient. They were controlled through the incidence angle of the elliptical hydrofoil, the free-stream velocity and the sphere diameter. A single helical breakdown of the vortex was systematically observed over a wide range of experimental parameters. The helical breakdown coiled around the sphere in the direction opposite to the vortex but rotated along the vortex direction. We have observed that the location of vortex breakdown moved upstream as the swirl number or the sphere diameter was increased. LDV measurements were corrected using a reconstruction procedure taking into account the so-called vortex wandering and the size of the LDV measurement volume. This allows us to investigate the spatio-temporal linear stability properties of the flow and demonstrate that the flow transition from columnar to single helical shape is due to a transition from convective to absolute instability.     

Some Observations of Vortex Breakdown in a Confined Flow with Solid Body Rotation

The occurrence of breakdown in slender vortex flows as a ``bubble'' or ``spiral'' pattern depends on the degree of radial deflection of the vortex core from its original axis as shown in [1]. A smooth transition from a bubble to a spiral-type ``mode'' can be forced by inducing a small asymmetric disturbance which led to the conclusion, that the patterns do not represent different fundamental modes of breakdown. The subject presented herein addresses the following question: how does breakdown evolve in a swirling flow in which the vortex core is forced on a straight axis? In addition, what is the effect of turbulent inflow conditions? This type of vortex conditions is achieved in a spinning tube flow. The swirl is introduced at the entrance of the rotating tube with a honeycomb package and maintained by the viscous action in the boundary layer of the spinning tube. A diffuser at the end induces an adverse pressure gradient to force the breakdown. Flow visualization experiments are carried out to characterize the nature of breakdown over a range of different flow conditions. For some selected characteristic stages, detailed velocity fields were obtained using the method of Digital Particle-Image-Velocimetry (DPIV). The results show, that for the range of parameters investigated, breakdown is initiated at Rossby-numbers below a critical value of Ro ≈ 0.6 similar to those observed in other experiments. The bursted part of the vortex has a near axi-symmetric slender conical shape containing approximately stagnant flow. Its downstream end is characterized by a jump-like contraction where the flow evolves into a jet with enhanced swirl on the axis. It is only in this region downstream of the jump-like contraction that asymmetric instabilities and wavy flow patterns could be observed. Perturbations caused by them travel upstream but do not change the near-axisymmetric shape of the bursted part of the vortex.     

Flow in a simple swirl chamber with and without controlled inlet forcing

Results are presented from a swirl chamber with and without controlled inlet forcing. The controlled inlet forcing is induced using arrays of vortex generators placed along one wall of the swirl chamber inlet duct. Flow visualization results are given, along with surveys of circumferential mean velocity, static pressure, and total pressure, at Reynolds numbers (based on inlet duct characteristics) as high as 8000. The controlled inlet forcing provides means to alter and control: (i) the spacing and number of Görtler vortices across the span of the swirl chamber, (ii) the amount of vortex development at a particular Reynolds number and circumferential location, (iii) the circumferential location and Reynolds number of initial Görtler vortex development, and (iv) the circumferential location and Reynolds number of Görtler vortex breakup into more chaotic flow.     

Dynamics of annular gas–liquid two-phase swirling jets

The dynamics of annular gas–liquid two-phase swirling jets have been examined by means of direct numerical simulation and proper orthogonal decomposition. An Eulerian approach with mixed-fluid treatment, combined with an adapted volume of fluid and a continuum surface force model, was used to describe the two-phase flow system. The unsteady, compressible, three-dimensional Navier–Stokes equations have been solved by using highly accurate numerical methods. Two computational cases have been performed to examine the effects of liquid-to-gas density ratio on the flow development. It was found that the higher density ratio case is more vortical with larger spatial distribution of the liquid, in agreement with linear theories. Proper orthogonal decomposition analysis revealed that more modes are of importance at the higher density ratio, indicating a more unstable flow field. In the lower density ratio case, both a central and a geometrical recirculation zone are captured while only one central recirculation zone is evident at the higher density ratio. The results also indicate the formation of a precessing vortex core at the high density ratio, indicating that the precessing vortex core development is dependent on the liquid-to-gas density ratio of the two-phase flow, apart from the swirl number alone.     

The effects of inlet conditions and expansion ratio on the onset of flow reversal in swirling flow in a sudden expansion

Experiments are reported in which the minimum swirl intensity required to produce a central recirculation zone in a swirling sudden expansion flow is determined as a function of expansion ratio and inlet conditions. Using a swirl generator which allows for independent variation of velocity profile shape and swirl number, it is shown that an inlet tangential velocity distribution with a large solid body vortex core or an axial velocity profile with a maximum on the axis will lead to a higher critical swirl.