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.

Global nonlinear dynamics of thin aerofoil wakes

In the present investigation of thin aerofoil wakes we compare the global nonlinear dynamics, obtained by direct numerical simulations, to the associated local instability features, derived from linear stability analyses. A given configuration depends on two control parameters: the Reynolds number Re and the adverse pressure gradient m (with m < 0) prevailing at the aerofoil trailing edge. Global instability is found to occur for large enough Re and |m|; the naturally selected frequency is determined by the local absolute frequency prevailing at the trailing edge. © 2010 Springer Science+Business Media B.V.

Phase synchronization between stratospheric and tropospheric quasi-biennial and semi-annual oscillations

A combination of singular systems analysis and analytic phase techniques are used to investigate the possible occurrence in observations of coherent synchronization between quasi-biennial and semi-annual oscillations (QBOs; SAOs) in the stratosphere and troposphere. Time series of zonal mean zonal winds near the Equator are analysed from the ERA-40 and ERA-interim reanalysis datasets over a ∼ 50-year period. In the stratosphere, the QBO is found to synchronize with the SAO almost all the time, but with a frequency ratio that changes erratically between 4:1, 5:1 and 6:1. A similar variable synchronization is also evident in the tropical troposphere between semi-annual and quasi-biennial cycles (known as TBOs). Mean zonal winds from ERA-40 and ERA-interim, and also time series of indices for the Indian and West Pacific monsoons, are commonly found to exhibit synchronization, with SAO/TBO ratios that vary between 4:1 and 7:1. Coherent synchronization between the QBO and tropical TBO does not appear to persist for long intervals, however. This suggests that both the QBO and tropical TBOs may be separately synchronized to SAOs that are themselves enslaved to the seasonal cycle, or to the annual cycle itself. However, the QBO and TBOs are evidently only weakly coupled between themselves and are frequently found to lose mutual coherence when each changes its frequency ratio to its respective SAO. This suggests a need to revise a commonly cited paradigm that advocates the use of stratospheric QBO indices as a predictor for tropospheric phenomena such as monsoons and hurricanes. © 2012 Royal Meteorological Society.

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

Oscillations in a decentralized economic network

The study of exchange markets dates back to LeonWalras's general equilibrium theory. Since then the economic market has been studied for its' equilibrium properties, fairness of allocations of private and public goods, and even the psychological incentives of participants. This paper studies the dynamics of an exchange economy built on a network of markets where consumers trade with suppliers to optimize utility. Viewing the market in as a decentralized network we study the system from the usual control theory point of view, evaluating the system's dynamic performance, stability and robustness. It is shown that certain consumer demand dynamics can lead to oscillations while others can converge to optimal allocations. © 2011 IFAC.

Linear and nonlinear model large-eddy simulations of a plane jet

The nonlinear Kosovic, and mixed Leray and α subgrid scale models are contrasted with linear Smagorinsky and Yoshizawa Large Eddy Simulations for a Re = 4000 plane jet simulation. Comparisons are made with Direct Numerical Simulation data and measurements. Global properties of the jet such as the spreading and centreline velocity decay rates are investigated. The mean-flow and turbulence parameters in the self-similar region are also studied. All models generally give encouraging agreement with the Direct Numerical Simulation data and reliable measurements. Solution differences for the models are relatively minor, none giving clear improvements for all data comparisons.

Mixed nonlinear LES for DES suitable flows

It is shown that a new mixed nonlinear/eddy viscosity LES model reproduces profiles better than a number of competing nonlinear and mixed models for plane channel flow. The objective is an LES method that produces a fully resolved turbulent boundary layer and could be applied to a variety of aerospace problems that are currently studied with RANS, RANS-LES, or DES methods that lack a true turbulent boundary layer. There are two components to the new model. One an eddy viscosity based upon the advected subgrid scale energy and a relatively small coefficient. Second, filtered nonlinear terms based upon the Leray regularization. Coefficients for the eddy viscosity and nonlinear terms come from LES tests in decaying, isotropic turbulence. Using these coefficients, the velocity profile matches measurements data at Reτ ≈ 1000 exactly. Profiles of the components of kinetic energy have the same shape as in the experiment, but the magnitudes differ by about 25%. None of the competing LES gets the shape correct. This method does not require extra operations at the transition between the boundary layer and the interior flow.

Nonlinear optics in Silicon photonic crystal cavities

Silicon is known to be a very good material for the realization of high-Q, low-volume photonic cavities, but at the same it is usually considered as a poor material for nonlinear optical functionalities like second-harmonic generation, because its second-order nonlinear susceptibility vanishes in the dipole approximation. In this work we demonstrate that nonlinear optical effects in silicon nanocavities can be strongly enhanced and even become macroscopically observable. We employ photonic crystal nanocavities in silicon membranes that are optimized simultaneously for high quality factor and efficient coupling to an incoming beam in the far field. Using a low-power, continuous-wave laser at telecommunication wavelengths as a pump beam, we demonstrate simultaneous generation of second- and third harmonics in the visible region, which can be observed with a simple camera. The results are in good agreement with a theoretical model that treats third-harmonic generation as a bulk effect in the cavity region, and second-harmonic generation as a surface effect arising from the vertical hole sidewalls. Optical bistability is also observed in the silicon nanocavities and its physical mechanisms (optical, due to two-photon generation of free carriers, as well as thermal) are investigated. © 2011 IEEE.

Modelling and control of nonlinear systems using Gaussian processes with partial model information

Gaussian processes are gaining increasing popularity among the control community, in particular for the modelling of discrete time state space systems. However, it has not been clear how to incorporate model information, in the form of known state relationships, when using a Gaussian process as a predictive model. An obvious example of known prior information is position and velocity related states. Incorporation of such information would be beneficial both computationally and for faster dynamics learning. This paper introduces a method of achieving this, yielding faster dynamics learning and a reduction in computational effort from O(Dn2) to O((D - F)n2) in the prediction stage for a system with D states, F known state relationships and n observations. The effectiveness of the method is demonstrated through its inclusion in the PILCO learning algorithm with application to the swing-up and balance of a torque-limited pendulum and the balancing of a robotic unicycle in simulation. © 2012 IEEE.

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.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Wave propagation in nonlinear one-dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one-dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross-correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross-correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loop under steady-state condition. Copyright © 2008 John Wiley & Sons, Ltd.