Numerical simulation of the upward propagation of a flame in a vertical tube filled with a very lean mixture.

Upwardpropagation of a premixed flame in averticaltubefilled with a very leanmixture is simulated numerically using a single irreversible Arrhenius reaction model with infinitely high activation energy. In the absence of heat losses and preferential diffusion effects, a curved flame with stationary shape and velocity close to those of an open bubble ascending in the same tube is found for values of the fuel mass fraction above a certain minimum that increases with the radius of the tube, while the numerical computations cease to converge to a stationary solution below this minimum mass fraction. The vortical flow of the gas behind the flame and in its transport region is described for tubes of different radii. It is argued that this flow may become unstable when the fuel mass fraction is decreased, and that this instability, together with the flame stretch due to the strong curvature of the flame tip in narrow tubes, may be responsible for the minimum fuel mass fraction. Radiation losses and a Lewis number of the fuel slightly above unity decrease the final combustion temperature at the flame tip and increase the minimum fuel mass fraction, while a Lewis number slightly below unity has the opposite effect.

Advances in global instability computations: from incompressible to hypersonic flow

Esta tesis constituye un gran avance en el conocimiento del estudio y análisis de inestabilidades hidrodinámicas desde un punto de vista físico y teórico, como consecuencia de haber desarrollado innovadoras técnicas para la resolución computacional eficiente y precisa de la parte principal del espectro correspondiente a los problemas de autovalores (EVP) multidimensionales que gobiernan la inestabilidad de flujos con dos o tres direcciones espaciales inhomogéneas, denominados problemas de estabilidad global lineal. En el contexto del trabajo de desarrollo de herramientas computacionales presentado en la tesis, la discretización mediante métodos de diferencias finitas estables de alto orden de los EVP bidimensionales y tridimensionales que se derivan de las ecuaciones de Navier-Stokes linealizadas sobre flujos con dos o tres direcciones espaciales inhomogéneas, ha permitido una aceleración de cuatro órdenes de magnitud en su resolución. Esta mejora de eficiencia numérica se ha conseguido gracias al hecho de que usando estos esquemas de diferencias finitas, técnicas eficientes de resolución de problemas lineales son utilizables, explotando el alto nivel de dispersión o alto número de elementos nulos en las matrices involucradas en los problemas tratados. Como más notable consecuencia cabe destacar que la resolución de EVPs multidimensionales de inestabilidad global, que hasta la fecha necesitaban de superordenadores, se ha podido realizar en ordenadores de sobremesa. Además de la solución de problemas de estabilidad global lineal, el mencionado desarrollo numérico facilitó la extensión de las ecuaciones de estabilidad parabolizadas (PSE) lineales y no lineales para analizar la inestabilidad de flujos que dependen fuertemente en dos direcciones espaciales y suavemente en la tercera con las ecuaciones de estabilidad parabolizadas tridimensionales (PSE-3D). Precisamente la capacidad de extensión del novedoso algoritmo PSE-3D para el estudio de interacciones no lineales de los modos de estabilidad, desarrollado íntegramente en esta tesis, permite la predicción de transición en flujos complejos de gran interés industrial y por lo tanto extiende el concepto clásico de PSE, el cuál ha sido empleado exitosamente durante las pasadas tres décadas en el mismo contexto para problemas de capa límite bidimensional. Típicos ejemplos de flujos incompresibles se han analizado en este trabajo sin la necesidad de recurrir a restrictivas presuposiciones usadas en el pasado. Se han estudiado problemas vorticales como es el caso de un vórtice aislado o sistemas de vórtices simulando la estela de alas, en los que la homogeneidad axial no se impone y así se puede considerar la difusión viscosa del flujo. Además, se ha estudiado el chorro giratorio turbulento, cuya inestabilidad se utiliza para mejorar las características de funcionamiento de combustores. En la tesis se abarcan adicionalmente problemas de flujos compresibles. Se presenta el estudio de inestabilidad de flujos de borde de ataque a diferentes velocidades de vuelo. También se analiza la estela formada por un elemento rugoso aislado en capa límite supersónica e hipersónica, mostrando excelentes comparaciones con resultados obtenidos mediante simulación numérica directa. Finalmente, nuevas inestabilidades se han identificado en el flujo hipersónico a Mach 7 alrededor de un cono elíptico que modela el vehículo de pruebas en vuelo HIFiRE-5. Los resultados comparan favorablemente con experimentos en vuelo, lo que subraya aún más el potencial de las metodologías de análisis de estabilidad desarrolladas en esta tesis. ABSTRACT The present thesis constitutes a step forward in advancing the frontiers of knowledge of fluid flow instability from a physical point of view, as a consequence of having been successful in developing groundbreaking methodologies for the efficient and accurate computation of the leading part of the spectrum pertinent to multi-dimensional eigenvalue problems (EVP) governing instability of flows with two or three inhomogeneous spatial directions. In the context of the numerical work presented in this thesis, the discretization of the spatial operator resulting from linearization of the Navier-Stokes equations around flows with two or three inhomogeneous spatial directions by variable-high-order stable finite-difference methods has permitted a speedup of four orders of magnitude in the solution of the corresponding two- and three-dimensional EVPs. This improvement of numerical performance has been achieved thanks to the high-sparsity level offered by the high-order finite-difference schemes employed for the discretization of the operators. This permitted use of efficient sparse linear algebra techniques without sacrificing accuracy and, consequently, solutions being obtained on typical workstations, as opposed to the previously employed supercomputers. Besides solution of the two- and three-dimensional EVPs of global linear instability, this development paved the way for the extension of the (linear and nonlinear) Parabolized Stability Equations (PSE) to analyze instability of flows which depend in a strongly-coupled inhomogeneous manner on two spatial directions and weakly on the third. Precisely the extensibility of the novel PSE-3D algorithm developed in the framework of the present thesis to study nonlinear flow instability permits transition prediction in flows of industrial interest, thus extending the classic PSE concept which has been successfully employed in the same context to boundary-layer type of flows over the last three decades. Typical examples of incompressible flows, the instability of which was analyzed in the present thesis without the need to resort to the restrictive assumptions used in the past, range from isolated vortices, and systems thereof, in which axial homogeneity is relaxed to consider viscous diffusion, as well as turbulent swirling jets, the instability of which is exploited in order to improve flame-holding properties of combustors. The instability of compressible subsonic and supersonic leading edge flows has been solved, and the wake of an isolated roughness element in a supersonic and hypersonic boundary-layer has also been analyzed with respect to its instability: excellent agreement with direct numerical simulation results has been obtained in all cases. Finally, instability analysis of Mach number 7 ow around an elliptic cone modeling the HIFiRE-5 flight test vehicle has unraveled flow instabilities near the minor-axis centerline, results comparing favorably with flight test predictions.

Influence of Strouhal number on pulsating methane–air coflow jet diffusion flames

Four periodically time-varying methane–air laminar coflow jet diffusion flames, each forced by pulsating the fuel jet's exit velocity Uj sinusoidally with a different modulation frequency wj and with a 50% amplitude variation, have been computed. Combustion of methane has been modeled by using a chemical mechanism with 15 species and 42 reactions, and the solution of the unsteady Navier–Stokes equations has been obtained numerically by using a modified vorticity-velocity formulation in the limit of low Mach number. The effect of wj on temperature and chemistry has been studied in detail. Three different regimes are found depending on the flame's Strouhal number S=awj/Uj, with a denoting the fuel jet radius. For small Strouhal number (S=0.1), the modulation introduces a perturbation that travels very far downstream, and certain variables oscillate at the frequency imposed by the fuel jet modulation. As the Strouhal number grows, the nondimensional frequency approaches the natural frequency of oscillation of the flickering flame (S≃0.2). A coupling with the pulsation frequency enhances the effect of the imposed modulation and a vigorous pinch-off is observed for S=0.25 and S=0.5. Larger values of S confine the oscillation to the jet's near-exit region, and the effects of the pulsation are reduced to small wiggles in the temperature and concentration values. Temperature and species mass fractions change appreciably near the jet centerline, where variations of over 2% for the temperature and 15% and 40% for the CO and OH mass fractions, respectively, are found. Transverse to the jet movement, however, the variations almost disappear at radial distances on the order of the fuel jet radius, indicating a fast damping of the oscillation in the spanwise direction.

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.

The structure of lean hydrogen-air flame balls

An analysis of the structure of flame balls encountered under microgravity conditions, which are stable due to radiant energy losses from H₂O, is carried out for fuel-lean hydrogen-air mixtures. It is seen that, because of radiation losses, in stable flame balls the maximum flame temperature remains close to the crossover temperature, at which the rate of the branching step H + O₂ -> OH + O equals that of the recombination step H + O₂ + M -> HO₂ + M. Under those conditions, all chemical intermediates have very small concentrations and follow the steady-state approximation, while the main species react according to the overall step 2H₂ + O₂-> 2H₂O; so that a one-step chemical-kinetic description, recently derived by asymptotic analysis for near-limit fuel-lean deflagrations, can be used with excellent accuracy to describe the whole branch of stable flame balls. Besides molecular diffusion in a binary-diffusion approximation, Soret diffusion is included, since this exerts a nonnegligible effect to extend the flammability range. When the large value of the activation energy of the overall reaction is taken into account, the leading-order analysis in the reaction-sheet approximation is seen to determine the flame ball radius as that required for radiant heat losses to remove enough of the heat released by chemical reaction at the flame to keep the flame temperature at a value close to crossover. The results are relevant to burning velocities at lean equivalent ratios and may influence fire-safety issues associated with hydrogen utilization.

Flammability conditions for ultra-lean hydrogen premixed combustion based on flame-ball analyses

It has been reasoned that the structures of strongly cellular flames in very lean mixtures approach an array of flame balls, each burning as if it were isolated, thereby indicating a connection between the critical conditions required for existence of steady flame balls and those necessary for occurrence of self-sustained premixed combustion. This is the starting assumption of the present study, in which structures of near-limit steady sphericosym-metrical flame balls are investigated with the objective of providing analytic expressions for critical combustion conditions in ultra-lean hydrogen-oxygen mixtures diluted with N2 and water vapor. If attention were restricted to planar premixed flames, then the lean-limit mole fraction of H2 would be found to be roughly ten percent, more than twice the observed flammability limits, thereby emphasizing the relevance of the flame-ball phenomena. Numerical integrations using detailed models for chemistry and radiation show that a onestep chemical-kinetic reduced mechanism based on steady-state assumptions for all chemical intermediates, together with a simple, optically thin approximation for water-vapor radiation, can be used to compute near-limit fuel-lean flame balls with excellent accuracy. The previously developed one-step reaction rate includes a crossover temperature that determines in the first approximation a chemical-kinetic lean limit below which combustión cannot occur, with critical conditions achieved when the diffusion-controlled radiation-free peak temperature, computed with account taken of hydrogen Soret diffusion, is equal to the crossover temperature. First-order corrections are found by activation-energy asymptotics in a solution that involves a near-field radiation-free zone surrounding a spherical flame sheet, together with a far-field radiation-conduction balance for the temperature profile. Different scalings are found depending on whether or not the surrounding atmosphere contains wáter vapor, leading to different analytic expressions for the critical conditions for flame-ball existence, which give results in very good agreement with those obtained by detailed numerical computations.

Flamelet structures in spray ignition

In typical liquid-fueled burners the fuel is injected as a high-velocity liquid jet that breaks up to form the spray. The initial heating and vaporization of the liquid fuel rely on the relatively large temperatures of the sourrounding gas, which may include hot combustion products and preheated air. The heat exchange between the liquid and the gas phases is enhanced by droplet dispersion arising from the turbulent motion. Chemical reaction takes place once molecular mixing between the fuel vapor and the oxidizer has occurred in mixing layers separating the spray flow from the hot air stream. Since in most applications the injection velocities are much larger than the premixed-flame propagation velocity, combustion stabilization relies on autoignition of the fuel-oxygen mixture, with the combustion stand-off distance being controlled by the interaction of turbulent transport, droplet heating and vaporization, and gas-phase chemical reactions. In this study, conditions are identified under which analyses of laminar flamelets canshed light on aspects of turbulent spray ignition. This study extends earlier fundamental work by Liñan & Crespo (1976) on ignition in gaseous mixing layers to ignition of sprays. Studies of laminar mixing layers have been found to be instrumental in developing un-derstanding of turbulent combustion (Peters 2000), including the ignition of turbulent gaseous diffusion flames (Mastorakos 2009). For the spray problem at hand, the configuration selected, shown in Figure 1, involves a coflow mixing layer formed between a stream of hot air moving at velocity UA and a monodisperse spray moving at velocity USUA. The boundary-layer approximation will be used below to describe the resulting sl ender flow, which exhibits different igniting behaviors depending on the characteristics of t he fuel. In this approximation, consideration of the case U A = U S enables laminar ignition distances to be related to ignition times of unstrained spray flamelets, thereby pro viding quantitative information of direct applicability in regions of low scala r dissipation-rate in turbulent reactive flows (see the discussion in pp. 181–186 of Peters (2000)) . This report is organized as follows. Effects of droplet dispersion dynamics on ignition of sprays in turbulent mixing layers are discussed in Section 2. The formulation f or ignition in laminar mixing layers is outlined in Sections 3 and 4. The results are presented in Section 5. In Section 6, the mixture-fraction field and associated scalar dissipat ion rates for spray ignition are discussed. Finally, some brief conclusions are drawn in Section 7.

Generalized flows, intrinsic stochasticity, and turbulent transport

The study of passive scalar transport in a turbulent velocity field leads naturally to the notion of generalized flows, which are families of probability distributions on the space of solutions to the associated ordinary differential equations which no longer satisfy the uniqueness theorem for ordinary differential equations. Two most natural regularizations of this problem, namely the regularization via adding small molecular diffusion and the regularization via smoothing out the velocity field, are considered. White-in-time random velocity fields are used as an example to examine the variety of phenomena that take place when the velocity field is not spatially regular. Three different regimes, characterized by their degrees of compressibility, are isolated in the parameter space. In the regime of intermediate compressibility, the two different regularizations give rise to two different scaling behaviors for the structure functions of the passive scalar. Physically, this means that the scaling depends on Prandtl number. In the other two regimes, the two different regularizations give rise to the same generalized flows even though the sense of convergence can be very different. The “one force, one solution” principle is established for the scalar field in the weakly compressible regime, and for the difference of the scalar in the strongly compressible regime, which is the regime of inverse cascade. Existence and uniqueness of an invariant measure are also proved in these regimes when the transport equation is suitably forced. Finally incomplete self similarity in the sense of Barenblatt and Chorin is established.

On the inverse parabolicity of pdf equations in turbulent flows

The focus of the present work is the well-known feature of the probability density function (PDF) transport equations in turbulent flows-the inverse parabolicity of the equations. While it is quite common in fluid mechanics to interpret equations with direct (forward-time) parabolicity as diffusive (or as a combination of diffusion, convection and reaction), the possibility of a similar interpretation for equations with inverse parabolicity is not clear. According to Einstein's point of view, a diffusion process is associated with the random walk of some physical or imaginary particles, which can be modelled by a Markov diffusion process. In the present paper it is shown that the Markov diffusion process directly associated with the PDF equation represents a reasonable model for dealing with the PDFs of scalars but it significantly underestimates the diffusion rate required to simulate turbulent dispersion when the velocity components are considered.

Differential diffusion: Often a finite-mising length effect

Many instances of differential diffusion, i e, different species having different turbulent diffusion coefficients in the same flow, can be explained as a finite mixing length effect. That is, in a simple mixing length scenario, the turbulent diffusion coefficient has the form 1 ( m )2 m m c l K w l OL   =  +    where, wm is the mixing velocity, lm the mixing length and Lc the overall distribution scale for a particular species. The first term represents the familiar gradient diffusion while the second term becomes important when lm/Lc is finite. This second term shows that different species will have different diffusion coefficients if they have different overall distribution scales. Such different Lcs may come about due to different boundary conditions and different intrinsic properties (molecular diffusivity, settling velocity etc) for different species. For momentum transfer in turbulent oscillatory boundary layers the second term is imaginary and explains observed phase leads of shear stresses ahead of velocity gradients.

Soot measurements in steady and pulsed ethylene/air diffusion flames using Laser-Induced Incandescence

Combustion-generated carbon black nano particles, or soot, have both positive and negative effects depending on the application. From a positive point of view, it is used as a reinforcing agent in tires, black pigment in inks, and surface coatings. From a negative point of view, it affects performance and durability of many combustion systems, it is a major contributor of global warming, and it is linked to respiratory illness and cancer. Laser-Induced Incandescence (LII) was used in this study to measure soot volume fractions in four steady and twenty-eight pulsed ethylene diffusion flames burning at atmospheric pressure. A laminar coflow diffusion burner combined with a very-high-speed solenoid valve and control circuit provided unsteady flows by forcing the fuel flow with frequencies between 10 Hz and 200 Hz. Periodic flame oscillations were captured by two-dimensional phase-locked LII images and broadband luminosity images for eight phases (0° – 360°) covering each period. A comparison between the steady and pulsed flames and the effect of the pulsation frequency on soot volume fraction in the flame region and the post flame region are presented. The most significant effect of pulsing frequency was observed at 10 Hz. At this frequency, the flame with the lowest mean flow rate had 1.77 times enhancement in peak soot volume fraction and 1.2 times enhancement in total soot volume fraction; whereas the flame with the highest mean flow rate had no significant change in the peak soot volume fraction and 1.4 times reduction in the total soot volume fraction. A correlation (fvRe-1 = a + b·Str) for the total soot volume fraction in the flame region for the unsteady laminar ethylene flames was obtained for the pulsation frequency between 10 Hz and 200 Hz, and the Reynolds number between 37 and 55. The soot primary particle size in steady and unsteady flames was measured using the Time-Resolved Laser-Induced Incandescence (TIRE-LII) and the double-exponential fit method. At maximum frequency (200 Hz), the soot particles were smaller in size by 15% compared to the steady case in the flame with the highest mean flow rate.

Soot Measurements in Steady and Pulsed Ethylene/Air Diffusion Flames Using Laser-Induced Incandescence

Combustion-generated carbon black nano particles, or soot, have both positive and negative effects depending on the application. From a positive point of view, it is used as a reinforcing agent in tires, black pigment in inks, and surface coatings. From a negative point of view, it affects performance and durability of many combustion systems, it is a major contributor of global warming, and it is linked to respiratory illness and cancer. Laser-Induced Incandescence (LII) was used in this study to measure soot volume fractions in four steady and twenty-eight pulsed ethylene diffusion flames burning at atmospheric pressure. A laminar coflow diffusion burner combined with a very-high-speed solenoid valve and control circuit provided unsteady flows by forcing the fuel flow with frequencies between 10 Hz and 200 Hz. Periodic flame oscillations were captured by two-dimensional phase-locked LII images and broadband luminosity images for eight phases (0°- 360°) covering each period. A comparison between the steady and pulsed flames and the effect of the pulsation frequency on soot volume fraction in the flame region and the post flame region are presented. The most significant effect of pulsing frequency was observed at 10 Hz. At this frequency, the flame with the lowest mean flow rate had 1.77 times enhancement in peak soot volume fraction and 1.2 times enhancement in total soot volume fraction; whereas the flame with the highest mean flow rate had no significant change in the peak soot volume fraction and 1.4 times reduction in the total soot volume fraction. A correlation (ƒv Reˉ1 = a+b· Str) for the total soot volume fraction in the flame region for the unsteady laminar ethylene flames was obtained for the pulsation frequency between 10 Hz and 200 Hz, and the Reynolds number between 37 and 55. The soot primary particle size in steady and unsteady flames was measured using the Time-Resolved Laser-Induced Incandescence (TIRE-LII) and the double-exponential fit method. At maximum frequency (200 Hz), the soot particles were smaller in size by 15% compared to the steady case in the flame with the highest mean flow rate.

Asymmetric hydrogen flame in a heated micro-channel: role of Darrieus-Landau and thermal-diffusive instabilities

Present work examines numerically the asymmetric behavior of hydrogen/air flame in a micro-channel subjected to a non-uniform wall temperature distribution. A high resolution (with cell size of 25 μm × 25 μm) of two-dimensional transient Navier–Stokes simulation is conducted in the low-Mach number formulation using detailed chemistry evolving 9 chemical species and 21 elementary reactions. Firstly, effects of hydrodynamic and diffusive-thermal instabilities are studied by performing the computations for different Lewis numbers. Then, the effects of preferential diffusion of heat and mass transfer on the asymmetric behavior of the hydrogen flame are analyzed for different inlet velocities and equivalence ratios. Results show that for the flames in micro-channels, interactions between thermal diffusion and molecular diffusion play major role in evolution of a symmetric flame into an asymmetric one. Furthermore, the role of Darrieus–Landau instability found to be minor. It is also found that in symmetric flames, the Lewis number decreases behind the flame front. This is related to the curvature of flame which leads to the inclination of thermal and mass fluxes. The mass diffusion vectors point toward the walls and the thermal diffusion vectors point toward the centerline. Asymmetric flame is observed when the length of flame front is about 1.1–1.15 times of the channel width.