Study of vortex breakdown by particle tracking velocimetry (PTV) Part 3: Time-dependent structure and development of breakdown-modes

This third part of the study deals with the time-dependent nature of vortex breakdown. The results show the unsteady velocity and vorticity field of the initiation and development of breakdown and transition between both predominant breakdown modes, the bubble and the spiral. During the development to breakdown, the generated amount of circumferential vorticity follows the theoretical prediction by Brown and Lopez (1990). This confirms the idea of positive feedback as the key-mechanism leading to vortex breakdown. We regard the bubble-type as the fundamental breakdown type, that is stationary and nearly axisymmetric. The circumferential vorticity is distributed in a form of an elliptical vortex-ring-like structure. Starting from this stage, an increase of volume flux to a higher Reynolds number leads to the transition to the spiral-type with an initial stretching of the vortex ring-like structure and a subsequent change to an asymmetric circumferential vorticity distribution. This in combination with the inductive effect causes the front stagnation point to be deflected radially away and later to rotate around the centerline. Consequently the approaching vortex core is radially deflected in opposite direction and evolves in a spiral path. The idea of a second positive feedback-mechanism gives a possible explanation for the transition. Following this theory the asymmetry of circumferential vorticity will trigger itself at a certain degree by the interaction with the inductively affected stagnation point and its influence on the approaching vortex core. This self-enhancing process will finally lead to the spiral-type breakdown in which the radial distance between rotating stagnation point and deflected vortex core is of the order of the characteristic vortex core radius. The reversed transition from the spiral to a stable bubble-type can be regenerated by decreasing the Reynolds number down to the value that corresponds to the stable bubble state. The flow structure evolves nearly in the time-reversed way as during transition from bubble towards the spiral.     

Downstream boundary effects on the spectral characteristics of a swirling flowfield

The spectral characteristics and the structural response of a swirling flowfield are investigated subject to a non-axisymmetric disturbance and a contraction imposed downstream. Two natural frequencies are noted in different regions of the undisturbed swirling flowfield, one is due to a precessing vortex core and the other to the most amplified downstream azimuthal instability. The downstream contraction usually causes compression of the central recirculation zone into two side-lobes, increases the dominant frequencies and forms a straight central vortex core with a high axial velocity. The dominant downstream instability frequency depends linearly on the inlet Reynolds number and on the mode of the breakdown. For the downstream non-axisymmetric disturbance, such as the passing of the turbine blades, the fundamental frequency is not altered by the disturbance and the oscillation strength of the downstream instability is greatly reduced as the excitation frequency remains unmatched with the dominant downstream natural frequency. Downstream azimuthal instability promotes the breakdown recirculation.A version of this paper was presented at the 26th AIAA Aerospace Sciences Meeting, Reno, Nevada, 11–14, Jan. 1988     

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.     

Axial plus tangential entry swirling jet

This paper presents an experimental investigation on swirling jets with well-defined initial conditions. The axial, radial, and azimuthal velocity components, with their respective fluctuations were measured using high spatial–resolution particle image velocimetry. These detailed measurements allow the initial conditions of the swirling jets to be established and the jets to be characterized using various swirl number definitions. The significance of each term in the swirl number calculations are quantified, and the effect of the common assumptions and simplifications are examined. The characteristics of the jets in relation to the initial conditions are then investigated and compared with the previous studies using similar characterization parameters. Jets with Reynolds number of approximately 5700 and swirl conditions ranging from a non-swirling reference case to high swirl are studied. General properties of swirling jets such as higher spreading rate, higher centerline velocity decay, and higher turbulence level are observed. When the degree of swirl is sufficiently high, vortex breakdown occurs. A swirl number of 0.94 is recorded for a high swirl case prior to vortex breakdown, much higher than the critical swirl number reported in the literature. This behavior is attributed to the effect of the initial conditions on the swirl number calculation.     

三角翼前缘涡的某些破裂形式及特性研究

基于流动显示和PIV技术测量的实验结果，对三角翼前缘涡破裂的一些形式和破裂特性进行了分析和讨论。通过PIV测量所得到的涡量分布证实了在螺旋破裂的情况下，涡核的螺旋方向与前缘涡的旋转方向相反，及双螺旋破裂形式的存在等。进而对螺旋波的形成机理提出了与有关文献不同的看法。     

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.     

Large Eddy Simulation of Low Swirl Flames Under External Flow Excitations

Low swirl flame characteristics under external flow excitations are numerically investigated using large eddy simulations with a dynamically thickened flame combustion model. A finite volume scheme on a Cartesian grid with a dynamic one equation eddy viscosity subgrid scale model is used for large eddy simulations. The excitations are imposed on inlet velocity profiles by a sinusoidal forcing function over a wide range of amplitudes and frequencies. Present investigation shows that although, the swirling motion of the low swirl flame is not intense enough to induce a recirculation zone in ensemble averaged results, external flow excitations increase the local swirl number upstream of the flame front. Such increase in the local swirl number is at maximum value when the low swirl flame is excited at the dominant frequency of the flow field, which in turn induces a vortex breakdown and hence a central recirculation zone. The strength and size of the time averaged recirculation zone depend on both the amplitude and frequency of the excitations. Moreover, phase-locked results indicate that external flow excitations induce local swirl fluctuations ahead of the flame front which alter the strength of the recirculation zone at different phase angles of the excitations.     

Formation and evolution of vortex rings induced by interactions between shock waves and a low-density bubble

The deformation and instability of a low-density spherical bubble induced by an incident and its reflected shock waves are studied by using the large eddy simulation method. The computational model is firstly validated by experimental results from the literature and is further used to examine the effect of incident shock wave strength on the formations and three-dimensional evolutions of the vortex rings. For the weak shock wave case (Ma?=?1.24), the baroclinic effect induced by the reflected shock wave is the key mechanism for the formation of new vortex rings. The vortex rings not only move due to the self-induced effect and the flow field velocity, but also generate azimuthal instability due to the pressure disturbance. For the strong shock wave case (Ma?=?2.2), a boundary layer is formed adjacent to the end wall owing to the approach of vortex ring, and unsteady separation of the boundary layer near the wall results in the ejection and formation of new vortex rings. These vortex rings interact in the vicinity of the end wall and finally collapse to a complicated vortex structure via azimuthal instability. For both shock wave strength cases, the evolutions of vortex rings due to the instability lead to the formation of the complicated structure dominated by the small-scale streamwise vortices.     

Interaction of heat release and vortex breakdown during flame flashback driven by combustion induced vortex breakdown

The interaction of heat release by chemical reaction and the flow dominates flame transition in swirling flows caused by combustion induced vortex breakdown (CIVB). The simultaneous application of 1 kHz high-speed particle imaging velocimetry (PIV) for the analysis of the flow field and OH planar laser-induced fluorescence for the detection of the flame front is particularly useful for the improvement of the understanding of the observed fast CIVB driven flame propagation. For the first time, the combination of both techniques was successfully applied to confined swirling flows. In the study, the flow field characteristics of an aerodynamically stabilized burner system with CIVB are analyzed in great depth. The influence of geometric parameters of the swirl generator was investigated and conclusions concerning the proper burner design of vortex breakdown premix burners are drawn from the experimental results. In particular, the effect of the vortex core with respect to the stability of the swirl stabilized burner is analyzed. The contribution of combustion to vortex breakdown is shown comparing isothermal and reacting flows. The presented data reveals that at the onset of CIVB driven flame transition, the azimuthal vorticity leads to the formation of a closed recirculation bubble at the tip of the internal recirculation zone. Once this bubble propagates upstream, the flame is able to follow and propagate relative to the bulk flow velocity with a velocity far beyond the turbulent flame speed. The interaction of reaction and flow was observed for different volumetric heat releases. The experiments confirm the CIVB theory of the authors, which was initially developed on the basis of a CFD study alone. Both the volume expansion and the baroclinic torque have an effect on whether fast flame propagation occurs. Whereas the volume expansion caused by the heat release stabilizes the flow field and the reaction, the baroclinic torque stimulates flame transition. For upstream propagation the flame tip has to have a position downstream of the stagnation point of the bubble. Else, the required transition inducing force is insufficient and the flame remains stable. In case the flame reaches positions too close or even upstream of the stagnation point, the fast propagation is interrupted or even prohibited. The key finding that the relative position of flame and stagnation bubble governs CIVB is discussed on the basis of high-speed LIF/PIV data as well as chemiluminescence. Since essentially the same behavior has been observed before in tests of a totally different swirler design and flow field, the conclusion can be made that the root cause for CIVB independent of the special geometry has been found.     

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.     

Absolute inviscid instability of a ring jet with back-flow and swirl

The absolute instability of a ring jet with back-flow and swirl has been investigated on the basis of the inviscid linearized theory for incompressible flow. An axisymmetric disturbance mode is found to be only convectively unstable. The first azimuthal mode can become absolutely unstable, if the ring jet has a back-flow on the jet axis, and an additional swirl can increase the instability. However, for large swirl the absolute instability is suppressed. A ring jet without back-flow becomes absolutely unstable only in the presence of swirl.     

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

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

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.     

Study of vortex breakdown by particle tracking velocimetry (PTV)

Using a quasi three-dimensional instantaneous measurement technique, which combines particle tracking velocimetry (PTV) with volume scanning, first quantitative experimental results of the unsteady and asymmetric interior region of vortex breakdown were obtained. The study was carried out in a low speed flow through a cylindrical tube. A vortex was generated by a set of guidevanes and subjected to an adverse pressure gradient causing its breakdown. By scanning a pulsed illuminated planar laser light sheet, a set of meridional and azimuthal cuts of the flow was obtained. With PTV the recorded particle paths in the cuts were processed in order to obtain the instantaneous two-dimensional velocity field, mean streamlines and vorticity distribution. Moreover, the three-dimensional shape of the appearing breakdown, visualized with fluorescent dye, was reconstructed from the cuts. The results revealed that the shape of the bubble nearly equals the streamsurface of the stagnation point. According to the conditions in the water tunnel a single tilted vortex ring at the open rear part of the bubble dominates the interior flow structure of the bubble as first noted by Sarpkaya (1971). The vortical flow is bulged over the bubble, restored and intensified at the lower end. The gathered data lead to the conclusion that the vortex axis remains parallel to the centerline.     

Shape Oscillations of Rising Bubbles

The paper details results from an experimental study on bubbles rising in still tap water. Shape and motion parameters of the bubbles were measured using a combination of high speed cinematography and digital image processing. The Reynolds numbers of the bubbles studied ranged from about 700 to 1300, with the bubbles exhibiting all the familiar shape and motion characteristics: oblate spheroids becoming wobbly, and spiralling or zig-zagging motion becoming rocking as the bubble size increased. Time series of the bubble major axes revealed regular oscillations in the bubble shape. In most cases three frequencies could be readily identified, corresponding to those of vortex shedding from the bubble and two modes of ellipsoidal harmonics (modes 2,0 and 2,2). Comparison of time series of bubble shape and motion indicated a strong interaction between the shape oscillations of mode 2,0 and bubble motion. As the bubble size increased the frequency of both shape oscillation modes approached that of the vortex shedding, which remained constant at about 12 Hz for all of our experiments. The frequencies become equal for bubbles larger than in our study, at a Reynolds number of about 3000. Using data from the literature we found that the vortex shedding appears to become locked-in on the mode 2,0 shape oscillation.     

Dynamic Response of Swirl Stabilized Turbulent Premixed Flames Based on the Helmholtz-Hodge Velocity Decomposition

The dynamic response of fully premixed flames stabilized in strongly swirled flows undergoing vortex breakdown is investigated with axisymmetric unsteady RANS simulations. The analysis relies on the well known Helmholtz-Hodge decomposition of the velocity field into its irrotational and rotational components. A novel methodology based on the linearization of the progress variable transport equation is developed to determine the separate contribution of these velocity components to the Flame Transfer Function (FTF). Due to the phase delay between the convected tangential velocity and instantaneously propagating axial velocity perturbations, a non-monotonic frequency dependence of the swirl number amplitude downstream the swirl generator is detected. In line with experimental observations, such non-monotonic frequency dependence is found also for the amplitude and phase of the FTF. This behaviour is associated here with rotational velocity perturbations generated by the Central Recirculation Zone (CRZ) generated by the phenomenon of vortex breakdown which, responding in a fashion totally similar to the swirl number perturbation, produces flame surface area fluctuations with the same distribution versus frequency.     

超声速尾涡的稳定性分析

从Ｎ－Ｓ方程出发,通过正则模方法,研究了超声速尾涡的绝对／对流不稳定性性质．计算了流动的稳定性特征随马赫数Ｍ,周向波数ｎ．,轴向自由流速度Ｗ０和旋转度ｑ等流动参数的变化规律,找到了绝对／对流不稳定区域的边界．通过比较发现,马赫数的增加使流动由绝对不稳定向对流不稳定乃至稳定转化．在所计算的参数范围,周向波数的增加加速了这一转化过程,而且,轴向速度的增加,同样使流动向着稳定的方向转化．同时还分析了不同旋拧程度的流动受可压缩影响的不同．这些结果对于了解旋拧流动稳定性的物理机理以及进行流动控制都有着重要意义．     

Vortex dynamics of turbulence–coherent structure interaction

We study the interaction between a coherent structure (CS) and imposed external turbulence by employing direct numerical simulations (DNS) designed for unbounded flows with compact vorticity distribution. Flow evolution comprises (i) the reorganization of turbulence into finer-scale spiral filaments, (ii) the growth of wave-like perturbations within the vortex core, and (iii) the eventual arrest of production, leading to the decay of ambient turbulence. The filaments, preferentially aligned in the azimuthal direction, undergo two types of interactions: parallel filaments pair to form higher-circulation “threads”, and anti-parallel threads form dipoles that self-advect radially outwards. The consequent radial transport of angular momentum manifests as an overshoot of the mean circulation profile—a theoretically known consequence of faster-than-viscous vortex decay. It is found that while the resulting centrifugal instability can enhance turbulence production, vortex decay is arrested by the dampening of the instability due to the “turbulent mixing” caused by instability-generated threads. Ensemble-averaged turbulence statistics show strong fluctuations within the core; these are triggered by the external turbulence, and grow even as the turbulence decays. This surprising growth on a normal-mode-stable vortex results from algebraic amplification through “linear transient growth”. Transient growth is examined by initializing DNS with the “optimal” modes obtained from linear analysis. The simulations show that the growth of transient modes reproduces the prominent dynamics of CS-turbulence interaction: formation of thread-dipoles, growth of core fluctuations, and appearance of bending waves on the column’s core. At the larger Reynolds numbers prevailing in practical flows, transient growth may enable accelerated vortex decay through vortex column breakdown.     

圆环旋转粘性液膜射流三维不稳定性研究

利用线性稳定性理论进行了射流液体粘性对圆环旋转液膜射流稳定性影响的研究,推导出了三维扰动下具有固体旋涡型速度分布的圆环旋转粘性液膜射流的色散方程;在此基础上进行了类反对称模式与类对称模式下的圆环旋转粘性液膜射流的三维不稳定性分析。研究结果表明,在类反对称模式下,液体粘性超过一定值后,射流最大扰动增长率随液体粘性的增加而迅速减小;轴对称模态的射流特征频率产生一个突降变化;随液体粘性增加,轴对称模态不稳定波数范围减小,非轴对称模态不稳定波数范围呈现出先减小后增大趋势。在类对称模式下,液体粘性对射流最大扰动增长率的影响主要体现在对非轴对称模态的影响上;液体粘性只在粘性较大时才会对非轴对称模态射流特征频率产生一定影响;液体粘性超过一定值后,轴对称模态与非轴对称模态的不稳定波数范围都会快速下降。     

脉冲激光等离子体与正激波相互作用的PIV实验研究

脉冲激光等离子体与超声速流场相互作用在飞行器减阻隔热、点火助燃等方面具有重要的应用价值.纹影实验方法只能定性或半定量地反映流动状态.为定量研究速度分布和旋涡结构,针对激光等离子体及其与正激波相互作用过程开展粒子图像测速PIV实验研究.在激波管实验平台上建立了纳秒脉冲激光能量沉积系统和PIV测量系统,通过定量测量,探明了激光等离子体引致的激光空气泡以及热核的流动特性,揭示了激光等离子体在正激波冲击下的流动特性与演化规律,并给出了激光能量大小和位置对相互作用过程的影响.结果表明:激光空气泡内的速度分布在激光入射方向上并不关于击穿点对称,而是在靠近激光入射方向一侧的流速略大于远离激光入射方向一侧;斜压导致热核在演化初期产生涡环,后期则由剪切主导;正激波与激光空气泡界面、热核界面相互作用时,产生斜压涡量,当激光能量为87.8 mJ、正激波马赫数1.4时,热核在正激波作用下产生的涡量比在静止空气中演化时大1个数量级;激光与正激波相互作用的关键过程是热核在正激波冲击下演化成涡环,在激波波前注入激光能量能够获得更加显著的涡环.