La modélisation par simulation des grandes échelles (Large-Eddy Simulation - LES) des processus physiques régissant la couche limite atmosphérique (CLA) demeure complexe de part la difficulté des modèles à capter l'évolution de la turbulence entre différentes conditions de stratification. De ce fait, l'étude LES du cycle diurne complet de la CLA comprenant des situations convectives la journée et des conditions stables la nuit est très peu documenté. La simulation de la couche limite stable où la turbulence est faible, intermittente et qui est caractérisée par des structures turbulentes de petite taille est tout particulièrement compliquée. En conséquence, la capacité de la LES à bien reproduire les conditions météorologiques de la CLA, notamment en situation stable, est étudiée à l'aide du code de mécanique des fluides développé par EDF R&D, Code_Saturne. Dans une première étude, le modèle LES est validé sur un cas de couche limite convective quasi stationnaire sur terrain homogène. L'influence des modèles sous-maille de Smagorinsky, Germano-Lilly, Wong-Lilly et WALE (Wall-Adapting Local Eddy-viscosity) ainsi que la sensibilité aux méthodes de paramétrisation sur les champs moyens, les flux et les variances est discutées. Dans une seconde étude le cycle diurne complet de la CLA pendant la campagne de mesure Wangara est modélisé. L'écart aux mesures étant faible le jour, ce travail se concentre sur les difficultés rencontrées la nuit à bien modéliser la couche limite stable. L'impact de différents modèles sous-maille ainsi que la sensibilité au coefficient de Smagorinsky ont été analysés. Par l'intermédiaire d'un couplage radiatif réalisé en LES, les répercussions du rayonnement infrarouge et solaire sur le jet de basse couche nocturne et le gradient thermique près de la surface sont exposées. De plus l'adaptation de la résolution du domaine à l'intensité de la turbulence et la forte stabilité atmosphérique durant l'expérience Wangara sont commentées. Enfin un examen des oscillations numériques inhérentes à Code_Saturne est réalisé afin d'en limiter les effets

This thesis focused on investigating flow, coherent structures, and mechanisms leading to scour around a series of three submerged groynes in a high curvature (radius of curvature (R)/channel width (B)=1.5) channel bend using a Large Eddy Simulation Numerical (LES) model. Flow was investigated during both an initial and a later stage of scour. The results showed that the groynes appeared effective in keeping the main core of high streamwise velocity away from the outer bank wall in the region where they were installed, although high potential still existed for local scour around the groynes. During the initial stage of scour, horseshoe vortices (HVs) showed the greatest propensity to induce scour immediately upstream of the groyne tips. During the later stage of scour, the HV in front of the first upstream groyne (G1) induced very high mean pressure fluctuations on the outer bank wall. Scour was also of very great concern around the tip of G1 due to severe mean bed pressure fluctuations. Downstream of the groyne field, the presence of a counter-rotating outer bank cell was capable of endangering the stability of the outer bank. The second focus of this thesis was to investigate flow in a 135° channel bend using both Reynolds Averaged Navier Stokes (RANS) and LES numerical models. The RANS study examined the effects of curvature ratio (R/B), and aspect ratio (B/H, where H is the inlet flow depth), on secondary circulation strength, and bed shear stresses. The study revealed that a decrease in R/B was associated with an increase in secondary circulation strength and peak bed shear stress. A change in B/H also substantially affected cross stream circulation strength. The LES study was conducted in a 135° (R/B = 1.5) bend flume with a fixed bed corresponding to near equilibrium scour conditions, and the results were compared to a similar high curvature 193° bend numerical study. Inner bank vortices and shear layers were present in both cases although their characteristics were substantially different. Distributions of boundary friction velocities, and turbulence were also quite different for each case.

Le contexte général de cette thèse s'inscrit dans le cadre de la modélisation eulérienne instationnaire des écoulements turbulents anisothermes gaz - particules. La modélisation de ces écoulements est cruciale pour de nombreuses applications industrielles et pour la prédiction de certains phénomènes naturels. Par exemple, la combustion diphasique dans les moteurs automobiles et aéronautiques est précédée par l'injection et la dispersion de carburant liquide dans la chambre de combustion. Les phénomènes mis en jeu exigent alors une prédiction locale tenant compte du caractère instationnaire de l'écoulement turbulent et de la présence de géométries complexes. De plus, de nombreuses études expérimentales et numériques récentes ont mis en évidence le rôle prépondérant de l'inertie des particules sur les mécanismes de dispersion et de concentration préférentielle en écoulement turbulent. Ceci rend donc indispensable la prise en compte de ces mécanismes dans la modélisation diphasique. Au cours de ce travail de thèse, une approche eulérienne locale et instantanée a été développée pour prédire les écoulements gaz-particules anisothermes et turbulents. Elle est basée sur l'approche statistique du Formalisme Eulérien Mésoscopique (MEF) introduite par Février et al. (JFM, 2005). Cette approche a été ici étendue aux variables thermiques pour la prise en compte du caractère anisotherme de l'écoulement. Cette approche a été ensuite utilisée dans le cadre de la méthode des moments (Kaufmann et al., JCP, 2008), et un système d'équations locales et instantanées pour la phase dispersée a été proposé. La modélisation au premier ordre exige la fermeture des moments de second ordre apparaissant dans les équations de la quantité de mouvement et de l'énergie. La proposition de telles relations constitutives fait l'objet d'une partie de la thèse. Afin de fournir une méthode capable de prédire le comportement local, instantané et anisotherme de la phase dispersée dans des configurations `a une échelle réaliste, les équations pour la phase dispersée ont été filtrées et une modélisation aux grandes échelles (LES) est effectuée. Cette modélisation étends, par la prise en compte des variables thermiques, le travail de Moreau et al. (FtaC, 2010) sur l'approche LES Euler-Euler en conditions isothermes. L'approche complète est enfin appliquée aux résultats de simulation numérique d'un jet plan turbulent gazeux froid, chargé en particules, dans une turbulence homogène isotrope chaude monophasique. / The aim of this thesis is to provide an Eulerian modeling for the dispersed phase interacting with unsteady non-isothermal turbulent flows. The modeling of these flows is crucial for several industrial applications and for predictions of natural events. Examples are the combustion chambers of areo engines where the combustion is preceded by the injection and dispersion of liquid fuel. The prediction of such phenomena involves a local modeling of the mixture for taking into account the unsteady behavior of the turbulent flow and the presence of complex geometries. Moreover, many experimental and numerical studies have recently highlighted the significant role of the particle inertia on the mechanisms of dispersion and preferential concentration. Accounting for such mechanisms is therefore essential for modeling the particle-laden turbulent flows. In this thesis, a local and instantaneous Eulerian approach able to describe and to predict the local behavior of inertial particles interacting with non-isothermal turbulent flows has been developed. It is based on the statistical approach known as Mesoscopic Eulerian formalism (MEF) introduced by Février et al. (JFM, 2005). The statistical approach has been extended to the thermal quantities in order to account for the non-isothermal conditions into the modeling. This formalism is then used in the framework of the moment approach (Kaufmann et al., JCP, 2008) and a system of local and instantaneous equations for the non-isothermal dispersed phase has been suggested. The first order modeling requires to close second-order moments appearing in momentum and energy equations. The proposal of such constitutive relations makes the object of a part of this study. In order to provide an Eulerian approach usable in real configurations at industrial scale, the equations of the dispersed phase are filtered and the approach developed in the framework of the Large-Eddy Simulations. From the work of Moreau et al. (FTaC, 2010), the Eulerian-Eulerian LES approach is then extended to non-isothermal conditions. The whole modeling is then a priori tested against numerical simulations of a cold planar turbulent particle-laden jet crossing a homogeneous isotropic decaying hot turbulence.

La réduction des émissions de polluants et la diminution de la consommation sont deux challenges fortement liés auxquels les constructeurs automobiles doivent faire face tout en maintenant les performances des moteurs. Les nouvelles stratégies telle que la réduction de la cylindrée associée à une optimisation de la boucle d'air (forte suralimentation et recirculation de gaz brûlés) possèdent ce potentiel. Cependant elles affectent la stabilité du moteur en favorisant les variations cycle à cycle (VCC) qui correspondent à une fluctuation de la combustion d'un cycle sur l'autre. L'objectif de cette thèse est de développer une méthodologie s'appuyant sur la Simulation aux Grandes Echelles (SGE) capable de prédire les VCC d'un moteur à allumage commandé. Les prédictions sont validées grâce à une base de données expérimentale conçue à l'IFP qui leur est dédiée. Une approche graduelle est employée : l'outil numérique est tout d'abord évalué sur une configuration simplifiée de moteur à piston sans combustion, puis appliqué à un moteur à allumage commandé entraîné pour valider la prédiction de l'aérodynamique interne. Sur cette dernière configuration le couplage avec le modèle de combustion DTFLES est rajouté pour simuler deux points de fonctionnement réactifs. Chacune de ces simulations intègre un ou plusieurs points de modélisation (les tétraèdres en maillage mobile, les modèles de choc et d'allumage, et la cinétique chimique) au préalable testés sur des configurations académiques. Ce travail de recherche montre que l'approche SGE, dans un contexte de calcul massivement parallèle, est un outil prometteur dans l'étude des VCC d'un moteur à allumage commandé de conception récente. / A major challenge for the development of internal combustion engines is to improve fuel economy and to reduce pollutant emissions while maintaining or enhancing engine performances. New strategies using downsizing with high levels of exhaust gas recirculation have this potential, but can impact on the combustion stability and trigger high cycle-to-cycle variations (CCV). The objective of this thesis is to set a methodology based on Large Eddy Simulation (LES) to study CCV of a Spark-Ignition (SI) engine. A gradual approach is used : the numerical tool is first evaluated on a motored axisymmetric piston-cylinder assembly, and then applied on a motored SI engine to validate the in-cylinder aerodynamic predictions. On this last configuration, the coupling with the turbulent combustion model DTFLES is added to simulate two operating points of the indirect injection engine mode. Each of these simulations takes into account one or several modeling key points (tetrahedra with moving mesh, the modelings of shock and ignition, and chemical kinetics) previously tested in academic configurations. This research work shows that LES approach, in the context of massively parallel computing, can be used to study the CCV of a realistic SI engine.

Les turbines à gaz aéronautiques doivent satisfaire des normes d'émissions polluantes toujours en baisse. La formation de polluants est directement liée à la qualité du mélange d’air et de carburant en amont du front de flamme. Ainsi, leur réduction implique une meilleure prédiction de la formation du spray et de son interaction avec l'écoulement turbulent gazeux. La simulation aux grandes échelles (SGE) semble un outil numérique approprié pour étudier ces mécanismes. Le but de cette thèse est d’évaluer l'impact de modèles d'injection simplifiés sur la SGE de l'écoulement diphasique évaporant d’une configuration complexe. La configuration cible choisie est un brûleur aéronautique installé sur le banc expérimental MERCATO. Le banc expérimental est equipé d’un système d’injection d'air vrillé et d’un système d'injection liquide avec un atomiseur pressurisé swirlé représentatif de foyers aéronautiques réels. Dans un premier temps, un modèle d'injection simplifié pour atomiseurs pressurisés swirlés négligeant les effets de l'atomisation sur la dynamique du spray est présenté. L'objectif principal de ce modèle réside dans la reproduction de conditions d’injection similaires pour des traitements Eulériens et Lagrangiens de la phase particulaire. Dans un second temps, la composante Lagrangienne de ce modèle d'injection est combinée à un modèle d'atomisation secondaire de la litérature pour permettre une prise en compte partielle des phénomènes de pulvérisation liquide. Les SGE de l'écoulement diphasique évaporant de la configuration MERCATO présentées comportent deux aspects. Premièrement, différents modèles d’injection sont évalués pour quantifier leur impact sur la dynamique de la phase particulaire. Deuxièmement, une comparaison de simulations Euler-Euler et Euler-Lagrange reposant sur un modèle d'injection unifié est effectuée. / Aeroautical gas turbines need to satisfy growingly stringent demands on pollutant emission. Pollutant emissions are directly related to the quality of fuel air mixing prior to combustion. Therefore, their reduction relies on a more accurate prediction of spray formation and interaction of the spray with the gaseous turbulent flowfield. Large-Eddy Simulation (LES) seems an adequate numerical tool to predict these mechanisms. The objective of this thesis is to evaluate the impact of simplified injection methods on the LES of the evaporating two-phase flow inside a complex geometry. The chosen target configuration is an aeronautical combustor installed on the MERCATO test-rig. The experimental setup includes an air-swirler injection system and a pressureswirl atomizer typical of realistic aeronautic combustors. In a first step, a simplified injection model for pressure swirl atomizers neglecting the impact of liquid disintegration on spray dynamics is presented. The main objective of this model lies in the reproduction of similar injection conditions for Eulerian and Lagrangian representations of the dispersed phase. In a second step, the Lagrangian injection method is combined to a secondary breakup model of the literature to partly account for the liquid disintegration process. The presented LES’s of the evaporating two-phase flow inside the MERCATO geometry consider two different aspects. First, the impact of injection modeling on spray dynamics is assessed. Second, Euler-Euler and Euler-Lagrange simulations relying on the common simplified injection model are compared.

Une part importante du bruit généré par les moteurs d'avion est liée à la combustion. Afin de réduire cette source de bruit, une compréhension fine des phénomènes associés est nécessaire. Deux mécanismes générateurs de bruit, et ayant pour origine la combustion, ont été identifié dans les moteurs d'avions dans les années 1970: un premier mécanisme dit direct, qui est lié directement à un dégagement de chaleur instationnaire, et l'autre dit indirect qui est lié aux interactions entre les étages de turbine et les fluctuations de température en sortie de chambre, également produites par la combustion. Des méthodes analytiques et des simulations numériques sont utilisées ici à la fois pour montrer l'importance du bruit de combustion indirect par rapport au bruit direct, et pour donner des limites de validité des approches analytiques qui sont basées sur l'hypothèse de tuyère compacte. Trois configurations différentes sont étudiées dans un premier temps: une tuyère quasi-1-D, une tuyère axi-symétrique 2-D, ainsi qu'une aube de turbine 2-D. Finalement, un secteur de chambre de combustion 3-D réelle (SNECMA) est calculé à l'aide de la Simulation aux Grandes Echelles. Les fluctuations en sortie du brûleur sont utilisées pour évaluer le bruit total généré par la combustion (direct et indirect) à l'aide des approches analytiques précédemment étudiées. / An important part of the noise generated by aero-engines is caused by the combustion. To decrease this source of noise, a precise comprehension of the underlying phenomenon is required. Two different mechanisms generating noise in aero-engines because of the combustion have been identified in the 1970’s: the direct mechanism that is directly related to the unsteady heat release, and the indirect one that is caused by the interactions between the turbine stages and the temperature fluctuations also produced by the combustion. Analytical methods and numerical simulations are used here both to show the importance of the indirect combustion noise compared with the direct one, and to provide some validity limits of compact nozzles analytical approaches. Three different configurations dealing with indirect noise are investigated: quasi-1- D nozzle, axisymmetric 2-D nozzle and 2-D turbine blade. Finally, an actual 3-D combustion chamber sector (SNECMA) is addressed with Large-Eddy Simulation. Fluctuations at the outlet of the combustor are used to compute the overall noise caused by the combustion (direct and indirect), by means of the investigated analytical models.

Comme le groupe de symétries de Lie des équations aux dérivées partielles représentent les propriétés physiques intrinsèques contenues dans les équations, il offre un outil efficace pour étudier et modéliser les phénomènes physiques. Ainsi, dans cette thèse, on se propose d’appliquer la théorie du groupe de symétries de Lie à la modélisation des écoulements anisothermes.On calcule alors des lois de paroi, et, plus généralement des lois d’échelle, pour la vitesse et la température dans le cas d’un écoulement parallèle. En fait, ces lois d’échelle se révèlent être simplement des solutions auto-similaires des équations de Navier-Stokes moyennées par rapport aux symétries des équations.Ensuite, par l’approche de la théorie des groupes de Lie, on construit une classe de modèles de sous-maille qui sont invariants par les symétries des équations de Navier-Stokes anisothermes.Ces modèles ont l’avantage de respecter les propriétés physiques des équations qui sont contenues dans les symétries. De plus, par cette approche, le modèle de flux de chaleur apparaît naturellement,sans qu’on ait besoin de faire appel à la notion de nombre de Prandtl de sous-maille,ce qui augmente la portée de ces modèles par rapport à la plupart des modèles existants. Par ailleurs, le comportement proche de la paroi de certains des modèles proposés est étudié. Enfin,des tests numériques en convection naturelle et en convection mixtes sont effectués. / Since the Lie group of a given partial differential equation, represent the intrinsic physical propertiesof the equation, it gives a strong tool for modeling its physical phenomenas. The mainpurpose of this thesis, is to apply the Lie group theory, in modeling non-isothermal flows. Therefore,we calculate wall laws and more generally scaling laws for the velocity and the temperatureof a parallel flow. In fact, these scaling laws are simply self-similar solutions of the Navier-Stokesequations averaged with respect to their symmetry.The approach of the Lie group theory, leads to a class of sub-grade models which are invariantvia the symmetries of the non-isothermal Navier-Stokes equations. These models respectthe physical properties contained in these symmetries. Moreover, via this approach the heat flowmodel appears naturally in this class, without introducing the notion of the Prandlt number,which is not the case for any other existing model. From the other side, the behavior near thewall of particular models in this class, is studied. Finally, numerical tests are done in both casesof the natural convection and the mixed one.

Le développement de nouvelles technologies pour le transport aérien moins polluant est de plus en plus basé sur la simulation numérique, qui nécessite alors une description fiable de la chimie. Pour la plupart des carburants, la description de la combustion nécessite des mécanismes détaillés mais leur utilisation dans une simulation numérique de combustion turbulente est limitée par le coût calcul. Des mécanismes cinétiques réduits et des méthodes de tabulation ont été proposés pour surmonter ce problème. Ces descriptions chimiques simplifiées ayant été développées dans le cadre de configurations laminaires, cette thèse propose de les évaluer dans des configurations turbulentes: une DNS de flamme prémélangée méthane/air de type Bunsen et une LES d’un brûleur expérimental. Les mécanismes sont analysés en termes de structure de flamme, paramètres de flamme globaux, longuer de flamme, prediction des concentrations en espèces majoritaires et des émissions polluantes. Une méthodologie pour évaluer a priori la capacité d’un mécanisme à prédire correctement des phénomènes chimiques tridimensionnels est proposée en se basant sur les résultats de flammes laminaires monodimensionnelles non étirées et étirées. Il ressort que, d’une part, pour construire un mécanisme réduit, il est nécessaire de faire un compromis entre coût calcul, robustesse et qualité des résultats. D’autre part, la qualité des résultats de DNS et LES de configurations tridimensionnelles turbulentes peut être anticipée par une analyse du comportement des schémas réduits dans des configurations simplifiées de flammes monodimensionnelles laminaires non étirées et étirées. / A growing need for numerical simulations based on reliable chemistries has been observed in the last years in order to develop new technologies which could guarantee the reduction of the enviromental impact on air transport. The description of combustion requires the use of detailed kinetic mechanisms for most hydro-carbons. Their use in turbulent combustion simulation is still prohibitive because of their high computational cost. Reduced chemistries and tabulation methods have been proposed to over-come this problem. Since all these reductions have been developed for laminar configurations, this thesis proposes to evaluate their performances in simulations of turbulent configurations such as a DNS of a premixed Bunsen methane/air flame and a LES of an experimental PREC-CINSTA burner. The mechanisms are analysed in terms of flame structure, global burning parameters, flame length, prediction of major species concentrations and pollutant emissions. An a priori methodology based on one-dimensional unstrained and strained laminar flames to evaluate the mechanism capability to predict three-dimensional turbulent flame features is therefore proposed. On the one hand when building a new reduced scheme, its requirements should be fixed compromising the computational cost, the robustness of the chemical description and the desired quality of results. On the other hand, the quality of DNS or LES results in three-dimensional configurations could be anticipated testing the reduced mechanism on laminar one-dimensional premixed unstrained and strained flames.

Ces cinquantes dernières années, la majorité des paramètres de conception des moteurs cryotechniques ont été ajustés en l'absence d'une compréhension détaillée des phénomènes de combustion, en raison des limites des diagnostiques expérimentaux et des capacités de calcul. L'objectif de cette thèse est de réaliser des simulations numériques instationnaires d'écoulements réactifs transcritiques de haute fidélité, pour permettre une meilleure compréhension de la dynamique de flamme dans les moteurs cryotechniques et finalement guider leur amélioration. Dans un premier temps, la thermodynamique gaz-réel et son impact sur les schémas numériques sont présentés. Comme la Simulation aux Grandes Echelles (SGE) comporte des équations filtrées, les effets de filtrages induits par la thermodynamique gaz-réel sont ensuite mis en évidence dans une configuration transcritique type et un opérateur de diffusion artificiel, spécifique au gaz réel, est proposé pour lisser les gradients transcritiques en SGE. Dans un deuxième temps, une étude fondamentale du mélange turbulent et de la combustion dans la zone proche-injecteur des moteurs cryotechniques est menée grâce à la Simulation Numérique Directe (SND). Dans le cas non-réactif, les lâchers tourbillonnaires dans le sillage de la lèvre de l’injecteur jouent un rôle majeur dans le mélange turbulent et provoquent la formation de structures en peigne déjà observées expérimentalement dans des conditions similaires. Dans le cas réactif, la flamme reste attachée à la lèvre de l'injecteur, sans extinction locale, et les structures en peigne disparaissent. La structure de flamme est analysée et différents modes de combustion sont identifiés. Enfin, une étude de flamme-jet transcritique H2/O2, accrochée à un injecteur coaxial avec et sans retrait interne, est menée. Les résultats numériques sont d'abord validés par des données expérimentales pour l'injecteur sans retrait. Ensuite, la configuration avec retrait est comparée à la solution de référence sans retrait et à des données experimentales pour observer les effets de ce paramètre de conception sur l'efficacité de combustion. / In the past fifty years, most design parameters of the combustion chamber of Liquid Rocket Engines (LREs) have been adjusted without a detailed understanding of combustion phenomena, because of both limited experimental diagnostics and numerical capabilities. The objective of the present thesis work is to conduct high-fidelity unsteady numerical simulations of transcritical reacting flows, in order to improve the understanding of flame dynamics in LRE, and eventually provide guidelines for their improvement. First real-gas thermodynamics and its impact on numerical schemes are presented. As Large-Eddy Simulation (LES) involves filtered equations, the filtering effects induced by real-gas thermodynamics are then highlighted in a typical 1D transcritical configuration and a specific real-gas artificial dissipation is proposed to smooth transcritical density gradients in LES. Then, a Direct Numerical Simulation (DNS) study of turbulent mixing and combustion in the near-injector region of LREs is conducted. In the non-reacting case, vortex shedding in the wake of the lip of the injector is shown to play a major role in turbulent mixing, and induces the formation of finger-like structures as observed experimentally in similar operating conditions. In the reacting case, the flame is attached to the injector rim without local extinction and the finger-like structures disappear. The flame structure is analyzed and various combustion modes are identified. Finally, a LES study of a transcritical H2/O2 jet flame, issuing from a coaxial injector with and without inner recess, is conducted. Numerical results are first validated against experimental data for the injector without recess. Then, the recessed configuration is compared to the reference solution and to experimental results, to scrutinize the effects of this design parameter on combustion efficiency.

La Simulation aux Grandes Echelles (SGE) est de plus en plus utilisée dans les processus de développement et la conception des réacteurs aéronautiques industriels. L'une des raisons pour ce besoin résulte dans la capacité de la SGE à fournir des informations instantanées d'un écoulement turbulent augmentant la quantité des prédictions de la composition des gaz d'échappement. Ce manuscrit de thèse aborde deux sujets récurrents de la SGE. D'une part, les schémas numériques pour la SGE nécessitent certaines propriétés, notamment une précision élevée avec une diffusivité faible pour ne pas nuire aux modèles de turbulence. Afin de répondre à ce pré requis, une famille de schémas d'intégration temporelle d'ordre élevée est proposée, permettant de modifier la diffusion numérique du schéma. D'autre part, la SGE étant intrinsèquement instationnaire, elle est très consommatrice en temps CPU. De plus, une géométrie complexe prend beaucoup de temps de simulation même avec les super calculateurs d'aujourd'hui. Dans le cas particulier d'intérêt et souvent rencontré dans les applications industrielles, l'approche bas-Mach est constitue une alternative intéressante permettant de réduire le coût et le temps de retour d'une simulation LES. L'impact et la comparaison des formalismes compressible et incompressible sont toutefois rarement quantifiés, ce qui est proposé dans ce travail pour une configuration représentative d'un brûleur swirlé industriel mesuré au CORIA / Large-Eddy Simulation (LES) becomes a more and more demanded tool to improve the design of aero-engines. The main reason for this request stems from the constraints imposed on the next generation low-emission engines at the industrial development level and the ability for LES to provide information on the instantaneous turbulent flow field which greatly contributes to improving the prediction of mixing and combustion thereby offering an improved prediction of the exhaust emission. The work presented in this thesis discusses two recurring issues of LES. For one, numerical schemes for LES require certain properties, i.e. low-diffusion schemes of high order of accuracy so as not to interfere with the turbulence models. To meet this purpose in the context of fully unstructured solvers, a new family of high-order time-integration schemes is proposed. With this class of schemes, the diffusion implied by the numerical scheme become adjustable and built-in. Second, since fully unsteady by nature, LES is very consuming in terms of CPU time. Even with today's supercomputers complex problems require long simulation times. Due to the low flow velocities often occurring in industrial applications, the use of a low-Mach number solver seems suitable and can lead to large reductions in CPU time if comparable to fully compressible solvers. The impact of the incompressibility assumption and the different nature of the numerical algorithms are rarely discussed. To partly answer the question, detailed comparisons are proposed for an experimental swirled configuration representative of a real burner that is simulated by LES using a fully explicit compressible solver and an incompressible solution developed at CORIA