Effects of nonuniform reactant stoichiometry on combustion instability

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to acoustic forcing in a laboratory-scale stratified burner. The double-swirler, double-channel annular burner was specially designed to generate acoustic velocity oscillations and radial fuel stratification at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the flame response are not considered. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a hot wire anemometer and photomultiplier tubes with narrowband OH*/CH* interference filters. Time-averaged CH* chemiluminescence intensities were measured using an intensified CCD camera. Results show that flame stabilization mechanisms vary depending on stratification ratio for a constant global equivalence ratio. For a uniformly premixed condition, an enveloped M-shaped flame is observed. For stratified conditions, however, a dihedral V-flame and a detached flame are developed for outer stream and inner stream fuel enrichment cases, respectively. Flame transfer function (FTF) measurement results indicate that a V-shaped flame tends to damp incident flow oscillations, while a detached flame acts as a strong amplifier relative to the uniformly premixed condition. The phase difference of FTF increases in the presence of stratification. More importantly, the dynamic characteristics obtained from the forced stratified flame measurements are well correlated with unsteady flame behavior under limit-cycle pressure oscillations. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, which has not been well explored to date. Copyright © 2011 by ASME.

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.

Calibration of the spring constant of cantilevers of arbitrary shape using the phase signal in an atomic force microscope.

The measurement of cantilever parameters is an essential part of performing a calibrated measurement with an atomic force microscope (AFM). The thermal motion method is a widely used technique for calibrating the spring constant of an AFM cantilever, which can be applied to non-rectangular cantilevers. Given the trend towards high frequency scanning, calibration of non-rectangular cantilevers is of increasing importance. This paper presents two results relevant to cantilever calibration via the thermal motion method. We demonstrate the possibility of using the AFM's phase signal to acquire the thermal motion. This avoids the challenges associated with connecting the raw photodiode signal to a separate spectrum analyser. We also describe how numerical calculations may be used to calculate the parameters needed in a thermal motion calibration of a non-rectangular cantilever. Only accurate knowledge of the relative size of the in-plane dimensions of the cantilever is needed in this computation. We use this pair of results in the calibration of a variety of rectangular and non-rectangular cantilevers. We observe an average difference between the Sader and thermal motion values of cantilever stiffness of 10%.

Finite element modelling of the sismic behaviour of water front structures

Water front structures have suffered significant damage in many of the recent earthquakes. These include gravity type quay walls, vertically composite walls, cantilever retaining walls, anchored bulkheads and similar structures. One of the primary causes for the poor performance of these classes of structures is the liquefaction of the foundation soil and in some instances liquefaction of the backfill soil. The liquefaction of the soil in-front of the quay wall tends to cause large lateral displacements and rotation of the wall. Often such gravity walls are placed on rubble mound deposited onto the sea bed.This paper presents finite element analyses of such a problem in which strength degradation of the foundation soil and the backfill material will be modelled using PZ mark III constitutive model. The performance of the wall in terms of its lateral displacement, vertical settlement and/or the rotation suffered by the wall will be presented. In addition, the contours of the horizontal and vertical effective stresses and the excess pore pressure ratio will be presented at different time instants together with hyrdraulic gradients. Immediately after the earthquake, the hydraulic gradients indicate migration of pore water into the region below the wall, suggesting further softening of the foundation soil below the wall.

Sensitivity analysis of a time-delayed thermo-acoustic system via an adjoint-based approach

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat-release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot-wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates. © Cambridge University Press 2013.

Parallels between wind and crowd loading of bridges

Parallels between the dynamic response of flexible bridges under the action of wind and under the forces induced by crowds allow each field to inform the other.Wind-induced behaviour has been traditionally classified into categories such as flutter, galloping, vortex-induced vibration and buffeting. However, computational advances such as the vortex particle method have led to a more general picture where effects may occur simultaneously and interact, such that the simple semantic demarcations break down. Similarly, the modelling of individual pedestrians has progressed the understanding of human–structure interaction, particularly for large amplitude lateral oscillations under crowd loading. In this paper, guided by the interaction of flutter and vortexinduced vibration in wind engineering, a framework is presented, which allows various human–structure interaction effects to coexist and interact, thereby providing a possible synthesis of previously disparate experimental and theoretical results.

Controlled growth orientation of carbon nanotube pillars by catalyst patterning in microtrenches

We present a novel method for controlling the growth orientation of individual carbon nanotube (CNT) microstructures on a silicon wafer substrate. Our method controls the CNT forest orientation by patterning the catalyst layer used in the CNTs growth on slanted KOH edges. The overlap of catalyst area on the horizontal bottom and sloped sidewall surfaces of the KOH-etched substrate enables precise variation of the growth direction. These inclined structures can profit from the outstanding mechanical, electrical, thermal, and optical properties of CNTs and can therefore improve the performance of several MEMS devices. Inclined CNT microstructures could for instance be used as cantilever springs in probe card arrays, as tips in dip-pen lithography, and as sensing element in advanced transducers. ©2009 IEEE.

Controlling the phase of an oscillator: A phase response curve approach

The paper discusses elementary control strategies to control the phase of an oscillator. Both feedforward and feedback (P and PI) control laws are designed based on the phase response curve (PRC) calculated from the linearized model. The performance is evaluated on a popular model of circadian oscillations. ©2009 IEEE.

Alternating spatial patterns for coordinated group motion

Motivated by recent observations of fish schools, we study coordinated group motion for individuals with oscillatory speed. Neighbors that have speed oscillations with common frequency, amplitude and average but different phases, move together in alternating spatial patterns, taking turns being towards the front, sides and back of the group. We propose a model and control laws to investigate the connections between these spatial dynamics, communication when sensing is range or direction limited, and convergence of coordinated group motions. ©2007 IEEE.

Analysis of interconnected oscillators by dissipativity theory

This paper employs dissipativity theory for the global analysis of limit cycles in particular dynamical systems of possibly high dimension. Oscillators are regarded as open systems that satisfy a particular dissipation inequality. It is shown that this characterization has implications for the global stability analysis of limit cycle oscillations: i) in isolated oscillators, ii) in interconnections of oscillators, and iii) for the global synchrony analysis in interconnections of identical oscillators. © 2007 IEEE.

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.

Fluid-structure interaction simulations of a tension-cone inflatable aerodynamic decelerator

A series of fluid-structure interaction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structure interaction computations combine large deformation analysis of thin shells with large-eddy simulation of compressible turbulent flows using a loosely coupled approach to enable quantification of the dynamics of the vehicle. The simulation results are compared with experiments carried out at the NASA Glenn Research Center. Reasonably good agreement between the simulations and the experiment is observed throughout a deflation cycle. The simulations help to illuminate the details of the dynamic progressive buckling of the tension-cone decelerator that ultimately results in the collapse of the structure as the inflation pressure is decreased. Furthermore, the tension-cone decelerator exhibits a transient oscillatory behavior under impulsive loading that ultimately dies out. The frequency of these oscillations was determined to be related to the acoustic time scale in the compressed subsonic region between the bow shock and the structure. As shown, when the natural frequency of the structure and the frequency of the compressed subsonic region approximately match, the decelerator exhibits relatively large nonaxisymetric oscillations. The observed response appears to be a fluid-structure interaction resonance resulting from an acoustic chamber (pistonlike) mode exciting the structure. Copyright © 2013 by Christopher Porter, R. Mark Rennie, Eric J. Jumper.

Modeling mode i fracture of bitumen films

The fracture behavior of thin films of bitumen in double cantilever beam (DCB) specimens was investigated over a wide range of temperature and loading rate conditions using finite-element analysis. The model includes a phenomenological model for the mechanical behavior of bitumen, implemented into a special-purpose finite-element user material subroutine, combined with a cohesive zone model (CZM) for simulating the fracture process. The finite-element model is validated against experimental results from laboratory tests of DCB specimens by comparing measured and predicted load-line deflection histories and fracture energy release rates. Computer simulation results agreed well with experimental data of DCB joints containing bitumen films in terms of peak stress, fracture toughness, and stress-strain history response. The predicted "normalized toughness," G=2h, was found to increase in a power-law manner with effective temperaturecompensated strain rate in the ductile region as previously observed experimentally. In the brittle regime, G=2h is virtually constant. The model successfully captured the ductile and brittle failure behavior of bitumen films in opening mode (tension) for stable crack growth conditions. © 2013 American Society of Civil Engineers.

Azimuthal instabilities in annular combustors: Standing and spinning modes

This theoretical study investigates spinning and standing modes in azimuthally symmetric annular combustion chambers. Both modes are observed in experiments and simulations, and an existing model predicts that spinning modes are the only stable state of the system. We extend this model to take into account the effect that the acoustic azimuthal velocity, u, has on the flames, and propose a phenomenological model based on experiments performed on transversely forced flames. This model contains a parameter, ä, that quantifies the influence that the transversal excitation has on the fluctuating heat release. For small values of ä, spinning modes are the only stable state of the system. In an intermediate range of ä, both spinning and standing modes are stable states. For large values of ä, standing modes are the only stable state. This study shows that a flame's response to azimuthal velocity fluctuations plays an important role in determining the type of thermoacoustic oscillations found in annular combustors. © 2013 The Authors.

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.