Le travail s'intéresse aux effets d'un pré-écrouissage cyclique axial ou en torsion sur le comportement cyclique et la durée de vie en fatigue sous sollicitation axiale à température ambiante d'un acier inoxydable austénitique 304L. Les essais cycliques séquentiels à amplitude de déformation croissante ou décroissante montrent que la courbe cyclique du 304L est non-unique. Cette caractéristique est liée à la persistance de l'état microstructural généré pendant les cycles d'amplitude de déformation maximale de la première séquence. En augmentant le nombre des séquences, l'acier 304L montre une tendance vers une courbe cyclique asymptotique, l'écrouissage semble se stabiliser. Des essais de fatigue sous chargement axial ont été réalisés sur des éprouvettes vierges ou pré-écrouies en traction-compression ou en torsion. Les durées de vie ont été sensiblement réduites pour les éprouvettes pré-écrouies. Ce phénomène est lié à la formation de structures de dislocations denses héritées de la phase de pré-écrouissage. Cependant, l'augmentation de l'amplitude de déformation en fatigue réduit l'effet du pré-écrouissage.

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.

Les écoulements turbulents de gaz denses, qui sont d’un grand intérêt pour un large éventail d'applications, sont le siège de phénomènes physiques encore peu connus et difficiles à étudier par des approches expérimentale. Dans ce travail, nous étudions pour la première fois l’influence des effets de gaz denses sur la structure de la turbulence compressible à l’aide de simulations numériques. Le fluide considéré est le PP11, un fluorocarbure lourd, dont le comportement thermodynamique a été représenté à l’aide de différentes lois d’état, afin de quantifier la sensibilité des solutions aux choix de modélisation. Nous avons considéré d’abord la décroissance d’une turbulence homogène isotrope compressible. Les fluctuations de température sont négligeables, alors que celles de la vitesse du son sont importantes à cause de leur forte dépendance de la densité. Le comportement particulier de la vitesse du son modifie de manière significative la structure de la turbulence, conduisant à la formation de shocklets de détente. L’analyse de la contribution des différentes structures à la dissipation d’énergie et à la génération d’enstrophie montre que, pour un gaz dense, les régions de forte dilatation jouent un rôle similaire à celles de forte compression, contrairement aux gaz parfaits, dans lesquels le comportement est fortement dissymétrique. Ensuite, nous avons mené des simulations numériques pour une configuration de canal plan en régime supersonique, pour plusieurs valeurs des nombres de Mach et de Reynolds. Les résultats confirment la validité de l’hypothèse de Morkovin. L’introduction d’une loi d’échelle semi-locale prenant en compte le variations de densité et viscosité, permet de comparer les profils des grandeurs turbulentes (contraintes de Reynolds, anisotropie, budgets d’énergie) avec ces observés en gaz parfait. Les variables thermodynamiques, quant à elles, présentent une évolution très différente pour un gaz parfait et pour un gaz dense, la chaleur spécifique élevée de ce dernier conduisant à un découplage des effets dynamiques et thermiques et à un comportement proche de celui d’un fluide incompressible avec des propriétés variables. / Dense gas turbulent flows, of great interest for a wide range of engineering applications, exhibit physical phenomena that are still poorly understood and difficult to reproduce experimentally. In this work, we study for the first time the influence of dense gas effects on the structure of compressible turbulence by means of numerical simulations. The fluid considered is PP11, a heavy fluorocarbon, whose thermodynamic behavior is described by means of different equations of state to quantify the sensitivity of solutions to modelling choices. First, we considered the decay of compressible homogeneous isotropic turbulence. Temperature fluctuations are found to be negligible, whereas those of the speed of sound are large because of the strong dependence on density. The peculiar behavior of the speed of sound significantly modifies the structure of the turbulence, leading to the occurrence of expansion shocklets. The analysis of the contribution of the different structures to energy dissipation and enstrophy generation shows that, for a dense gas, high expansion regions play a role similar to high compression ones, unlike perfect gases, in which the observed behaviour is highly asymmetric. Then, we carried out numerical simulations of a supersonic turbulent channel flow for several values of Mach and Reynolds numbers. The results confirm the validity of the Morkovin' hypothesis. The introduction of a semi-local scaling, taking into account density and viscosity variations across the channel, allow to compare the wall-normal profiles of turbulent quantities (Reynolds stresses, anisotropy, energy budgets) with those observed in ideal gases. Nevertheless, the thermodynamic variables exhibit a different evolution between perfect and dense gases, since the high specific heats of the latter lead to a decoupling of dynamic and thermal effects, and to a behavior close to that of variable property incompressible fluids.

Le travail s’intéresse aux effets de l’histoire de chargement sur le comportement et la durée de vie en fatigue de deux nuances (THYSSEN et CLI)d’un acier inoxydable austénitique 304L à la température ambiante. Les essais ont été réalisés en utilisant deux catégories d’éprouvettes. Les éprouvettes de la première catégorie (vierges) ont été soumises à des essais classiques de fatigue,alors que celles de la deuxième ont subi, avant les essais de fatigue, un pré-écrouissage monotone ou cyclique en déformation imposée. Les éprouvettes vierges manifestent un adoucissement cyclique suivi d’un durcissement cyclique alors que les éprouvettes pré-écrouies ne présentent qu’un durcissement cyclique. Les résultats montrent une grande influence du pré-écrouissage qui semble bénéfique en contrainte imposée, mais néfaste en déformation imposée,même en présence d'une contrainte moyenne de compression. Ces résultats sont discutés en termes d'évolution cyclique du module d'élasticité, des écrouissages isotropes et cinématiques, et de la densité d’énergie absorbées par cycle, dans différentes configurations : avec ou sans pré-écrouissage, en contrainte ou déformation imposées... / This study deals with the effect of the loading history on the cyclic behavior and the fatigue life of two kinds (THYSSEN and CLI) of 304L stainless steel at room temperature. The experiments have been performed using two specimens’ categories. The first one (virgin) has been submitted to only classical fatigue tests while in the second category, prior to the fatigue test, the specimen is subjected to a pre-hardening process under either monotonic or cyclic strain control. Cyclic softening followed by cyclic hardening are observed for the virgin specimens while only cyclic softening is exhibited by the pre-hardened specimens. The obtained results show that fatigue life is strongly influenced by the pre-hardening: it seems beneficial under stress control but detrimental under strain control, even in the presence of a compressive mean stress. The results are discussed regarding the cyclic evolution of the elastic modulus aswell as the isotropic and kinematic parts of the strain hardening, and strain energy density per cycle, in different configurations: with or without prehardening,stress or strain control...

Le développement de méthodes d’identification de paramètres de lois de comportement de matériaux est devenu primordial pour avoir accès à la connaissance complète du comportement. En effet, les méthodes de mesure optiques, comme la Corrélation d’Images Numériques, permettent d’obtenir les quantités cinématiques de la relation de comportement sous forme de champs de vecteurs. En revanche, les contraintes ne sont généralement pas mesurables et il est nécessaire d’identifier les paramètres de la loi de comportement du matériau considéré pour y avoir accès. Plusieurs méthodes ont vu le jour et permettent de répondre à cette problématique mais la plupart d’entre elles supposent une homogénéité du matériau. Ce mémoire traite de l’application de certaines de ces méthodes, notamment la méthode de l’écart à l’équilibre (MEQ) et la méthode de recalage de modèle éléments finis (MREF), dans des matériaux hétérogènes à microstructure complexe où les propriétés mécaniques évoluent spatialement dans le volume. L’objectif est d’identifier ces propriétés mécaniques locales qui régissent la cinématique mesurée de tels matériaux dans le cadre de l’élasticité linéaire isotrope. Dans un premier temps, les deux méthodes citées sont décrites, implémentées et comparées sur des cas simulés en 2D. La MREF est préférée à la MEQ car plus robuste vis-à-vis des incertitudes de mesure. Basée sur un formalisme itératif, une parallélisation de l’algorithme a été opérée pour diminuer le coût en temps de la méthode. Des expérimentations dans le plan sur des éprouvettes en polyuréthane où les hétérogénéités sont maîtrisées ont permis de valider la méthode. Enfin, deux applications en 3D sur un matériau en mousse polyuréthane et un composite à base de fibres de bois démontrent l’intérêt d’une telle méthode pour l’identification de propriétés mécaniques locales. La mise en évidence d’une relation entre les propriétés locales identifiées et les propriétés locales de la microstructure de ces matériaux est effectuée. / The development of identification methods of material constitutive equation parameters has become fundamental to completely know the mechanical behavior. Indeed, optical methods, such as Digital Image Correlation, allows to get kinematics quantities of the constitutive equation as vectors fields. But, stresses are usually not available experimentally and one has to identify constitutive equation parameters to compute them. Several methods have been developed and answer to that problematic but most of them suppose the materials as homogeneous. This memoir is about the application of some of these methods, such as the equilibrium gap method (EGM) and the finite element model updating method (FEMU), in the case of heterogeneous materials with complex structures where mechanical properties vary spatially in the volume. The objective is to identify these local mechanical properties which rule the measured kinematics of such materials considering the isotropic linear elasticity. Firstly, both methods are detailed, implemented and compared on 2D simulated cases. The FEMU method is preferred because it is more robust in the presence of noisy data. Based on an iterative process, a parallelisation of the algorithm is achieved in order to reduce the cost of the method. In-plane experiments on polyurethane samples where heterogeneities are controlled have validated the method. Finally, two 3D applications on a polyurethane foam material and a wood-based fibrous composite have demonstrated the interest of this approach to identify local mechanical properties. The highlighting of a relationship between identified local properties and microstructural properties of these materials is made.

La modélisation numérique des structures bois dans des conditions d’incendie nécessite la connaissance : de la variation des propriétés physiques du bois telles que la conductivité thermique, la chaleur spécifique et la densité en fonction de la température ; de la dégradation thermique du bois au cours des phases de séchage, de pyrolyse et de combustion. En particulier, nous nous sommes intéressés à l’étude du comportement thermomécanique du matériau bois. La loi thermique est décrite par l’équation de la chaleur. Le modèle choisi intègre les trois modes du transfert de chaleur : la conduction, le rayonnement et la convection. La loi mécanique est modélisée dans le cadre de la thermodynamique des processus irréversibles utilisant la notion des variables d’état. Elle tient compte du couplage entre le comportement élastique orthotrope, plastique anisotrope à écrouissage non linéaire isotrope et un endommagement isotrope. L’intégration numérique de la loi mécanique par un schéma implicite itératif combinant la technique du retour radial avec la réduction du nombre des équations est présentée. Le couplage thermomécanique est réalisé, selon l’approche réglementaire de l’Eurocode 5 relatif à la résistance au feu des structures en bois, en appliquant le facteur de réduction Kθ sur la résistance mécanique d’un résineux. Les aspects théoriques et les conditions aux limites associés au modèle thermomécanique sont abordés. L’identification des paramètres du modèle est réalisée sur des données expérimentales obtenues sur des tests réels d’incendie disponibles dans la littérature. À ce titre, plusieurs comparaisons avec différentes applications sont réalisées. Le modèle éléments finis reproduit avec précision la distribution du champ de température dans l’épaisseur des panneaux en bois, la formation du charbon ainsi que l’évolution de la résistance mécanique au cours de l’exposition au feu / Numerical modelling of timber structures in fire conditions requires the knowledge of the variation with temperature of the physical properties of the wood material (the thermal conductivity, the specific heat and the density) in order to take into account the thermal degradation of wood under high temperatures during the drying, pyrolysis and combustion phases, as well as the temperature profiles in the thickness of the surfaces exposed to fire. In particular, this work focusses on the thermomechanical behaviour of timber. The heat transfer analysis is described by the standard equations of heat conduction. It includes the three modes of heat transfer: conduction, radiation and convection. The structural response is modelled within the framework of thermodynamics of irreversible processes using the notion of state variables. It takes into account the coupling between the orthotropic elastic behaviour, the anisotropic plastic behaviour with isotropic nonlinear hardening, and isotropic damage. The numerical integration of the equilibrium equations is carried out with an iterative implicit scheme combining the technique of radial re- turn with the reduction of the number of equations. The thermomechanical coupling is carried out according to the approach recommended by Eurocode 5 for the fire resistance of timber structures by applying the reduction factor Kθ to the strength of a softwood. The theoretical aspects and boundary conditions associated with the thermomechanical model are also discussed. The parameters of the model are identified with experimental data obtained from actual fire tests available in the literature. Several comparative applications are carried out. The finite element model accurately reproduces the distribution of the temperature profile in the thickness of timber planks, the formation of the charred layer, and the evolution of the mechanical resistance during exposure to fire

Die Bestrahlungstherapie von Aderhautmelanomen mit der Protonentherapie stellt ein sehr präzises Behandlungsverfahren dar, welches hohe Ansprüche an die zugrundliegende Bildgebung stellt. Die Magnetresonanztomographie (MRT) und Computertomographie (CT) scheinen prinzipiell in der Lage, mit geeigneten Oberflächenspulen bzw. angepaßter Untersuchungstechnik diese Ansprüche zu erfüllen. Ziel der Arbeit ist die Bereitstellung und Überprüfung eines bildgebenden Systems, das einerseits die diagnostischen und differentialdiagnostischen Aufgaben für die Differenzierung und das Staging von Aderhautmelanomen erfüllt und andererseits als Grundlage der Bestrahlungsplanung dienen kann. Erstrebenswert ist dabei die Erstellung von Schnittbildern bestehend aus isotropen Voxeln, d.h. Voxeln mit identischer Kantenlänge, mit der Magnetresonanztomographie (MRT) und der Computertomographie (CT). Zur Differenzierung der unterschiedlichen Tumorentitäten wurden verschiedene Untersuchungen mit zum Teil großen Patientenkollektiven durchgeführt und histologisch und klinisch korreliert bzw. verifiziert. Dabei konnte das MR-morphologische Erscheinungsbild der relevanten Erkrankungen studiert werden. Die Genauigkeit der MR-gestützten Tumorvolumetrie wurde in vitro im Tierversuch und in vivo am Patientenauge überprüft. Zur Überprüfung der Genauigkeit der MR-Bildgebung mit isotropen Voxeln wurden unterschiedliche Spulen getestet und verschiedene Phantomuntersuchungen durchgeführt. Entsprechend wurde auch die CT evaluiert. Mit der digitalen Bildfusion wurden die MRT- und die CT-Daten in einen gemeinsamen Datensatz überführt, der sowohl zur Bestimmung der Genauigkeit beider Verfashren als auch als Basis zur Bestrahlungsplanung verwendet werden konnte. Mit einer optimierten Untersuchungstechnik konnten Schnittbilder auf der Basis von CT und MRT zur Verfügung gestellt werden, die eine präzise Planung der Protonenbestrahlung beim Aderhautmelanom ermöglichen. / The radiation therapy of uveal melanomas using proton therapy is a very precise therapy procedure, that makes high demands on the underlying imaging modalities. In general magnetic resonance imaging (MRI) and computed tomography (CT) have the potential to fulfill these criteria using suited surface coils and optimized imaging techniques, respectively. Purpose of the investigations was the preparation and evaluation of an imaging system, that allows diagnosis and differetial diagnosis of uveal melanomas including staging on the one hand. On the other hand it should be the base for high precision radiation therapy planning. The preparation of images with isotropic voxels, i.e. the edges of the voxels have identical size, is desirable with both magnetic resonance imaging (MRI) and computed tomography (CT). Several examinations on large patient groups were performed for the differentiation of different tumor entities and compared with clinical and histopathological outcome. The precision of the MR-based tumor volumetry was evaluated with animal studies in vitro and patient eyes in vivo. The MR-imaging with isotropic voxels was evaluated using several surface coils and phantom examinations. Comparable examinations were performed regarding CT. Digital image fusion of CT and MRI data sets was used to evaluate the precision of both modalities and for planning of radiation therapy. Using optimized examination techniques both CT and especially MRI are suitable for a precise proton therapy planning of uveal melanoma.

Une nouvelle approche appelée “Approche Pseudo Inverse” (API) est développée pour la modélisation du procédé de forgeage à froid des pièces axisymétriques. L'API est basée sur la connaissance de la forme de la pièce finale. Certaines configurations intermédiaires « réalistes » ont été introduites dans l'API pour considérer le chemin de déformations. Elles sont créées géométriquement sans traitement de contact et ensuite corrigées par la méthode de surface libre afin de respecter l'équilibre, les conditions aux limites et la condition d'incompressibilité. Un nouvel algorithme direct de plasticité est développé, conduisant à une méthode d'intégration plastique très rapide, précise et robuste même dans le cas de très grands incréments de déformations. Un modèle d'endommagement en déformation, est couplé à la plasticité et implémenté dans l'API. Les validations numériques montrent que l'API donne des résultats très proches des résultats de l'approche incrémentale mais en utilisant beaucoup moins de temps de calcul.L'API est adoptée comme solveur du forgeage pour la conception et l'optimisation des préformes du forgeage multi-passes. La rapidité et la robustesse de l'API rendent la procédure d'optimisation très performante. Une nouvelle technique est développée pour générer automatiquement le contour initial d'un outil de préforme pour la procédure d'optimisation. Les variables de conception sont les positions verticales des points de contrôle des courbes B-spline définissant les formes des outils de préforme. Notre optimisation multi-objectif consiste à minimiser la variation de la déformation plastique équivalente dans la pièce finale et la force du poinçon au cours du forgeage. Un algorithme génétique et un algorithme de recuit simulé sont utilisés pour trouver les points d'optimum de Pareto. Pour réduire le nombre de simulations de forgeage, un méta-modèle de substitution basé sur la méthode de krigeage est adopté pour établir une surface de réponse approximative. Les résultats obtenus par l'API en utilisant les outils de préforme optimaux issues de l'optimisation sont comparés à ceux obtenus par les approches incrémentales classiques pour montrer l'efficacité et les limites de l'API. La procédure d'optimisation combinée avec l'API peut être un outil numérique rapide et performant pour la conception et l'optimisation des outillages de préforme. / A new method called “Pseudo Inverse Approach” (PIA) is developed for the axi-symmetrical cold forging modelling. The PIA is based on the knowledge of the final part shape. Some « realistic » intermediate configurations are introduced in the PIA to consider the deformation path. They are created geometrically without contact treatment, and then corrected by using a free surface method in order to satisfy the equilibrium, the boundary conditions and the metal incompressibility. A new direct algorithm of plasticity is proposed, leading to a very fast, accurate and robust plastic integration method even in the case of very large strain increments. An isotropic damage model in deformation is coupled with the plasticity and implemented in the PIA. Numerical tests have shown that the Pseudo Inverse Approach gives very close results to those obtained by the incremental approach, but using much less calculation time.The PIA is adopted as forging solver for the design and optimization of preform tools in the multi-stage forging process. The rapidity and robustness of the PIA make the optimization procedure very powerful. A new method is developed to automatically generate the initial preform tool shape for the optimization procedure. The design variables are the vertical positions of the control points of B-spline curves describing the preform tool shape. Our multi-objective optimization is to minimize the equivalent plastic strain in the final part and the punch force during the forging process. The Genetic algorithm and Simulated Annealing algorithm are used to find optimal Pareto points. To reduce the number of forging simulations, a surrogate meta-model based on the kriging method is adopted to build an approximate response surface. The results obtained by the PIA using the optimal preform tools issued from the optimization procedure are compared to those obtained by using the classical incremental approaches to show the effectiveness and limitations of the PIA. The optimization procedure combined with the PIA can be a rapid and powerful tool for the design and optimization of the preform tools.

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.