Cette thèse s’intéresse à la turbulence incompressible homogène et anisotrope, et plus particulièrement à l’effet d’une rotation solide et d’un champ magnétique uniforme et stationnaire. En plus des simulations numériques directes (DNS), nous présentons également un modèle synthétique de turbulence (Kinematic Simulation ou KS) construit à partir d’une superposition de modes de Fourier, et au sein duquel la dynamique linéaire (basée sur la Rapid Distorsion Theory) peut être incluse. Dans un premier temps, l’effet de la rotation solide est étudié avec un effort particulier porté sur les corrélations en deux-temps. Une comparaison entre DNS et KS est proposée dans le cas isotrope comme dans le cas en rotation. Dans le contexte aéroacoustique, on montre dans quelle mesure la rotation modifie l’émission acoustique d’une turbulence homogène. L’effet d’un champ magnétique est ensuite considéré et l’anisotropie de l’écoulement est étudiée en fonction du nombre de Reynolds magnétique. Il est montré dans quelle mesure la turbulence magnétohydrodynamique quasi-statique est similaire à la turbulence “deux-dimensions trois-composantes” du fait de la dissipation Joule anisotrope et comment l’écoulement restaure son isotropie lorsque le nombre de Reynolds magnétique croît. Enfin, l’effet couplé d’une rotation et d’un champ magnétique est considéré. Les propriétés énergétiques ainsi que l’anisotropie sont étudiés et une étude paramétrique en fonction du nombre d’Elsasser est proposée. / This thesis focuses on incompressible homogeneous anisotropic turbulence, with a particuliar interest on the effect of a solid-body rotation and a uniform steady magnetic field. In addition to Direct Numerical Simulations (DNS), we also propose a stochastic model of anisotropic turbulence (Kinematic Simulation or KS) based on a superposition of Fourier modes, in which linear dynamics (based on Rapid Distorsion Theory) is included. First, the effect of a solid body rotation is studied and we focus on two-time velocity correlations. Comparisons between KS and DNS are presented in the isotropic and rotating cases. We then apply these results to aeroacoutics and show how the rotation modifies the acoustic emission of homogeneous turbulence. Secondly, we focus on the effect of an imposed magnetic field and we describe the anisotropy of the flow in function of the magnetic Reynolds number. The quasi-static magnetohydrodynamic turbulence is very similar to the so-called 2D-3C turbulence due to anisotropic Joule dissipation and we show how the isotropy is restored as the magnetic Reynolds numbers increases. Finally, we consider the coupled effect of rotation and magnetic field. Energetic properties and anisotropy are studied and a dependance of the results on the Elsasser number is also proposed.

La finalité de cette thèse est de fournir une meilleure compréhension de la turbulence homogène anisotrope soumise à un forçage barocline. À cette fin, nous utilisons une approche numérique pseudo-spectrale basée sur la transformation de Rogallo. L'utilisation d'un tel algorithme nous permet de considérer une asymétrie des fonctions de probabilité en faveur des évènements négatifs est observée. Le lien entre la distribution de vorticité potentielle et celle d'un scalaire passif est également étudié. Il est montré qu'à faible nombre de Richardson, c'est le mode vortex (à vorticité potentielle nulle) qui contient les plus importantes fuctuations de scalaire. Un écoulement homogène dans les trois directions de l'espace. Plusieurs simulations numériques directes (DNS) sont effectuées dans un contexte assez proche des écoulements géophysiques que l'on retrouve entre autre dans la stratosphère, où un gradient constant de vitesse zonale vient se coupler à un gradient constant de densité dans un repère tournant. Les résultats obtenus s'articulent autour de quatre axes principaux. Tout d'abord, une étude linéaire à temps fini est présentée en vue de compléter les résultats existants sur la dynamique linéaire asymptotique. La solution linéaire est décomposée en une partie 'onde' (qui se propage) et une partie dite 'vortex'(stationnaire). L'étude analytique est complétée par un modèle synthétique de turbulence (Kinematic Simulation ou KS) basé sur la théorie de la distorsion rapide(RDT). Nous montrons qu'une distribution initiale non nulle de vorticité potentielle linéarisée peut conduire à d'importantes croissances transitoires. Ce résultat pourrait s'étendre à des modélisations du climat ou météorologique, où la distribution initiale de vorticité potentielle semble avoir autant d'importance que la distribution initiale de température ou de vitesse. Ensuite, nous consacrons une partie de notre étude à l'analyse paramétrique et à la stabilité de l'écoulement. Plusieurs DNS sont effectuées pour différents taux de rotation et stratification. Le diagramme de stabilité obtenu montre que pour de faibles taux de rotation, la limite de stabilité est identique à celle connue des écoulements sans rotation. À plus faible nombre de Rossby -- lorsque la baroclinicité devient importante -- la limite linéaire de stabilité Ri = 1 relative à l'instabilité symétrique est confirmée. La coexistance de l'instabilité barocline avec l'instabilité symétrique est également clarifiée. Une analyse énergétique détaillée mène à la conclusion suivante : la stratification doit être suffisamment importante (Ri ' 1) pour que l'instabilité barocline soit dominante i.e. que la conversion d'énergie potentielle soit la source principale d'énergie cinétique turbulente. Dans le cas contraire, l'instabilité symétrique -- qui tire son énergie de l'énergie cinétique de l'écoulement moyen et non de son énergie potentielle -- domine la dynamique de l'écoulement. Le troisième axe d'étude concerne la turbulence à proprement parler. En conséquence de l'ajustement géostrophique, le vent thermique force la turbulence d'une manière naturelle, en opposition à d'autres méthodes de forçage stochastique. L'émergence de structures dans le contexte barocline est approfondie. Des statistiques Euleriennes sont présentées afin de fournir une caractérisation fine de l'anisotropie de l'écoulement. Enfin, nous étendons notre étude à la caractérisation de la vorticité potentielle turbulente. Les fonctions de probabilité de la vorticité potentielle d'Ertel montrent que des anomalies sont présentes dans les configurations instables. En particulier, une asymétrie des fonctions de probabilité en faveur des évènements négatifs est observée. Le lien entre la distribution de vorticité potentielle et celle d'un scalaire passif est également étudié. Il est montré qu'à faible nombre de Richardson, c'est le mode vortex (à vorticité potentielle nulle) qui contient les plus importantes fuctuations de scalaire.

Cette thèse s'intéresse à la turbulence incompressible homogène et anisotrope, et plus particulièrement à l'effet d'une rotation solide et d'un champ magnétique uniforme et stationnaire. En plus des simulations numériques directes (DNS), nous présentons également un modèle synthétique de turbulence (Kinematic Simulation ou KS) construit à partir d'une superposition de modes de Fourier, et au sein duquel la dynamique linéaire (basée sur la Rapid Distorsion Theory) peut être incluse. Dans un premier temps, l'effet de la rotation solide est étudié avec un effort particulier porté sur les corrélations en deux-temps. Une comparaison entre DNS et KS est proposée dans le cas isotrope comme dans le cas en rotation. Dans le contexte aéroacoustique, on montre dans quelle mesure la rotation modifie l'émission acoustique d'une turbulence homogène. L'effet d'un champ magnétique est ensuite considéré et l'anisotropie de l'écoulement est étudiée en fonction du nombre de Reynolds magnétique. Il est montré dans quelle mesure la turbulence magnétohydrodynamique quasi-statique est similaire à la turbulence "deux-dimensions trois-composantes" du fait de la dissipation Joule anisotrope et comment l'écoulement restaure son isotropie lorsque le nombre de Reynolds magnétique croît. Enfin, l'effet couplé d'une rotation et d'un champ magnétique est considéré. Les propriétés énergétiques ainsi que l'anisotropie sont étudiés et une étude paramétrique en fonction du nombre d'Elsasser est proposée.

La finalité de cette thèse est de fournir une meilleure compréhension de la turbulence homogène anisotrope soumise à un forçage barocline. À cette fin, nous utilisons une approche numérique pseudo-spectrale basée sur la transformation de Rogallo. L’utilisation d’un tel algorithme nous permet de considérer une asymétrie des fonctions de probabilité en faveur des évènements négatifs est observée. Le lien entre la distribution de vorticité potentielle et celle d’un scalaire passif est également étudié. Il est montré qu’à faible nombre de Richardson, c’est le mode vortex (à vorticité potentielle nulle) qui contient les plus importantes fuctuations de scalaire. Un écoulement homogène dans les trois directions de l’espace. Plusieurs simulations numériques directes (DNS) sont effectuées dans un contexte assez proche des écoulements géophysiques que l’on retrouve entre autre dans la stratosphère, où un gradient constant de vitesse zonale vient se coupler à un gradient constant de densité dans un repère tournant. Les résultats obtenus s’articulent autour de quatre axes principaux. Tout d’abord, une étude linéaire à temps fini est présentée en vue de compléter les résultats existants sur la dynamique linéaire asymptotique. La solution linéaire est décomposée en une partie ‘onde’ (qui se propage) et une partie dite ‘vortex’(stationnaire). L’étude analytique est complétée par un modèle synthétique de turbulence (Kinematic Simulation ou KS) basé sur la théorie de la distorsion rapide(RDT). Nous montrons qu’une distribution initiale non nulle de vorticité potentielle linéarisée peut conduire à d’importantes croissances transitoires. Ce résultat pourrait s’étendre à des modélisations du climat ou météorologique, où la distribution initiale de vorticité potentielle semble avoir autant d’importance que la distribution initiale de température ou de vitesse. Ensuite, nous consacrons une partie de notre étude à l’analyse paramétrique et à la stabilité de l’écoulement. Plusieurs DNS sont effectuées pour différents taux de rotation et stratification. Le diagramme de stabilité obtenu montre que pour de faibles taux de rotation, la limite de stabilité est identique à celle connue des écoulements sans rotation. À plus faible nombre de Rossby — lorsque la baroclinicité devient importante — la limite linéaire de stabilité Ri = 1 relative à l’instabilité symétrique est confirmée. La coexistance de l’instabilité barocline avec l’instabilité symétrique est également clarifiée. Une analyse énergétique détaillée mène à la conclusion suivante : la stratification doit être suffisamment importante (Ri ' 1) pour que l’instabilité barocline soit dominante i.e. que la conversion d’énergie potentielle soit la source principale d’énergie cinétique turbulente. Dans le cas contraire, l’instabilité symétrique — qui tire son énergie de l’énergie cinétique de l’écoulement moyen et non de son énergie potentielle — domine la dynamique de l’écoulement. Le troisième axe d’étude concerne la turbulence à proprement parler. En conséquence de l’ajustement géostrophique, le vent thermique force la turbulence d’une manière naturelle, en opposition à d’autres méthodes de forçage stochastique. L’émergence de structures dans le contexte barocline est approfondie. Des statistiques Euleriennes sont présentées afin de fournir une caractérisation fine de l’anisotropie de l’écoulement. Enfin, nous étendons notre étude à la caractérisation de la vorticité potentielle turbulente. Les fonctions de probabilité de la vorticité potentielle d’Ertel montrent que des anomalies sont présentes dans les configurations instables. En particulier, une asymétrie des fonctions de probabilité en faveur des évènements négatifs est observée. Le lien entre la distribution de vorticité potentielle et celle d’un scalaire passif est également étudié. Il est montré qu’à faible nombre de Richardson, c’est le mode vortex (à vorticité potentielle nulle) qui contient les plus importantes fuctuations de scalaire. / The main objective of this thesis is to provide a better understandingof homogeneous turbulence dynamics under an external baroclinic forcing.To achieve this goal, we use a pseudo-spectral code based on the Rogallo transformation.The use of such an algorithm allows to assume homegeneity in the threespatial directions. Direct Numerical Simulations (DNS) are done in a context representativeof geophysical baroclinic flows in the middle atmosphere: superpositionof a uniform mean zonal flow with stable vertical stratification and frame rotation.The results we obtained are then presented along four axes.First, a finite-time linear analysis is done to complete previous asymptotic results.The linearized flow is decomposed into a propagating (wave) and stationary (vortex)part. The analytical work is completed by a Kinematic Simulation (KS) modelbased on Rapid Distortion Theory (RDT). It is shown that the linearized potentialvorticity mode can produce dramatic transient growth of the kinetic energy if nonzeroinitially. The consequence of such a result is then of capital interest in climatemodelling, where the initial distribution of potential vorticity seems to have moreimportance than other eulerian quantities (temperature or velocity).The second axis is dedicated to a parametric analysis of the flow stability. SeveralDNS are done for different rotation and stratification rates. The derived stabilitydiagram shows that at low rotation rates, the stability bound for purely shearedstratifiedflows is recovered. At higher rotation rates — when baroclinicity is dominant— the linear bound for the so-called symmetric instability is confirmed. Thecoexistence of baroclinic and symmetric instabilies is also clarified. A complete energeticanalysis leads to the conclusion that stratification must be sufficiently highto enhance potential energy release through baroclinic instability. If not, symmetricinstability — driving its energy from the kinetic energy of the mean flow and notfrom the potential energy of the mean flow — is found to dominate the dynamics.The third axis is devoted to a characterisation of homogeneous turbulence submittedto an external baroclinic forcing. As coming from the geostrophic adjustment,the thermal shear allows an organic forcing of turbulence, in opposition to ad-hocarticificial forcing. The structures associated with the simultaneous presence of rotation,stratification and shear are investigated. Eulerian statistics are gathered togive a sharp characterisation of the spatial anisotropy of the flow.Finally, we open our work to the study of turbulent potential vorticity. Probabilitydensity functions of Ertel’s potential vorticity show that potential vorticityanomalies are present in unstable configurations. In particular, an asymmetry ofthe probability density functions toward negative events is observed. An attemptto link potential vorticity dynamics with scalar mixing in baroclinic flows is donethrough joint probability functions analysis.

Le but de cette étude expérimentale est de caractériser la dynamique de particules solides, à flottabilité nulle, incluse dans un écoulement turbulent isotrope en décroissance libre. Les particules utilisées sont de forme sphérique et ont un diamètre de 4 à 5 fois plus grand que l'échelle spatial de Kolmogorov de l'écoulement. De part leur taille, les particules ont également un nombre de Stokes proche de l'unité. On s'attend alors à ce que ces particules aient une dynamique différente de celle du fluide environnant. Dans cette étude, ont se propose de quantifier les différences de vitesses entre les deux phases à l'aide d'une technique de vélocimétrie simultanée / In this experimental study, the focus is made on the characterization of the dynamics of solid neutrally buoyant particles embedded in a freely decaying, nearly isotropic turbulence, with a weak mean flow. The particles are spherical with diameters several times larger than the Kolmogorov scale. The study of this flow configuration is still challenging both theoretically and numerically. Due to large particle sizes, the local flow around particles can not be considered as uniform and due to fluid-particle density ratio of around unity, the history and Basset forces cannot be neglected in comparison with the viscous drag force. Particle equation of motion is then fully non-linear, in contrast to the equation for heavy particles with diameters smaller then the Kolmogorov scale, for which only the Stokes drag is considered. In several experimental and numerical studies, the effect of particle size on velocity and acceleration statistics has been investigated (Homann and Bec 2010 ; Qureshi et al. 2008 ; Ouellette et al. 2008 ; Xu and Bodenschatz 2008). In the case of isotropic turbulence, Homann and Bec (2010) show that while the PDF of the particle velocity normalized by the square root of its variance does not vary with particle size, the variance itself is size dependent. A scaling relation for particle velocity variance has been proposed by using the Faxen correction (Gatignol 1983) which takes into account the non uniformity of the fluid flow at the scale of the particle. The aim of our research is to further study the dependence of particle dynamics on particle size. To that purpose, a turbulence generator has been set-up and the resulting turbulence is characterized. Then the flow was seeded with millimeter sized, neutrally-buoyant particles and the velocity of the two phases have been measured simultaneously. Simultaneous measurements of particle and surrounding fluid velocities show that although the global velocity statistics of the two phases have comparable values, the particles may have different local velocity from the velocity of the neighboring fluid

Dans cette thèse, l'application de l'assimilation de données (AD) à la MFN est étudiée, avec comme objectif global de contribuer à l'amélioration de la prévision numérique d'écoulements complexes. L'AD consiste à fusionner les outils de prévision numérique avec des données expérimentales afin d'améliorer l'estimation des paramètres d'entrée du code MFN. Les aspects méthodologiques de l'AD et son application pour des études physiques sont tous deux examinés dans cette thèse. Dans un premier temps, l'AD est utilisée pour une étude théorique de la turbulence de grille. Un modèle spectral pour les écoulements turbulents homogènes et anisotropes est également proposé. Plusieurs méthodes d'AD sont ensuite implémentées pour un code MFN et appliquées à la reconstruction d'écoulements instationnaires et compressibles en présence d'incertitudes sur des paramètres d'entrée de grandes dimensions afin d'évaluer les forces et faiblesses respectives de ces techniques. Des stratégies pour le placement optimal de réseaux de capteurs sont élaborées afin d'améliorer les performances du processus d'AD. Enfin, l'AD est appliquée à l'identification de sources de polluants et à la reconstruction de conditions météorologiques pour des écoulements en milieu urbain prédits par Simulation des Grandes Echelles. / In this thesis, we investigate the use of various data assimilation (DA) techniques in the context of CFD, with the ultimate goal of enhancing the prediction of real-world flows. DA consists in merging numerical predictions and experimental observations in order to improve the estimation of the CFD solver inputs. Both methodological aspects of DA and its potential application to physics investigations are explored for various flow configurations. First, DA is considered for the theoretical analysis of grid turbulence decay. Fundamental aspects of anisotropic homogeneous turbulence are also investigated through spectral modelling. Various DA methodologies are deployed in conjunction with a Navier-Stokes solver and are assessed for the reconstruction of unsteady compressible flows with large control vectors. Sensor placement strategies are developed to enhance the performances of the DA process. Finally, a first application of DA to Large Eddy Simulations of full-scale urban flows is proposed with the aim of identifying source and wind parameters from concentration measurements.

Le moteur downsizé à allumage commandé constitue l'une des voies principales explorées par les constructeurs automobiles pour améliorer le rendement et réduire les émissions de dioxyde de carbone des motorisations essence. Il s'agit de combiner une réduction de la cylindrée avec une forte suralimentation afin d'améliorer le rendement du moteur, en particulier à faibles et moyennes charges. Leur mise au point est limitée par l'augmentation des combustions anormales, dont le contrôle par forte dilution peut également entraîner l'apparition de variabilités cycliques importantes. Actuellement, la compréhension des nombreux paramètres intervenant dans l'apparition de ces phénomènes et de leurs interactions, reste encore imparfaite. Dans ce contexte, l'objectif de ce travail est de contribuer à la compréhension des mécanismes impliqués dans les processus de propagation des flammes turbulentes. Cette étude est réalisée dans une enceinte de combustion sphérique haute pression haute température, équipée de ventilateurs générant une turbulence homogène et isotrope. La première partie de ce travail est consacrée à l'étude de la combustion prémélangée laminaire isooctane/air. Dans un deuxième temps, l'aérodynamique de l'écoulement dans l'enceinte est finement caractérisée par Vélocimétrie Laser Doppler et Vélocimétrie par Images de Particules. Enfin, la propagation des flammes turbulentes est étudiée en termes de vitesse à partir de visualisations par ombroscopie. Une loi unifiée, permettant de décrire la propagation des flammes turbulentes indépendamment des conditions thermodynamiques initiales, de l'intensité de la turbulence et de la nature du mélange réactif est notamment proposée.

Le moteur downsizé à allumage commandé constitue l’une des voies principales explorées par les constructeurs automobiles pour améliorer le rendement et réduire les émissions de dioxyde de carbone des motorisations essence. Il s’agit de combiner une réduction de la cylindrée avec une forte suralimentation afin d’améliorer le rendement du moteur, en particulier à faibles et moyennes charges. Leur mise au point est limitée par l’augmentation des combustions anormales, dont le contrôle par forte dilution peut également entraîner l'apparition de variabilités cycliques importantes. Actuellement, la compréhension des nombreux paramètres intervenant dans l'apparition de ces phénomènes et de leurs interactions, reste encore imparfaite. Dans ce contexte, l’objectif de ce travail est de contribuer à la compréhension des mécanismes impliqués dans les processus de propagation des flammes turbulentes. Cette étude est réalisée dans une enceinte de combustion sphérique haute pression haute température, équipée de ventilateurs générant une turbulence homogène et isotrope. La première partie de ce travail est consacrée à l’étude de la combustion prémélangée laminaire isooctane/air. Dans un deuxième temps, l’aérodynamique de l’écoulement dans l’enceinte est finement caractérisée par Vélocimétrie Laser Doppler et Vélocimétrie par Images de Particules. Enfin, la propagation des flammes turbulentes est étudiée en termes de vitesse à partir de visualisations par ombroscopie. Une loi unifiée, permettant de décrire la propagation des flammes turbulentes indépendamment des conditions thermodynamiques initiales, de l’intensité de la turbulence et de la nature du mélange réactif est notamment proposée. / Engine downsizing is a promising way explored for reducing carbon dioxide and pollutant emissions of spark ignition engines. Its principle is to reduce the engine size and to increase its specific power with the use of a turbocharger, especially at middle and high loads. Its development leads to the appearance of abnormal combustions, whose control by a high dilution rate can also lead to important cyclic variabilities. The understanding of the multiple parameters controlling these phenomena and their interactions, is still very limited. In this context, the present thesis addresses the issue related to the understanding of the mechanisms implied in turbulent flame propagation processes. This study is carried out in a fan-stirred high pressure high temperature spherical combustion vessel where turbulence is homogeneous and isotropic. The first part of this study concerns the study of laminar premixed flames of isooctane/air mixtures. Second, the turbulent flow in the vessel is accurately characterized using Laser Doppler Velocimetry and Particle Image Velocimetry measurements. Then, turbulent flame propagation is investigated in terms of flame velocity using shadowgraph visualizations. A unified scaling law for the turbulent propagation speed is especially proposed, regardless of the thermodynamic initial conditions, the intensity of the turbulence and the nature of the air/fuel mixture.

Une des problématiques majeures dans l’industrie des bioprocédés réside dans l’extrapolation des procédés biologiques à grande échelle. On observe généralement à l’échelle industrielle des écarts de rendement de croissance de la biomasse, ainsi que la formation de sous-produits comparativement à l’échelle du laboratoire. La formation de gradients de concentration à l’échelle des bioréacteurs est souvent évoquée. Dans ce travail, les interactions entre micromélange et assimilation du substrat sont abordées à l’échelle du microorganisme. Un modèle couplant transport et assimilation à l’échelle d’un microorganisme est proposé. L’existence de régimes physique et biologique, limitant l’assimilation du substrat est mise en lumière. Une approche basée sur le suivi Lagrangien de particules dans un champ de turbulence homogène isotrope est ensuite retenue. Les effets des hétérogénéités de concentration vues par les microorganismes, sont traduits à l’échelle de la population entière. Une loi analytique permettant de construire la distribution de flux reçus par les microorganismes à partir de la distribution de concentration en substrat dans le fluide, est proposée. Partant de cette distribution de concentrations vues, l’adjonction d’un modèle métabolique simplifié permet d’expliquer les baisses de vitesse spécifiques de croissance et la formation de sous-produits observées expérimentalement. Enfin, de premiers résultats sur le couplage inverse biologique sont présentés. L’effet des microorganismes sur le champ de concentration est caractérisé et une étude paramétrique sur les propriétés dynamiques et biologiques est réalisée. / The scale-up of biological process is a critical issue in the bioprocess industry. When passing from a laboratory to an industrial scale, the conversion yield of substrate into biomass is often overestimated and by-products are formed. Different existing works attempt to predict the effect of mixing on biomass growth and the emergence of substrate concentration gradients at the reactor scale are a first explanation of the degraded performances. In this work the interactions between micro-mixing and substrate assimilation are addressed at the microorganism scale. A coupled transport-assimilation model is proposed for an isolated microorganism. The emergence of physical and biological regimes limiting the substrate assimilation is enlightened. An approach based on the Lagrangian tracking of microorganisms in a homogeneous isotropic turbulent field is then chosen. The effects of local concentration heterogeneities seen by microorganisms are observed at the population scale. An analytical expression is proposed for the assimilated substrate flux distribution by the microorganisms, based on the substrate concentration distribution in the fluid. From these concentrations encountered by microorganisms, we coupled a simplified metabolic model that explains the decreased specific growth rate, and the by-products formation often observed in many experiments. Finally, first results on the biological two-way coupling are proposed. The effect of microorganisms on the substrate field is characterised and a parametric study on the dynamics as well as biological parameters is realised.

La propagation de flammes turbulentes dans des milieux réactifs inhomogènes concerne un grand nombre d’applications pratiques, y compris celles qui reposent sur des cycles de combustion à volume constant. Les hétérogénéités de composition (richesse, température,dilution par des gaz brûlés, etc.) sont issues de plusieurs facteurs distincts tels que la dispersion du spray de gouttelettes de combustible et son évaporation, la topologie de l’écoulement ainsi que la présence éventuelle de gaz brûlés résiduels issus du cycle précédent. La structure des flammes partiellement prémélangées qui en résultent est significativement plus complexe que celles des flammes plus classiques de diffusion ou de prémélange. L’objectif de ce travail de thèse est donc de contribuer à l’amélioration de leur connaissance, en s’appuyant sur la génération et l’analyse de base de données de simulations numériques directes ou DNS (Direct Numerical Simulation). Celles-ci sont conduites avec le code de calcul Asphodele qui est basé sur l’approximation de faible nombre de Mach. Le combustible de référence retenu est l’iso-octane.La base de données est structurée suivant cinq paramètres qui permettent de caractériser l’écoulement turbulent ainsi que l’hétérogénéité de composition du milieu réactif. Dans un premier temps, des configurations bidimensionnelles ont été considérées en raison du coût élevé induit par la description détaillée de la cinétique chimique. L’étude des ces différents cas de calcul a permis de mettre en lumière plusieurs mécanismes fondamentaux de propagation dans les milieux hétérogènes en composition. Une réduction significative des coûts de calcula pu ensuite être obtenue grâce au développement d’un modèle chimique simplifié optimisé.Son utilisation a permis d’étendre les analyses à de / The propagation of turbulent flames in non-homogeneous reactive mixtures of reactants concerns a large number of practical applications, including those based on constant volume combustion cycles. The composition heterogeneities (equivalence ratio, temperature, dilution by burnt gases, etc.) result from several distinct factors such as the dispersion of the spray of fuel droplets and its evaporation, the flow field topology as well as the possible presence of residual burnt gases issued from the previous cycle. The resulting partially premixed flames structure is significantly more complex than the one of more conventional diffusion or premixed flames.The aim of this thesis work is therefore to contribute to the improvement of their understanding, by proceeding to the generation and analysis of a new set of direct numerical simulations (DNS) databases. The present computations are performed with the low-Mach number DNS solver Asphodele. The database is structured according to five parameters that characterize the turbulent flow as well as the composition heterogeneity of the reactive mixture. First, because of the high numerical costs induced by the detailed description of chemical kinetics, two-dimensional configurations were considered. The study of these various simulations highlights several fundamental mechanisms of flame propagation in heterogeneous mixtures. Then, a significant computational cost saving has been achieved through the development of an optimized simplified chemistry model. The use of the latter allowed to overcome the major bottleneck of high CPU costs related to chemical kinetics description and thus to extend the analysis to three-dimensional configurations. Some of the conclusions obtained previously were reinforced.