京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第21100号 / 工博第4464号 / 新制||工||1694(附属図書館) / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 黒瀬 良一, 教授 中部 主敬, 教授 吉田 英生 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

An experimental study of the effects of turbulent flow regime on the flame structure is conducted by using perforated-plate-stabilized hydrogen-piloted lean premixed methane/air turbulent flames. The underlying non-reacting turbulent flow field was investigated using two-dimensional three-components particle imaging velocimetry (2D3C-PIV) with and without three perforated plates. The non-reacting flow data allowed a separation of the turbulent flow regime into axial velocity dominated and vortex dominated flows. A plate with 62\% blockage ratio was used to represent the stream-dominant flow regime and another with 86\% blockage ratio was used to represent the vortex-dominant flow regime. OH laser-induced fluorescence was used to study the effects of the turbulent flow regime on the mean progress variable, flame brush thickness, flame surface density, and global consumption speed. In comparison with the stream-dominant flow, the vortex-dominant flow makes a wider and shorter flame. Also, the vortex-dominant flow has a thicker horizontal flame brush thickness and a thinner longitudinal flame brush thickness. Especially, the horizontal flame brush thickness for the vortex-dominant flow does not follow the turbulence diffusion theory. Then, the vortex-dominant flow shows a relatively constant flame surface density along the stream-wise direction, while the stream-dominant flow shows a decreasing flame surface density. Lastly, the vortex-dominant turbulent flow improves the consumption speed in comparison to the stream-dominant turbulent flow regime with the same velocity fluctuation level.

La combustion pré mélangée a été depuis longtemps un domaine vaste d’étude au niveau appliqué et fondamental. Bien que la plupart des applications industrielles en combustion aient lieu dans des régimes turbulents, le passage par l’étude laminaire est indispensable pour comprendre les mécanismes fondamentaux des flammes turbulentes. Ce travail de thèse porte essentiellement sur l’étude des différentes instabilités agissant sur un front de flamme laminaire de pré mélange pour des mélanges de propane-air et méthane-air, enrichis ou non en azote. L’étude consiste à mesurer les taux de croissance des perturbations dans un brûleur d’Hele-Shaw formé par deux plaques de verre ($150\times50cm$) très rapprochées (brûleur 2D). Grâce à un système de forçage constitué par des plaques modulées spatialement avec différentes longueurs d’onde, le taux de croissance peut être mesuré même en présence du développement spontané de l’instabilité avec la longueur d’onde la plus instable. A richesse constante et pour des valeurs croissantes de la dilution en oxygène le vecteur d'onde de coupure augmente avec la vitesse de flamme qui devient plus instable. Pour des mélanges de propane ce nombre d’onde augmente également lorsqu'on augmente la richesse à vitesse constante et il diminue dans le cas des mélanges de méthane, cela est en relation avec l'évolution des effets thermodiffusifs dans chacun des cas. Une augmentation de l’épaisseur de la cellule d’Hele-Shaw, aboutit à une augmentation du taux de croissance pour les petits nombres d’onde ainsi que du nombre de Markstein, et à une diminution du taux de croissance pour les grands nombres, du fait des effets des pertes thermiques. / Premixed combustion has been the subject of extensive work, concerning both applications and fundamental aspects. Although in most practical applications combustion occurs in a turbulent environment, the study of laminar flames is important to understand the fundamental mechanisms of turbulent flame propagation. The objective of this work is to study the various instabilities acting on a laminar premixed flame for mixtures of propane-air and methane-air, enriched or not with nitrogen. it consists in measuring the linear growth rates of disturbances in a Hele-Shaw burner formed by two glass plates ($150\times50cm$) separated by a thin gap width (2D burner). Using spatially modulated plates with different wavelengths, the linear growth rate of perturbations can be measured even in the presence of the most unstable wavelength. The experimental values of the linear growth rate as a function of wavenumber are fitted by a linear dispersion relation to estimate the Markstein number and the cutoff wavenumber. For a constant equivalence ratio with increasing values of the dilution in oxygen, the cutoff wavenumber grows with the flame velocity and it is becoming more unstable. The cutoff wave number rises also when the equivalence ratio increases for propane-air mixture and decreases for methane-air mixture, in relation to the evolution of thermal diffusive effects. An enlargement in the thickness of the Hele-Shaw cell results in an increase of the growth rate for small wavenumbers thus in the Markstein number, and in a decline in the growth rate for the large wavenumbers, in relation with the effects of heat losses.

A full Strouhal number thermo-acoustic model is purposed for the feedback control of self excited noise driven combustors. The inclusion of time delays in the volumetric heat release perturbation models create unique behavioral characteristics which are not properly reproduced within current low Strouhal number thermo acoustic models. New analysis tools using probability density functions are introduced which enable exact expressions for the statistics of a time delayed system. Additionally, preexisting tools from applied mathematics and control theory for spectral analysis of time delay systems are introduced to the combustion community. These new analysis tools can be used to extend sensitivity function analysis used in control theory to explain limits to control effectiveness in self-excited combustors. The control effectiveness of self-excited combustors with actuator constraints are found to be most sensitive to the location of non-minimum phase zeros. Modeling the non-minimum phase zeros correctly require accurate volumetric heat release perturbation models. Designs that removes non-minimum phase zeros are more likely to have poles in the right hand complex plane. As a result, unstable combustors are inherently more responsive to feedback control.

Le monoxyde d’azote (NO) est un polluant atmosphérique responsable d’effets nuisibles sur l’environnement et la santé. Afin de mieux contrôler ces émissions, il est indispensable de comprendre et de maîtriser leur formation,en particulier lors de la combustion à haute pression, domaine d’application industrielle (cas des turbines à gaz,des moteurs…). On distingue quatre voies principales de formation de NO : la voie thermique, la voie du NO précoce, la voie NNH et la voie N2O. L’objectif de cette thèse à caractère expérimentale est de compléter la base de données expérimentale déjà existante nécessaire à la compréhension et à l’identification de la contribution de chaque voie à la formation du NO à haute pression.Dans cette thèse, un dispositif de brûleurs à contre-courants a été utilisé pour étudier la structure de flammes laminaires, prémélangées à haute pression. Les profils de concentration de NO dans les flammes CH4/O2/N2 à différentes richesses (Фc =0,7-1,2) et différentes pressions (P=0,1-0,7 MPa) ont été mesurés par Fluorescence Induite par Laser. L’effet de l’ajout d’hydrogène (80%CH4/20%H2 : Application Hythane®) sur la formation de NO a également été étudié dans les flammes pauvres CH4/O2/N2. Le mécanisme cinétique GDF-Kin®3.0_NCN a été comparé aux mesures de NO disponibles dans la littérature ainsi qu’aux simulations des mécanismes cinétiques du Gaz Research Institute (version 2.11 et 3.0). Ces trois mécanismes ont été ensuite comparés aux mesures expérimentales réalisées dans ces travaux de thèse. / The nitric oxide (NO) is a pollutant responsible of detrimental effects on the environment and health. To better control these emissions, it’s crucial to understand and to control their formation, in particular during the combustion process at high pressure, area of industrial applications (gas turbines, engines…).There are four major routes of the NO formation: the thermal route, the prompt-NO route, the NNH route and theN2O route. The aim of this experimental thesis is to complete the existing experimental database which isnecessary to the understanding and the identification of the contribution from each route to the NO formation at high pressure.In this thesis, a facility of two twin counter-flow burners was used to study the structure of the laminar, premixed flames at high pressure. Experimental NO concentration profiles have been measured in CH4/O2/N2 flames for arange of equivalence ratio (from 0.7 to 1.2) and pressures (from 0.1 to 0.7 MPa) by Laser Induced Fluorescence.The effect of adding hydrogen (80%CH4/20%H2: Hythane® application) on the NO formation has been also studied in lean CH4/O2/N2 flames. The GDF-Kin®3.0_NCN kinetic mechanism has been compared to experimental data from the literature and also compared to the simulations from the Gas Research Institute mechanisms (version 2.11 and 3.0). These three mechanisms have been finally compared to the experimental data from this thesis.

<div>This research has been studied a laser-induced spark in flowing gases. The relationship between the minimum ignition energy (MIE), the turbulence intensity, and the flame kernel propagation speed is considered. Plasma emission, produced by the laser-induced spark, and flame kernel generation by the plasma are investigated. The energy balance equation between an ignition energy and energy losses by heat transfer is studied at laminar flows and turbulent flows. Hydrogen and air mixtures were used in a premixed jet burner for ignition experiments. Particle image velocimetry (PIV) examined the velocity and the turbulence intensity under the turbulent flows. The flame kernel development was visualized using Schlieren imaging and infrared images (IR camera). Flame kernel temperatures were measured through Rayleigh scattering and infrared images (IR camera). Plasma evaluations were captured through an intensified CCD camera (ICCD camera). Minimum ignition energies were measured at the laminar flows and the turbulent flows. The MIE decreases with an increase in the turbulence intensity which changed by ignition locations and perforated plates at the constant bulk velocity. Improved mixing rates due to the ignition locations or the geometry of the perforated plates decrease the MIE at the constant bulk velocity. The turbulence intensity increases wrinkles in the flame kernel surface, thus the contact between the flame kernel and reactants increases due to the wrinkles. Therefore, the flame kernel propagation speed increases as the turbulence intensity is higher since the increased reaction by the wrinkles and the contact. Thus, the MIE decreases as the turbulence intensity increases at the constant ignition condition, including bulk velocities and ignition heights, since the high turbulence intensity increases the flame kernel propagation speed. Laser energy differences affect the plasma expansions by the laser absorption. Laser-supported radiation (LSR) wave speeds were measured and calculated using energy balance equations. Velocity does not affect the flame kernel temperature distribution during the early reaction steps because the plasma generates a flame kernel and determines the flame kernel temperature distribution. The MIE increases with increasing the bulk velocity. The energy losses considering convection, conduction, and radiation were calculated using the flame kernel radius, the flame kernel temperature, mixture properties, and the flame speed. The energy balance equation in the ignition of flowing gases is newly written at the laminar flows and the turbulent flows.</div>

<p>Combustion instability, a complex phenomenon observed in combustion chambers is due to the coupling between heat release and other unsteady flow processes. Combustion instability has long been a topic of interest to rocket scientists and has been extensively investigated experimentally and computationally. However, to date, there is no computational tool that can accurately predict the combustion instabilities in full-size combustors because of the amount of computational power required to perform a high-fidelity simulation of a multi-element chamber. Hence, the focus is shifted to reduced fidelity computational tools which may accurately predict the instability by using the information available from the high-fidelity simulations or experiments of single or few-element combustors. One way of developing reduced fidelity computational tools involves using a reduced fidelity solver together with the flame transfer functions that carry important information about the flame behavior from a high-fidelity simulation or experiment to a reduced fidelity simulation.</p> <p> </p> <p>To date, research has been focused mainly on premixed flames and using acoustic solvers together with the global flame transfer functions that were obtained by integrating over a region. However, in the case of rockets, the flame is non-premixed and distributed in space and time. Further, the mixing of propellants is impacted by the level of flow fluctuations and can lead to non-uniform mean properties and hence, there is a need for reduced fidelity solver that can capture the gas dynamics, nonlinearities and steep-fronted waves accurately. Nonlinear Euler equations have all the required capabilities and are at the bottom of the list in terms of the computational cost among the solvers that can solve for mean flow and allow multi-dimensional modeling of combustion instabilities. Hence, in the current work, nonlinear Euler solver together with the spatially distributed local flame transfer functions that capture the coupling between flame, acoustics, and hydrodynamics is explored.</p> <p> </p> <p>In this thesis, the approach to extract flame transfer functions from high-fidelity simulations and their integration with nonlinear Euler solver is presented. The dynamic mode decomposition (DMD) was used to extract spatially distributed flame transfer function (FTF) from high fidelity simulation of a single element non-premixed flame. Once extracted, the FTF was integrated with nonlinear Euler equations as a fluctuating source term of the energy equation. The time-averaged species destruction rates from the high-fidelity simulation were used as the mean source terms of the species equations. Following a variable gain approach, the local species destruction rates were modified to account for local cell constituents and maintain correct mean conditions at every time step of the nonlinear Euler simulation. The proposed reduced fidelity model was verified using a Rijke tube test case and to further assess the capabilities of the proposed model it was applied to a single element model rocket combustor, the Continuously Variable Resonance Combustor (CVRC), that exhibited self-excited combustion instabilities that are on the order of 10% of the mean pressure. The results showed that the proposed model could reproduce the unsteady behavior of the CVRC predicted by the high-fidelity simulation reasonably well. The effects of control parameters such as the number of modes included in the FTF, the number of sampling points used in the Fourier transform of the unsteady heat release, and mesh size are also studied. The reduced fidelity model could reproduce the limit cycle amplitude within a few percent of the mean pressure. The successful constraints on the model include good spatial resolution and FTF with all modes up to at least one dominant frequency higher than the frequencies of interest. Furthermore, the reduced fidelity model reproduced consistent mode shapes and linear growth rates that reasonably matched the experimental observations, although the apparent ability to match growth rates needs to be better understood. However, the presence of significant heat release near a pressure node of a higher harmonic mode was found to be an issue. This issue was rectified by expanding the pressure node of the higher frequency mode. Analysis of two-dimensional effects and coupling between the local pressure and heat release fluctuations showed that it may be necessary to use two dimensional spatially distributed local FTFs for accurate prediction of combustion instabilities in high energy devices such as rocket combustors. Hybrid RANS/LES-FTF simulation of the CVRC revealed that it might be necessary to use Flame Describing Function (FDF) to capture the growth of pressure fluctuations to limit cycle when Navier-Stokes solver is used.</p> <p> </p> <p>The main objectives of this thesis are:</p> <p>1. Extraction of spatially distributed local flame transfer function from the high fidelity simulation using dynamic mode decomposition and its integration with nonlinear Euler solver</p> <p>2. Verification of the proposed approach and its application to the Continuously Variable Resonance Combustor (CVRC).</p> <p>3. Sensitivity analysis of the reduced fidelity model to control parameters such as the number of modes included in the FTF, the number of sampling points used in the Fourier transform of the unsteady heat release, and mesh size.</p> <p> </p> <p>The goal of this thesis is to contribute towards a reduced fidelity computational tool which can accurately predict the combustion instabilities in practical systems using flame transfer functions, by providing a path way for reduced fidelity multi-element simulation, and by defining the limitations associated with using flame transfer functions and nonlinear Euler equations for non-premixed flames.</p> <p> </p><br>

L'épuisement des réserves pétrolières et l'augmentation de la concentration du gaz à effet de serre CO2 sont les deux principaux problèmes connus liés à l'utilisation des carburants fossiles. Les biocarburants apparaissent comme un des moyens permettant à la fois une diminution de la dépendance au pétrole et une réduction de l'impact néfaste des moteurs automobiles sur l'environnement. Les biocarburants sont en effet considérés comme une source d'énergie renouvelable. L'objectif de cette thèse était de développer et valider les modèles cinétiques de combustion des composés oxygénés de biocarburants : l'éthanol, les biocarburants de deuxième-génération des familles du furane (furane, 2-méthylfurane, 2,5-diméthylfurane), du tétrahydrofurane (tétrahydrofurane, 2-méthyltétrahydrofurane) et le tétrahydropyrane, en utilisant les nouvelles données obtenues en flamme laminaire pré-mélangée à basse pression. De 20 à 60 produits ont été quantifiés par chromatographie en phase gazeuse et identifiés par couplage avec la spectrométrie de masse. Les résultats obtenus ont ensuite été utilisés pour analyser les voies de consommation des réactifs et de formation des produits, surtout pour les polluants, dans le but de mieux comprendre la chimie de la combustion de ces biocarburants. Ce rapport comprend 5 chapitres et une conclusion. Le premier chapitre présente une revue bibliographique des travaux antérieurs sur l'oxydation de l'éthanol et des éthers cycliques. Dans le second chapitre, le dispositif expérimental est décrit, en détaillant en particulier les nouveaux développements. Enfin les chapitres 3, 4, 5 présentent les résultats de l'étude de la combustion des composés étudiés / The decrease of petroleum reserves and the increase of concentration of greenhouse gas CO2 are the two major known problems related to the use of fossil fuels. Bio-fuels appear as a means allowing a decrease of the dependence on fossil fuels and a reduction of the harmful impact of engine on the environment. Bio fuels are considered as a source of renewable energy. The aim of this thesis was to develop and validate experimentally the high temperature kinetic models for the combustion of oxygenated compounds of bio-fuels: ethanol, second-generation bio-fuels of families of furan (furan, 2-methylfuran, 2,5-dimethylfuran), of tetrahydrofuran (tetrahydrofuran, 2 methyltetrahydrofuran), and tetrahydropyran, using new data obtained in laminar premixed low-pressure flame. About 20-60 products were quantified by gas chromatography and identified using mass spectrometry. The results obtained were then used to analyze the consumption pathways of fuels and the formation pathways of products, especially for pollutants, in order to better understand the combustion chemistry of these bio-fuels. This thesis report includes 5 chapters and a conclusion. The first chapter presents a review of the major works already published in the literature for the oxidation of ethanol and cyclic ethers. In the second chapter, the experimental setup of laminar premixed flame with the analytical techniques is described, detailing in particular new developments. Eventually, chapters 3, 4, 5 present the experimental and modeling results of the study of the combustion chemistry of the compounds studied

This thesis aims at improving our knowledge on swirl combustors. The work presented here is based on Large Eddy Simulations (LES) coupled to an advanced combustion model: the Conditional Moment Closure (CMC). Numerical predictions have been systematically compared and validated with detailed experimental datasets. In order to analyze further the physics underlying the large numerical datasets, Proper Orthogonal Decomposition (POD) has also been used throughout the thesis. Various aspects of the aerodynamics of swirling flames are investigated, such as precession or vortex formation caused by flow oscillations, as well as various combustion aspects such as localized extinctions and flame lift-off. All the above affect flame stabilization in different ways and are explored through focused simulations. The first study investigates isothermal air flows behind an enclosed bluff body, with the incoming flow being pulsated. These flows have strong similarities to flows found in combustors experiencing self-excited oscillations and can therefore be considered as canonical problems. At high enough forcing frequencies, double ring vortices are shed from the air pipe exit. Various harmonics of the pulsating frequency are observed in the spectra and their relation with the vortex shedding is investigated through POD. The second study explores the structure of the Delft III piloted turbulent non-premixed flame. The simple configuration allows to analyze further key combustion aspects of combustors, with further insights provided on the dynamics of localized extinctions and re-ignition, as well as the pollutants emissions. The third study presents a comprehensive analysis of the aerodynamics of swirl flows based on the TECFLAM confined non-premixed S09c configuration. A periodic component inside the air inlet pipe and around the central bluff body is observed, for both the inert and reactive flows. POD shows that these flow oscillations are due to single and double helical vortices, similar to Precessing Vortex Cores (PVC), that develop inside the air inlet pipe and whose axes rotate around the burner. The combustion process is found to affect the swirl flow aerodynamics. Finally, the fourth study investigates the TECFLAM configuration again, but here attention is given to the flame lift-off evident in experiments and reproduced by the LES-CMC formulation. The stabilization process and the pollutants emission of the flame are investigated in detail.