<p><br></p> <p>A  combustion simulation involves various physiochemical processes, such as molecular and turbulent diffusion, smoke and soot formation, thermal radiation, chemical reaction mechanisms, and kinetics. In the last decade, computational fluid dynamics (CFD) has been increasingly used in combustion modeling. It is critically important to improve and enhance the predictive capabilities of combustion models. This work presents an analysis of two types of diffusion flames: the momentum-dominant jet flames and buoyancy-controlled pool fires. The gap between the existing knowledge of differential molecular diffusion in turbulent high momentum jet flow and the practical applications has been reduced. The importance of mixing modeling in pool fire simulations has been revealed, and enhancement for predicting fire extinction limits has been proposed.</p> <p><br></p> <p>Modeling differential molecular diffusion in turbulent non-premixed combustion remains a great challenge for flamelet models. The laminar flamelet is a key component of a flamelet model for turbulent combustion. One significant challenge that has not been well addressed is the representativity of laminar flamelet for the characteristics of differential molecular diffusion in turbulent combustion problems. Laminar flamelet is generated typically based on two conceptual burner configurations, the opposed jet burner, and the Tsuji burner. They are commonly considered equivalent when dealing with the description of laminar flamelet structures. A difference between them is revealed in this work for the first time when they are used to represent differential molecular diffusion. The traditionally opposed jet burner yields an almost fixed equal diffusion location in the mixture fraction space for the transport of different elements. The Tsuji burner can produce a continuous variation of the equal diffusion location in the mixture fraction space with a slight extension. This variation of the equal diffusion location is shown to be an essential characteristic of turbulent non-premixed combustion, as demonstrated in a laminar jet mixing layer problem, a turbulent jet mixing layer problem, and a turbulent jet non-premixed flame. The Tsuji burner is thus potentially a more suitable choice than the opposed jet burner for laminar flamelet generation that can be consequently used in flamelet modeling of differential molecular diffusion for turbulent non-premixed combustion.</p> <p><br></p> <p>Capturing fire extinction limits in simulations is essential for developing predictive capabilities for fire. In this work, the combined large-eddy simulation (LES) and transported probability density function (PDF) methods are assessed for the predictions of fire extinction. The University of Maryland line burner is adopted as a validation test case. The NIST Fire Dynamics Simulator (FDS) code for LES is combined with an in-house PDF code called HPDF for the fire simulations. The simulation results were verified by using the available experimental data. The combustion efficiency under the different oxygen depletion levels in the oxidizer is analyzed. Fire extinction occurs when the oxygen depletion level reduces to a certain level. The model’s capability to capture this extinction limit is assessed by using the experimental data. Different mixing models and model parameters are examined. It is found that the fire extinction limit is very sensitive to the different mixing models and mixing parameters. The level of sensitivity is higher than in momentum-driven turbulent flames, which suggests the importance of mixing modeling in fire simulations. The existing mixing models need further enhancement for predicting fire extinction. </p> <p><br></p>

A numerical study of turbulent combustion systems was pursued to examine different computational modeling techniques, namely computational fluid dynamics (CFD) and chemical reactor network (CRN) methods. Both methods have been studied and analyzed as individual techniques as well as a coupled approach to pursue better understandings of the mechanisms and interactions between turbulent flow and mixing, ignition behavior and pollutant formation. A thorough analysis and comparison of both turbulence models and chemistry representation methods was executed and simulations were compared and validated with experimental works. An extensive study of turbulence modeling methods, and the optimization of modeling techniques including turbulence intensity and computational domain size have been conducted. The final CFD model has demonstrated good predictive performance for different turbulent bluff-body flames. The NOx formation and the effects of fuel mixtures indicated that the addition of hydrogen to the fuel and non-flammable diluents like CO2 and H2O contribute to the reduction of NOx. The second part of the study focused on developing chemical models and methods that include the detailed gaseous reaction mechanism of GRI-Mech 3.0 but cost less computational time. A new chemical reactor network has been created based on the CFD results of combustion characteristics and flow fields. The proposed CRN has been validated with the temperature and species emission for different bluff-body flames and has shown the capability of being applied to general bluff-body systems. Specifically, the rate of production of NOx and the sensitivity analysis based on the CRN results helped to summarize the reduced reaction mechanism, which not only provided a promising method to generate representative reactions from hundreds of species and reactions in gaseous mechanism but also presented valuable information of the combustion mechanisms and NOx formation. Finally, the proposed reduced reaction mechanism from the sensitivity analysis was applied to the CFD simulations, which created a fully coupled process between CFD and CRN, and the results from the reduced reaction mechanism have shown good predictions compared with the probability density function method. / Ph. D.

This thesis is concerned with investigations on two related issues of turbulent jet diffusion flame, namely (a) stabilization at liftoff and (b) autoignition in a turbulent jet diffusion flame. The approach of Conditional Moment Closure (CMC) has been taken. Fully elliptic first order CMC equations are solved with detailed chemistry to simulate lifted H2/N2 flame in vitiated coflow. The same approach is further used to simulate transient autoignition process in inhomogeneous mixing layers. In Chapter 1, difficulties involved in numerical simulation of turbulent combustion problems are explained. Different numerical tools used to simulate turbulent combustion are briefly discussed. Previous experimental, theoretical and numerical studies of lifted jet diffusion flames and autoignition are reviewed. Various research issues related to objectives of the thesis are discussed. In Chapter 2, the first order CMC transport equations for the reacting flows are presented. Various closure models that are required for solving the governing equations are given. Calculation of mean reaction rate term for detailed chemistry is given with special focus on the reaction rates for pressure dependent reactions. In Chapter 3, starting with the laminar flow code, further extension is carried to include kε turbulence model and PDF model. The code is validated at each stage of inclusion of different model. In this chapter, the code is first validated for the test problem of constant density, 2D, axisymmetric turbulent jet. Further, validation of PDF model is carried out by simulating the problem of nonreacting jet of cold air issuing into a vitiated coflow. The results are compared with the published data from experiments as well as numerical simulations. It is shown that the results compare well with the data. In Chapter 4, numerical results of lifted jet diffusion flame are presented. Detailed chemistry is modelled using Mueller mechanism for H2/O2 system with 9 species and 21 reversible reactions. Simulations are carried out for different jet velocities and coflow stream temperatures. The predicted liftoff generally agrees with experimental data, as well as joint PDF results. Profiles of mean scalar fluxes in the mixture fraction space, for different coflow temperatures reveal that (1) Inside the flamezone, the chemical term balances the molecular diffusion term, and hence the structure is of a diffusion flamelet for both cases. (2) In the preflame zone, the structure depends on the coflow temperature: for low coflow temperatures, the chemical term being small, the advective term balances the axial diffusion term. However, for the high coflow temperature case, the chemical term is large and balances the advective term, the axial diffusion term being small. It is concluded that, liftoff is controlled (a) by turbulent premixed flame propagation for low cofflow temperature while (b) by autoignition for high coflow temperature. In Chapter 5, the numerical results of autoignition in inhomogeneous mixing layer are presented. The configuration consists of a fuel jet issued into hot air for which transient simulations are performed. It is found that the constants assumed in various modelling terms can severely influence the results, particularly the flame temperature. Hence, modifications to these constants are suggested to obtain improved predictions. Preliminary work is carried out to predict autoignition lengths (which may be defined by Tign × Ujet incase of jet- and coflowvelocities being equal) by varying the coflow temperature. The autoignition lengths show a reasonable agreement with the experimental data and LES results. In Chapter 6, main conclusions of this thesis are summarized. Possible future studies on this problem are suggested.

Les travaux réalisés durant la thèse s'inscrivent dans une problématique scientifique liée à l'étude des transferts radiatifs. Plus particulièrement, l'application de cette étude est la télédétection infrarouge d'un écoulement diphasique réactif et turbulent à haute température. Cette étude a pour objectif la modélisation et la simulation du rayonnement infrarouge émis par cet écoulement et re¸cu par un détecteur. Pour développer un outil de simulation numérique de la signature infrarouge d'un jet de gaz et de particules à haute température, les espèces majoritaires qui sont responsables de l'émission du rayonnement ont été identifiées lors d'expérimentations. Les campagnes expérimentales ont permis la construction de bases de données concernant les gaz (CO2-CO-H2O) et les particules (oxydes métalliques) présents dans le jet. Connaissant la nature des gaz et des particules, le calcul de leurs propriétés radiatives doit être réalisé. Cette étape est nécessaire puisque ces propriétés caractérisent l'émission de rayonnement par le jet et elles doivent être connues pour résoudre l'équation de transfert radiatif. Pour les gaz, un code de calcul raie par raie de spectres synthétiques a été développé. De plus, pour diminuer le temps de calcul d'une signature infrarouge, il est préférable d'utiliser des modèles spectraux de bandes étroites. Le modèle de télédétection infrarouge est un modèle spectral utilisant des k(coefficient d'absorption)-distributions sous l'hypothèse des k-corrélés avec l'approximation d'un gaz unique pour le mélange associée à l'hypothèse des gaz fictifs. Les paramètres de ce modèle (CKFG-SMG), ont été tabulés et validés dans l'étude. En ce qui concerne les propriétés radiatives des nuages de particules sphériques, le modèle de Mie est utilisé car il est valable pour les gammes de fractions volumiques rencontrées. Pour tester l'influence de la diffusion, une étude de sensibilité à la diffusion a été réalisée. En effet, nous avons quantifié l'erreur commise sur le flux émis par différentes couches si les processus de diffusion du rayonnement sont négligés. Cette étude a montré que l'influence de la diffusion peut être négligée dans le cadre de notre étude. La modélisation de la signature infrarouge du jet diphasique réactif issu de la combustion du matériau pyrotechnique, nécessite la connaissance des températures et des concentrations en gaz et particules, en tous les points du jet. Ce jet diphasique réactif a été simulé à l'aide du logiciel Fluent. De plus, une interface graphique a été développée qui recrée la scène optronique en se servant des profils aérothermochimiques du jet diphasique et des données concernant la position du détecteur. De cette fa¸con, un outil de simulation numérique de la signature infrarouge du jet (SIRJET) a été développé qui inclue un modèle de transfert radiatif (lancer de rayon) ainsi que les paramètres tabulés (gaz et particules) du modèle spectral de télédétection infrarouge (CK, CKFG, CK-SMG, CKFG-SMG). Enfin, une confrontation est présentée entre une mesure et le résultat d'une simulation de la signature infrarouge d'un jet diphasique à haute température.

Turbulent premixed propane/air and methane/air flames were studied using planar Rayleigh scattering and particle image velocimetry on a stabilized Bunsen type burner. The fuel-air equivalence ratio was varied from Φ=0.7 to 1.0 for propane flames, and from Φ=0.6 to 1.0 for methane flames. The non-dimensional turbulence intensity, u'/SL (ratio of fluctuation velocity to laminar burning velocity), covered the range from 3 to 24, equivalent to conditions of corrugated flamelets and thin reaction zones regimes. Temperature gradients decreased with the increasing u'/SL and levelled off beyond u'/SL > 10 for both propane and methane flames. Flame front thickness increased slightly as u'/SL increased for both mixtures, although the thickness increase was more noticeable for propane flames, which meant the thermal flame front structure was being thickened. A zone of higher temperature was observed on the average temperature profile in the preheat zone of the flame front as well as some instantaneous temperature profiles at the highest u'/SL. Curvature probability density functions were similar to the Gaussian distribution at all u'/SL for both mixtures and for all the flame sections. The mean curvature values decreased as a function of u'/SL and approached zero. Flame front thickness was smaller when evaluated at flame front locations with zero curvature than that with curvature. Temperature gradients and FSD were larger when the flame curvature was zero. The combined thickness and FSD data suggest that the curvature effect is more dominant than that of the stretch by turbulent eddies during flame propagation. Integrated flame surface density for both propane and methane flames exhibited no dependance on u'/SL regardless of the FSD method used for evaluation. This observation implies that flame surface area may not be the dominant factor in increasing the turbulent burning velocity and the flamelet assumption may not be valid under the conditions studied. Dκ term, the product of diffusivity evaluated at conditions studied and the flame front curvature, was a magnitude smaller than or the same magnitude as the laminar burning velocity.

The use of Computational Fluid Dynamics (CFD) tools is crucial for the development of novel and cost-effective combustion technologies and the minimization of environmental concerns at industrial scale. CFD simulations facilitate scaling-up procedures that otherwise would be complicated by strong interactions between reaction kinetics, turbulence and heat transfer. CFD calculations can be applied directly at the industrial scale of interest, thus avoiding scaling-up from lab-scale experiments. However, this advantage can only be obtained if CFD tools are quantitatively predictive and trusted as so. Despite the improvements in the computational capability, the implementation of detailed physical and chemical models in CFD simulations can still be prohibitive for real combustors, which require large computational grids and therefore significant computational efforts. Advanced simulation approaches like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) guarantee higher fidelity in computational modeling of combustion at, unfortunately, increased computational cost. However, with adequate, reduced, and cost-effective modeling of physical phenomena, such as chemical kinetics and turbulence-chemistry interactions, and state of the art computing, LES will be the tool of choice to describe combustion processes at industrial scale accurately. Therefore, the development of reduced physics and chemistry models with quantified model-form uncertainty is needed to overcome the challenges of performing LES of industrial systems. Reduced-order models must reproduce the main features of the corresponding detailed models. They feature predictivity and capability of bridging scales when validated against a broad range of experiments and targeted by Validation and Uncertainty Quantification (V/UQ) procedures. In this work, V/UQ approaches are applied for reduced-order modeling of pulverized coal devolatilization and subsequent char oxidation, and furthermore for modeling NOx emissions in combustion systems.For coal devolatilization, a benchmark of the Single First-Order Reaction (SFOR) model was performed concerning the accuracy of the prediction of volatile yield. Different SFOR models were implemented and validated against experimental data coming from tests performed in an entrained flow reactor at oxy-conditions, to shed light on their drawbacks and benefits. SFOR models were chosen because of their simplicity: they can be easily included in CFD codes and are very appealing in the perspective of LES of pulverized coal combustion burners. The calibration of kinetic parameters was required to allow the investigated SFOR model to be predictive and reliable for different heating rates, hold temperatures and coal types. A comparison of several calibration approaches was performed to determine if one-step models can be adaptive and able to bridge scales, without losing accuracy, and to select the calibration method to employ for wider ranges of coal rank and operating conditions. The analysis pointed out that the main drawback of the SFOR models is the assumption of a constant ultimate volatile yield, equal to the value from the coal proximate analysis. To overcome this drawback, a yield model, i.e. a simple functional form that relates the ultimate volatile yield to the particle temperature, was proposed. The model depends on two parameters that have a certain degree of uncertainty. The performances of the yield model were assessed using a collaboration of experiments and simulations of a pilot-scale entrained flow reactor. A consistency analysis, based on the Bound-to-Bound Data Collaboration (B2B-DC) approach, and a Bayesian method, based on Gaussian Process Regression (GPR), were employed for the investigation of experiments and simulations. In Bound-to- Bound Data Collaboration the model output, evaluated at specified values of the model parameters, is compared with the experimental data: if the prediction of the model falls within the experimental uncertainty, the corresponding parameter values would be included in the so-called feasible set. The existence of a non-empty feasible set signifies consistency between the experiments and the simulations, i.e. model-data agreement. Consistency was indeed found when a relative error of 19% for all the experimental data was applied. Hence, a feasible set of the two SFOR model parameters was provided. A posterior state of knowledge, indicating potential model forms that could be explored in yield modeling, was obtained by Gaussian Process Regression. The model form evaluated through the consistency analysis is included within the posterior derived from GPR, indicating that it can satisfactorily match the experimental data and provide reliable estimation in almost every range of temperatures. CFD simulations were carried out using the proposed yield model with first-order kinetics, as in the SFOR model. Results showed promising agreement between predicted and experimental conversion for all the investigated cases.Regarding char combustion modeling, the consistency analysis has been applied to validate a reduced-order model and quantify the uncertainty in the prediction of char conversion. The model capability to address heterogeneous reaction between char carbon and O2, CO2 and H2O reagents, mass transport of species in the particle boundary layer, pore diffusion, and internal surface area changes was assessed by comparison with a large number of experiments performed in air and oxy-coal conditions. Different model forms had been considered, with an increasing degree of complexity, until consistency between model outputs and experimental results was reached. Rather than performing forward propagation of the model-form uncertainty on the predictions, the reduction of the parameter uncertainty of a selected model form was pursued and eventually achieved. The resulting 11-dimensional feasible set of model parameters allows the model to predict the experimental data within almost ±10% uncertainty. Due to the high dimensionality of the problem, the employed surrogate models resulted in considerable fitting errors, which led to a spoiled UQ inverse problem. Different strategies were taken to reduce the discrepancy between the surrogate outputs and the corresponding predictions of the simulation model, in the frameworks of constrained optimization and Bayesian inference. Both strategies succeeded in reducing the fitting errors and also resulted in a least-squares estimate for the simulation model. The variety of experimental gas environments ensured the validity of the consistent reduced model for both conventional and oxy-conditions, overcoming the differences in mass transport and kinetics observed in several experimental campaigns.The V/UQ-aided modeling of coal devolatilization and char combustion was done in the framework of the Predictive Science Academic Alliance Program II (PSAAP-II) funded by the US Department of Energy. One of the final goals of PSAAP-II is to develop high-fidelity simulation tools that ensure 5% uncertainty in the incident heat flux predictions inside a 1.2GW Ultra-Super-Critical (USC) coal-fired boiler. The 5% target refers to the expected predictivity of the full-scale simulation without considering the uncertainty in the scenario parameters. The data-driven approaches used in this Thesis helped to improve the predictivity of the investigated models and made them suitable for LES of the 1.2GW USC coal-fired boiler. Moreover, they are suitable for scale-bridging modeling of similar multi-phase processes involved in the conversion of solid renewable sources, such as biomass.In the final part of the Thesis, the sensitivity to finite-rate chemistry combustion models and kinetic mechanisms on the prediction of NO emissions was assessed. Moreover, the forward propagation of the uncertainty in the kinetics of the NNH route (included in the NOx chemistry) on the predictions of NO was investigated to reveal the current state of the art of kinetic modeling of NOx formation. The analysis was carried out on a case where NOx formation comes from various formation routes, both conventional (thermal and prompt) and unconventional ones. To this end, a lab-scale combustion system working in Moderate and Intense Low-oxygen Dilution (MILD) conditions was selected. The results showed considerable sensitivity of the NO emissions to the uncertain kinetic parameters of the rate-limiting reactions of the NNH pathway when a detailed kinetic mechanism is used. The analysis also pointed out that the use of one-step global rate schemes for the NO formation pathways, necessary when a skeletal kinetic mechanism is employed, lacks the required chemical accuracy and dims the importance of the NNH pathway in this combustion regime. An engineering modification of the finite-rate combustion model was proposed to account for the different chemical time scales of the fuel-oxidizer reactions and NOx formation pathways. It showed an equivalent impact on the emissions of NO than the uncertainty in the kinetics of the NNH route. At the cost of introducing a small mass imbalance (of the order of ppm), the adjustment led to improved predictions of NO. The investigation established a possibility for the engineering modeling of NO formation in MILD combustion with a finite-rate chemistry combustion model that can incorporate a detailed mechanism at affordable computational costs. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

Two new models to calculate the species instantaneous and filtered reaction rates for multi-step, multi-species chemical kinetics mechanisms are developed based on the artificial neural networks (ANN) approach. The proposed methodologies depend on training the ANNs off-line on a thermo-chemical database representative of the actual composition and turbulence level of interest. The thermo-chemical database is constructed by stand-alone linear eddy mixing (LEM) model simulations under both premixed and non-premixed conditions, where the unsteady interaction of turbulence with chemical kinetics is included as a part of the training database. In this approach, the information regarding the actual geometry of interest is not needed within the LEM computations. The developed models are validated extensively on the large eddy simulations (LES) of (i) premixed laminar-flame-vortex-turbulence interaction, (ii) temporally mixing non-premixed flame with extinction-reignition characteristics, and (iii) stagnation point reverse flow combustor, which utilizes exhaust gas re-circulation technique. Results in general are satisfactory, and it is shown that the ANN provides considerable amount of memory saving and speed-up with reasonable and reliable accuracy. The speed-up is strongly affected by the stiffness of the reduced mechanism used for the computations, whereas the memory saving is considerable regardless.

Turbulent premixed propane/air and methane/air flames were studied using planar Rayleigh scattering and particle image velocimetry on a stabilized Bunsen type burner. The fuel-air equivalence ratio was varied from Φ=0.7 to 1.0 for propane flames, and from Φ=0.6 to 1.0 for methane flames. The non-dimensional turbulence intensity, u'/SL (ratio of fluctuation velocity to laminar burning velocity), covered the range from 3 to 24, equivalent to conditions of corrugated flamelets and thin reaction zones regimes. Temperature gradients decreased with the increasing u'/SL and levelled off beyond u'/SL > 10 for both propane and methane flames. Flame front thickness increased slightly as u'/SL increased for both mixtures, although the thickness increase was more noticeable for propane flames, which meant the thermal flame front structure was being thickened. A zone of higher temperature was observed on the average temperature profile in the preheat zone of the flame front as well as some instantaneous temperature profiles at the highest u'/SL. Curvature probability density functions were similar to the Gaussian distribution at all u'/SL for both mixtures and for all the flame sections. The mean curvature values decreased as a function of u'/SL and approached zero. Flame front thickness was smaller when evaluated at flame front locations with zero curvature than that with curvature. Temperature gradients and FSD were larger when the flame curvature was zero. The combined thickness and FSD data suggest that the curvature effect is more dominant than that of the stretch by turbulent eddies during flame propagation. Integrated flame surface density for both propane and methane flames exhibited no dependance on u'/SL regardless of the FSD method used for evaluation. This observation implies that flame surface area may not be the dominant factor in increasing the turbulent burning velocity and the flamelet assumption may not be valid under the conditions studied. Dκ term, the product of diffusivity evaluated at conditions studied and the flame front curvature, was a magnitude smaller than or the same magnitude as the laminar burning velocity.

Le syngas est un mélange gazeux de CO et H2 qui constitue un intermédiaire important dans l’industrie pétrochimique. Plusieurs approches sont utilisées pour le produire. L’oxydation partielle non catalytique (POX) et le reformage à la vapeur (SMR) en font partie. Le reformage auto thermique du méthane (ATR) combine quant à lui ces deux procédés au sein d’un même réacteur. L’amélioration du rendement global du procédé ATR requiert une meilleure caractérisation du comportement des gaz au sein de la chambre. La simulation numérique apparaît comme un outil efficace pour y parvenir. Pour réduire le coût CPU, c'est généralement l'approche RANS (Reynolds Average Numerical Simulation) qui est privilégiée pour la simulation complète de la chambre. Cette approche repose sur l'utilisation de modèles, parmi lesquels le modèle de combustion turbulente, qui a pour objectif de représenter les interactions entre la turbulence et la réaction chimique au sein du mélange. Plusieurs stratégies ont été proposées pour le calculer, qui bénéficient globalement d'une large expérience pour les systèmes classiques mettant en jeu la combustion. Cependant, les flammes observées dans les réacteurs ATR présentent des propriétés assez différentes de ces configurations classiques. La validité des modèles de combustion turbulente classiques doit donc y être vérifiée. L'objectif de cette thèse est de répondre à ce besoin, en testant la validité de différents modèles de combustion turbulente. La première partie du travail a consisté à analyser les propriétés des flammes CH4/O2 enrichies en vapeur d'eau à haute pression, et a notamment permis le développement d’une méthode d’évaluation des temps caractéristiques d’un système chimique. Dans un deuxième temps, une expérience numérique à l’aide d’un code DNS a été réalisée, afin de servir de référence pour tester a priori sur des configurations ATR plusieurs modèles RANS de combustion turbulente couramment utilisés dans le milieu industriel. / Syngas is a gaseous mixture mainly composed of CO and H2, which constitutes a major feedstock in petrochemical industry. Several industrial approaches are commonly used to produce it. Non catalytic Partial Oxidation (POX) and Steam Methane Reforming (SMR) are two of them. Autothermal Reforming (ATR) is a third process that combines both POX and SMR in the same reactor. A better knowledge of the reactive flow properties inside the chamber is required in order to improve the ATR process efficiency. Numerical simulation appears as an efficient tool to reach this goal. Because of the high CPU cost required for these simulations, RANS (Reynolds Average Numerical Simulation) formulation is usually preferred for the simulation of the whole chamber. This approach relies on the use of models, like the turbulent combustion model that aims at describing the interactions between turbulence and chemical reactions. Several approaches have been proposed to compute it, which benefit from a relatively wide experience for the simulation of classical combustion systems. However, ATR flames have some specific properties that make them quite different from these classical configurations, especially because of high pressure, reactants dilution with water and high global equivalence ratio. The validity of classical turbulent combustion models therefore requires to be assessed in ATR configurations. The objective of this thesis is to meet this need by testing the validity of several turbulent combustion models. The first part of this work has been to analyze water-enriched CH4/O2 flames properties at high pressure. In particular, a strategy for evaluating characteristic chemical time scales of a reactive system has been proposed within this context. In a second part, a DNS numerical experiment has been performed. Its results are then used as a benchmark for a priori testing several turbulent combustion models in the context of ATR reactor RANS simulations.

El modelado CFD se ha convertido en una herramienta aceptada y ampliamente utilizada en el ámbito del diseño de motores de combustión interna alternativos. Los modelos de combustión avanzados ayudan a comprender los fenómenos complejos químicos y físicos del proceso de combustión y aportan información detallada que no se puede obtener con experimentos. Indudablemente, el modelado del proceso de combustión turbulenta parcialmente premezclada característico de los chorros Diesel es particularmente difícil y por lo tanto es un tema de gran interés para la comunidad científica. Los retos más importantes del modelado de este tipo de llamas son la predicción del proceso del auto-encendido, caracterizado por el tiempo de retraso, y la estructura de la llama cuasi-estacionaria con su característica longitud de lift-off. Estos dos parámetros globales de los chorros Diesel son importantes por varios aspectos. Primero, es relativamente sencillo medir estos dos parámetros y por lo tanto utilizarlos para la validación de modelos y segundo, son factores determinantes en el proceso de la combustión en un motor. El auto-encendido marca el inicio de la tasa de liberación de calor y la longitud de lift-off desempeña un papel fundamental en la formación de hollín. El mecanismo de estabilización de la llama en la zona del lift-off todavía no es bien conocido aunque existen diferentes teorías en la literatura, por lo que su modelado es en la actualidad un reto no resuelto. De acuerdo con el contexto descrito previamente, en este trabajo se pretende implementar un modelo de combustión integrado en un solver RANS utilizando la plataforma CFD OpenFOAM de código abierto. El modelo propuesto está basado en el concepto de flamelets usando una química detallada combinado con funciones de probabilidad determinadas a priori (presumed-PDF) para considerar el efecto de interacción entre la química y las características del flujo turbulento, que implica hipótesis importantes. En primer lugar, con el concepto flamelet se asume que una llama Diesel turbulenta quema localmente como un conjunto de llamas laminares de difusión de flujos opuestos. En segundo lugar se asume que las fluctuaciones de las propiedades introducidas por el flujo turbulento, que son las responsables de los fenómenos de interacción entre la química y la turbulencia durante la combustión, siguen un comportamiento estadístico en el tiempo de acuerdo a una distribución de probabilidad conocida a priori. Los fenómenos complejos del auto-encendido de hidrocarburos exigen el uso de mecanismo químicos detallados para recuperar satisfactoriamente los tiempos de retraso del auto-encendido en un rango amplio de condiciones termoquímicas. Una estrategia de interés para mantener los costes computacionales dentro de límites aceptables consiste en pre-tabular los resultados del cálculo de la química en tablas. Los parámetros independientes de estas tablas son la fracción de mezcla, la variable de progreso y la tasa de disipación escalar. Además, la hipótesis de que la distribuciones de probabilidad de las fluctuaciones generadas por la turbulencia sobre las propiedades del flujo son conocidas permite generar una tabla con la información química del problema apta para su aplicación en cálculo CFD en un entorno RANS. Esta aproximación basada en la pre-tabulación de los resultados químicos presenta dos ventajas fundamentales, siendo la primera de ellas la posibilidad de considerar modelos avanzados de interacción química-turbulencia y la segunda el relevante ahorro de tiempo de cálculo. Sin embargo, estas tablas representan un gran espacio de datos cuya gestión eficiente no es trivial. El desarrollo de un almacenamiento adecuado para un acceso de datos rápido y directo así como un esquema de interpolación multidimensional también forma parte del presente trabajo. / Winklinger, JF. (2014). Implementation of a combustion model based on the flamelet concept and its application to turbulent reactive sprays [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/48488 / TESIS