Study of spike and modal stall phenomena in a low-speed axial compressor

This paper presents a study of stall inception mechanisms a in low-speed axial compressor. Previous work has identified two common flow breakdown sequences, the first associated with a short lengthscale disturbance known as a `spike', and the second with a longer lengthscale disturbance known as a `modal oscillation'. In this paper the physical differences between these two mechanisms are illustrated with detailed measurements. Experimental results are also presented which relate the occurrence of the two stalling mechanisms to the operating conditions of the compressor. It is shown that the stability criteria for the two disturbances are different: long lengthscale disturbances are related to a two-dimensional instability of the whole compression system, while short lengthscale disturbances indicate a three-dimensional breakdown of the flow-field associated with high rotor incidence angles. Based on the experimental measurements, a simple model is proposed which explains the type of stall inception pattern observed in a particular compressor. Measurements from a single stage low-speed compressor and from a multistage high-speed compressor are presented in support of the model.

Measurements of triggering and transient growth in a model lean-premixed gas turbine combustor

Instability triggering and transient growth of thermoacoustic oscillations were experimentally investigated in combination with linear/nonlinear flame transfer function (FTF) methodology in a model lean-premixed gas turbine combustor operated with CH 4 and air at atmospheric pressure. A fully premixed flame with 10kW thermal power and an equivalence ratio of 0.60 was chosen for detailed characterization of the nonlinear transient behaviors. Flame transfer functions were experimentally determined by simultaneous measurements of inlet velocity fluctuations and heat release rate oscillations using a constant temperature anemometer and OH */CH * chemiluminescence emissions, respectively. The phase-resolved variation of the local flame structure at a limit cycle was measured by planar laser-induced fluorescence of OH. Simultaneous measurements of inlet velocity, OH */CH * emission, and acoustic pressure were performed to investigate the temporal evolution of the system from a stable to a limit cycle operation. This measurement allows us to describe an unsteady instability triggering event in terms of several distinct stages: (i) initiation of a small perturbation, (ii) exponential amplification, (iii) saturation, (iv) nonlinear evolution of the perturbations towards a new unstable periodic state, (v) quasi-steady low-amplitude periodic oscillation, and (vi) fully-developed high-amplitude limit cycle oscillation. Phase-plane portraits of instantaneous inlet velocity and heat release rate clearly show the presence of two different attractors. Depending on its initial position in phase space at infinitesimally small amplitude, the system evolves towards either a high-amplitude oscillatory state or a low-amplitude oscillatory state. This transient phenomenon was analyzed using frequency- and amplitude-dependent damping mechanisms, and compared to subcritical and supercritical bifurcation theories. The results presented in this paper experimentally demonstrate the hypothesis proposed by Preetham et al. based on analytical and computational solutions of the nonlinear G-equation [J. Propul. Power 24 (2008) 1390-1402]. Good quantitative agreement was obtained between measurements and predictions in terms of the conditions for the onset of triggering and the amplitude of triggered combustion instabilities. © 2011 The Combustion Institute.

SIMPLE REMESHING SCHEME FOR THE FINITE ELEMENT BASED SIMULATION OF DOPANT DIFFUSION IN SILICON.

Process simulation programs are valuable in generating accurate impurity profiles. Apart from accuracy the programs should also be efficient so as not to consume vast computer memory. This is especially true for devices and circuits of VLSI complexity. In this paper a remeshing scheme to make the finite element based solution of the non-linear diffusion equation more efficient is proposed. A remeshing scheme based on comparing the concentration values of adjacent node was then implemented and found to remove the problems of oscillation.

Effects of nonuniform reactant stoichiometry on thermoacoustic instability in a lean-premixed gas turbine combustor

Detailed experimental investigations of the amplitude dependence of flame describing functions (FDF) were performed using a stratified swirl-stabilized combustor, in order to understand the combustion-acoustic interactions of CH4/air flames propagating into nonhomogeneous reactant stoichiometry. Phase-synchronized OH planar laser induced fluorescence (OH PLIF) measurements were used to investigate local reaction zone structures of forced flames. To determine the amplitude-and frequency-dependent forced flame response, simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. The measurements were made over a wide range of stratification ratios, including inner stream enrichment ( θ o>θ i) and outer stream enrichment ( θ o>θ i)) conditions, and compared to the baseline condition of spatially and temporally homogeneous cases ( θ o=θ i)). Results show that for the inlet conditions investigated, fuel stratification has a significant influence on local and global flame structures of unforced and forced flames. Under stratified conditions, length scales of local contours were found to be much larger than the homogeneous case due to high kinematic viscosities associated with high temperature. Stratification has a remarkable effect on flame-vortex interactions when the flame is subjected to high-amplitude acoustic forcing, leading to different evolution patterns of FDF (amplitude and disturbance convective time) in response to the amplitude of the imposed inlet velocity oscillation. The present experimental investigation reveals that intentional stratification has the potential to eliminate or suppress the occurrence of detrimental combustion instability problems in lean-premixed gas turbine combustion systems. © 2012 Copyright Taylor and Francis Group, LLC.

Surface tension-induced global instability of planar jets and wakes

The effect of surface tension on global stability of co-flow jets and wakes at a moderate Reynolds number is studied. The linear temporal two-dimensional global modes are computed without approximations. All but one of the flow cases under study are globally stable without surface tension. It is found that surface tension can cause the flow to be globally unstable if the inlet shear (or equivalently, the inlet velocity ratio) is strong enough. For even stronger surface tension, the flow is re-stabilized. As long as there is no change of the most unstable mode, increasing surface tension decreases the oscillation frequency. Short waves appear in the high-shear region close to the nozzle, and their wavelength increases with increasing surface tension. The critical shear (the weakest inlet shear at which a global instability is found) gives rise to antisymmetric disturbances for the wakes and symmetric disturbances for the jets. However, at stronger shear, the opposite symmetry can be the most unstable one, in particular for wakes at high surface tension. The results show strong effects of surface tension that should be possible to reproduce experimentally as well as numerically.

Combustion rumble prediction with integrated computational-fluid-dynamics/ low-order-model methods

The pressure oscillation within combustion chambers of aeroengines and industrial gas turbines is a major technical challenge to the development of high-performance and low-emission propulsion systems. In this paper, an approach integrating computational fluid dynamics and one-dimensional linear stability analysis is developed to predict the modes of oscillation in a combustor and their frequencies and growth rates. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone, which enables the computational fluid dynamics calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in agreement with experimental measurements. Furthermore, results from the computational fluid dynamics calculation provide the flame transfer function to describe unsteady heat release rate. Departures from ideal one-dimensional flows are described by shape factors. Combined with this information, low-order models can work out the possible oscillation modes and their initial growth rates. The approach developed here can be used in more general situations for the analysis of combustion oscillations. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Internal representations of temporal statistics and feedback calibrate motor-sensory interval timing.

Humans have been shown to adapt to the temporal statistics of timing tasks so as to optimize the accuracy of their responses, in agreement with the predictions of Bayesian integration. This suggests that they build an internal representation of both the experimentally imposed distribution of time intervals (the prior) and of the error (the loss function). The responses of a Bayesian ideal observer depend crucially on these internal representations, which have only been previously studied for simple distributions. To study the nature of these representations we asked subjects to reproduce time intervals drawn from underlying temporal distributions of varying complexity, from uniform to highly skewed or bimodal while also varying the error mapping that determined the performance feedback. Interval reproduction times were affected by both the distribution and feedback, in good agreement with a performance-optimizing Bayesian observer and actor model. Bayesian model comparison highlighted that subjects were integrating the provided feedback and represented the experimental distribution with a smoothed approximation. A nonparametric reconstruction of the subjective priors from the data shows that they are generally in agreement with the true distributions up to third-order moments, but with systematically heavier tails. In particular, higher-order statistical features (kurtosis, multimodality) seem much harder to acquire. Our findings suggest that humans have only minor constraints on learning lower-order statistical properties of unimodal (including peaked and skewed) distributions of time intervals under the guidance of corrective feedback, and that their behavior is well explained by Bayesian decision theory.

Response of two acoustically excited turbulent premixed flames to an imposed phase lag

This paper describes an experimental investigation into the interactions that occur between two acoustically forced lean turbulent premixed flames for an induced phase lag. Phase-averaged FSD from cinematographic OH-PLIF measurements and global heat release measurements were obtained for a range phase lags (ψ S) and amplitudes (A) as a function of flame separation distance, S. The effect of bringing two flames closer together causes jet merging, which alters the vortex flame interactions that drive the thermo-acoustic response. To simulate circumferential modes a phase lag was introduced, which affected the flame dynamics in the region of flame-flame interaction, with these changes dependent on S. For moderate separation distances, the flame structure becomes increasingly asymmetric inducing a very small transverse oscillation. However, for moderate phase lags (φ s ≤ 20) the magnitude of these changes and their subsequent influence on the thermo-acoustic response was found to be slight in comparison with changes in S. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

Studying adsorbent dynamics on a quartz crystal resonator using its nonlinear electrical response

The quartz crystal resonator has been traditionally employed in studying surface-confined physisorbed films and particles by measuring dissipation and frequency shifts. However, theoretical interpretation of the experimental observations is often challenged due to limited understanding of physical interaction mechanisms at the interfaces involved. Here we model a physisorbed interaction between particles and gold electrode surface of a quartz crystal and demonstrate how the nonlinear modulation of the electric response of the crystal due to the nonlinear interaction forces may be used to study the dynamics of the particles. In particular, we show that the graphs of the deviation in the third Fourier harmonic response versus oscillation amplitude provide important information about the onset, progress and nature of sliding of the particles. The graphs also present a signature of the surface-particle interaction and could be used to estimate the interaction energy profile. Interestingly, the insights gained from the model help to explain some of the experimental observations with physisorbed streptavidin-coated polystyrene microbeads on quartz resonators. © 2012 Elsevier B.V. All rights reserved.

A novel measurement of dynamic friction

A tribometer, based on a pin-on-disc machine, uses a PZT drive to produce small sinusoidal fluctuations of sliding speed. The frequency and amplitude of these fluctuations can be controlled, and the dynamic response measured. Preliminary test results show that the dynamic friction variation is influenced by the contact materials, normal force, oscillation frequency and steady sliding speed. The variation of friction force amplitude and phase with frequency gives clues about the underlying state variables determining the friction. Modelling studies illustrate the expected behaviour for idealized friction laws governed by, for example, sliding speed, contact temperature, and "rate-state" laws. © 2008 SAE International.

Dissipativity characterization of a class of oscillators and networks of oscillators

This paper uses dissipativity theory to provide the system-theoretic description of a basic oscillation mechanism. Elementary input-output tools are then used to prove the existence and stability of limit cycles in these "oscillators". The main benefit of the proposed approach is that it is well suited for the analysis and design of interconnections, thus providing a valuable mathematical tool for the study of networks of coupled oscillators.

Enhanced infra-red emission from sub-millimeter microelectromechanical systems micro hotplates via inkjet deposited carbon nanoparticles and fullerenes

In this paper, we demonstrate a micro-inkjet printing technique as a reproducible post-process for the deposition of carbon nanoparticles and fullerene adlayers onto fully CMOS compatible micro-electro-mechanical silicon-on-insulator infrared (IR) light sources to enhance their infrared emission. We show experimentally a significant increase in the infrared emission efficiency of the coated emitters. We numerically validate these findings with models suggesting a dominant performance increase for wavelengths <5.5 μm. Here, the bimodal size distribution in the diameter of the carbon nanoparticles, relative to the fullerenes, is an effective mediator towards topologically enhanced emittance of our miniaturised emitters. A 90% improvement in IR emission power density has been shown which we have rationalised with an increase in the mean thickness of the deposited carbon nanoparticle adlayer. © 2013 AIP Publishing LLC.

Nonlinear phenomena in thermoacoustic systems with premixed flames

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.

Statistics of reaction progress variable and mixture fraction gradients from DNS of turbulent partially premixed flames

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios = 0.7 and 1.0 using three-dimensional, simplified chemistry-based direct numerical simulations (DNS). The simulation parameters are chosen such that the flames represent the thin reaction zones regime combustion. A random bimodal distribution of equivalence ratio is introduced in the unburned gas ahead of the flame to account for the mixture inhomogeneity. The results suggest that the probability density functions (PDFs) of the mixture fraction gradient magnitude |Δξ| (i.e., P(|Δξ|)) can be reasonably approximated using a log-normal distribution. However, this presumed PDF distribution captures only the qualitative nature of the PDF of the reaction progress variable gradient magnitude |Δc| (i.e., P(|Δc|)). It has been found that a bivariate log-normal distribution does not sufficiently capture the quantitative behavior of the joint PDF of |Δξ| and |Δc| (i.e., P(|Δξ|, |Δc|)), and the agreement with the DNS data has been found to be poor in certain regions of the flame brush, particularly toward the burned gas side of the flame brush. Moreover, the variables |Δξ| and |Δc| show appreciable correlation toward the burned gas side of the flame brush. These findings are corroborated further using a DNS data of a lifted jet flame to study the flame geometry dependence of these statistics. © 2013 Copyright Taylor and Francis Group, LLC.

Nonlinear thermoacoustics of ducted premixed flames: The influence of perturbation convection speed

When a premixed flame is placed within a duct, acoustic waves induce velocity perturbations at the flame's base. These travel down the flame, distorting its surface and modulating its heat release. This can induce self-sustained thermoacoustic oscillations. Although the phase speed of these perturbations is often assumed to equal the mean flow speed, experiments conducted in other studies and Direct Numerical Simulation (DNS) conducted in this study show that it varies with the acoustic frequency. In this paper, we examine how these variations affect the nonlinear thermoacoustic behaviour. We model the heat release with a nonlinear kinematic G-equation, in which the velocity perturbation is modelled on DNS results. The acoustics are governed by linearised momentum and energy equations. We calculate the flame describing function (FDF) using harmonic forcing at several frequencies and amplitudes. Then we calculate thermoacoustic limit cycles and explain their existence and stability by examining the amplitude-dependence of the gain and phase of the FDF. We find that, when the phase speed equals the mean flow speed, the system has only one stable state. When the phase speed does not equal the mean flow speed, however, the system supports multiple limit cycles because the phase of the FDF changes significantly with oscillation amplitude. This shows that the phase speed of velocity perturbations has a strong influence on the nonlinear thermoacoustic behaviour of ducted premixed flames. © 2013 The Combustion Institute.