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.

Electrical detection of trapped charge displacements using an ultra-thin SOI percolation transistor

Electrical detection of solid-state charge qubits requires ultrasensitive charge measurement, typically using a quantum point contact or single-electron-transistor, which imposes strict limits on operating temperature, voltage and current. A conventional FET offers relaxed operating conditions, but the back-action of the channel charge is a problem for such small quantum systems. Here, we discuss the use of a percolation transistor as a measurement device, with regard to charge sensing and backaction. The transistor is based on a 10nm thick SOI channel layer and is designed to measure the displacement of trapped charges in a nearby dielectric. At cryogenic temperatures, the trapped charges result in strong disorder in the channel layer, so that current is constrained to a percolation pathway in sub-threshold conditions. A microwave driven spatial Rabi oscillation of the trapped charge causes a change in the percolation pathway, which results in a measurable change in channel current. © The Electrochemical Society.

Model-based control of thermoacoustic instabilities in partially premixed lean combustion - A design case study

Self-excited oscillation is becoming a major issue in low-emission, lean partially premixed combustion systems, and active control has been shown to be a feasible method to suppress such instabilities. A number of robust control methods are employed to obtain a feedback controller and it is observed that the robustness to system uncertainty is significantly better for a low complexity controller in spite of the norms being similar. Moreover, we demonstrate that closed-loop stability for such a complex system can be proved via use of the integral quadratic constraint method. Open- and closed-loop nonlinear simulations are provided. © 2013 Copyright Taylor and Francis Group, LLC.

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.

Adaptive combustion instability control with saturation

The control of a class of combustion systems, suceptible to damage from self-excited combustion oscillations, is considered. An adaptive stable controller, called Self-Tuning Regulator (STR), has recently been developed, which meets the apparently contradictory challenge of relying as little as possible on a particular combustion model while providing some guarantee that the controller will cause no harm. The controller injects some fuel unsteadily into the burning region, thereby altering the heat release, in response to an input signal detecting the oscillation. This paper focuses on an extension of the STR design, when, due to stringent emission requirements and to the danger of flame extension, the amount of fuel used for control is limited in amplitude. A Lyapunov stability analysis is used to prove the stability of the modified STR when the saturation constraint is imposed. The practical implementation of the modified STR remains straightforward, and simulation results, based on the nonlinear premixed flame model developed by Dowling, show that in the presence of a saturation constraint, the self-excited oscillations are damped more rapidly with the modified STR than with the original STR. © 2001 by S. Evesque. Published by the American Institute of Aeronautics and Astronautics, Inc.

Minimalistic control of biped walking in rough terrain

Toward our comprehensive understanding of legged locomotion in animals and machines, the compass gait model has been intensively studied for a systematic investigation of complex biped locomotion dynamics. While most of the previous studies focused only on the locomotion on flat surfaces, in this article, we tackle with the problem of bipedal locomotion in rough terrains by using a minimalistic control architecture for the compass gait walking model. This controller utilizes an open-loop sinusoidal oscillation of hip motor, which induces basic walking stability without sensory feedback. A set of simulation analyses show that the underlying mechanism lies in the "phase locking" mechanism that compensates phase delays between mechanical dynamics and the open-loop motor oscillation resulting in a relatively large basin of attraction in dynamic bipedal walking. By exploiting this mechanism, we also explain how the basin of attraction can be controlled by manipulating the parameters of oscillator not only on a flat terrain but also in various inclined slopes. Based on the simulation analysis, the proposed controller is implemented in a real-world robotic platform to confirm the plausibility of the approach. In addition, by using these basic principles of self-stability and gait variability, we demonstrate how the proposed controller can be extended with a simple sensory feedback such that the robot is able to control gait patterns autonomously for traversing a rough terrain. © 2010 Springer Science+Business Media, LLC.

Evidence of print gap airflow affecting web printing quality

The current study extends our earlier investigation on the real-time dynamics of print gap airflow around a single jetted drop over a moving substrate. In the present work, simulated web press printing was performed using a stationary grey-scale commercial inkjet print head to print full-width block of solid colour images onto a paper substrate with extended print gaps. The resultant printed images exhibit patterns or 'wood-graining' effects which become more prevalent as the relevant Reynolds number (Re) increases. The high-resolution scans of the printed images revealed that the patterns are created by oscillation and coalescence of neighboring printed tracks across the web. The phenomenon could be a result of drop stream perturbations caused by unsteady print gap airflow of the type similar to that observed in the previous study. ©2013; Society for Imaging Science and Technology.

Determination of some reactive oxygen species scavengers based on the peroxidase-oxidase oscillator

This paper reports a new method for detection of ROS scavengers including superoxide dismutase, ascorbic acid and glutathione based on a 'probe' of peroxidase-oxidase biochemical oscillator. The oscillation period and amplitude change with different concentrations of scavengers. The linear ranges of superoxide dismutase, ascorbic acid and glutathione are respectively 1.56 x 10(-4)-1.56 x 10(-3) mg mL(-1), 1.75 x 10(-7) -1.75 x 10(-5) mol L-1 and 9.38 x 10(-7) -7.5 x 10(-5) mol L-1. The selectivity, linearity and precision for superoxide dismutase, ascorbic acid, and glutathione are presented and discussed. The results compared well with other standard methods for determination of superoxide dismutase, ascorbic acid and glutathione. Some possible steps in the overall reaction mechanisms are discussed.

HYDROGRAPHIC VARIATIONS IN SIZE STRUCTURE OF PLANKTON BIOMASS IN A CHANGJIANG FLOODPLAIN LAKE

The size structure of the planktonic community in a Changjiang floodplain lake (Lake Chenhu, Hubei, P. R. China) was described for the inundation period of May through September 1983. The modality of the Sheldon-type size distributions changed hydrographically with the spectral profiles being bimodal during low, rising, mid-high and falling water phases, and trimodal soon after filling and shortly before falling. The modal peaks corresponded respectively to the dominant organisms of chlorophytes and nauplii, while the troughs centered on the bacteria and macrocrustacean size classes in the lake. The slope of the normalized biomass spectrum (an index of plankton size distribution) was less than -1.0 for the filling and falling phases or close to -1.0 for the high water period, indicating that the planktonic biomass tended to decrease or evenly distributes across logarithmically ordered size classes, respectively. This observed variation in the size distribution of the plankton community mainly resulted from changes in water levels and contents of particulate inorganic matter (PIM) in the lake.

Interference in transport through double barriers in interacting quantum wires

We investigate interference effects of the backscattering current through a double-barrier structure in an interacting quantum wire attached to noninteracting leads. Depending on the interaction strength and the location of the barriers, the backscattering current exhibits different oscillation and scaling characteristics with the applied voltage in the strong and weak interaction cases. However, in both cases, the oscillation behaviors of the backscattering current are mainly determined by the quantum mechanical interference due to the existence of the double barriers.

Directional emission InP/GaInAsP square-resonator microlasers

InP/GaInAsP square-resonator microlasers with an output waveguide connected to the midpoint of one side of the square are fabricated by standard photolithography and inductively-coupled-plasma etching technique. For a 20-mu m-side square microlaser with a 2-mu m-wide output waveguide, cw threshold current is 11 mA at room temperature, and the highest mode Q factor is 1.0 X 10(4) measured from the mode linewidth at the injection current of 10 mA. Multimode oscillation is observed with the lasing mode wavelength 1546 nm and the side-mode suppression ratio of 20 dB at the injection current of 15 mA. (C) 2008 Optical Society of America

Coherence loss and recovery of an electron spin coupled inhomogeneously to a one-dimensional interacting spin bath: An adaptive time-dependent density-matrix renormalization group study

Coherence evolution and echo effect of an electron spin, which is coupled inhomogeneously to an interacting one-dimensional finite spin bath via hyperfine-type interaction, are studied using the adaptive time-dependent density-matrix renormalization group method. It is found that the interplay of the coupling inhomogeneity and the transverse intrabath interactions results in two qualitatively different coherence evolutions, namely, a coherence-preserving evolution characterized by periodic oscillation and a complete decoherence evolution. Correspondingly, the echo effects induced by an electron-spin flip at time tau exhibit stable recoherence pulse sequence for the periodic evolution and a single peak at root 2 tau for the decoherence evolution, respectively. With the diagonal intrabath interaction included, the specific feature of the periodic regime is kept, while the root 2 tau-type echo effect in the decoherence regime is significantly affected. To render the experimental verifications possible, the Hahn echo envelope as a function of tau is calculated, which eliminates the inhomogeneous broadening effect and serves for the identification of the different status of the dynamic coherence evolution, periodic versus decoherence.

High-performance 2 mm gate width GaNHEMTs on 6H-SiC with output power of 22.4 W @ 8 GHz

Optimized AlGaN/AlN/GaN high electron mobility transistors (HEMTs) structures were grown on 2-in semi-insulating (SI) 6H-SiC substrate by metal-organic chemical vapor deposition (MOCVD). The 2-in. HEMT wafer exhibited a low average sheet resistance of 305.3 Omega/sq with a uniformity of 3.85%. The fabricated large periphery device with a dimension of 0.35 pm x 2 nun demonstrated high performance, with a maximum DC current density of 1360 mA/mm, a transconductance of 460 mS/mm, a breakdown voltage larger than 80 V, a current gain cut-off frequency of 24 GHz and a maximum oscillation frequency of 34 GHz. Under the condition of continuous-wave (CW) at 9 GHz, the device achieved 18.1 W output power with a power density of 9.05 W/mm and power-added-efficiency (PAE) of 36.4%. While the corresponding results of pulse condition at 8 GHz are 22.4 W output power with 11.2 W/mm power density and 45.3% PAE. These are the state-of-the-art power performance ever reported for this physical dimension of GaN HEMTs based on SiC substrate at 8 GHz. (c) 2008 Elsevier Ltd. All rights reserved.