Simulation of vortex shedding behind a bluff body flame stabilizer using a hybrid U-RANS/PDF method

The present study aims at the investigation of the effects of turbulence-chemistry interaction on combustion instabilities using a probability density function(PDF) method.The instantaneous quantities in the flow field were decomposed into the Favre-averaged variables and the stochastic fluctuations,which were calculated by unsteady Reynolds averaged Navier-Stokes(U-RANS) equations and the PDF model,respectively.A joint fluctuating velocityfrequency-composition PDF was used.The governing equations are solved by a consistent hybrid finite volume/MonteCarlo algorithm on triangular unstructured meshes.A nonreacting flow behind a triangular-shaped bluff body flame stabilizer in a rectilinear combustor was simulated by the present method.The results demonstrate the capability of the present method to capture the large-scale coherent structures.The triple decomposition was performed,by dividing the coherent Favre-averaged velocity into time-averaged value and periodical coherent part,to analyze the coherent and incoherent contributions to Reynolds stresses.A simple modification to the coefficients in the turbulent frequency model will help to improve the simulation results.Unsteady flow fields were depicted by streamlines and vorticity contours.Moreover,the association between turbulence production and vorticity saddle points is illustrated.     

On the Flame-generated Vorticity Dynamics of Bluff-body-stabilized Premixed Flames

This investigation considers the dynamics of flame-generated vorticity for a premixed, submerged bluff-body stabilized flame. Experimentation characterizes the far-field region in particular with a level of detail not previously afforded to this type of flow. Simultaneous particle imaging velocimetry (PIV), Mie scattering and CH ^? chemiluminescence are used to obtain velocity fields and flame location. Mean static pressure measurements at the combustion chamber wall capture the pressure field. Analysis of the flame fronts in relation to the mean velocity and vorticity fields provides useful insight into the interaction of the flame and the flow. The unique nature of the velocity and vorticity fields and their effect on downstream flame structures are explained by the baroclinic torque generation of vorticity. The coupling that exists among pressure, heat release, and baroclinic generation is acknowledged and will influence strategies for control of the baroclinic mechanism.     

Large Eddy Simulation and Extended Dynamic Mode Decomposition of Flow-Flame Interaction in a Lean Premixed Low Swirl Stabilized Flame

In this paper, numerical studies are reported on the effect of flow-flame interaction at large and medium scales and its impact on flame stabilization in a lean premixed low swirl stabilized methane/air flame. The numerical simulations are based on a large eddy simulation (LES) approach with a three-scalar flamelet model with equations for mixture fraction and fuel mass fraction and the level-set G-equation to account respectively for stratification of the mixture, fuel leakage at the trailing edge of the flame, and tracking of the flame front. Distinct frequencies, associated with the flame stabilization process, are identified from point data of LES in the outer and inner shear layers of the burner induced flow field. To understand the effect of the spatial structures related to the observed flow frequencies, a dynamic mode decomposition (DMD) is performed. Based on the analysis of LES data, frequency specific coherent flow structures are extracted along with associated flame structures through an extended version of DMD. The inner shear layer generated vortices are associated with recurring frequency specific coherent structures of both flow and flame and contribute to the flame stabilization in the outer regions of the flow.     

Numerical Study of Turbulent Jet Ignition in a Lean Premixed Configuration

Direct numerical simulations (DNS) of a hot combustion product jet interacting with a lean premixed hydrogen-air coflow are conducted to fundamentally investigate turbulent jet ignition (TJI) in a three-dimensional configuration. TJI is an efficient method for initiating and controlling combustion in ultra-lean combustion systems. Fully compressible gas dynamics and species equations are solved with high order finite difference methods. The hydrogen-air reaction is simulated with a reliable detailed chemical kinetics mechanism. The physical processes involved in the TJI-assisted combustion are investigated by considering the flame heat release, temperature, species concentrations, vorticity, and Baroclinc torque. The complex turbulent flame and flow structures are delineated in three main: i) hot product jet, ii) burned-mixed, and iii) flame zones. In the TJI-assisted combustion, the flow structures and the flame features such as flame speed, temperature, and species distribution are found to be quite different than those in “standard” turbulent premixed combustion due to the existence of a high energy turbulent hot product jet. The flow structures and statistics are also found to be different than those normally seen in non-isothermal non-reacting jets.     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

Simultaneous OH and Schlieren visualization of premixed flames at the lean blow-out limit

The structure of an air-propane premixed flame was studied experimentally at the lean flammability limit, using Schlieren photography synchronized with OH-imaging done with the Planar Laser Induced Fluorescence (PLIF) technique. The flame was studied in a wide range of fuel equivalence ratios. Various steps in the process of the flame destabilization were investigated, including partial lift-off, stable lift-off, and final blow-out conditions. The flame structure was visualized for each stage showing the transition from a flame held at the nozzle to a flame held by the flow structures. In order to study the latter conditions in more detail the flame was acoustically excited at the preferred mode frequency generating large, stable, coherent structures in the core region. The modified flame structure was visualized to understand the interaction between the flame and vortical flow dynamics.It is shown that for the flow conditions when the flame cannot be stabilized at the nozzle, a new anchoring point is reached at the location of the initial vortex roll-up in the jet shear layer. At this point the flow reversal and transition to turbulence produce stagnation points with relatively low local velocities and velocity gradients where the flame can be stabilized. When the flame jet is being forced at the jet most unstable frequency, large coherent structures are formed and the flame is stabilized intermittently on these vortices.     

Large Eddy Simulations of Unconfined Non-reacting and Reacting Turbulent Low Swirl Jets

The low swirl flow is a novel method for stabilizing lean premixed combustion to achieve low emissions of nitrogen oxides. Understanding the characteristics of low swirl flows is of both practical and fundamental interest. In this paper, in order to gain better insight into low swirl stabilized combustion, large eddy simulation and dynamically thickened flame combustion modeling are used to characterize various features of non-reacting and reacting low swirl flows including vortex breakdown, shear layers’ instability, and coherent structures. Furthermore, four test cases with different equivalence ratios are studied to evaluate the effects of equivalence ratio on the flame and flow characteristics. A finite volume scheme on a Cartesian grid with a dynamic one equation eddy viscosity subgrid model is used for large eddy simulations. The obtained results show that the combustion heat release and increase in equivalence ratio toward the stoichiometric value decrease the local swirl number of the flow field, while increasing the flow spreading at the burner outlet. Results show that the flame becomes W shaped as the equivalence ratio increases. Moreover, the combination of the swirling motion and combustion heat release temporally imposes a vortex breakdown in the post-flame region, which leads to occurrence of a transient recirculation zone. The temporal recirculation zone disappears downstream of the burner outlet due to merging of the inner shear layer from all sides at the centerline. Also, various analyses of shear layers’ wavy and vortical structures show that combustion heat release has the effect of decreasing the instability amplitude and vortex shedding frequency.     

Instability of buoyant diffusion flames

Buoyant jet diffusion flames are known to exhibit large scale vortical flow structures strongly interacting with flame structures. In the present work, the formation and evolution of coherent flow structures is studied in a methane/ air coflow arrangement. This is accomplished by utilizing visualization techniques (planar laser induced hydroxyl fluorescence and Mie-scattering) and Laser Doppler Velocimetry. A striking repeatability and correlation of evolving coherent structures of the air co-flow and the reaction zone is observed. In the transitional region, flow and flame structures oscillate at very pure frequencies ranging from 10–15 Hz. A local absolutely unstable velocity profile close to the burner rim seems to be responsible. Self-excited axisymmetric wavelike structures propagate up- and downstream of this location. We study the influence of the exit velocities and the type of coflowing oxidizer (air or oxygen) on the location of transition to periodic flow structures and related frequencies. Conditional averages of image and velocity data are employed to describe the evolution of coherent flow structures and their interaction with flame structures.The authors wish to thank the Deutsche Forschungsgemeinschaft for financial support under contract Kn 118/22-2.     

Noncircular jets in combustion systems

Combustion dynamics of burners with corners were studied using Planar Laser Induced Fluorescence (PLIF) imaging. The effect of sharp corners on the air flow dynamics, shown earlier in cold flow tests, was also found in the reacting flow of a flame. The sharp corners interrupted the coherent structures generated in an axisymmetric shear flow. The combustion at the flat sections of the flame occurred in periodic, coherent large scale structures but was continuous and homogeneous in the vertices sections. The azimuthal structure of the noncircular flame changed in a pattern similar to that found in nonreacting flows. Combined regions of small- and large-scale mixing in the same flow, a unique feature of burners having sharp corners, is beneficial for combustion applications.     

A study on coherent structures and drag-reduction in the wall turbulence with polymer additives by TRPIV

An experimental measurement was performed using time-resolved particle image velocimetry (TRPIV) to investigate the spatial topological character of coherent structures in wall-bounded turbulence of polymer additive solution. The fully developed near-wall turbulent flow fields with and without polymer additives at the same Reynolds number were measured by TRPIV in a water channel. The comparisons of turbulent statistics confirm that due to viscoelastic structure of long-chain polymers, the wall-normal velocity fluctuation and Reynolds shear stress in the near-wall region are suppressed significantly. Furthermore, it is noted that such a behavior of polymers is closely related to the decease of the motion of the second and forth quadrants, i.e., the ejection and sweep events, in the near-wall region. The spatial topological mode of coherent structures during bursts has been extracted by the new mu-level criteria based on locally averaged velocity structure function. Although the general shapes of coherent structures are unchanged by polymer additives, the fluctuating velocity, velocity gradient, velocity strain rate and vorticity of coherent structures during burst events are suppressed in the polymer additive solution compared with that in water. The results show that due to the polymer additives the occurrence and intensity of coherent structures are suppressed, leading to drag reduction.     

Dynamics of interaction modes in excited annular jets

Evolution of coherent structures and their interaction dynamics are educed in the near field of an acoustically excited basic annular jet using conditional sampling technique based on a multiple triggering criterion to detect the two dominating modes of structure pattern. Acoustic excitation is applied with an aim to better organize the phase alignment of initial rolling and pairing process in the outer shear layer. Negligible modification of the time-averaged flow field results from the excitation. The educed coherent vorticities show that the two modes of evolution are due to the corresponding two modes of shedding pattern of the wake structures from the centerbody, namely the mode one wake and the mode zero wake. In both modes, the shear-layer mode jet vortex rings in the outer layer are perturbed by the shedding of wake structures in the inner region and interaction involving primary merging of three successive jet vortex rings or their partial circumferential sections is found. This results in the formation of wake-induced structures of the corresponding mode pattern, which possesses concentration of coherent vorticity and fluid circulation over a large spatial extent at 1 < x/D < 2. Secondary interactions, such as vortex tearing, are also observed.     

Coherent vortical and straining structures in the finite wall-mounted square cylinder wake

Topological aspects of the turbulent wake of a finite, surface-mounted, square-cross-section cylinder of h/d = 4 are addressed by decomposing the velocity field into a quasi-periodic coherent part and the unresolved incoherent fluctuations. The three-dimensional large scale structure is educed through a reconstruction of planar phase-averaged PIV measurements using the simultaneously sampled surface pressure difference on opposing sides of the obstacle as a phase reference. A topological model for the vortex structure is educed and mean streamwise wake vorticity is explained in terms of the connections between initially vertical structures shed alternately from either side of the obstacle, rather than previously proposed ‘tip’ vortex structures generated at the obstacle free-end. The coherent structure educed accounts for a significant portion of the fluctuating energy in the wake. The turbulent field is further analyzed by finding Lagrangian straining structures that form by induction of the coherent vorticity field, and these structures are related to the energy transfer from the base phase-averaged flow since they act to stretch incoherent vorticity fluctuations in their neighbourhood.     

湍流的相干结构:推演方法和结果

本文讨论壁面湍流发展的相干结构的观点.在简要的历史文献综述后,我们回顾一些基本观点, 并且介绍相干结构的思想.基于大量主要是由实验所得的结果,本文通过广泛运用的事件检 测技术,探讨湍流边界层内部和外部区域发生的现象.我们从边界层内部区域发生的现象、 边界层外部大尺度运动的发展和涡结构动力学的角度来描述流动的现象.在文章的 第2部分,介绍从背景流动中推演出湍流相干结构的各种方法以及在各种方法框架下所得到的结果, 讨论速度梯度张量不变量、压力的Hessian矩阵分析和本征正交分解等方法.每一个过程 都有``相干结构'的特定的定义,满足恰当的数学构架,并可以对湍流数据做相干结构动力 学分析.这一工作可能会对当前流体动力学家在湍流研究中用到的最新理论和技术的传播有 所贡献.     

Coherent structures and riblets

This work deals with the effect of the riblets on the coherent structures near the wall. The emphasis is put on the genesis of the quasi-streamwise vortices in the presence of the riblets. The quasi-streamwise vortices regenerate by the tilting of wall normal vorticity induced by prevailing structures. This requires a mechanism which leads to a temporal streamwise dependence near the elongated flow structures and to a subsequent formation of new wall normal vorticity. It is suggested here that the action of existing quasi-streamwise vortices on the sidewalls of wall normal vorticity may create a local, streamwise dependent spanwise velocity and therefore, a secondary wall normal vorticity field. A preliminary analysis of the set-up and the time and space development of this secondary three-dimensional flow associated with the regeneration mechanism, is given. An attempt is made, in order to explain the drag reduction performed by the riblets through an intermittent model, based on the protrusion height. Logical estimates of the amount of drag reduction are obtained. The differences between the mechanism suggested here and those based on forced control experiments are also discussed.     

Effect on flow field characteristics in methane–air non-premixed flame with steam addition

The present paper reports on an experimental study to determine the effect of humidity on the flow field and the flame stability limit in turbulent non-premixed flame, and examines the dynamical behavior of the unsteady aerodynamic flow structures observed on a bluff-body burner at both humid and non-humid air combustion states. Particle image velocimetry (PIV) is used to capture the instantaneous appearance of vortex structures and obtain the quantitative velocity field. Streamlines and velocity contours analysis are used to identify specific flame structures and reveal the effect of steam added on the vortex structure. The results show both central fuel penetration limit and partially quenching limit in the humid air case reduce. The decrease in the critical penetration limit is primarily attributed to a reduction in momentum of the humid air. The flamelet concepts are applied to discuss the partially quenching limit in the blue neck region. The analysis reveals that the large decrease in the partially quenching limit is due to the increase in chemical reaction time of the humid air combustion.     

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.     

Role of Buoyancy on Instabilities and Structure of Transitional Gas Jet Diffusion Flames

Transitional jet diffusion flames provide the link between dynamics of laminar and turbulent flames. In this study, instabilities and their interaction with the flow structure are explored in a transitional jet diffusion flame, with focus on isolating buoyancy effects. Experiments are conducted in hydrogen flames with fuel jet Reynolds number of up to 2,200 and average jet velocity of up to 54 m/s. Since the fuel jet is laminar at the injector exit, the transition from laminar to turbulent flame occurs by the hydrodynamic instabilities in the shear layer of fuel jet. The instabilities and the flow structures are visualized and quantified by the rainbow schlieren deflectometry technique coupled with a high-speed imaging system. The schlieren images acquired at 2,000 frames per second allowed exposure time of 23 μs with spatial resolution of 0.4 mm. Results identify a hitherto unknown secondary instability in the flame surface, provide explanation for the observed intermittency in the breakpoint length, show coherent vortical structures downstream of the flame breakpoint, and illustrate gradual breakdown of coherent structures into small-scale random structures in the far field turbulent region.     

吹吸扰动对壁湍流相干结构的影响

采用热线风速仪测量受吹吸扰动的壁湍流边界层的流向速度,用傅里叶变换和子波变换研究 吹吸扰动对壁湍流能谱的影响,结果显示施加的低频扰动使边界层内层大尺度结构的能量减 少,小尺度结构的能量有所增强,远离壁面时扰动强度逐步衰减直到在外层中消失;通过VITA 法和子波变换法检测猝发事件,表明该扰动降低了猝发强度,使猝发周期延长,条件平均速 度波形的幅值降低、持续时间变短,说明扰动明显抑制了相干结构的猝发过程. 利用子波变 换可以实现湍谱分析,能有效检测猝发中的湍流结构,是一种客观的分析工具.     

Flow structures in initial region of two interacting parallel plane jets

The initial region of two interacting parallel plane jets with a truncated bluff body at the plane of symmetry is studied with velocity measurements in the frequency and time domains. The mean flows of the two plane jets are nearly parallel to the plane of symmetry, and their inner and outer mixing regions are found. Within these mixing regions their respective trains of coherent structures, which are of vortical form, are established. The inner vortices are inwardly rotational, while the outer vortices rotate outwardly. The successive initial vortices in the inner mixing region seem to undergo two processes of either pairing or amalgamation. The former results in the formation of a nearly circular coherent structure, while the latter results in the elongated structure. Both the inner and other structures experience fairly rapid decay. The process of decay also seems to be either the usual decay or the division of the coalesced structure back into the individual vortical structures followed by their decay.     

Turbulence in a skewed three‐dimensional wall‐bounded shear flow: effect of mean vorticity on structure modification

Mean‐flow three‐dimensionalities affect both the turbulence level and the coherent flow structures in wall‐bounded shear flows. A tailor‐made flow configuration was designed to enable a thorough investigation of moderately and severely skewed channel flows. A unidirectional shear‐driven plane Couette flow was skewed by means of an imposed spanwise pressure gradient. Three different cases with 8°, 34°and 52°skewing were simulated numerically and the results compared with data from a purely two‐dimensional plane Couette flow. The resulting three‐dimensional flow field became statistically stationary and homogeneous in the streamwise and spanwise directions while the mean velocity vector V and the mean vorticity vector Ω remained parallel with the walls. Mean flow profiles were presented together with all components of the Reynolds stress tensor. The mean shear rate in the core region gradually increased with increasing skewing whereas the velocity fluctuations were enhanced in the spanwise direction and reduced in the streamwise direction. The Reynolds shear stress is known to be closely related to the coherent flow structures in the near‐wall region. The instantaneous and ensemble‐averaged flow structures were turned by the skewed mean flow. We demonstrated for the medium‐skewed case that the coherent structures should be examined in a coordinate system aligned with V to enable a sound interpretation of 3D effects. The conventional symmetry between Case 1 and Case 2 vortices was broken and Case 1 vortices turned out to be stronger than Case 2. This observation is in conflict with the common understanding on the basis of the spanwise (secondary) mean shear rate. A refined model was proposed to interpret the structure modifications in three‐dimensional wall‐flows. What matters is the orientation of the mean vorticity vector Ω relative to the vortex vorticity vector ω _v, that is, the sign of Ω · ω _v. In the present situation, Ω · ω _v > 0 for the Case 1 vortices causing a strengthening relative to the Case 2 vortices. Copyright © 2011 John Wiley & Sons, Ltd.