La circulation atmosphérique en Afrique de l'Ouest est caractérisée par des vents de sud-ouest (mousson) pendant la saison humide et par des vents de nord-est (harmattan) pendant la saison sèche. Cette alternance des saisons est due aux variations de pression liée à l'état des surfaces (rugosité, albédo, végétation) en réaction au forçage solaire. Ces mêmes états de surface génèrent une variabilité de flux turbulents de surface et des circulations secondaires qui rendent complexes les analyses des mesures effectuées sur place en vue de documenter les interactions surface-atmosphère. La modélisation fine échelle (LES) couramment utilisée dans l'étude de la couche limite atmosphérique est requise pour pouvoir palier à ces difficultés en raison de sa capacité à prendre en compte les flux turbulents en 3D et sur topographie complexe. Notre site d'étude est le bassin versant d'ARA située au Nord du Bénin dans un contexte Soudanien avec des propriétés de surface variables. Une analyse climatique est effectuée sur la base des observations de radiosondage, de radar UHF et de stations au sol afin d'extraire des données composites représentatives des saisons sèche et humide. Ces données composites ont servi par la suite à forcer le modèle Méso-NH dans sa version LES. Pour pouvoir caractériser les échelles de longueur des flux turbulents de surface relatives aux saisons sèche et humide, les données standard de forçage de surface de Méso-NH que sont le relief GTOPO30 (1km de résolution) et la végétation ECOCLIMAP (1km de résolution) ont été respectivement remplacer par le SRTM (90m de résolution) et les données de SPOT/HRV (20m de résolution) reéchantillonné à 90m de résolution. A l'aide d'outils statistiques comme la variographie 2D et le suivi Lagrangien, il ressort que la variabilité spatiale de la chaleur sensible H est gouvernée par le couple vent-relief tandis que celle de la chaleur latente E est difficile à mettre en lien sur végétation hétérogène (SPOT/HRV) en saison sèche. En saison humide, la variabilité spatiale du champ H dépend du vent tandis que celle du champ E dépend de la végétation. Cette étude révèle dans tous les cas que les échelles caractéristiques de ces deux champs diffèrent dans les mêmes conditions de forçage de surface et atmosphérique. / West Africa atmosphere circulation is characterized by south-westerly wind (monsoon regime) during the wet season and north-easterly wind (harmattan regime) during the dry season. This alternation of wind regime is due to surface pressure variability linked to surface heterogeneities. Surface heterogeneities generate surface flux variability, secondary circulation and make complex analysis when trying to document surface-atmosphere feedbacks. LES modelling usually used for boundary-layer studies due to its potential to take into account 3D turbulence over complex topography, is used here to overcome these difficulties. Our site of interest is located in north of Benin characterized by Soudanian climate and heterogeneous surface properties. Climate analysis are first performed with radiosoundings, UHF radar, and EC station data in order to extract composite profile representing dry and wet season.. These composite profiles are then used to force atmosphere part of the Méso-NH LES model. To characterize turbulent fluxes length scales relative to dry and wet season, standard surface forcing data with Méso-NH like GTOPO30 orography (1km ) and ECOCLIMAP vegetation (1km) are respectively replaced by SRTM (90m) and SPOT/HRV vegetation data (20m) resampled to 90m. Along with statistical tools like 2D variography and Lagrangian, we notice that during dry season on heterogeneous vegetation, sensible heat flux H is more driven by wind and orography while we not able to discuss the latent heat flux E case. During wet season with the same surface forcing, it appears that H is driven by wind while E is more dependent to vegetation variability. Our study concludes in all case that H and E are not characterized by the same length scale.

Le contexte de cette étude est l’optimisation de la combustion turbulente prémélangée de syngaz pour la production propre d’énergie. Un brûleur CH4/air de type bec Bunsen avec forçage turbulent multi-échelles produit par un système de trois grilles, est simulé numériquement par différentes techniques de simulation des grandes échelles (SGE), et les résultats sont comparés à l’expérience. On a développé et appliqué une formulation bas-Mach du solveur Navier-Stokes basé sur différents schémas numériques, allant des différences finies centrées d’ordre 4 à des versions avancées des schémas WENO d’ordre 5. La méthodologie est évaluée sur une série de cas-tests classiques (flamme laminaire 1D prémélangée, turbulence homogène et isotrope en auto-amortissement), et sur des simulations 2D de la flamme turbulente prémélangée expérimentale. Les SGE implicites (ILES), i.e. sans aucune modélisation sous-maille, et explicites avec le modèle de flamme épaissie et un modèle de plissement sous-maille nouvellement élaboré (TFLES), sont appliquées à la simulation 3D du brûleur expérimental. Les résultats montrent que l’approche TFLES avec un schéma d’ordre élevé à faible dissipation numérique prédit correctement la longueur de la flamme et la densité de surface de flamme. La SGE implicite avec un schéma WENO avancé produit une flamme trop courte mais réaliste à condition que la taille de la maille soit de l’ordre de l’épaisseur de flamme laminaire. La représentation des interactions flamme/turbulence est néanmoins très différente entre TFLES et ILES. / The context of this study is the optimization of premixed turbulent combustion of syngas for clean energy production. A Bunsen-type CH4/air turbulent premixed burner with a multi-scale grid generator is simulated with different Large Eddy Simulation (LES) strategies and compared to experimental results. A low-Mach formulation of a compressible Navier-Stokes solver based on different numerical methods, ranging from 4th order central finite difference to 5th order advanced WENO schemes, is developed and applied. Classical test cases (1D laminar premixed flame, decaying HIT), and 2D simulations of the turbulent premixed flame are performed to assess the numerical methodology. Implicit LES (ILES), i.e. LES without any explicit subgrid modeling, and explicit LES with the Thickened Flame model and subgrid scale flame wrinkling modelling (TFLES) are applied to simulate numerically the 3D experimental burner. Results show that TFLES with a high-order low dissipation scheme predicts quite well the experimental flame length and flame surface density. ILES with advanced WENO schemes produces a slightly shorter although realistic flame provided the grid spacing is of order of the laminar flame thickness. The representation of flame/turbulence interactions in TFLES and ILES are however quite different.

O objetivo deste trabalho é estudar, dentro da Engenharia do Vento Computacional (EVC), a dispersão de poluentes em zonas urbanas, empregando-se um modelo numérico baseado em técnicas da Dinâmica dos Fluidos Computacional para escoamentos incompressíveis, não isotérmicos e com transporte de massa. Um esquema explícito de dois passos é usado para a discretização temporal das equações governantes, considerando expansões em séries de Taylor de segunda ordem para as derivadas no tempo. O processo de discretização espacial é realizado através da aplicação do Método dos Elementos Finitos (MEF), onde hexaedros de oito nós com um ponto de integração são utilizados. A turbulência é tratada numericamente através da Simulação de Grandes Escalas (LES) e os modelos clássico e dinâmico de Smagorinsky são empregados na modelagem das escalas inferiores à resolução da malha. Efeitos de temperatura sobre o escoamento são considerados na forma de forças de flutuação presentes na equação de balanço de momentum, as quais são calculadas a partir da aproximação de Boussinesq. Técnicas de paralelização em memória compartilhada (OpenMP) são também usadas a fim de melhorar a eficiência computacional do presente modelo para problemas com grande número de elementos. Exemplos clássicos de Dinâmica de Fluidos e Fenômenos de Transporte são inicialmente analisados para teste das ferramentas numéricas implementadas. Problemas de dispersão de poluentes com e sem a inclusão dos efeitos de temperatura são abordados para configurações geométricas bi e tridimensionais de street canyons, representando a unidade geométrica básica encontrada em centros urbanos de grandes cidades. / The main goal of the present work is to study the pollutant dispersion in urban areas using a numerical model based on techniques developed by Computational Fluid Dynamics, where applications of Computational Wind Engineering (CWE) are analyzed considering incompressible flows with heat and mass transport. A two-step explicit scheme is adopted for the time discretization of the governing equations considering second order Taylor series expansions of the time derivative terms. Spatial discretization is performed by applying the Finite Element Method (FEM), where eight-node hexahedral elements with one-point quadrature are utilized. Turbulence is numerically analyzed by using Large Eddy Simulation (LES) with the classical and dynamic Smagorinsky’s models for subgrid scale modeling. Thermal effects on the flow field are taken into account through buoyancy forces acting on the momentum balance equation, which are calculated considering the Boussinesq approximation. Shared memory parallelization techniques (OpenMP) are also employed in order to improve computational efficiency for problems with large number of elements. Classic examples of Fluid Dynamics and Transport Phenomena are first analyzed to verify the numerical tools implemented. Problems involving pollutant dispersion with and without the inclusion of thermal effects are investigated for two and three-dimensional geometric configurations of street canyons, which represent the basic geometric unit observed in urban centers of large cities.

O experimento Unsteady Aerodynamics Experiment Phase VI, realizado no ano de 2000 pelo Laboratório Nacional norte-americano para as Energias Renováveis (NREL) no túnel de vento Ames da NASA, foi reproduzido numericamente neste trabalho. O objetivo é o estudo das características da esteira aerodinâmica produzida pela turbina eólica de duas pás e 10 metros de diâmetro, operando à velocidade de rotação constante de 72 RPM, sujeita a uma velocidade de corrente livre do vento uniforme de 9 m/s, em um túnel de vento cuja seção de testes mede 36,6 m de largura por 24,4 m de altura e o comprimento mede 170 m. Para isso, foi utilizado o programa comercial ANSYS FLUENT versão 13.0, baseado no Método dos Volumes Finitos para a solução numérica das Equações de Navier-Stokes em regime transiente em conjunto com a Simulação de Grandes Escalas (SGE) para resolver a turbulência. As geometrias de todos os componentes da máquina foram criadas em software CAD. Um domínio móvel em forma de disco, contendo as pás do rotor e o hub da máquina, foi criado separadamente, e posteriormente inserido no domínio principal, estático, usando a ferramenta Moving Mesh disponível no software FLUENT. Ambos os domínios foram preenchidos por malhas compostas por tetraedros. Dados provenientes das simulações numéricas foram comparados aos dados experimentais de velocidade fornecidos por dois anemômetros sônicos instalados 5,8 m à jusante do rotor, ao que foi verificada boa concordância, com diferenças da ordem de 1% para o anemômetro 1 e 6% para o anemômetro 2. Resultados de velocidade na linha de centro do túnel e perfis de velocidade à jusante foram comparados com recente estudo numérico, e revelam diferenças importantes entre dados obtidos pela SGE, principalmente no que se refere à detecção de picos e flutuações relacionados às escalas turbulentas, e dados obtidos através da modelagem clássica da turbulência, RANS. As perturbações ultrapassaram a marca dos 10 diâmetros à jusante e atingiram o final do domínio localizado a 15 diâmetros. A esteira não apresentou simetria axial, e o ponto de maior redução na velocidade do escoamento foi detectado fora da linha de centro do rotor. / The Unsteady Aerodynamics Experiment Phase VI, which has been carried out in 2000 by the US National Renewable Energy Laboratory (NREL) at the NASA Ames wind tunnel, has been numerically reproduced. The purpose of this work is to study the characteristics of the wind wake produced by the 10 meter two bladed wind turbine, operating at a constant rotational speed of 72 RPM, subject to a free stream wind velocity of 9 m/s, inside a wind tunnel in which dimensions are 36.6 m in width, 24.4 m in height and length of 170 m. To achieve that, the ANSYS FLUENT version 13.0 commercial code, based in the Finite Volume Method to numerically solve the Navier-Stokes equations in transient state, has been used, together with the Large Eddy Simulation (LES) to characterize the turbulence. Geometries of all the machine components have been created in CAD software. A disc shaped moving domain, containing the blades and hub, has been created separately, and later inserted into the main, static domain, using the Moving Mesh tool available in the software. Both domains have been filled with meshes composed by tetrahedra. Data collected at the numerical simulations have been compared to experimental wind speed data provided by two sonic anemometers installed 5.8 m downstream from the rotor, for which a good agreement has been found, with differences of approximately 1% to the anemometer 1 and 6% to the anemometer 2. Results of wind velocity at the tunnel centerline and velocity profiles downstream have been compared with recent numerical study, and show important differences between data obtained by LES, especially with regard to the detection of peaks and fluctuations related to the turbulent scales, and data obtained by the classic turbulence modeling, RANS. Disturbances have passed the 10 diameter mark and reached the end at the domain located at 15 diameters. The wake did not show axial symmetry and the point of maximum reduction in the flow speed was detected outside the rotor centerline.

La conception des turbines à gaz est aujourd'hui contrainte par des normes d'émissions de plus en plus draconiennes, couplées à l'urgente nécessité d'économiser les ressources en carburant fossile. Les choix technologiques adoptés pour répondre à ces exigences entraînent parfois l'apparition d'instabilités de combustion. Dans les chambres de combustion annulaires, ces instabilités prennent souvent la forme de modes azimutaux. Prédire ces modes reste un défi à l'heure actuelle et impose de considérer la totalité de la géométrie annulaire, ce qui n'est rendu possible, dans le domaine de la simulation numérique en mécanique des fluides, que par l'avènement très récent des supercalculateurs massivement parallèles. Dans ce travail de thèse, les modes azimutaux pouvant apparaître dans les chambres de combustion annulaires sont abordés avec plusieurs approches: un modèle analytique 1D, un solveur acoustique de Helmholtz 3D et enfin des Simulations aux Grandes Echelles. Combiner ces méthodes permet une meilleure compréhension de la structure de ces modes et peut amener à considérer des solutions innovantes pour concevoir des chambres inconditionnellement stables. / Increasingly stringent regulations and the need to tackle rising fuel prices have placed great emphasis on the design of aeronautical gas turbines. This drive towards innovation has resulted sometimes in new concepts being prone to combustion instabilities. Combustion instabilities arise from the coupling of acoustics and combustion. In the particular field of annular combustion chambers, these instabilities often take the form of azimuthal modes. To predict these modes, one must consider the full combustion chamber, which, in the numerical simulation domain, remained out of reach until very recently and the development of massively parallel computers. In this work, azimuthal modes that may develop in annular combustors are studied with different numerical approaches: a low order model, a 3D Helmholtz solver and Large Eddy Simulations. Combining these methods allows a better understanding of the structure of the instabilities and may provide guidelines to build intrinsically stable combustion chambers.

Cette thèse se focalise sur la compréhension et la prédiction des instabilités de combustion dans les systèmes à haute pression. Elle s'oriente autour de la simulation numérique d’un banc d'essai, opéré à l'université de Purdue, comprenant un injecteur caractéristique des moteurs-fusées et dont les propriétés acoustiques peuvent varier à l'aide d’un tube d'injection mobile. Une méthode d'initialisation et d'allumage pour les calculs LES de chambres de combustions terminées par une tuyère est présentée. Un point de fonctionnement instable est choisi pour étudier le mécanisme de l'instabilité. Les simulations sont comparées aux résultats expérimentaux en terme de fréquence et structure du mode instable. La fonction de transfert de flamme est calculée à l'aide du modèle n − τ puis implémentée dans un solveur acoustique (ne résolvant que les perturbations acoustiques à partir de l'équation de Helmholtz en écoulement réactif). Différents modèles d'impédance de tuyère, extraits de la littérature, sont comparés et leurs impacts sur les résultats de stabilité sont analysés. Le théorème d’impédance translatée est implémenté dans le solveur acoustique pour analyser, à faible coût de calcul, l’influence de la variation de la longueur du tube d'injection. Des écarts entre les fréquences prédites et celles trouvées expérimentalement subsistent mais la carte de stabilité de l’expérience est bien reproduite. / This thesis focuses on the understanding and the prediction of combustion instability in high-pressure devices. A model rocket combustor, tested experimentally at Purdue University, with continuously variable acoustic properties, thanks to a variable-length injector tube, is simulated. A method to initialize and ignite Large-Eddy-Simulation (LES) calculation of combustion chamber surrounded by nozzle is proposed. An unstable operating point is then chosen to investigate the mechanism of the instability. The simulations are compared to experimental results in terms of frequency and mode structure. The flame transfer function is calculated using the n − τ model to feed an acoustic solver which solves only the acoustic perturbation using a Helmholtz equation in reacting flows. The importance of the modeling of the nozzles impedance is studied through the main theories in the literature. The impedance translation theorem is implemented in the acoustic solver to analyze at low cost the influence of the variation of the injector tube. Despite differences in frequency of the instability, the stability map of the experiment is well reproduced.

Pratiquement toutes les chambres de combustion présentent des instabilités. Par conséquent, il est nécessaire de mieux les comprendre afin de les contrôler. Une possibilité est de simuler l’écoulement réactif à l’intérieur d’une chambre de combustion grâce à la Simulation aux Grandes Echelles (SGE). Cependant la SGE est très coûteuse en terme de capacité de calcul. Une autre possibilité est de réduire la complexité du problème à une simple équation d’onde thermoacoustique (équation dite de Helmholtz), qui peut être résolue en fréquence comme un problème aux valeurs propres. Le couplage entre l’acoustique et la flamme est alors prise en compte au travers des modèles appropriés. Le principal problème de cette méthode est qu’elle repose sur l’hypothèse d’un nombre de Mach nul. Tous les phénomènes liés à l’écoulement moyen sont donc négligés. La présente thèse propose une nouvelle stratégie pour prendre en compte certains effets de l’écoulement dans un contexte à Mach nul. Dans une première partie, la manière la plus judicieuse d’imposer un élément présentant un écoulement très rapide est étudiée. La seconde partie se focalise sur le couplage entre l’acoustique et les hétérogénéités de température qui sont générées par la flamme et naturellement convectées par l’écoulement moyen. Ce phénomène est important car il est responsable du bruit indirect de combustion qui peut conduire à une instabilité thermoacoustique. Un nouveau type de condition limite (DECBC) est proposé afin de prendre en compte ce mécanisme dans un contexte de résolution de l’équation de Helmholtz à Mach nul. Dans la dernière partie, une chambre de combustion aéronautique présentant une instabilité mixte acoustique/entropique est étudiée. Le bénéfice des méthodes développées dans la présente thèse est testé et comparé à des calculs avec la SGE. Il est montré que les calculs avec un solveur de Helmholtz peuvent reproduire une instabilité de combustion complexe, et que cet outil s’avère avoir le potentiel pour prédire les instabilités afin de concevoir de nouvelles chambres de combustion. / Virtually all combustion chambers are subject to instabilities. Consequently there is a need to better understand them so as to control them. A possibility is to simulate the reactive flow within a combustor with the Large-Eddy Simulation (LES) method. However LES results come at a tremendous computational cost. Another route is to reduce the complexity of the problem to a simple thermoacoustic Helmholtz wave equation, which can be solved in the frequency domain as an eigenvalue problem. The coupling between the flame and the acoustics is then taken into account via proper models. The main drawback of this latter methodology is that it relies on the zero-Mach number assumption. Hence all phenomena inherent to mean flow effects are neglected. The present thesis aims to provide a novel strategy to introduce back some mean flow effects within the zero-Mach number framework. In a first part, the proper way to impose high-speed elements such as a turbine is investigated. The second part focuses on the coupling between acoustics and temperature heterogeneities that are naturally generated at the flame and convected downstream by the flow. Such phenomenon is important because it is responsible for indirect combustion noise that may drive a thermoacoustic instability. A Delayed Entropy Coupled Boundary Condition (DECBC) is then derived in order to account for this latter mechanism in the framework of a Helmholtz solver where the baseline flow is assumed at rest. In the last part, a realistic aero-engine combustor that features a mixed acoustic/entropy instability is studied. The methodology developed in the present thesis is tested and compared to LES computations. It is shown that computations with the Helmholtz solver can reproduce a complex combustion instability, and that this latter methodology is a potential tool to design new combustors so as to predict and avoid combustion instabilities.

L'augmentation globale du trafic aérien mondial est susceptible de fortement augmenter les nuisances sonores induites autour des aéroports. De fait, des initiatives à diverses échelles sont prises afin de réduire le bruit des aéronefs. Le "bruit de cellule" de l'avion (airframe noise) est un contributeur majeur du bruit total en phase d'approche, dont environ un tiers est émis par le train d'atterrissage. La simulation aux grandes échelles (LES) des équations de Navier-Stokes compressibles a été identifiée comme une candidate pour la prédiction du bruit de train d'atterrissage. Elle permet de résoudre l'écoulement intrinsèquement instationnaire autour du train, et de capturer les différentes échelles des structures turbulentes mise en jeu dans le sillage. En outre, une approche non-structurée est retenue afin de traiter des géométries réalistes et complexes de trains d'atterrissage. Répondant à ces critères, le logiciel AVBP est choisi pour de telles simulations. Une analogie de Ffowcs-Williams et Hawkings est ensuite utilisée pour propager les sources acoustiques en champ lointain. La détermination des paramètres de simulation optimaux ainsi qu'une validation du code sont effectuées sur le cas académique du "barreau-profil" réalisé à l'Ecole Centrale de Lyon. Un bon accord des résultats aérodynamiques et acoustiques est observé avec la base de donnée expérimentale ainsi qu'avec d'autres résultats numériques. Une analyse des sources acoustiques a également été mise en oeuvre, mettant en évidence une contribution significative du cylindre aux angles rasants ainsi que des phénomènes interférentiels entre les émissions du cylindre et du profil, constructifs ou destructifs selon l'angle d'observation. Le code AVBP est ensuite utilisé pour la simulation des écoulements autour des trois trains simplifiés de la base de donnée LAGOON, supportée par Airbus. Un bon accord est globalement obtenu pour les résultats aérodynamiques et acoustiques pour les trois géométries de train, et les tendances des résultats expérimentaux lors de le complexification de la géométrie sont retrouvés à l'aide de la LES. Une étude des sources acoustique est ensuite réalisée. Les contributions respectives de chaque élément des trains sont identifiées, et une attention particulière est portée sur les tons de cavité. La structure de ces derniers est clairement mise en évidence à l'aide d'un solveur des équations d'Helmholtz (AVSP), expliquant également l'absence de tons pour les micros en survol. Les effets de l'ajout d'éléments supplémentaires pour les deux dernières configurations de train sont également étudiés, et il est montré en particulier qu'un bruit d'interaction, similaire à celui observé dans le cas "barreau-profil", entre le sillage de la barre de traction et l'essieu du train explique la prédominance d'un ton supplémentaire dans les spectres en champ lointain sur les tons de cavité. / Global air traffic growth is most likely to increase the community noise around airports. Consequently, several initiatives have been taken to reduce aircraft noise. Airframe noise has been identified as a major contributor to aircraft noise at approach, and one third is emitted by the landing gear. Compressible Large Eddy Simulation (LES) is a promising candidate for landing gear noise prediction. It is able to solve the unsteady flow around the landing gear, and also to capture the broad range of turbulent structures in the wake. Unstructured meshing approach has been selected to handle industrial-like and complex landing gear geometries. The AVBP software fulfills these criteria, and thus has been chosen to conduct the simulations. A Ffowcs-Williams and Hawkings analogy is also used to propagate the noise sources to the far-field. The best practices determination and the code validation are carried out on the academic "rodairfoil" configuration, which has been experimentally studied at Ecole Centrale de Lyon. A good agreement is observed between the numerical aerodynamic and acoustic results and the experimental database, as well with the other numerical simulations from the literature. A noise source analysis is also carried out, showing a significant contribution of the rod at grazing angles, as well as interferences between the rod and airfoil contributions, either being constructives or destructives depending on the observation angle. Then, the AVBP code is used to simulate unsteady flows around the three landing gear configurations of the LAGOON database, supported by Airbus. A good agreement is found for the aerodynamic and acoustic results for the three geometries, and experimental results tendencies are retrieved with LES while increasing the landing gear complexity. A noise source analysis leads to the identification of the noise source contributions of each landing gear element. A special attention is also paid on the cavity modes. Their structure, exhibited with a Helmholtz equation solver (AVSP), explains the absence of tone in the flyover direction. For the second and third geometries, it is also shown that interactions between the tow bar wake and the axle, similarly to the ones observed in the rod-airfoil benchmark, lead to an additional tone in the far-field spectra, which prevails on the initial cavity modes.

Pour améliorer le rendement des moteurs essence, une méthode efficace est le downsizing qui consiste en la diminution de la cylindrée moteur compensée par l’ajout d’un compresseur pour maintenir la puissance. Lorsque le niveau de downsizing est trop important les fortes pression et températures rencontrées favorisent l’apparition de phénomènes d’auto-allumage de type cliquetis ou rumble néfastes pour l’intégrité du moteur. Ce type de phénomène, aujourd’hui encore mal compris, constitue une limite à l’utilisation du downsizing. Dans cette thèse la Simulation aux Grandes Echelles est utilisée pour étudier ce type de combustion dite anormale. L’objectif est de proposer une méthodologie numérique capable de reproduire leurs apparitions pour en étudier les mécanismes. L’auto-allumage est un mode de combustion sensible aux variations des conditions thermodynamiques locales. Des méthodes numériques précises et des modèles appropriés, en particulier pour la thermique paroi doivent donc être utilisés. La première partie de ce manuscrit présente la méthodologie numérique proposée et en particulier deux aspects développés lors de cette thèse: un modèle d’auto-allumage qui permet de reproduire le délai d’auto-allumage des gaz frais avec un schéma cinétique réduit et une méthodologie de couplage entre la chambre de combustion et la culasse permettant de définir des champs de températures paroi réalistes. La seconde partie de ce manuscrit présente les résultats de deux études numériques reproduisant certains points de fonctionnement d’un moteur expérimental. La première étude est réalisée à l’aide de modèles de combustion de la littérature et vise à reproduire le comportement expérimental pour diverses variations paramétriques influant sur la combustion. La seconde étude est réalisée à l’aide des modèles développés dans cette thèse afin d’étudier l’impact de la thermique paroi dans les mécanismes d’apparition des combustions anormales. / In order to improve the efficiency of gasoline engines, one efficient solution resides in engine downsizing which consists in the diminution of the engine size with the adjunction of a compressor to keep the power output. When the downsizing level is important, the high pressure and temperature levels promote auto-ignition phenomena such as knocking or rumble that can damage the engine. This kind of combustion, still misunderstood, is a limit to further use downsizing. In this thesis, Large Eddy Simulation is used to study this kind of abnormal combustions. The objective is to propose a numerical methodology able to reproduce its apparition and to understand its mechanisms. Auto-ignition is a combustion regime very sensitive to the variations of local thermodynamic conditions. Precise numerical methods and appropriate models, especially for thermal boundary conditions must be used. The first part of this manuscript presents the proposed numerical methodology and in particular two aspects implemented during this thesis: an auto-ignition model that permits to reproduce auto-ignition delays with reduced kinetic schemes and a coupling methodology between combustion chamber and cylinder head in order to obtain realistic temperature fields for the boundary conditions. The second part of this manuscript presents the results of two numerical studies that reproduce some operating points from an experimental engine database. The first study is performed using combustion models from the literature and aims at reproducing experimental behavior for various parametric variations impacting the combustion. The second study is performed thanks to the numerical models implanted in this thesis in order to evaluate the impact of the thermal boundary conditions on the mechanisms leading to abnormal combustions.

Nous nous intéressons dans cette thèse au développement et à la mise en oeuvre de schémas numériques Volumes Finis d’ordre élevé pour des maillages non-structurés. Il s’agit de mettre en place les ingrédients numériques pour réaliser des simulations aux grandes échelles avec le code numérique elsA. Les schémas numériques proposés sont basés sur une approche directionnelle, afin de limiter le coût CPU et de réduire la molécule de points. La partie convective du schéma numérique doit être d’ordre élevé. L’ordre élevé est obtenu en utilisant différents gradients sur un stencil prédéféni utilisant 4 cellules. Deux gradients sont utilisés pour la partie convective : le gradient GreenGauss et le gradient “UIG”. Pour la partie diffusive, le gradient “UIG” est utilisé. Ce gradient a été développé durant la thèse et permet d’avoir un gradient moyen d’ordre 2 sur chaque interface. Ce gradient a été étudié et validé sur différents cas-tests. Les schémas numériques d’ordre élevé ont été analysés théoriquement avec des analyses d’ordre et de stabilité. Il a été montré que ces schémas peuvent atteindre l’ordre 5 sur des hexaèdres et l’ordre 3 sur des triangles équilatéraux. Suite à cette analyse, les différents schémas ont été d’abord testés en 1D sur un cas classique d’advection, puis ont été validés sur le cas de convection du vortex isentropique. / This study will present the development and results of high-order Finite Volume schemes for unstructured grids. The goal is to prepare numerical tools to perform Large Eddy Simulations with the indutrial solver elsA. These numerical schemes are based on a directional approach in order to limitate the CPU cost and reduce the stencil. The convective part of the scheme needs to be high order and this is obtained by the use of gradients on a four-cell stencil. Two gradients are used for the convective part, the Green-Gauss gradient and the “UIG” gradient. For the diffusive part, the “UIG” gradient is used. It was developped during this study and allows to recover a secondorder accurate scheme. This gradient was validated theorically and numerically on some test cases. High order numerical schemes were studied theorically with order and frequency analysis. It was shown that these schemes are fifth-order accurate on regular hexaedral elements and third-order accurate on equilateral triangles. Following this analysis, these schemes were tested in 1D on an advection test case and were then validated on the convection of an isentropic vortex.