对流速度对燃料液滴火焰形态及燃烧的影响

考虑气相非稳态及液滴内部环流,建立运动液滴非稳态蒸发燃烧模型.模型采用动网格方法精确追踪液滴表面位置,采用守恒方程组更新液滴表面边界条件.根据单步全局化学反应机理,仿真研究正庚烷燃料液滴在不同对流速度下的火焰形态及燃烧.结果表明:运动液滴内部环流使液滴内部低温区向环流中心移动.当液滴运动速度大于某临界值后,火焰形态由包覆火焰转变为尾迹火焰.包覆火焰的富燃区范围、高温区范围及燃烧速率明显较尾迹火焰大;包覆火焰的液滴表面温度及表面蒸发流率分布也明显不同于尾迹火焰.     

On the thermo-acoustic Fant equation

A ‘reduced complexity’ equation is derived to investigate combustion instabilities of a Rijke burner. The equation is nonlinear and furnishes limit cycle solutions for finite amplitude burner modes. It is a generalisation to combustion flows of the Fant equation used to investigate the production of voiced speech by unsteady throttling of flow by the vocal folds [G. Fant, Acoustic Theory of Speech Production. Mouton, The Hague, 1960]. In the thermo-acoustic problem the throttling occurs at the flame holder. The Fant equation governs the unsteady volume flow past the flame holder which, in turn, determines the acoustics of the entire system. The equation includes a fully determinate part that depends on the geometry of the flame holder and the thermo-acoustic system, and terms defined by integrals involving thermo-aerodynamic sources, such as a flame and vortex sound sources. These integrals provide a clear indication of what must be known about the flow to obtain a proper understanding of the dynamics of the thermo-acoustic system. Illustrative numerical results are presented for the linearised equation. This governs the growth rates of the natural acoustic modes, determined by system geometry, boundary conditions and mean temperature distribution, which are excited into instability by unsteady heat release from the flame and damped by large scale vorticity production and radiation losses into the environment. In addition, the equation supplies information about the ‘combustion modes’ excited by the local time-delay feedback dynamics of the flame.     

A generalized flame surface density modelling approach for the auto-ignition of a turbulent non-premixed system

Auto-ignition of turbulent non-premixed systems is encountered in practical devices such as diesel internal combustion engines. It remains a challenge for modellers, as it exhibits specific features such as unsteadiness, flame propagation and combustion far from stoichiometric conditions. In this paper, a two-dimensional DNS database of an igniting H₂/O₂/N₂ mixing layer, including detailed chemistry and transport, is extensively post-processed in order to gain physical insight into the flame structure and dynamics during auto-ignition. The results are used as a framework for the development of a generalized flame surface density modelling approach by integrating the equations over all possible mixture fraction values. The mean reaction rate is split into two contributions: a generalized flame surface density and a mean reaction rate per unit generalized flame surface density. The unsteadiness of the ignition phenomenon is accounted for via a generalized progress variable. Closures for the generalized surface average of the reaction rate and for the generalized progress variable are proposed, and the modelling approach is tested a priori versus the DNS data. The use of a laminar database for the chemistry coupled to the mean turbulent field via the generalized progress variable shows very promising results, capturing the correct ignition delay and the premixed peak in the turbulent mean heat release rate evolution. This allows confidence in future inclusion and validation of this approach in a RANS-CFD code.     

Nonlinear hydrodynamic and thermoacoustic oscillations of a bluff-body stabilised turbulent premixed flame

Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.     

An analytical study of droplet combustion under microgravity: Quasi-steady transient approach

An analytical model based on an assumption of combined quasi-steady and transient behavior of the process is presented to exemplify the unsteady, sphero-symmetric single droplet combustion under microgravity. The model used in the present study includes an alternative approach of describing the droplet combustion as a process where the diffusion of fuel vapor residing inside the region between the droplet surface and the flame interface experiences quasi-steadiness while the diffusion of oxidizer inside the region between the flame interface and the ambient surrounding experiences unsteadiness. The modeling approach especially focuses on predicting; the variations of droplet and flame diameters with burning time, the effect of vaporization enthalpy on burning behavior, the average burning rates and the effect of change in ambient oxygen concentration on flame structure. The modeling results are compared with a wide range of experimental data available in the literature. It is shown that this simplified quasi-steady transient approach towards droplet combustion yields behavior similar to the classical droplet theory.     

Large eddy simulation of piloting effects on turbulent swirling flames

Pilot flames, created by additional injectors of pure fuel, are often used in turbulent burners to enhance flame stabilization and reduce combustion instabilities. The exact mechanisms through which these additional rich zones modify the flame anchoring location and the combustion dynamics are often difficult to identify, especially when they include unsteady hydrodynamic motion. This study presents Large Eddy Simulations (LES) of the reacting flow within a large-scale gas turbine burner for two different cases of piloting, where either 2 or 6% of the total methane used in the burner is injected through additional pilot flame lines. For each case, LES shows how the pilot fuel injection affects both flame stabilization and flame stability. The 6% case leads to a stable flame and limited hydrodynamic perturbations in the initial flame zone. The 2% case is less stable, with a small-lift-off of the flame and a Precessing Vortex Core (PVC) in the cold stabilization zone. This PVC traps some of the lean cold gases issuing from the pilot passage stream, changes the flame stabilization point and induces instability.     

Numerical investigation of laminar cross-flow non-premixed flames in the presence of a bluff-body

A knowledge of flame stability regimes in the presence of cylindrical bluff-bodies of various dimensions is essential to design non-premixed burners. The reacting flow field in such cases is reported to be three-dimensional and unsteady. In the literature, only a few experimental investigations with limited measurements are available. Therefore, in this work, a detailed numerical study of laminar cross-flow non-premixed methane–air flames in the presence of a square cylinder is presented. The flow, temperature, species and reaction fields have been predicted using a comprehensive transient three-dimensional reacting flow model with detailed chemical kinetics and variable thermo-physical properties, in order to get a good insight into the flame stabilisation phenomena. Further, analyses of quantities such as local equivalence ratio, cell Damköhler number, species velocity, net consumption rate of methane, which are not easily obtained through experiments even with detailed diagnostics, have been carried out. The influence of the flow field due to varying inlet velocity of the oxidiser, in the presence of the bluff-body, on flame anchoring location has been analysed in detail. Local equivalence ratio contours obtained from non-reacting flow calculations are seen to be quite useful in analysing the mixing process and in the prediction of flame anchoring locations when the flames are not separated. Cell Damköhler number has been calculated using cell size, species velocity of the fuel, which is a derived quantity, and the net reaction rate of the fuel. The flame zone, which is customarily inferred from the contours of temperature, CO and OH, is also shown to be predicted well by the contour line corresponding to a Damköhler number equal to unity. The net reaction rate of CH₄ and the net rates of two dominant reactions, which consume methane, show clearly the variation in the flame anchoring locations in these three cases. Further, the three-dimensionality of these flames are analysed by plotting the mean temperature contours in y–z planes. Finally, the unsteadiness in the separated flame case is analysed.     

Nonlinear hydrodynamic phenomena in Stokes flow regime

We investigate nonlinear phenomena in dispersed two-phase systems under creeping-flow conditions. We consider nonlinear evolution of a single deformed drop and collective dynamics of arrays of hydrodynamically coupled particles. To explore physical mechanisms of system instabilities, chaotic drop evolution, and structural transitions in particle arrays we use simple models, such as small-deformation equations and effective-medium theory. We find numerical and analytical solutions of the simplified governing equations. The small-deformation equations for drop dynamics are analyzed using results of dynamical systems theory. Our investigations shed new light on the dynamics of complex fluids, where the nonlinearity often stems from the evolving boundary conditions in Stokes flow.     

Controlled detonation initiation in hypersonic flow

Experimental evidence of controlled detonation initiation and propagation in a hypersonic flow of premixed hydrogen-air is presented. This controlled detonation initiation is created in a hypersonic facility capable of producing a Mach 5 flow of hydrogen-air. Flow diagnostics such as high-speed schlieren and OH* chemiluminescence results show that a flame deflagration-to-detonation transition occurs as a combined result of turbulent flame acceleration and shock-focusing. The experimental results define three new distinct regimes in a Mach 5 premixed flow: deflagration-to-detonation transition (DDT), unsteady compressible turbulent flames, and shock-induced combustion. A two-dimensional implicit-LES (ILES) simulation, which solves the compressible, reactive Navier-Stokes equations on an adapting grid is conducted to provide additional insight into the local physical mechanism of detonation transition and propagation.     

Ethanol turbulent spray flame response to gas velocity modulation

A numerical investigation of the interaction between a spray flame and an acoustic forcing of the velocity field is presented in this paper. In combustion systems, a thermoacoustic instability is the result of a process of coupling between oscillations in heat released and acoustic waves. When liquid fuels are used, the atomisation and the evaporation process also undergo the effects of such instabilities, and the computational fluid dynamics of these complex phenomena becomes a challenging task. In this paper, an acoustic perturbation is applied to the mass flow of the gas phase at the inlet and its effect on the evaporating fuel spray and on the flame front is investigated with unsteady Reynolds averaged Navier-Stokes numerical simulations. Two flames are simulated: a partially premixed ethanol/air spray flame and a premixed pre-vaporised ethanol/air flame, with and without acoustic forcing. The frequencies used to perturb the flames are 200 and 2500 Hz, which are representative for two different regimes. Those regimes are classified based on the Strouhal number St = (D/U)f_f: at 200 Hz, St = 0.07, and at 2500 Hz, St = 0.8. The exposure of the flame to a 200 Hz signal results in a stretching of the flame which causes gas field fluctuations, a delay of the evaporation and an increase of the reaction rate. The coupling between the flame and the flow excitation is such that the flame breaks up periodically. At 2500 Hz, the evaporation rate increases but the response of the gas field is weak and the flame is more stable. The presence of droplets does not play a crucial role at 2500 Hz, as shown by a comparison of the discrete flame function in the case of spray and pre-vaporised flame. At low Strouhal number, the forced response of the pre-vaporised flame is much higher compared to that of the spray flame.     

Consequence analysis of blast wave from accidental gas explosions

Recently, consequence analyses of accidental gas explosions are often carried out to assess the risk of chemical plants, hazardous-materials sites and new energy systems. In these consequence analyses, it is indispensable to adequately predict the blast-wave (pressure-wave) intensity from gas deflagrations. Some prediction models already exist; however, most of them are based on the theory for explosives and adjusting parameters are needed for evaluating gas deflagrations. In this study, new prediction methods for gas deflagrations were developed. From theoretical analysis of blast-wave generation by a gas deflagration, an evaluation equation of the blast-wave intensity was derived. As the scale of gas deflagration becomes larger, flame front instability (especially hydrodynamic instability) would be more effective and the flame propagating velocity starts to be accelerated. Therefore, the equation was modified considering the effect of flame instability. The evaluations by this modified equation agreed well with the results of large scale experiments. By this analysis, it was found that not only total energy release but also combustion reaction rate has to be introduced into the prediction of gas deflagrations. Using this concept, a modified scale model to predict the blast-wave intensity was developed by improving the previous scale model introducing the term of combustion reaction rate as burning velocity. Furthermore, scale analysis was performed to develop the new scaling law. The universal relationship between scaled distance and overpressure has been realized by this new scaling law for gas deflagrations. In summary, these results provide new methods for accurate prediction of the blast-wave intensity from gas deflagrations.     

Flame dynamics and consideration of deflagration-to-detonation transition in central gravitational field

Expanding reaction fronts are central to many terrestrial processes as well as such cosmic phenomena as the thermonuclear combustion in supernovae. While in terrestrial conditions the effect of intrinsic flamefront instabilities is generally believed to be supplementary to, say, external turbulence and chamber dynamics, at the astrophysical scale the role of flamefront instabilities in the flame acceleration is presumably dominant. Moreover, while in terrestrial systems we focus mainly on the hydrodynamic, Darrieus–Landau (DL) instability, the Rayleigh–Taylor (RT), body-force instability could be a key issue for astrophysical flames because of the enormous gravity and the associated large scales. Consequently, if RT effects dominate over DL effects, the globally-spherical flamefront can be replaced by an expanding bubble with a non-spherical shape, possibly that of digit “8”. In the present work we have developed a self-similar formulation describing a globally-spherical expanding flamefront corrugated due to the DL instability in a central gravitation field. The associated scenario of the flame acceleration, the evolution of the upstream flow, and the instant when a gas parcel ahead of the flamefront first explodes, hence initiating the deflagration-to-detonation transition (DDT), as well as the radial coordinate of this parcel, are determined. We have also compared the effects of DL and RT instabilities, estimating whether a globally-spherical, DL-corrugated flamefront is subsequently terminated by the formation of RT bubbles. It is shown how the instant of such a transition and the relevant global flame radius are coupled to various flame and flow parameters.     

Stability and dynamical properties of material flow systems on random networks

The theory of complex networks and of disordered systems is used to study the stability and dynamical properties of a simple model of material flow networks defined on random graphs. In particular we address instabilities that are characteristic of flow networks in economic, ecological and biological systems. Based on results from random matrix theory, we work out the phase diagram of such systems defined on extensively connected random graphs, and study in detail how the choice of control policies and the network structure affects stability. We also present results for more complex topologies of the underlying graph, focussing on finitely connected Erdös-Réyni graphs, Small-World Networks and Barabási-Albert scale-free networks. Results indicate that variability of input-output matrix elements, and random structures of the underlying graph tend to make the system less stable, while fast price dynamics or strong responsiveness to stock accumulation promote stability.     

An efficient approach of unsteady flamelet modeling of a cross-flow-jet combustion system using LES

Steady flamelet models have been widely used in turbulent combustion simulations because of their simplicity, efficiency, yet physics-based nature. They are, however, unable to handle slow chemical and physical processes such as pollutant formation. Unsteady flamelet models have been shown to be able to provide accurate predictions especially for pollutants, but their implementations are usually not as straightforward as for the steady models, and additional assumptions are involved. One relatively straightforward approach of implementing the unsteady flamelet model is to tabulate the time history of unsteady flamelet solutions. This often leads to flamelet libraries of large sizes because of increased dimensions for the new physics. The purpose of this paper is to introduce a new and efficient approach of tabulating unsteady flamelet solutions in the LES of complex systems, here demonstrated in simulations of a cross-flow-jet combustion system. This approach employs Taylor series expansions to represent the time history of unsteady flamelet solutions. Compared with other approaches, the new approach retains the efficiency and simplicity benefits of steady flamelet models but possesses the accuracy of unsteady flamelet models. Various issues associated with the formulation and implementation of this approach are discussed, which include the selection of the base solution, the order of accuracy of the expansion, and the treatment of simultaneous wall heat losses and heat transfer through thermal radiation. This approach is validated in large eddy simulations of a cross-flow-jet combustion system. Good agreement with experiments is obtained for both temperature and NO concentration, as well as for major species.     

Study of vortex core precession in combustion chambers

The article presents the results of experimental investigation of swirling flow of lean propane/air flame in a model combustion chamber at atmospheric pressure. To study the unsteady turbulent flow, the particle image velocimetry technique was used. It was concluded that dynamics of high swirl flows with and without combustion was determined by a global helical mode, complying with a precessing double-spiral coherent vortex structure. The studied low swirl flame had similar size and stability characteristics, but amplitude of the coherent helical structure substantially oscillated in time. The oscillations were associated with intermittently appearing central recirculation zone that was absent in the nonreacting flow. It is expected that the low swirl flow without the permanent central recirculation zone should be more sensitive to an external active control. In particular, this result may be useful for suppression of thermoacoustic resonance in combustion chambers.     

A physical relationship between consumption and displacement speed for premixed flames with finite thickness

The characterization of premixed flames by a flame speed has been a subject that has occupied much interest in the literature in many systematic studies on combustion phenomena. Consumption and displacement speeds are two such flame speeds that are understood to describe the flame dynamics under the effect of flame curvature, flow non-uniformities, Lewis number and turbulence effects along with heat transfer with flame holders and cold walls. As such, much work has been done in the past where either one of these two speeds has been employed along with a linear sensitivity coefficient (Markstein length) for describing different sensitivities to stretch effects. However, despite recent attempts using the asymptotic theory, the relationship between these two quantities has only been clarified in a limited manner for flames of finite thickness. In this study, we use flame stretch theory that takes into account changes of stretch, curvature, heat transfer and Lewis number effects throughout the pre-heat zone and its integral effect on the flame reaction zone. A sound mathematical and physical basis is provided for understanding the two speeds that is valid for weak as well as strong stretch effects. Understanding from theory is further demonstrated by analysing several example 1D stretched flames along with a 2D bluff body flame near extinction.     

Large eddy simulation of turbulent combustion systems

This paper reviews recent and ongoing work on numerical models for turbulent combustion systems based on a classical LES approach. The work is confined to single-phase reacting flows. First, important physico-chemical features of combustion-LES are discussed along with several aspects of overall LES models. Subsequently, some numerical issues, in particular questions associated with the reliability of LES results, are outlined. The details of chemistry, its reduction, and tabulation are not addressed here. Second, two illustrative applications dealing with non-premixed and premixed flame configurations are presented. The results show that combustion-LES is able to provide predictions very close to measured data for configurations where the flow is governed by large turbulent structures. To meet the future demands, new key challenges in specific modelling areas are suggested, and opportunities for advancements in combustion-LES techniques are highlighted. From a predictive point of view, the main target must be to provide a reliable method to aid combustion safety studies and the design of combustion systems of practical importance.     

Impact of pressure fluctuations on the dynamics of laminar premixed flames

Thermo-acoustic instabilities are problematic in the design of continuous-combustion propulsion systems such as gas turbine engines, rocket motors, jet engine afterburners, and ramjets. Conceptually, the coupling between acoustics and flame dynamics can be divided into two categories: flame area fluctuations and changes in the local flame speed. The latter can be caused by the thermodynamic fluctuations that accompany an acoustic wave. This coupling is the focus of the present work. In this paper, we are concerned with the dynamics of laminar premixed flames involving large hydrocarbon species. Through high-fidelity numerical simulations, we investigate the flame response for a wide range of fuels and acoustic frequencies. The combustion of hydrogen and methane is considered for verification purposes and as baseline cases for comparison with two large hydrocarbon fuels, n-heptane and n-dodecane. We extract the phase and gain of the unsteady heat release response, which are directly related to the Rayleigh criterion and thus the stability of the system. For all fuels, we observe a local peak in the heat release gain. At high frequencies, we find that the fluctuations of the different species mass fractions decrease with the inverse of the acoustic frequency, leading to chemistry being “frozen” in the high-frequency limit. This allows us to predict the flame behavior directly from the steady-state solution.     

A theoretical study on the unsteady aerothermodynamics for attached flow models

The principle of the unsteady aerothermodynamics was theoretically investigated for the attached flow. Firstly, two simplified models with analytic solutions to the N-S equations were selected for the research, namely the compressible unsteady flows on the infinite flat plate with both time-varying wall velocity and time-varying wall temperature boundary conditions. The unsteady temperature field and the unsteady wall heat flux (heat flow) were analytically solved for the second model. Then, the interaction characteristic of the unsteady temperature field and the unsteady velocity field in the simplified models and the effects of the interaction on the transient wall heat transfer were studied by these two analytic solutions. The unsteady heat flux, which is governed by the energy equation, is directly related to the unsteady compression work and viscous dissipation which originates from the velocity field governed by the momentum equation. The main parameters and their roles in how the unsteady velocity field affects the unsteady heat flux were discussed for the simplified models. Lastly, the similarity criteria of the unsteady aerothermodynamics were derived based on the compressible boundary layer equations. Along with the Strouhal number Stu, the unsteadiness criterion of the velocity field, StT number, the unsteadiness criterion of the temperature field was proposed for the first time. Different from the traditional method used in unsteady aerodynamics which measures the flow unsteadiness only by the Stu number, present results show that the flow unsteadiness in unsteady aerothermodynamics should be comprehensively estimated by comparing the relative magnitudes of the temperature field unsteadiness criterion StT number with the coefficients of other terms in the dimensionless energy equation.     

Axis switching in turbulent buoyant diffusion flames

Dynamics of buoyant diffusion flames from rectangular, square, and round fuel sources were investigated using direct numerical simulation (DNS). Fully three-dimensional simulations were performed employing high-order numerical methods and boundary conditions to solve governing equations for variable-density flow and finite-rate Arrhenius chemistry. Significant differences among the different cases were revealed in the vortex dynamics, entrainment rate, small-scale mixing, and consequently flame structures. Mixing and entrainment enhancement in non-circular flames in comparison with circular ones was explained using the Biot–Savart instability theory, which relates vortex dynamics to the local azimuthal curvature. An extension of the theory elucidated why rectangular flames entrain more efficiently and spread wider than square ones, although both configurations have corners. It also provided an explanation for the aspect ratio effects in the near field. In the far field, nonlinear effects were dominant and the general transport equations for vorticity were analyzed in detail. The corner effects and aspect ratio effects were shown to be augmented by the intricate interactions among vortex dynamics, combustion, and buoyancy through the various terms in the equations. The presence of corners in non-circular flames led to concentrated regions of fine-scale mixing and intense reactions centered around the corners. Moreover, the rectangular flames exhibited a different dynamic behavior from even the square one, by creating discrepancies in entrainment, mixing, and combustion between the minor and major axis directions. Increasing the aspect ratio exacerbated such directional discrepancies, and ultimately led to axis switching. It was the first time that axis switching was observed by DNS in a rectangular flame of aspect ratio 3, which raised further questions in combustion prediction and control. Finally, a unified explanation for corner and aspect ratio effects was given on the basis of the Biot–Savart instability theory and the vorticity transport equations.