The oxyfuel cutting method is still widely used nowadays, even though it is not a fully autonomous process. Thisthesis presents a computational model to study ion and electron transport and current-voltage characteristics inside a methane-oxygen flame. By finding the relationship between current-voltage characteristics and critical parameters,such as standoff, fuel oxygen ratio, and flow rate, a control algorithm could be implemented into the system and make it autonomous. Star CCM+ software is used to develop preheat phase computational models by splitting the simulations into the combustion and electrochemical transport parts. Both the laminar and turbulent flows are considered. Several laboratory experiments are used to compare test data with the numerical results generated using this model. The initial and boundary conditions used in the simulation were to the extent possible similar to the experimental conditions in the laboratory experiment. In the combustion part, the general GRI3.0 mechanism plus three additional ionization reactions are applied, and the combustion part results are then used as input into the electrochemical transport part. A particular inspection line inside the domain is created to analyze the results of the electrochemical transport part. Ions, electrons number density, and current density are studied in the interval from -40V to 40V electric potential. The ions are heavier and more challenging to move than electrons. The results show that at both the torch and work surfaces, charged sheaths are formed, which cause three different regions of current-voltage relations to form in a similar manner as observed in the tests. / Master of Science / Oxyfuel cutting is essential to numerous industries, such as shipbuilding, rail, earth moving equipment, commercial building construction, etc. Tuning the process parameters and diagnosing problems with the oxyfuel process still relies on experienced operators. The main obstacles to the automation of the oxyfuel process come from the limitations of the sensing suites currently in use. Since typical sensors are highly unreliable in the harsh environment near the high-temperature flame, an alternate method is proposed to find the co-dependence between the flame's electrical characteristics and critical parameters of the oxyfuel cutting system (standoff, flow rate, F/O ratio, etc.). The relevant electrical characteristics are the electrical potential and distribution of ions and electrons. Two-dimensional models are created to analyze the combustion of methane-oxygen flame and transport of ions and electrons. The models allow the derivation of the current-voltage characteristic between the torch and work surface. Also, the way sheath phenomena of ions and electrons on the surface affect the current-voltage relationship can be analyzed from ions and electrons distribution. The electric field is added to the model by applying a constant voltage to the torch tip surface. To validate the models, a laboratory experiment with a similar geometry arrangement is used as a comparison. The models' results reveal three different regimes in the current-voltage relationship.

The utilization of lean premixed combustors has become attractive to designers of industrial gas turbines as a means of meeting strict emissions standards without compromising efficiency.  Mixing the fuel and air prior to combustion allows for lower temperature flame zones, creating the potential for drastically reduced nitrous oxide emissions.  While effective, these systems are commonly plagued by combustion driven instabilities.  These instabilities produce large pressure and heat release rate fluctuations due to a resonant interaction between the combustor acoustics and the flame.  A primary feedback mechanism responsible for driving these systems is the propagation of Fuel/Air Ratio (FAR) fluctuations into the flame zone.  These fluctuations are formed inside of the premixing chamber when fuel is injected into and mixed with an oscillating air flow. The research presented here aimed to develop new technology for premixer designs, along with an application strategy, to avoid resonant thermo-acoustic events driven by FAR fluctuations.  A passive fuel control technique was selected for investigation and implementation. The selected technique utilized fuel injections at multiple, strategically placed axial locations to target and inhibit FAR fluctuations at the dominant resonant mode of the combustor.  The goal of this research was to provide an understanding of the mixing response inside a realistic premixer geometry and investigate the effectiveness of the proposed suppression technique. The mixing response was investigated under non-reacting flow conditions using a unique modular premixer.  The premixer incorporated variable axial fuel injection locations, as well as interchangeable mixing chamber geometries.  Two different chamber designs were tested: a simple annular chamber and one incorporating an axial swirler.  The mixing response of the simple annular geometry was well characterized, and it was found that multiple injections could be effectively configured to suppress the onset of an unstable event at very lean conditions. Energy dense flame zones produced at higher equivalence ratios, however, were found to be uncontrollable using this technique. Additionally, the mixing response of the swirl geometry was difficult to predict. This was found to be the result of large spatial gradients formed in the dynamic velocity field as acoustic waves passed through the swirl vanes. / Ph. D.

The predictive design of fuel injection hardware used for active combustion control is not well established in the gas turbine industry. The primary reason for this is that the underlying mechanisms governing the flow rate authority downstream of the nozzle are not well understood. A detailed investigation of two liquid fuel flow modulation configurations is performed in this thesis: a piston and a throttle-valve configuration. The two systems were successfully built with piezoelectric actuation to drive the prime movers proportionally up to 800 Hz. Discussed in this thesis are the important constituents of the fuel injection system that affect heat release authority: the method of fuel modulation, uncoupled dynamics of several components, and the compressibility of air trapped in the fuel line. Additionally, a novel technique to model these systems by way of one-dimensional, linear transmission line acoustic models was developed to successfully characterize the principle of operation of the two systems. Through these models, insight was gained on the modes through which modulation authority was dissipated and on methods through which successful amplitude scaling would be possible. At high amplitudes, it was found that the models were able to successfully predict the actual performance reasonably well for the piston device. A proportional phase shifting controller was used to test the authority on a 40-kW rig with natural longitudinal modes. Results show that, under limited operating conditions, the sound pressure level at the limit cycle frequency was reduced by about 26 dB and the broadband energy was reduced by 23 dB. Attenuation of the fuel pulse at several combustor settings was due to fluctuating vorticity and temporal droplet distribution effects. / Master of Science

Detonation cycles are identified as an efficient alternative to the Brayton cycles used in power and propulsion applications. Rotating Detonation Engine (RDE) operating on a detonation cycle works by compressing the working fluid across a detonation wave, thereby reducing the number of compressor stages required in the thermodynamic cycle. Numerical analyses of RDEs are flexible in understanding the flow field within the RDE, however, three-dimensional analyses are expensive due to the differences in time-scale required to resolve the combustion process and flow-field. The alternate two-dimensional analyses are generally modeled with perfectly premixed fuel injection and do not capture the effects of improper mixing arising due to discrete injection of fuel and oxidizer into the chamber. To model realistic injection in a 2-D analysis, the current work uses an approach in which, a Probability Density Function (PDF) of the fuel mass fraction at the chamber inlet is extracted from a 3-D, cold-flow simulation and is used as an inlet boundary condition for fuel mass fraction in the 2-D analysis. The 2-D simulation requires only 0.4% of the CPU hours for one revolution of the detonation compared to an equivalent 3-D simulation. Using this method, a perfectly premixed RDE is comparing with a non-premixed case. The performance is found to vary between the two cases. The mean detonation velocities, time-averaged static pressure profiles are found to be similar between the two cases, while the local detonation velocities and peak pressure values vary in the non-premixed case due to local pockets fuel rich/lean mixtures. The mean detonation cell sizes are similar, but the distribution in the non-premixed case is closer due to stronger shock structures. An analytical method is used to check the effects of fuel-product stratification and heat loss from the RDE and these effects adversely affect the local detonation velocity. Overall, this method of modeling captures the complex physics in an RDE with the advantage of reduced computational cost and therefore can be used for design and diagnostic purposes. / Master of Science / The conventional Brayton cycle used in power and propulsion applications is highly optimized, at cycle and component levels. In pursuit of higher thermodynamic efficiency, detonation cycles are identified as an efficient alternative and gained increased attention in the scientific community. In a Rotating Detonation Engine (RDE), which is based on the detonation cycle, the compression of gases occurs across a shock wave. This method of achieving high compression ratios reduces the number of compressor stages required for operation. In an RDE (where combustion occurs between two coaxial cylinders), the fuel and oxidizer are injected axially into the combustion chamber where the detonation is initiated. The resultant detonation wave spins continuously in the azimuthal direction, consuming fresh fuel mixture. The combustion products expand and exhaust axially providing thrust/mechanical energy when coupled with a turbine. Numerical analyses of RDEs are flexible over experimental analysis, in terms of understanding the flow physics and the physical/chemical processes occurring within the engine. However, three-dimensional numerical analyses are computationally expansive, and therefore demanding an equivalent, efficient two-dimensional analysis. In most RDEs, fuel and oxidizer are injected from separate plenums into the chamber. This type of injection leads to inhomogeneity of the fuel-air mixture within the RDE which adversely affects the performance of the engine. The current study uses a novel method to effectively capture these physics in a 2-D numerical analysis. Furthermore, the performance of the combustor is compared between perfectly premixed injection and discrete, non-premixed injection. The method used in this work can be used for any injector design and is a powerful/efficient way to numerically analyze a Rotating Detonation Engine.

A new combustor, referred to as Stagnation Point Reverse Flow (SPRF) combustor has been developed at Georgia Tech to meet increasingly stringent emission regulations. The combustor incorporates a novel design to meet the conflicting requirements of low pollution and high efficiency in both premixed and non-premixed modes. The objective of this thesis is to perform Large Eddy Simulations (LES) on this lab-scale combustor and explain the underlying physics. To achieve this, numerical simulations are performed in both the premixed and non-premixed combustion modes. The velocity field, species field, entrainment characteristics, flame structure, emissions and mixing characteristics are then analyzed. Simulations have been carried out first for a non-reactive case and the flow features in the combustor are analyzed. Next, the simulations have been extended for the premixed reactive case by employing different sub-grid scale combustion chemistry closures - Eddy Break Up (EBU), Artificially Thickened Flame (TF) and Linear Eddy Mixing (LEM) models. Only LEMLES which is an advanced scalar approach is able to accurately predict both the velocity and species field in the combustor. The results from LEM with LES (LEMLES) using a reduced chemical mechanism have been analyzed in the premixed mode. The results showed that mass entrainment occurs along the shear layer in the combustor. The entrained mass carried products into the reactant stream and provided preheating. The product entrainment enhances the reaction rates and stabilizes the flame even at very lean conditions. These products are shown to enter into the flame through local extinction zones present on the flame surface. The flame structure is further analyzed and the combustion mode is found to be primarily in thin reaction zones. The emissions in the combustor are studied using simple global mechanisms for NOx. Computations show extremely low NOx values comparable to the measured emissions. These low emissions are shown to be primarily due to the low temperatures in the combustor. LEMLES computations are also performed with detailed chemistry to capture more accurately the flame structure. The flame in the detailed chemistry case is more sensitive to strain effects and show more extinction zones very near to the injector. LEMLES approach is also used to resolve the combustion mode in the non-premixed case. The studies indicate that mixing of fuel and air close to the injector controls the combustion process. The predictions in the near field are shown to be very sensitive to the inflow conditions. Analysis shows that fuel and air mixing occurs to lean proportions in the combustor before any burning takes place. The flame structure in the non-premixed mode is very similar to the premixed mode. Along with fuel-air mixing, the products also mix with the reactants and provide the preheating effects to stabilize the flame in the downstream region of the combustor.

The search for clean and efficient combustors is motivated by the increasingly stringent emissions regulations. New gas turbine engines are designed to operate under lean conditions with inhomogeneous reactants to ensure cleanliness and stability of the combustion. This ushers in a new mode of combustion, called the inhomogeneous reactants premixed combustion. The present study investigates the effects of inhomogeneous reactants on premixed combustion, specifically on the interactions of an initially planar flame with field of inhomogeneous reactants. Unsteady and unstrained laminar methane-air flames are studied in one- and two-dimensional simulations to investigate the effects of normally and tangentially (to the flame surface) stratified reactants. A three-dimensional DNS of turbulent inhomogeneous reactants premixed combustion is performed to extend the investigation into turbulent flames. The methaneair combustion is represented by a complex chemical reaction mechanism with 18 species and 68 steps. The flame surface density (FSD) and displacement speed S_d have been used as the framework to analyse the inhomogeneous reactants premixed flame. The flames are characterised by an isosurface of reaction progress variable. The unsteady flames are compared to the steady laminar unstrained reference case. An equivalence ratio dip is observed in all simulations and it can serve as a marker for the premixed flame. The dip is attributed to the preferential diffusion of carbon- and hydrogen- containing species. Hysteresis of S_d is observed in the unsteady and unstrained laminar flames that propagate into normally stratified reactants. Stoichiometric flames propagating into lean mixture have a larger S_d than lean flames propagating into stoichiometric mixtures. The cross-dissipation term contribution to S_d is small (~~10%) but its contribution to the hysteresis of S_d is not (~~50%). Differential propagation of the flame surface is observed in the laminar flame that propagates into tangentially stratified reactants. Stretch on the flame surface is induced by the differential propagation, which in turn increases the flame surface area.

La combustion prévaporisée prémélangée pauvre est une piste de choix pour réduire les émissions polluantes des moteurs d'avions mais peut conduire à l'apparition d'instabilités thermo-acoustiques. Afin d'améliorer la stabilité de telles flammes, l'étagement du combustible consiste à contrôler la distribution spatiale du carburant. Une telle procédure s'accompagne cependant d'une complexité accrue du système pouvant déboucher sur des phénomènes inattendus.Un brûleur à l'échelle de laboratoire alimenté par du dodécane liquide est utilisé dans cette thèse. Le combustible est injecté dans deux étages séparés, permettant ainsi de contrôler sa distribution. Cette particularité permet l'observation de différentes formes de flammes et notamment de points bistables pour lesquels deux flammes différentes peuvent exister malgré des conditions opératoires identiques.L'utilisation de diagnostics optiques à haute cadence (diffusion de Mie des gouttes de combustible et émission spontanée de la flamme) est couplée à des méthodes de post-traitement avancées comme la Décomposition en Modes Dynamiques. Ainsi, des mécanismes pilotant la stabilisation des flammes ainsi que leurs changements de forme sont proposés. Ils mettent notamment en lumière les interactions entre l'écoulement gazeux, les gouttes de combustible et la flamme. / A promising way to reduce jet engines pollutant emissions is the use of lean premixed prevaporized combustion but it tends to trigger thermo-acoustic instabilities. To improve the stability of these flames, a procedure called staging consists in splitting the fuel injection to control its spatial distribution. This however leads to an increased complexity and unexpected phenomena can occur.In the present work, a model gas turbine combustor fed with liquid dodecane is used. It is equipped with two fuel injection stages to control the fuel distribution in the burner. Different flame stabilizations can be observed and a bistable case where two flame shapes can exist for the same operating conditions is highlighted.High-speed optical diagnostics (fuel droplets Mie scatering and chemiluminescence measurements) are coupled with advanced post-processing methods like Dynamic Mode Decomposition. The results enable to propose mechanisms leading to flame stabilization and flame shape transitions. They show a strong interplay between the gaseous flow, the fuel droplets and the flame itself.

Large Eddy Simulation (LES) of bluff body stabilized premixed and partially premixed combustion close to the flammability limit is carried out in this thesis. The LES algorithm has no ad-hoc adjustable model parameters and is able to respond automatically to variations in the inflow conditions. Algorithm validation is achieved by comparison with reactive and non-reactive experimental data. In the reactive flow, two scalar closure models, Eddy Break-Up (EBULES) and Linear Eddy Mixing (LEMLES), are used and compared. Over important regions, the flame lies in the Broken Reaction Zone regime. Here, the EBU model assumptions fail. The flame thickness predicted by LEMLES is smaller and the flame is faster to respond to turbulent fluctuations, resulting in a more significant wrinkling of the flame surface. As a result, LEMLES captures better the subtle effects of the flame-turbulence interaction. Three premixed (equivalence ratio = 0.6, 0.65, and 0.75) cases are simulated. For the leaner case, the flame temperature is lower, the heat release is reduced and vorticity is stronger. As a result, the flame in this case is found to be unstable. In the rich case, the flame temperature is higher, and the spreading rate of the wake is increased due to the higher amount of heat release Partially premixed combustion is simulated for cases where the transverse profile of the inflow equivalence ratio is variable. The simulations show that for mixtures leaner in the core the vortical pattern tends towards anti-symmetry and the heat release decreases, resulting also in instability of the flame. For mixtures richer in the core, the flame displays sinusoidal flapping resulting in larger wake spreading. More accurate predictions of flame stability will require the use of detailed chemistry, raising the computational cost of the simulation. To address this issue, a novel algorithm for training Artificial Neural Networks (ANN) for prediction of the chemical source terms has been implemented and tested. Compared to earlier methods, the main advantages of the ANN method are in CPU time and disk space and memory reduction.

This thesis describes an experimental investigation of the response of lean turbulent swirling flames with different degrees of premixedness (i.e. different mixture patterns) to acoustic forcing using the same burner configuration and varying only the fuel injection strategy. Special emphasis was placed on the amplitude dependence of their response. Also, the behaviour of self-excited fully premixed flames was examined. kHz OH* chemiluminescence was used to study qualitatively the heat release response of the flames, while kHz OH Planar Laser Induced Fluorescence (PLIF) was employed to understand the response of the flame structure and the behaviour of the various parts of the flame. The Proper Orthogonal Decomposition (POD) method was used to extract the dominant structures of the flame and their periodicity. In the first part of the thesis, self-excited oscillations were induced by extending the length of the duct downstream of the bluff body. It was found that the longer the duct length and the higher the equivalence ratio, the stronger the self-excited oscillations were, with the effect of duct length being much stronger. The dominant frequencies of the system were found to increase with equivalence ratio and bulk velocity and decrease with duct length. For some conditions, three simultaneous periodic motions were observed, where the third motion oscillated at a frequency equal to the difference of the other two frequencies. A novel application of the POD method was proposed to estimate the convection velocity from the most dominant reaction zone structures detected by OH* chemiluminescence imaging. For a range of conditions, the convection velocity was found to be in the range of 1.4-1.7 bulk flow velocities at the inlet of the combustor. In the second part, the response of fully premixed, non-premixed with radial fuel injection (NPR) and axial fuel injection (NPA) flames was investigated and compared. All systems exhibited a nonlinear response to acoustic forcing. The highest response was observed by the NPR flame, followed by the fully premixed and the non-premixed with axial fuel injection flame. The proximity of forced flames to blow-off was found to be critical in their heat release response, as close to blow-off the flame response was significantly lower than that farther from blow-off. In the NPR and NPA systems, it was shown that the acoustic forcing reduced the stability of the flame and the stability decreased with the increase in forcing amplitude. In the fully premixed system, the flame area modulations constituted an important mechanism of the system, while in the NPR system both flame area and equivalence ratio modulations were important mechanisms of the heat release modulations. The quantification of the local response of the various parts of the flame at the forcing frequency showed that the ratio RL (OH fluctuation at 160 Hz to the total variance of OH) was greater in the inner shear layer region than in the other parts in the case of NPR and NPA flames. In fully premixed flames, greater RL values were observed in large regions on the downstream side of the flame than those in the ISL region close to the bluff body. The ratio of the convection velocity to the bulk velocity was estimated to be 0.54 for the NPR flame, while it was found to be unity for the respective fully premixed flame. In the last part of the thesis, the response of ethanol spray flames to acoustic oscillations was investigated. The nonlinear response was very low, which was reduced closer to blow-off. The ratio RL was the highest in the spray outer cone region, downstream of the annular air passage, while RL values were very low in the inner cone region, downstream of the bluff body. Unlike NPR and fully premixed flames, in case of spray and NPA systems, it was found that forcing did not affect greatly the flame structure. The understanding of the nonlinear response of flames with different degrees of premixedness in a configuration relevant to industrial systems contributes to the development of reliable flame response models and lean-burn devices, because the degree of premixedness affects greatly the flame response. Also, the understanding of the behaviour of forced spray flames is of great interest for industrial applications, contributing to the development of thermoacoustic models for liquid fuelled combustors. Finally, the estimation of the convection velocity is of importance in the modelling of self-excited flames and flame response models, since the convection velocity affects the flame response significantly.

L'influence de l'ajout d'hydrogène sur les flammes de premelange pauvre methane-air est simulée dans cette étude. Le modèle de la chambre a haute pression Orleans - ICARE (France), a été développé. Les propriétés du front de flamme sont examinées par deux modèles de combustion turbulente prémélangée, à savoir Zimont et Flamme Cohérente Model (CFM) modèles.Toutes les études de modélisation sont effectués avec le logiciel Fluent et les résultats sont comparés aux expériences. En suite, l'influence de la pression sur les statistiques de la front de flamme prémélangée a été examinée. Les simulations montrent que l'augmentation du ratio d'équivalence a diminué la hauteur des flammes et l'épaisseur de la flamme du méthane/air flames. D'autre part, l’ajout d’hydrogene de mélange pauvre méthane-air a modifié les propriétés de la flamme prémélangée. Lorsque le pourcentage volumique de l'hydrogène dans le mélange est augmenté, la position en hauteur de la flamme est réduite et l'épaisseur de la flamme devient plus mince. En outre, il a été observé que les propriétés de la flamme prémélangée ont été modifiées avec l'opération à des conditions de pression plus élevée. / Hydrogenated premixed methane/air flames under lean conditions are simulated in this study. The model of the high pressure chamber setup of Orleans - ICARE (France) has been developed. The flame front properties are investigated by two turbulent premixed combustion models, Zimont and Coherent Flame Model (CFM) models. All modeling studies are performed with Fluent software and compared to experiments. The influence of the pressure on the premixed flame front statistics has been examined as well. The simulations show that increasing the equivalence ratio decreases the flame tip height and the flame brush thickness for methane/air flames. In addition, enriching the methane-air mixture with hydrogen modifies the premixed flame front properties. When the volumetric percentage of hydrogen in the mixture is increased, the flame-end position is reduced and flame brush thickness becomes thinner. It is also observed that the premixed flame properties have been modified with operation at higher pressure conditions.