Absolute and convective instability in gas turbine fuel injectors

Hydrodynamic instabilities in gas turbine fuel injectors help to mix the fuel and air but can sometimes lock into acoustic oscillations and contribute to thermoacoustic instability. This paper describes a linear stability analysis that predicts the frequencies and strengths of hydrodynamic instabilities and identifies the regions of the flow that cause them. It distinguishes between convective instabilities, which grow in time but are convected away by the flow, and absolute instabilities, which grow in time without being convected away. Convectively unstable flows amplify external perturbations, while absolutely unstable flows also oscillate at intrinsic frequencies. As an input, this analysis requires velocity and density fields, either from a steady but unstable solution to the Navier-Stokes equations, or from time-averaged numerical simulations. In the former case, the analysis is a predictive tool. In the latter case, it is a diagnostic tool. This technique is applied to three flows: a swirling wake at Re = 400, a single stream swirling fuel injector at Re - 106, and a lean premixed gas turbine injector with five swirling streams at Re - 106. Its application to the swirling wake demonstrates that this technique can correctly predict the frequency, growth rate and dominant wavemaker region of the flow. It also shows that the zone of absolute instability found from the spatio-temporal analysis is a good approximation to the wavemaker region, which is found by overlapping the direct and adjoint global modes. This approximation is used in the other two flows because it is difficult to calculate their adjoint global modes. Its application to the single stream fuel injector demonstrates that it can identify the regions of the flow that are responsible for generating the hydrodynamic oscillations seen in LES and experimental data. The frequencies predicted by this technique are within a few percent of the measured frequencies. The technique also explains why these oscillations become weaker when a central jet is injected along the centreline. This is because the absolutely unstable region that causes the oscillations becomes convectively unstable. Its application to the lean premixed gas turbine injector reveals that several regions of the flow are hydrodynamically unstable, each with a different frequency and a different strength. For example, it reveals that the central region of confined swirling flow is strongly absolutely unstable and sets up a precessing vortex core, which is likely to aid mixing throughout the injector. It also reveals that the region between the second and third streams is slightly absolutely unstable at a frequency that is likely to coincide with acoustic modes within the combustion chamber. This technique, coupled with knowledge of the acoustic modes in a combustion chamber, is likely to be a useful design tool for the passive control of mixing and combustion instability. Copyright © 2012 by ASME.

On Mixing and Shaking : Structure and Dynamics of Turbulent Stratified Flames

The drive for low emission combustion systems encourages applications using premixed flames. Yet in many applications, considerations of flame stability or mixing times lead to systems with neither premixed nor diffusion flames, which are often called technically premixed or stratified flames. In this talk we discuss the current state of understanding of the effect of mixing and extent of stratification on the structure, microstructure and dynamics of selected turbulent stratified flames. Over the past few years, a significant database of scalar and velocity data has been built to analyze the effects of unmixedness on local and global flame structure. Microscale studies of the flame structures show in detail how the effect of local stratification affects (or not!) the flame structure, flame surface density and scalar dissipation rates, and production of selected species. The experiments place exacting demands on current spectroscopic diagnostics, and reveal the progress and limits to our understanding of turbulent flames in general. The dynamics of stratified flames with respect to instabilities is also shown to be very rich, as the particular shape of the flames and the stabilization points are is significantly affected by the fuel distribution, modifying the rate and location of heat release, and thus the coupling with the surrounding acoustics and determining the onset of self-excitations.

Instability of the steady-state regime of generation in laser diodes with a broad gain spectrum

The laser-diode parameters at which the steady-state regime of generation becomes unstable are analyzed within the framework of the mode-locking model. The crucial role of the transverse inhomogeneity of the field, pumping intensity, and spectrum width in developing the instabilities of the steady-state regime of generation is demonstrated. The calculated values of the instability threshold are shown to be consistent with the experimental results. © 2008 Springer Science+Business Media, Inc.

CFD investigation of turbulent premixed flame response to transverse forcing

Screech is a high frequency oscillation that is usually characterized by instabilities caused by large-scale coherent flow structures in the wake of bluff-body flameholders and shear layers. Such oscillations can lead to changes in flame surface area which can cause the flame to burn unsteadily, but also couple with the acoustic modes and inherent fluid-mechanical instabilities that are present in the system. In this study, the flame response to hydrodynamic oscillations is analyzed in a controlled manner using high-fidelity Computational Fluid Dynamics (CFD) with an unsteady Reynolds-averaged Navier-Stokes approach. The response of a premixed flame with and without transverse velocity forcing is analyzed. When unforced, the flame is shown to exhibit a self-excitation that is attributed to the anti-symmetric shedding of vortices in the wake of the flameholder. The flame is also forced using two different kinds of low-amplitude out-of-phase inlet velocity forcing signals. The first forcing method is harmonic forcing with a single characteristic frequency, while the second forcing method involves a broadband forcing signal with frequencies in the range of 500 - 1000 Hz. For the harmonic forcing method, the flame is perturbed only lightly about its mean position and exhibits a limit cycle oscillation that is characteristic of the forcing frequency. For the broadband forcing method, larger changes in the flame surface area and detachment of the flame sheet can be seen. Transition to a complicated trajectory in the phase space is observed. When analyzed systematically with system identification methods, the CFD results, expressed in the form of the Flame Transfer Function (FTF) are capable of elucidating the flame response to the imposed perturbation. The FTF also serves to identify, both spatially and temporally, regions where the flame responds linearly and nonlinearly. Locking-in between the flame's natural self-excited frequency and the subharmonic frequencies of the broadband forcing signal is found to alter the dynamical behaviour of the flame. Copyright © 2013 by ASME.

Global linear stability of the boundary-layer flow over a rotating sphere

We consider the linear global stability of the boundary-layer flow over a rotating sphere. Our results suggest that a self-excited linear global mode can exist when the sphere rotates sufficiently fast, with properties fixed by the flow at latitudes between approximately 55°-65° from the pole (depending on the rotation rate). A neutral curve for global linear instabilities is presented with critical Reynolds number consistent with existing experimentally measured values for the appearance of turbulence. The existence of an unstable linear global mode is in contrast to the literature on the rotating disk, where it is expected that nonlinearity is required to prompt the transition to turbulence. Despite both being susceptible to local absolute instabilities, we conclude that the transition mechanism for the rotating-sphere flow may be different to that for the rotating disk. © 2014 Elsevier Masson SAS. All rights reserved.

Comparison of acoustic velocity perturbation measurements using PIV vs. two-microphone technique

Understanding combustion instabilities requires accurate measurements of the acoustic velocity perturbation into injectors. This is often accomplished via the use of the two microphone technique, as this only requires two pressure transducers. However, measurements of the actual velocities emerging from the injectors are not often taken, leaving questions regarding the assumptions about the acoustic velocity. A comparison of velocity measured at downstream of the injector with that of two-microphone technique can show the accuracy and limitations of two-microphone technique. In this paper, velocity measurements are taken using both particle image velocimetry (PIV) and the two-microphone technique in a high pressure facility designed for aeroengine injector measurements. The flow is excited using an area modulation device installed on the choked end of the combustion chamber, with PIV measurements enabled by optical access downstream of the injector through a quartz tube and windows. Acoustic velocity perturbations at the injector are determined by considering the Fourier transformed pressure fluctuations for two microphones installed at a known distance upstream of the injector. PIV measurements are realized by seeding the air flow with micrometric water particles under 2.5 bar pressure at ambient temperature. Phase locked velocity fields are realized by synchronizing the acquisition of PIV images with the revolution of the acoustic modulator using the pressure signal measured at the face of injector. The mean velocity fluctuation is calculated as the difference between maximum and minimum velocities, normalized by the mean velocity of the unforced case. Those values are compared with the peak-to-peak velocity fluctuation amplitude calculated by the two-microphone technique. Although the ranges of velocity fluctuations for both techniques are similar, the variation of fluctuation with forcing frequencies diverges significantly with frequency. The differences can be attributed to several limitations associated with of both techniques, such as the quality of the signal, the signal/noise ratio, the accuracy of PIV measurements and the assumption of isentropic flow of the particle velocity from the plenum through the injector. We conclude that two-microphone methods can be used as a reference value for the velocity fluctuation in low order applications such as flame transfer functions, but not for drawing conclusions regarding the absolute velocity fluctuations in the injector. Copyright © 2013 by ASME.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.

Static and dynamic electric field domain formation in a doped GaAs/AlAs superlattice

We have observed the transition from static to dynamic electric field domain formation induced by a transverse magnetic field and the sample temperature in a doped GaAs/AlAs superlattice. The observations can be very well explained by a general analysis of instabilities and oscillations of the sequential tunnelling current in superlattices based solely on the magnitude of the negative differential resistance region in the tunnelling characteristic of a single barrier. Both increasing magnetic field and sample temperature change the negative differential resistance and cause the transition between static and dynamic electric field domain formation. (C) 2000 Elsevier Science B.V. All rights reserved.

Anomaly of the current self-oscillation frequency in the sequential tunneling of a doped GaAs/AlAs superlattice

An anomalous behavior of the current self-oscillation frequency is observed in the dynamic de voltage bands, emerging from each sawtoothlike branch of the current-voltage characteristic of a doped GaAs/A1As superlattice in the transition process from static to dynamic electric field domain formations. Varying the applied de voltage at a fixed temperature, we find that the frequency increases while the averaged current decreases. Inside each voltage band, the frequency has a strong voltage dependence in the temperature range where the averaged current changes with the applied de voltage. This dependence can be understood in terms of motion of the system along a limit cycle.

Negative differential resistance and the transition to current self-oscillation in GaAs/AlAs superlattices

We investigate the transition from static to dynamic electric field domains (EFDs) in a doped GaAs/AlAs superlattice (SL). We show that a transverse magnetic field and/or the temperature can induce current self-oscillations. This observation can be attributed to the negative differential resistance (NDR) effect. Transverse magnetic field and the temperature can increase the NDR of a doped SL. A large NDR can lead to an unstable EFD in a certain range of d.c. bias. (C) 1999 Elsevier Science Ltd. All rights reserved.

Effect of channel surface wettability and temperature gradients on the boiling flow pattern in a single microchannel

The objective of this paper is to investigate the effects of channel surface wettability and temperature gradients on the boiling flow pattern in a single microchannel. The test section consists of a bottom silicon substrate bonded with a top glass cover. Three consecutive parts of an inlet fluid plenum, a central microchannel and an outlet fluid plenum were etched in the silicon substrate. The central microchannel had a width of 800 mu m and a depth of 30 mu m. Acetone liquid was used as the working fluid. High outlet vapor qualities were dealt with here. The flow pattern consists of a fluid triangle (shrinkage of the liquid films) and a connected long liquid rivulet, which is generated in the central microchannel in the timescale of milliseconds. The peculiar flow pattern is formed due to the following reasons: (1) the liquid rivulet tends to have a large contact area with the top hydrophilic channel surface of the glass cover, but a smaller contact area with the bottom silicon hydrophobic surface. (2) The temperature gradient in the chip width direction at the top channel surface of the glass cover not only causes the shrinkage of the liquid films in the central microchannel upstream, but also attracts the liquid rivulet populated near the microchannel centerline. (3) The zigzag pattern is formed due to the competition between the evaporation momentum forces at the vapor-liquid interfaces and the force due to the Marangoni effect. The former causes the rivulet to deviate from the channel centerline and the latter draws the rivulet toward the channel centerline. (4) The temperature gradient along the flow direction in the central microchannel downstream causes the breakup of the rivulet to form isolated droplets there. (5) Liquid stripes inside the upstream fluid triangle were caused by the small capillary number of the liquid film, at which the large surface tension force relative to the viscous force tends to populate the liquid film locally on the top glass cover surface.

Static and dynamic electric field domain formation in a doped GaAs/AlAs superlattice

We have observed the transition from static to dynamic electric field domain formation induced by a transverse magnetic field and the sample temperature in a doped GaAs/AlAs superlattice. The observations can be very well explained by a general analysis of instabilities and oscillations of the sequential tunnelling current in superlattices based solely on the magnitude of the negative differential resistance region in the tunnelling characteristic of a single barrier. Both increasing magnetic field and sample temperature change the negative differential resistance and cause the transition between static and dynamic electric field domain formation. (C) 2000 Elsevier Science B.V. All rights reserved.

Linear spatio-temporal instability analysis of ice growth under a falling water film

A linear spatio-temporal stability analysis is conducted for the ice growth under a falling water film along an inclined ice plane. The full system of linear stability equations is solved by using the Chebyshev collocation method. By plotting the boundary curve between the linear absolute and convective instabilities (AI/CI) of the ice mode in the parameter plane of the Reynolds number and incline angle, it is found that the linear absolute instability exists and occurs above a minimum Reynolds number and below a maximum inclined angle. Furthermore, by plotting the critical Reynolds number curves with respect to the inclined angle for the downstream and upstream branches, the convectively unstable region is determined and divided into three parts, one of which has both downstream and upstream convectively unstable wavepackets and the other two have only downstream or upstream convectively unstable wavepacket. Finally, the effect of the Stefan number and the thickness of the ice layer on the AI/CI boundary curve is investigated.

HIRFL-CSR频谱现象：中心频率突变分析研究

In commissioning the HIRFL-CSR (Cooler Storage Ring of Heavy Ion Research Facility in Lanzhou),with the spectrum analyzer monitoring the beams on-line, we measured many parameters of the beam status for improving quality of the heavy ion beams. In this paper, the sudden change of the beam center frequency is analyzed.A theoretical analysis is done mainly on instabilities of the electron voltage and the neutralization factor. Simulations were made for the problem. According to the simulation result, we can conclude that it is the instability of neutralization factor that caused the problem. 中文摘要：HIRFL-CSR调试过程中，利用频谱分析仪实时监测束流，根据其上测的数据进行相关参数调整，改善束流的品质。本文对分析仪上观测到的现象一束流的中心频率发生突变进行了研究，主要就电子束的高压不稳定性和中性化因子的突变做了理论分析，并进行了相关的模拟计算 ，根据模拟结果初步断定该现象的出现由电子束的中性化因子导致。

Investigation of the dewetting inhibition mechanism of thin polymer films

The mechanism of inhibition of polymer film dewetting is investigated by adding a star comb-like polymer, four-arm P(S-ran-VB-g-PMMA), to PS film and PMMA film on different substrates. It is found that the mechanism of inhibition of polymer film dewetting is kinetic in nature, and is related to the miscibility between the additional compound and the polymer film. On addition to the miscible system [four-arm P(S-ran-VB-g-PMMA) and PMMA], the star comb-like polymers can increase the resistant force of dewetting with hole growth and inhibit the dewetting process of the thin polymer film by enrichment in the rim.