Mitigating instability by actuating the swirler in a combustor

Thermoacoustic instability in a lean premixed combustor is a major impediment towards reliable operation of gas turbine engines for both aerospace and land based applications. In this communication, we investigate the following concept: in a laboratory combustor, could the otherwise static swirler be actuated to a rotary motion, such that the higher intensity turbulence and higher swirl number generated in the flame stabilization region might alter the flame position, structure and thereby assist in mitigating thermoacoustic instabilities? Results obtained using microphone and high speed imaging, show prominent reductions in the amplitudes of the first mode of the thermoacoustically unstable flame, with increased rotation rate of the swirler. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Depletion of nonlinearity in magnetohydrodynamic turbulence: Insights from analysis and simulations

It is shown how suitably scaled, order-m moments, D-m(+/-), of the Elsasser vorticity fields in three-dimensional magnetohydrodynamics (MHD) can be used to identify three possible regimes for solutions of the MHD equations with magnetic Prandtl number P-M = 1. These vorticity fields are defined by omega(+/-) = curl z(+/-) = omega +/- j, where z(+/-) are Elsasser variables, and where omega and j are, respectively, the fluid vorticity and current density. This study follows recent developments in the study of three-dimensional Navier-Stokes fluid turbulence Gibbon et al., Nonlinearity 27, 2605 (2014)]. Our mathematical results are then compared with those from a variety of direct numerical simulations, which demonstrate that all solutions that have been investigated remain in only one of these regimes which has depleted nonlinearity. The exponents q(+/-) that characterize the inertial range power-law dependencies of the z(+/-) energy spectra, epsilon(+/-)(k), are then examined, and bounds are obtained. Comments are also made on (a) the generalization of our results to the case P-M not equal 1 and (b) the relation between D-m(+/-) and the order-m moments of gradients of magnetohydrodynamic fields, which are used to characterize intermittency in turbulent flows.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Evaluation of the Cumulative Formation of High-Speed Liquid Jets

This paper describes the generation of pulsed, high-speed liquid jets using the cumulation method. This work mainly includes (1) the design of the nozzle assembly, (2) the measurement of the jet velocity and (3) flow visualization of the injection sequences. The cumulation method can be briefly described as the liquid being accelerated first by the impact of a moving projectile and then further after it enters a converging section. The experimental results show that the cumulation method is useful in obtaining a liquid jet with high velocity. The flow visualization shows the roles of the Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the breakup of the liquid depend on the jet diameter and the downstream distance. When the liquid jet front is far downstream from the nozzle exit, the jet is decelerated by air drag. Meanwhile, large coherent vortex structures are formed surrounding the jet. The liquid will break up totally by the action of these vortices. Experimental results showing the effect of the liquid volume on the jet velocity are also included in this paper. Finally, a method for measuring the jet velocity by cutting two carbon rods is examined.

Surface-waves, stability, and scattering for a lined duct with flow

We consider a straight cylindrical duct with a steady subsonic axial flow and a reacting boundary (e.g. an acoustic lining). The wave modes are separated into ordinary acoustic duct modes, and surface modes confined to a small neighbourhood of the boundary. Many researchers have used a mass-spring-damper boundary model, for which one surface mode has previously been identified as a convective instability; however, we show the stability analysis used in such cases to be questionable. We investigate instead the stability of the surface modes using the Briggs-Bers criterion for a Flügge thin-shell boundary model. For modest frequencies and wavenumbers the thin-shell has an impedance which is effectively that of a mass-spring-damper, although for the large wavenumbers needed for the stability analysis the thin-shell and mass-spring-damper impedances diverge, owing to the thin shell's bending stiffness. The thin shell model may therefore be viewed as a regularization of the mass-spring-damper model which accounts for nonlocally-reacting effects. We find all modes to be stable for realistic thin-shell parameters, while absolute instabilities are demonstrated for extremely thin boundary thicknesses. The limit of vanishing bending stiffness is found to be a singular limit, yielding absolute instabilities of arbitrarily large temporal growth rate. We propose that the problems with previous stability analyses are due to the neglect of something akin to bending stiffness in the boundary model. Our conclusion is that the surface mode previously identified as a convective instability may well be stable in reality. Finally, inspired by Rienstra's recent analysis, we investigate the scattering of an acoustic mode as it encounters a sudden change from a hard-wall to a thin-shell boundary, using a Wiener-Hopf technique. The thin-shell is considered to be clamped to the hard-wall. The acoustic mode is found to scatter into transmitted and reflected acoustic modes, and surface modes strongly linked to the solid waves in the boundary, although no longitudinal or transverse waves within the boundary are excited. Examples are provided that demonstrate total transmission, total reflection, and a combination of the two. This thin-shell scattering problem is preferable to the mass-spring-damper scattering problem presented by Rienstra, since the thin-shell problem is fully determined and does not need to appeal to a Kutta-like condition or the inclusion of an instability in order to avoid a surface-streamline cusp at the boundary change.

Continuous spectrum growth and modal instability in swirling duct flow

We present results on the stability of compressible inviscid swirling flows in an annular duct. Such flows are present in aeroengines, for example in the by-pass duct, and there are also similar flows in many aeroacoustic or aeronautical applications. The linearised Euler equations have a ('critical layer') singularity associated with pure convection of the unsteady disturbance by the mean flow, and we focus our attention on this region of the spectrum. By considering the critical layer singularity, we identify the continuous spectrum of the problem and describe how it contributes to the unsteady field. We find a very generic family of instability modes near to the continuous spectrum, whose eigenvalue wavenumbers form an infinite set and accumulate to a point in the complex plane. We study this accumulation process asymptotically, and find conditions on the flow to support such instabilities. It is also found that the continuous spectrum can cause a new type of instability, leading to algebraic growth with an exponent determined by the mean flow, given in the analysis. The exponent of algebraic growth can be arbitrarily large. Numerical demonstrations of the continuous spectrum instability, and also the modal instabilities are presented.

Linear Stability Analysis of Convection in Two-Layer System with an Evaporating Vapor-Liquid Interface

Classical theories have successfully provided an explanation for convection in a liquid layer heated from below without evaporation. However, these theories are inadequate to account for the convective instabilities in an evaporating liquid layer, especially in the case when it is cooled from below. In the present paper, we study the onset of Marangoni convection in a liquid layer being overlain by a vapor layer.A new two-sided model is put forward instead of the one-sided model in previous studies. Marangoni-Bénard instabilities in evaporating liquid thin layers are investigated with a linear instability analysis. We define a new evaporation Biot number, which is different from that in previous studies and discuss the influences of reference evaporating velocity and evaporation Biot number on the vapor-liquid system. At the end, we explain why the instability occurs even when an evaporating liquid layer is cooled from below.

Influence of liquid bridge volume on instability of floating half zone convection

Thermocapillary instabilities on floating half zone convection in microgravity environment were investigated by linear instability analysis method. The critical Marangoni numbers were obtained and compared with the experimental ones. The influences of the liquid bridge volume and the aspect ratio on the critical Marangoni number were analyzed. It is found that the liquid bridge volume and the aspect ratio have great influence on the critical Marangoni number. There was a gap region where the oscillatory convection will not be observed in present analyses and in experiments in the curve of the critical Marangoni number vs the liquid bridge volume for the case of large Prandtl number and small aspect ratio.

The Hydrothermal Wave Of Large-Prandtl-Number Fluid In A Shallow Cavity

The hydrothermal wave was investigated numerically for large-Prandtl-number fluid (Pr = 105.6) in a shallow cavity with different heated sidewalls. The traveling wave appears and propagates in the direction opposite to the surface flow (upstream) in the case of zero gravity when the applied temperature difference grows and over the critical value. The phase relationships of the disturbed velocity, temperature and pressure demonstrate that the traveling wave is driven by the disturbed temperature, which is named hydrothermal wave. The hydrothermal wave is so weak that the oscillatory flow field and temperature distribution can hardly be observed in the liquid layer. The exciting mechanism of hydrothermal wave is analyzed and discussed in the present paper.

Progress In Modeling Of Fluid Flows In Crystal Growth Processes

Modeling of fluid flows in crystal growth processes has become an important research area in theoretical and applied mechanics. Most crystal growth processes involve fluid flows, such as flows in the melt, solution or vapor. Theoretical modeling has played an important role in developing technologies used for growing semiconductor crystals for high performance electronic and optoelectronic devices. The application of devices requires large diameter crystals with a high degree of crystallographic perfection, low defect density and uniform dopant distribution. In this article, the flow models developed in modeling of the crystal growth processes such as Czochralski, ammonothermal and physical vapor transport methods are reviewed. In the Czochralski growth modeling, the flow models for thermocapillary flow, turbulent flow and MHD flow have been developed. In the ammonothermal growth modeling, the buoyancy and porous media flow models have been developed based on a single-domain and continuum approach for the composite fluid-porous layer systems. In the physical vapor transport growth modeling, the Stefan flow model has been proposed based on the flow-kinetics theory for the vapor growth. In addition, perspectives for future studies on crystal growth modeling are proposed. (c) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Thermovibrational Instability Of Rayleigh-Marangoni-Benard Convection In Two-Layer Fluid Systems

The thermovibrational instability of Rayleigh-Marangoni-Benard convection in a two-layer system under the high-frequency vibration has been investigated by linear instability analysis in the present paper. General equations for the description of the convective flow and within this framework, the generalized Boussinesq approximation are formulated. These equations are dealt with using the averaging method. The theoretical analysis results show that the high-frequency thermovibrations can change the Marangoni-Benard convection instabilities as well as the oscillatory gaps of the Rayleigh-Marangoni-Benard convection in two-layer liquid systems. It is found that vertical high-frequency vibrations can delay convective instability of this system, and damp the convective flow down. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Experimental Study On Liquid Free Surface In Buoyant-Thermocapillary Convection

We investigate the surface deformations of buoyant-thermocapillary convection in a rectangular cavity clue to gravity and temperature gradient between the two sidewalls. The cavity is 52mm x 42mm in horizontal cross section, the thickness of liquid layer h is changed from 2.5mm to 6.5mm. Surface deformations of h = 3.5mm and 6.0mm are discussed and compared. Temperature difference is increased gradually, and the flow in the liquid layer will change from stable convection to unstable convection. Two kinds of optical diagnostic system with image processor are developed for study of the kinetics of buoyant-thermocapillary convection, they give out the information of liquid free surface. The quantitative results are calculated by Fourier transform and correlation analysis, respectively. With the increasing temperature gradient, surface deformations calculated are more declining. It is interesting phenomenon that the inclining directions of the convections in thin and thick liquid layers are different. For a thin layer, the convection is mainly controlled by thermocapillary effect. However, for a thick layer, the convection is mainly controlled by buoyancy effect. The surface deformation theoretically analysed is consistent with our experimental results. The present experiment proves that surface deformation is related to temperature gradient and thickness of the liquid layer. In other words, surface deformation lies on capillary convection and buoyancy convection.

Thermocapillary Convection In An Inhomogeneous Porous Layer

A new model consisting of an inhomogeneous porous medium saturated by incompressible fluid is investigated. We focus on the effects of inhomogeneity for the streamline patterns and instabilities of the system. Influences of the 'mean porosity' and gradient of distributions of porosity are also emphasized. The results cannot be obtained by studying the media with constant porosity as carried out by other researchers, and have not been discussed before.

On the mechanism of void growth and the effect of straining mode in ductile materials

Finite element analysis is employed to investigate void growth embedded in elastic-plastic matrix material. Axisymmetric and plane stress conditions are considered. The simulation of void growth in a unit cell model is carried out over a wide range of triaxial tensile stressing or large plastic straining for various strain hardening materials to study the mechanism of void growth in ductile materials. Triaxial tension and large plastic strain encircling around the void are found to be of most importance for driving void growth. The straining mode of incremental loading which favors the necessary strain concentration around void for its growth can be characterized by the vanishing condition of a parameter called "the third invariant of generalized strain rate". Under this condition, it accentuates the internal strain concentration and the strain energy stored/dissipated within the material layer surrounding the void. Experimental results are cited to justify the effect of this loading parameter. (C) 2000 Elsevier Science Ltd. All rights reserved.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.