Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Analysis of basic processes inside the keyhole during deep penetration ND-YAG CW laser welding

During laser welding, the keyhole is generated by the recoil pressure induced by the evaporation processes occurring mainly on the front keyhole wall (KW). In order to characterize the evaporation process, we have measured this recoil pressure by using a plume deflection technique, where the plume generated for static conditions (i. e. with no sample displacement) is deflected by a transverse side gas jet. From the measurement of the plume deflection angle, the recoil pressure can be determined as a function of incident intensity and sample material. From these data one can estimate the pressure generated on the front KW, during laser welding. Therefore, the corresponding dynamic pressure exerted by the vapor plume expansion on the rear KW, in contact with the melt pool, can be also estimated. These pressures appear to be in close agreement with those generated by an additional side jet that has been used in previous experiments, for stabilizing the observed melt pool oscillations or fluctuations.

Adaptive control of combustion instabilities for unknown sign of the high frequency gain

A type of adaptive, closed-loop controllers known as self-tuning regulators present a robust method of eliminating thermoacoustic oscillations in modern gas turbines. These controllers are able to adapt to changes in operating conditions, and require very little pre-characterisation of the system. One piece of information that is required, however, is the sign of the system's high frequency gain (or its 'instantaneous gain'). This poses a problem: combustion systems are infinite-dimensional, and so this information is never known a priori. A possible solution is to use a Nussbaum gain, which guarantees closed-loop stability without knowledge of the sign of the high frequency gain. Despite the theory for such a controller having been developed in the 1980s, it has never, to the authors' knowledge, been demonstrated experimentally. In this paper, a Nussbaum gain is used to stabilise thermoacoustic instability in a Rijke tube. The sign of the high frequency gain of the system is not required, and the controller is robust to large changes in operating conditions - demonstrated by varying the length of the Rijke tube with time. Copyright © 2008 by Simon J. Illingworth & Aimee S. Morgans.

VGs for a normal SBLI with a downstream diffuser

New types of vortex generators for boundary layer control were investigated experimentally in a flow field which contains a Mach 1.4 normal Shockwave followed by a subsonic diffuser. A parametric study of device height and distance upstream of the normal shock was undertaken with two novel devices: ramped-vanes and split-ramps. Flowfield diagnostics included high-speed Schlieren, oil flow visualization, and pitot-static pressure measurements. A number of flowfield parameters including flow separation, pressure recovery, centerline incompressible boundary layer shape factor, and shock stability were analyzed and compared to the baseline. All configurations tested yielded an elimination of centerline flow separation with the presence of the vortex generators. However, the devices also tended to increase the three-dimensionality of the flow with increased side-wall interaction. When located 25δo upstream of the normal shock, the largest ramped-vane device (whose height was about 0.75 the incoming uncontrolled boundary layer thickness, δo) yielded the smallest centerline incompressible shape factor and the least streamwise oscillations of the normal shock. However, additional studies are needed to better understand the corner interaction effects, which are substantial. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Model-based thermoacoustic oscillation control via ℋ ∞ loop-shaping and integral quadratic constraints

Lean premixed prevaporized (LPP) technology has been widely used in the new generation of gas turbines in which reduced emissions are a priority. However, such combustion systems are susceptible to the damage of self-excited oscillations. Feedback control provide a way of preventing such dynamic stabilities. A flame dynamics assumption is proposed for a recently developed unsteady heat release model, the robust design technique, ℋ ∞ loop-shaping, is applied for the controller design and the performance of the controller is confirmed by simulations of the closed-loop system. The Integral Quadratic Constraints(IQC) method is employed to prove the stability of the closed-loop system. ©2010 IEEE.

Triggering in a thermoacoustic system with stochastic noise

This paper explores the mechanism of triggering in a simple thermoacoustic system, the Rijke tube. It is demonstrated that additive stochastic perturbations can cause triggering before the linear stability limit of a thermoacoustic system. When triggering from low noise amplitudes, the system is seen to evolve to self-sustained oscillations via an unstable periodic solution of the governing equations. Practical stability is introduced as a measure of the stability of a linearly stable state when finite perturbations are present. The concept of a stochastic stability map is used to demonstrate the change in practical stability limits for a system with a subcritical bifurcation, once stochastic terms are included. The practical stability limits are found to be strongly dependent on the strength of noise.

Transient growth and triggering in the horizontal Rijke tube Matthew

This theoretical paper examines a non-normal and non-linear model of a horizontal Rijke tube. Linear and non-linear optimal initial states, which maximize acoustic energy growth over a given time from a given energy, are calculated. It is found that non-linearity and non-normality both contribute to transient growth and that, for this model, linear optimal states are only a good predictor of non-linear optimal states for low initial energies. Two types of non-linear optimal initial state are found. The first has strong energy growth during the first period of the fundamental mode but loses energy thereafter. The second has weaker energy growth during the first period but retains high energy for longer. The second type causes triggering to self-sustained oscillations from lower energy than the first and has higher energy in the fundamental mode. This suggests, for instance, that low frequency noise will be more effective at causing triggering than high frequency noise.

Late Pleistocene divergence and subsequent population expansion of two closely related fish species, Japanese anchovy (Engraulis japonicus) and Australian anchovy (Engraulis australis)

Climatic oscillations during the Pleistocene ice ages produced great changes in species' geographical distribution and abundance, which could be expected to have genetic consequences. Living in the temperate upwelling zones of the northwestern Pacific, Ja

Synchronization of modulated traveling baroclinic waves in a periodically forced, rotating fluid annulus.

Frequency entrainment and nonlinear synchronization are commonly observed between simple oscillatory systems, but their occurrence and behavior in continuum fluid systems are much less well understood. Motivated by possible applications to geophysical fluid systems, such as in atmospheric circulation and climate dynamics, we have carried out an experimental study of the interaction of fully developed baroclinic instability in a differentially heated, rotating fluid annulus with an externally imposed periodic modulation of the thermal boundary conditions. In quasiperiodic and chaotic amplitude-modulated traveling wave regimes, the results demonstrate a strong interaction between the natural periodic modulation of the wave amplitude and the externally imposed forcing. This leads to partial or complete phase synchronization. Synchronization effects were observed even with very weak amplitudes of forcing, and were found with both 1:1 and 1:2 frequency ratios between forcing and natural oscillations.

Synchronization in baroclinic systems

Synchronization of periodic and chaotic oscillations between two coupled rotating baroclinic fluid systems will be presented. The numerical part of the study involves a pair of coupled two-layer quasigeostrophic models, and the experimental part comprises two thermally coupled baroclinic fluid annuli, rotating one above the other on the same turntable. Phase synchronization and imperfect synchronization (phase slips) have been found in both model and experiments, and model simulations also exhibit chaos-destroying synchronization. © 2008 IOP Publishing Ltd.

Synchronization in climate dynamics and other extended systems

Synchronization is now well established as representing coherent behaviour between two or more otherwise autonomous nonlinear systems subject to some degree of coupling. Such behaviour has mainly been studied to date, however, in relatively low-dimensional discrete systems or networks. But the possibility of similar kinds of behaviour in continuous or extended spatiotemporal systems has many potential practical implications, especially in various areas of geophysics. We review here a range of cyclically varying phenomena within the Earth's climate system for which there may be some evidence or indication of the possibility of synchronized behaviour, albeit perhaps imperfect or highly intermittent. The exploitation of this approach is still at a relatively early stage within climate science and dynamics, in which the climate system is regarded as a hierarchy of many coupled sub-systems with complex nonlinear feedbacks and forcings. The possibility of synchronization between climate oscillations (global or local) and a predictable external forcing raises important questions of how models of such phenomena can be validated and verified, since the resulting response may be relatively insensitive to the details of the model being synchronized. The use of laboratory analogues may therefore have an important role to play in the study of natural systems that can only be observed and for which controlled experiments are impossible. We go on to demonstrate that synchronization can be observed in the laboratory, even in weakly coupled fluid dynamical systems that may serve as direct analogues of the behaviour of major components of the Earth's climate system. The potential implications and observability of these effects in the long-term climate variability of the Earth is further discussed. © 2010 Springer-Verlag Berlin Heidelberg.

Molecular phylogeography and population structure of a mid-elevation montane frog Leptobrachium ailaonicum in a fragmented habitat of southwest China

Leptobrachium ailaonicum is a vulnerable anuran restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes at mid-elevations (approximately 2000-2600 m) in the subtropical Mount Wuliang and Mount Ailao ranges in southwest China (Yunnan Province) and northern Vietnam. Given high habitat specificity and lack of suitable habitat in lower elevations between these ranges, we hypothesized limited gene flow between populations throughout its range. We used two mitochondrial genes to construct a phylogeographic pattern within this species in order to test our hypothesis. We also examined whether this phylogeographic pattern is a response to past geological events and/or climatic oscillations. A total of 1989 base pairs were obtained from 81 individuals of nine populations yielding 51 unique haplotypes. Both Bayesian and maximum parsimony phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographical regions separated by deep river valleys. These results suggest a long history of allopatric separation by vicariance. The distinct geographic distributions of four major clades and the estimated divergence time suggest spatial and temporal separations that coincide with climatic and paleogeographic changes following the orogeny and uplift of Mount Ailao during the late Miocene to mid Pliocene in southwest China. At the southern distribution, the presence of two sympatric yet differentiated clades in two areas are interpreted as a result of secondary contact between previously allopatric populations during cooler Pleistocene glacial cycles. Analysis of molecular variance indicates that most of the observed genetic variation occurs among the four regions implying long-term interruption of maternal gene flow, suggesting that L ailaonicum may represent more than one distinct species and should at least be separated into four management units corresponding to these four geographic lineages for conservation. (C) 2009 Elsevier Inc. All rights reserved.

Investigation of self excited and acoustically forced combustion instabilities with laser diagnostics

Combustion oscillations in gas turbines can result in serious damage. One method used to predict such oscillations is to analyze the combustor acoustics using a simple linear model. Such a model requires a flame transfer function to describe the response of the heat release to flow perturbations inside the combustor. This paper reports on the application of Planar Laser Induced Fluorescence (PLIF) of OH radicals to analyze the response of a lean premixed flame to oncoming flow perturbations. Both self-excited oscillations and low amplitude forced oscillations at various frequencies are investigated in an atmospheric pressure model combustor rig. In order to visualize fluctuations of local fuel distribution, acetone-PLIF was also applied in non-reacting and acoustically forced flows at oscillation frequencies of 200 Hz and 510 Hz, respectively. OH-PLIF images were acquired over a range of operating parameters. The results presented in this paper originate from data sets acquired at fixed phase angles during the oscillation cycle. Comparative experiments in self excited and forced acoustic oscillations show that the flame and the combustion intensity develop similarly throughout the pressure cycle in both cases. Although the peak fluorescence intensities differ between self excited and the forced instabilities, there is a clear correspondence in the observed frequency and phase information from the two cases. This result encourages a comparison of the OH-PLIF and the acetone-PLIF results. Quantitative measurements of the equivalence ratio in specific areas of the measurement plane offer insight on the complex phenomena coupling acoustic perturbations, i.e. flow velocity fluctuations, to fluctuations in fuel distribution and combustion intensity, ultimately resulting in self excited combustion oscillations.

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

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.