Dans le contexte aéronautique, cette thèse, financée par Safran Helicopter Engine, s’intéresse à l’application de l’approche Simulations aux Grandes Échelles (SGE) pour les systèmes de refroidissement de turbine. Le système de refroidissement industriel complexe est divisé en cas académiques plus simples donnant accès à des caractérisations expérimentales de la dynamique et de la thermique. Le jet impactant est traité en tant que système interne et l’écoulement de protection au bord de fuite en tant que système externe. Après une brève introduction du contexte lié au refroidissement de turbine et des objectifs scientifiques, ce manuscrit est divisé en 3 parties. La 1ère partie traite d’un écoulement de jet impactant sur plaque plane représentatif de l’impact à mi-corde. Elle se concentre sur la validation et la qualification des outils et modèles ainsi que sur l’analyse physique de l’écoulement. Les différentes instationarités de l’écoulement sont reliées à la thermique de paroi à l’aide de diagnostics statistiques et d’analyses modales. La 2ème partie s’intéresse à l’impact sur paroi concave représentatif de l’impact au bord d’attaque. Cette étude se concentre principalement sur la caractérisation de l’effet de courbure pour le jet impactant. Contrairement, au consensus actuel sur l’effet de courbure, la réduction des transferts thermiques est observée pour le cas d’étude de cette thèse. Au vu de ces résultats, une discussion est proposée pour tenter d’expliquer cet écart. Finalement, la 3ème partie de ce manuscrit contient une application de la SGE à un système de protection du bord de fuite par film isolant. Dans ce dispositif, des effets de groupe sont mis en évidence. L’impact des choix de modélisation tels que l’hypothèse de périodicité dans la direction de l’envergure est alors évalué. Il est montré que cette hypothèse de périodicité influe sur la prédiction locale de l’efficacité en forçant l’écoulement. En revanche, la prédiction de l’efficacité globale du système de protection n’est pas impactée.

La recherche en bruit de combustion est de nos jours majoritairement consacrée au développement d'outils de calcul du bruit rayonné par les flammes. Les méthodes actuelles de CFD telles que la LES ou la DNS sont capables de fournir le champ acoustique rayonné par des sources de bruit, mais elles sont cependant limitées à des domaines de faible taille, ceci dû à leur fort coût de calcul. Pour surmonter cette limitation, on a vu l'émergence de méthodes hybrides. Dans cette approche, les sources de bruit sont découplées du son rayonné. Les sources sont alors calculées par les méthodes de DNS et de LES tandis que l'analogie acoustique permet de calculer le son rayonné par des codes acoustiques, moins coûteux en temps de calcul.L'objet de cette étude est de développer un outil numérique sur la base de l'analogie acoustique de Phillips pour de faibles nombres de Mach. Il permet de prendre en compte l'impact des conditions limites sur le champ acoustique résultant. La LES et le code de calcul acoustique développé ont été utilisés pour évaluer le bruit produit par une flamme turbulente confinée. Les deux techniques donnent des résultats en accord tant que les bonnes quantités sont comparées: il a été observé que le signal de pression obtenu directement à partir de la LES contient une quantité non négligeable d'hydrodynamique, laquelle doit être négligée si on cherche à comparer seulement les champs acoustiques issus des deux codes.L'hypothèse d'un nombre de Mach faible est totalement réaliste si l'on considère l'écoulement présent dans une chambre de combustion. Elle conduit à des simplifications significatives lorsque les analogies acoustiques sont considérées. Cependant, cette hypothèse ne peut pas être utilisée pour l'écoulement en amont (entrée d'air, compresseur) ni en aval (turbine, tuyère) des chambres de combustion aéronautiques. Un outil numérique a été développé pour pallier ce problème. Il est basé sur les équations d'Euler Quasi-1D, qui prennent en compte des écoulements convectifs, non isentropiques et non isenthalpiques. Cet outil permet d'estimer les conditions limites acoustiques qui doivent être imposées sur les entrées/sorties d'une chambre de combustion pour prendre en compte la présence d'un écoulement de nombre de Mach non négligeable, alors que les calculs acoustiques sont eux effectués sous cette hypothèse fortement restrictive. / Today, much of the current effort in combustion noise is the development of efficient numerical tools to calculate the noise radiated by flames. Although unsteady CFD methods such as LES or DNS can directly provide the acoustic field radiated by noise sources, this evaluation is limited to small domains due to high computational costs. Hybrid methods have been developed to overcome this limitation. In these schemes, the noise sources are decoupled from the radiated sound. The sources are still calculated by DNS or LES codes whereas the radiated sound is evaluated by acoustic codes using an acoustic analogy.In the present study, a numerical tool based on the Phllips' analogy for low Mach numbers flows has been developed. This tool accounts for the role of the boundary conditions in the resulting acoustic field. Both LES and the acoustic code developed here are used to assess the noise produced by a turbulent confined flame of a turbulent swirled-stabilized staged combustor. Good agreements are obtained between both techniques as long as the good quantities are compared: the pressure signal obtained directly from LES contains a non negligible amount of hydrodynamics that must be removed when a suitable acoustics-acoustics comparison is sought. The low Mach number assumption is completely realistic when considering the flow within a combustion chamber; it also conducts to considerable simplifications when leading with acoustic analogies. However, it cannot be used for the upstream (air-intake, compressors) and downstream (turbines, nozzle) of an aeronautical combustion chamber. A numerical tool is developed based on the quasi-1D Linearized Euler Equations in order to account for convective, non-isentropic and non-isenthalpic flows. By means of this tool, it is possible to estimate the acoustic boundary conditions that should be imposed at the inlet/oultlet of a given combustion chamber when performing low-Mach number acoustic computations.

Cette thèse porte sur la simulation numérique d'un écoulement interne compressible, ins- tationnaire et turbulent à l'aide de la Simulation des Grandes Echelles (SGE). Dans cette approche, les grandes échelles énergétiques et instationnaires de l'écoulement sont calcu- lées, tandis que les petites échelles sont modélisées. L'objectif de ce travail, réalisé en collaboration avec le Consortium Industrie-Recherche dans les Turbomachines (CIRT) et le LEMFI, est d'analyser les potentialités de la SGE pour le calcul d'écoulements confinés en géométrie 3D en vue d'applications aux turboma- chines. Le cas d'étude retenu pour la validation de la méthodologie est l'écoulement dans une tuyère 3D transsonique pour lequel il existe de nombreux résultats statistiques (LEMFI) et expérimentaux (ONERA). La tuyère est caractérisée par une bosse en flèche sur la paroi basse génératrice d'effets 3D. L'écoulement, sans rotation, présente cependant les phénomènes physiques complexes d'interaction onde de choc/couche limite et de large dé- collement, comme ceux rencontrés au sein des turbomachines. La SGE est obtenue à partir d'un modèle de sous maille d'échelles mixtes développé au LIMSI-CNRS. La discrétisa- tion temporelle des équations est réalisée par un schéma de Runge-Kutta explicite d'ordre deux. Les flux Euler sont discrétisés par un schéma TVD d'Harten-Yee d'ordre deux tandis que les flux visqueux le sont par un schéma centré. La SGE a permis d'obtenir des informations instationnaires sur l'écoulement, et de mettre en évidence la formation et le lâcher de tourbillons dans la tuyère. La solution insta- tionnaire est différente de la solution stationnaire RANS obtenue avec une modélisation statistique classique et montre l'oscillation du choc et la déstabilisation du décollement au cours du temps. Les résultats de la SGE obtenus sur le nombre de Mach isentropique, les profils de la vitesse moyenne et les tensions de Reynolds sont discutés et comparés aux résultats expérimentaux (ONERA) et statistiques, obtenus avec le modèle de Launder- Shima (LEMFI). Le bilan de l'énergie cinétique turbulente (k) est analysé et comparé à celui donné par la modélisation statistique, dans le décollement. Les résultats de la SGE montrent des différences notables avec les résultats statistiques dans le c÷ur de la tuyère : les termes de fluctuations de pression importants mesurés dans le choc sont pris en compte par la SGE.

La distribution de la matière dans l'Univers est supposée homogène et isotrope à très grande échelle mais l'observation de la distribution des galaxies lors<br />de grandes campagnes de recensements comme le SDSS nous montre un véritable réseau d'amas et de filaments sur des échelles de plusieurs centaines de<br />mégaparsecs.\\<br /><br />De nombreuses méthodes ont été développées dans le but de caractériser cette distribution et nous nous proposons dans cette thèse de présenter l'adaptation en trois dimensions d'un nouvel outil: le squelette.<br />Cette méthode vise à donner une définition mathématique claire des filaments ainsi qu'un algorithme numérique robuste permettant leur identification ainsi que le calcul de leurs propriétés. \\<br /><br />Afin de pouvoir comparer les résultats obtenus à partir des simulations N-corps de matière noire aux observations, une nouvelle méthode, baptisée MoLUSC, spécialisée dans la création de catalogues virtuels de<br />galaxies a aussi été élaborée. Elle se base sur les modèles semi-analytiques et est particulièrement efficace pour la fabrication de catalogues de grande<br />taille simulant de manière suffisamment réaliste les propriétés galactiques.\\<br /><br />Les utilisations de ces deux outils sont nombreuses et nous montrons par exemple qu'il est possible en mesurant la densité de longueur des filaments à<br />une échelle donnée de contraindre la quantité de matière dans l'univers<br />$\Omega_m$. Ces méthodes peuvent aussi être appliquées avec succès à la mesure statistique des propriétés du flux de matière noire le long des filaments, une mesure inédite. Nous présentons enfin de nombreuses applications possibles dont les résultats préliminaires sont très encourageants.

Cette thèse se situe dans le cadre du projet CANOE, piloté par GDF SUEZ et l'ADEME, qui vise à étudier la faisabilité et la viabilité du régime de combustion sans flamme dans les chaudières industrielles. Il est maintenant établi que le préchauffage des réactifs permet d'améliorer le rendement thermique et de diminuer la consommation de combustible d'une configuration. Pour contourner la formation d'oxydes d'azote résultant de l'augmentation de température des réactifs, ces derniers peuvent être massivement dilués par des produits de combustion. Cela permet d'éviter la formation de points chauds et d'homogénéiser les gradients de température: c'est le principe de la combustion sans flamme. L'objectif de cette thèse est de développer un modèle de combustion turbulente adapté à ce type de régime. La cinétique chimique complexe et le contrôle des pertes thermiques est d'une importance capitale dans l'établissement et la stabilisation du processus de combustion sans flamme. Ici, ces effets sont considérés dans une approche de tabulation de la cinétique chimique de type FPV. Pour discriminer les effets associés aux évolutions suivant les paramètres de contrôle de nos bases de données, on effectue une analyse des réponses de flammes laminaires à différents niveaux de dilution et de pertes thermiques. De plus, nous évaluons l'importance de l'utilisation de tabulations de dimensions élevées, et les capacités prédictives des méthodes de tabulation développées sont mises en exergue. Puis, des simulations aux grandes échelles de la turbulence de configurations adiabatique et à parois refroidies sont effectuées. On compare des tabulations de nombre de degrés de libertés variés aux données expérimentales. Les résultats obtenus sont en très bon accord avec ces dernières pour les tabulations les plus complexes, alors que des limitations significatives apparaissent pour des tabulations de dimensions inférieures. Les simulations proposées indiquent la capacité de nos modèles à reproduire des structures de flammes réalistes.

Les chantiers de maintenance et d’assainissement dans les centrales nucléaires nécessitent la mise en place d’enceintes ventilées autour des zones contaminées afin de limiter la propagation de la contamination à l’environnement extérieur. L’air rentre dans l’enceinte aux travers d’ouvertures sous la forme d’un écoulement directionnel, orienté de l’extérieur vers l’intérieur, assurant ainsi le confinement dynamique. En raison des opérations qui se déroulent à l’intérieur de l’enceinte et des perturbations externes, il est possible que l’écoulement de confinement dynamique aux ouvertures soit perturbé et subisse des inversions locales et instationnaires, conduisant ainsi à transporter la contamination à l’extérieur de l’enceinte. La présente étude s’intéresse aux petites ouvertures de type fentes minces rectangulaires où l’écoulement au droit de celles-ci est généralement turbulent. Les principaux objectifs de la thèse sont d’une part d’identifier les conditions aérodynamiques susceptibles de produire le phénomène de rétrodiffusion aux ouvertures, d’autre part d’évaluer la capacité des approches de modélisation de la turbulence URANS et LES à reproduire les instabilités liées à ce type d’écoulement. Il a été montré que l’apparition du phénomène de rétrodiffusion est principalement liée à la présence d’une perturbation aéraulique additionnelle, de type jet turbulent ou sillage, en compétition avec l’écoulement initial de confinement dynamique. Des expériences de traçage gazeux ont été mises en place sur une maquette expérimentale dans le but de quantifier la rétrodiffusion en fonction des différentes conditions aérauliques à l’ouverture et des caractéristiques de celle-ci. Des visualisations des écoulements à l’ouverture ont également été réalisées à l’aide d’un dispositif de tomographie laser. Enfin, l’analyse des résultats des simulations CFD a démontré que les approches de type RANS ou URANS ne permettaient pas de reproduire les instabilités de l’écoulement conduisant au phénomène de rétrodiffusion, contrairement aux simulations des grandes échelles de la turbulence (LES) qui reproduisent fidèlement les structures locales et instantanées à l’origine du phénomène.

L’objet de cette thèse est de modéliser et de simuler des écoulements turbulents diphasiques incompressibles à phases non miscibles. La modélisation et la simulation de ce type d’écoulements sont traitées dans le cadre des méthodes de Simulation des Grandes Echelles (SGE) ou Large Eddy Simulation (LES) en anglais qui consistent à calculer directement les plus grandes structures de l’écoulement et à modéliser les plus petites. Ces méthodes adaptées aux écoulements turbulents monophasiques sont étendues au cadre des écoulements turbulents diphasiques. Pour cela, elles sont couplées avec une méthode eulérienne de type ’ Volume Of Fluid’ (VOF) spécifique au caractère diphasique de l’écoulement. La pertinence du couplage entre les modélisations SGE et VOF est testée sur la configuration industrielle proposée par le CEA-CESTA: l’impact d’un jet rond turbulent sur une surface libre eau/air dans une cavité. Des mesures expérimentales de vitesse (Particle Image Velocimetry PIV) réalisées au CEA-CESTA sont disponibles pour valider les résultats numériques issus des simulations. / The scope of this dissertation is to model and simulate non-miscible two-phase turbulent and incompressible flows. The modeling and the simulation of this kind of flows are carried out in the framework of the Large Eddy Simulation (LES) which consists in calculating directly the largest structures of the flow and in modeling the finest ones. These numerical methods, applied usually to the simulation of single-phase turbulent flows, are extended to the simulation of two-phase turbulent flows in this work. Hence, the LES methods are coupled with an Eulerian ’Volume of Fluid’ (VOF) approach which is particularly adapted to interfacial flows. The relevance of this numerical coupling bewtween LES and VOF methods is validated in the following industrial configuration of the CEA-CESTA: the impact of a turbulent round jet on a free water/air surface in a cavity. Some experimental velocity measurements (Particle Image Velocimetry PIV), carried out at the CEA-CESTA, are available to validate the numerical results.

Le moteur à allumage commandé downsizé, couplé à une stratégie d’injection directe, est l’une des solutions privilégiées par les constructeurs automobiles afin de réduire les émissions polluantes et d’augmenter le rendement. Toutefois, l’augmentation de la pression d’injection visant à favoriser l’atomisation du spray et donc l’homogénéité du mélange peut engendrer une forte interaction entre le spray et les parois de la chambre de combustion. Cette interaction est à l’origine d’hétérogénéités locales susceptibles d’altérer la combustion. Du fait de son caractère instationnaire, l’interaction spray/paroi (formation et évaporation d’un film liquide) et plus généralement la préparation du mélange en moteur à injection directe essence sont des phénomènes difficiles à analyser expérimentalement. En effet, un moteur muni d’accès optiques ne peut pas fonctionner dans les conditions thermodynamiques réelles (pression, température...). Dans ce contexte, la modélisation et plus particulièrement la Simulation aux Grandes Echelles (“Large Eddy Simulation” LES) est un moyen d’analyse complémentaire et indispensable. L’objectif de cette thèse est de développer les modèles physiques nécessaires à la description de la phase liquide avec une approche Euler-Lagrange pour la simulation dans les moteurs à piston. Dans un premier temps, une modélisation des caractéristiques physiques du spray en sortie d’injecteur, nommée GDI, est proposée et validée par comparaison avec des mesures expérimentales. Les résultats montrent la capacité du modèle GDI à reproduire la dynamique générale d’un spray pour deux types d’injecteurs multi-trous. Dans un deuxième temps, deux modèles sont développés pour traiter respectivement l’interaction entre le spray et les parois et l’évaporation du film liquide. Les premières validations de ces modèles sont faites sur des expériences académiques dédiées, permettant des comparaisons précises avec les mesures. Finalement deux configurations moteur sont simulées. La première, sans combustion, permet d’évaluer l’impact d’une modélisation fine de l’interaction spray/paroi par rapport à une approche simplifiée. Les résultats montrent que la prise en compte de la formation et de l’évaporation du film liquide modifie significativement la formation du mélange, notamment le champ de richesse au Point Mort Haut. La seconde est utilisée pour analyser l’impact de la phase liquide sur le mélange et la combustion. Ces calculs sont comparés à des calculs réalisés sans injection liquide et à des mesures expérimentales. Les résultats mettent en évidence que les stratifications de richesse et de température, causées par l’évaporation du liquide, ont un effet de plissement sur la flamme et diminuent sa vitesse de propagation. / Downsized spark ignition engines coupled with a direct injection strategy, are more and more attractive for car manufacturers in order to reduce pollutant emissions and increase efficiency. However, the high pressure levels used to promote spray atomization and consequently mixing can generate a strong interaction between the spray and the combustion chamber walls. The combustion process may be affected by local heterogeneities caused by this interaction. Spray/walls interaction (formation and evaporation of the liquid film) and mixture preparation are unsteady phenomena, explaining why their experimental studies are limited. In fact, it is difficult to reproduce the thermodynamic conditions (pressure, temperature...) representative of an engine with optical accesses. In this context, numerical simulation, and in particular Large Eddy Simulation (LES) is a complementary mean of analysis. This work aims at developing the necessary models for the two-phase combustion simulation for engines, using an Euler-Lagrange approach. First, a modeling of the spray physics downstream to the injector exit is proposed and validated by comparison with experimental data. Second, two models are proposed and implemented to adress respectively the spray/wall interaction and the liquid film evaporation. These models are first validated on dedicated academics experiments, allowing an accurate comparison with experimental data. Then, two engine configurations are simulated. The first one, without combustion, allows the evaluation of an accurate spray/wall interaction modeling in comparison with a simplified approach. Results show that accounting for the formation and evaporation of the liquid film has a significant impact on the fuel/air mixing, especially on the equivalence ratio distribution at the Top Dead Center. The second one is used to analyze the impact of liquid on the mixing and the combustion. The simulations are compared to experiments data and to simulations assuming a perfect gaseous mixing (without liquid injection). Results show that the temperature and equivalence ratio heterogeneities, created by the liquid evaporation, have a wrinkling effect on the flame and reduce its propagation speed.

Dans le contexte aéronautique, cette thèse, financée par Safran Helicopter Engine, s’intéresse à l’application de l’approche Simulations aux Grandes Échelles (SGE) pour les systèmes de refroidissement de turbine. Le système de refroidissement industriel complexe est divisé en cas académiques plus simples donnant accès à des caractérisations expérimentales de la dynamique et de la thermique. Le jet impactant est traité en tant que système interne et l’écoulement de protection au bord de fuite en tant que système externe. Après une brève introduction du contexte lié au refroidissement de turbine et des objectifs scientifiques, ce manuscrit est divisé en 3 parties. La 1ère partie traite d’un écoulement de jet impactant sur plaque plane représentatif de l’impact à mi-corde. Elle se concentre sur la validation et la qualification des outils et modèles ainsi que sur l’analyse physique de l’écoulement. Les différentes instationarités de l’écoulement sont reliées à la thermique de paroi à l’aide de diagnostics statistiques et d’analyses modales. La 2ème partie s’intéresse à l’impact sur paroi concave représentatif de l’impact au bord d’attaque. Cette étude se concentre principalement sur la caractérisation de l’effet de courbure pour le jet impactant. Contrairement, au consensus actuel sur l’effet de courbure, la réduction des transferts thermiques est observée pour le cas d’étude de cette thèse. Au vu de ces résultats, une discussion est proposée pour tenter d’expliquer cet écart. Finalement, la 3ème partie de ce manuscrit contient une application de la SGE à un système de protection du bord de fuite par film isolant. Dans ce dispositif, des effets de groupe sont mis en évidence. L’impact des choix de modélisation tels que l’hypothèse de périodicité dans la direction de l’envergure est alors évalué. Il est montré que cette hypothèse de périodicité influe sur la prédiction locale de l’efficacité en forçant l’écoulement. En revanche, la prédiction de l’efficacité globale du système de protection n’est pas impactée. / This PhD thesis, funded by Safran Helicopter Engines, focuses on the application of the Large Eddy Simulation (LES) formalism to cooling systems present in high pressure turbine. The complex industrial problem is simplified into academic test cases for which experimental data are available for the validation process. The manuscript is divided into 3 parts dealing respectively with the impinging jet system on flat and concave plates and with the film cooling at the trailing edge equipped with a cutback on the pressure side. The 1st part deals with a jet impinging on a flat plate representing the impingement at mid-chord. This part focuses on the validation and qualification of the tools and models as well as on the physical analysis of the flow. The unsteadiness present in such an impinging jet flow are linked to the thermal behavior of the wall through the use of statistical analysis and modal decomposition of the flow field. The 2nd part is dedicated to the study of a jet impinging on a concave surface. This study aims at characterizing the effect of curvature for an impinging jet flow. The results found in this study disagreed with the current consensus attributing heat transfer enhancement on concave surface to Gortler instability. Hence, a discussion is proposed in an attempt to explain this discrepancy. Finally, the 3rd part reports an LES of the film cooling at the trailing edge. Group effects, due to the presence of internal ribs, are highlighted for the configuration studied here. These simulations use a spatial periodicity assumption to reduce the size of the computational domain. It is shown that this specific assumption is not suited as it forces the flow and modifies the group effect. The local results, in terms of adiabatic effectiveness, are found to be sensitive to such a forcing. However, the global behavior of the effectiveness is not impacted by this periodic boundary condition.

Cette thèse se situe dans le cadre du projet CANOE, piloté par GDF SUEZ et l'ADEME, qui vise à étudier la faisabilité et la viabilité du régime de combustion sans flamme dans les chaudières industrielles. Il est maintenant établi que le préchauffage des réactifs permet d'améliorer le rendement thermique et de diminuer la consommation de combustible d'une configuration. Pour contourner la formation d'oxydes d'azote résultant de l'augmentation de température des réactifs, ces derniers peuvent être massivement dilués par des produits de combustion. Cela permet d'éviter la formation de points chauds et d'homogénéiser les gradients de température: c'est le principe de la combustion sans flamme. L'objectif de cette thèse est de développer un modèle de combustion turbulente adapté à ce type de régime. La cinétique chimique complexe et le contrôle des pertes thermiques est d'une importance capitale dans l'établissement et la stabilisation du processus de combustion sans flamme. Ici, ces effets sont considérés dans une approche de tabulation de la cinétique chimique de type FPV. Pour discriminer les effets associés aux évolutions suivant les paramètres de contrôle de nos bases de données, on effectue une analyse des réponses de flammes laminaires à différents niveaux de dilution et de pertes thermiques. De plus, nous évaluons l'importance de l'utilisation de tabulations de dimensions élevées, et les capacités prédictives des méthodes de tabulation développées sont mises en exergue. Puis, des simulations aux grandes échelles de la turbulence de configurations adiabatique et à parois refroidies sont effectuées. On compare des tabulations de nombre de degrés de libertés variés aux données expérimentales. Les résultats obtenus sont en très bon accord avec ces dernières pour les tabulations les plus complexes, alors que des limitations significatives apparaissent pour des tabulations de dimensions inférieures. Les simulations proposées indiquent la capacité de nos modèles à reproduire des structures de flammes réalistes. / This thesis is within a framework of the CANOE project, under the responsibility of GDF SUEZ and the ADEME, and aims at studying the feasibility of the flameless combustion regime in industrial boilers. It is now well established that reactants preheating allows an improvement on thermal efficiency as well as fuel savings. To avoid an increase in nitrogen oxides emissions arising from reactants temperature augmentation, massive dilution of the reactants by burnt gases can be used. While doing so, hot spots are averted and temperature gradients are smoothed: it’s the principle of flameless combustion. Even though this combustion regime is a subject of increasing interest to the industry, its mechanisms are not yet fully understood. The objective of this work is to develop and validate a new turbulent combustion model adapted to this kind of regimes. Complex chemistry as well as heat losses control is of paramount importance in the establishment and stabilization of the flameless combustion process. Here, these effects are taken into account in an FPV-type chemistry tabulation approach. To discriminate the effects associated with evolutions along the database control parameters, we analyze laminar flame responses to different levels of dilution and heat losses. Moreover, we assess the importance of high-order tabulations, and predictive capabilities of our approaches are highlighted. Then, large-eddy simulations of laboratory scale experiments (an adiabatic and a cooled wall configuration) are carried out. To this end, we compare databases featuring different number of dimensions to experimental data on temperature and species distributions available in the literature. Simulation results are in very good agreement with experimental data for complex tabulations, while discrepancies arise for lower order tabulations. Simulation results show that flameless combustion features a wide variety of flame structures, and that our models are able to tackle realistic flame structures.