Sensitivity analysis was performed on a soot model using a partially stirred reactor (PaSR) in order to determine the effects of mixing model parameters on soot scalar values. The sensitivities of the mixture fraction zeta and progress variable C to the mixing model constant C_phi were calculated; these values were used to compute the sensitivity of water mass fraction Y_H2O to C_phi and several soot quantities to soot moments. Results were validated by evaluating the mean mixture fraction sensitivity and a long simulation time case. From the baseline case, it was noted that soot moment sensitivities tended to peak on the rich side of the stoichiometric mixture fraction zeta_st. Timestep, number of notional particles, mixing timescale tau_mix, and residence time tau_res were varied independently. Choices for timestep and notional particle count were shown to be sufficient to capture relevant scalar profiles, and did not greatly affect sensitivity calculations. Altering tau_mix or tau_res was shown to affect sensitivity to mixing, and it was concluded that the soot model is more heavily influenced by the chemistry than mixing. / text

Soot in coal combustion simulations is often ignored due to its computational complexity, despite significant effects on flame temperature and radiation. In this research, a 40 kW oxy-coal combustion system is modeled using Large Eddy Simulations (LES) and a semi-empirical monodisperse coal soot model. Simulation results are compared to experimental measurements of temperature, species concentrations, and soot concentration. Cases where soot is modeled are compared with cases where soot is neglected to determine the accuracy benefits of modeling soot. The simulations were able to replicate experimental results within an acceptable level of error. Including soot in the simulations did not consistently increase accuracy for the simulation setup and modeling assumptions used in this research.

[ES] En los últimos años, las emisiones de hollín provenientes de los motores de combustión interna han recibido más atención debido al impacto negativo que éstas tienen no solo en el ambiente, sino también en la salud del ser humano. Como respuesta, leyes cada vez más estrictas han sido aplicadas impulsando así a la comunidad científica al desarrollo de motores más eficientes en el uso del combustible y por supuesto más limpios en términos de emisiones contaminantes. En este contexto, el modelado computacional ha sido la herramienta utilizada en numerosos esfuerzos que buscan contribuir a mejorar el entendimiento que se tiene sobre los altamente complejos fenómenos que componen el proceso de producción de hollín. El principal objetivo de esta tesis es simular la producción de hollín en chorros Diesel en condiciones de operación típicas de un motor de combustión interna utilizando CFD. La consecución del objetivo de la tesis comprende una evaluación preliminar de la configuración de los distintos modelos para el caso de chorros inertes. En segundo lugar, el estudio detallado de la hipótesis utilizada para caracterizar la estructura de la llama a nivel sub-grid (tomando como base los conceptos well-mixed o flamelet) y del enfoque para tener en cuenta la interacción entre turbulencia y química. Por último, se presentan resultados del modelado de la combustión y producción de hollín para diferentes condiciones de contorno de reactividad y mezcla del chorro utilizando un modelo de hollín de dos ecuaciones. En resumen, el lector encontrará a lo largo de este documento un estudio exhaustivo sobre la combustión y producción de hollín en chorros inyectados con toberas mono-orificio en ambientes quiescentes. De este tipo de chorros, el Spray A y Spray D de la Engine Combustion Network son utilizados como casos de referencia. / [CA] En els últims anys, les emissions de sutge provinents dels motors de combustió interna han rebut més atenció a causa de l'impacte negatiu que aquestes tenen no sols en l'ambient, sinó també en la salut de l'ésser humà. Com a resposta, lleis cada vegada més estrictes han sigut aplicades impulsant així a la comunitat científica al desenvolupament de motors més eficients en l'ús del combustible i per descomptat més nets en termes d'emissions contaminants. En aquest context, el modelatge computacional ha sigut l'eina utilitzada en nombrosos esforços que busquen contribuir a millorar l'enteniment que es té sobre els altament complexos fenòmens que componen el procés de producció de sutge. El principal objectiu d'aquesta tesi és simular la producció de sutge en rolls dièsel en condicions d'operació típiques d'un motor de combustió interna utilitzant CFD. La consecució de l'objectiu de la tesi comprèn una avaluació preliminar de la configuració dels diferents models per al cas de rolls inerts. En segon lloc, l'estudi detallat de la hipòtesi utilitzada per a caracteritzar l'estructura de la flama a nivell sub-grid (prenent com a base els conceptes well-mixed o flamelet) i de l'enfocament per a tindre en compte la interacció entre turbulència i química. Finalment, es presenten resultats del modelatge de la combustió i producció de sutge per a diferents condicions de contorn de reactivitat i mescla del doll utilitzant un model de sutge de dues equacions. En resum, el lector trobarà al llarg d'aquest document un estudi exhaustiu sobre la combustió i producció de sutge en dolls injectats amb toveres mono-orifici en ambients immòbils. D'aquesta mena de dolls, l'Spray A i Spray D de la Engine Combustion Network són utilitzats com a casos de referència. / [EN] Over the past few years, soot emissions from internal combustion engines have gained attention due to its impact on the environment and human health. In response, ever-stricter legislation has been enforced driving the research community toward more fuel-efficient and cleaner engines. Within this context, soot modeling has been the subject of many efforts seeking to contribute to the understanding of the highly complex phenomena that composes the soot production process. This thesis main objective aims at simulating soot production in Diesel sprays under engine-like conditions using computational fluid dynamics (CFD). The fulfillment of the thesis main objective entails a preliminary assessment of the inert spray computational setup for validation purposes. Then, a detailed study on the sub-grid flame structure and handling of turbulence-chemistry interaction is reported focusing on well-mixed and flamelet assumptions. Lastly, the study of reactivity and mixing boundary condition variations on combustion and soot production are assessed with a two-equation soot model. In summary, throughout this document the reader will find a comprehensive study of combustion and soot modeling in single-hole nozzle sprays in quiescent environments from which the Spray A and Spray D target conditions from the Engine Combustion Network are the main reference cases. / The respondent wishes to acknowledge the financial support received through Programa de Ayudas de Investigación y Desarrollo (PAID-01-16) and Ayudas para movilidad dentro del Programa para la Formación de Personal investigador 2017 of Universitat Politècnica de València and the Government of Spain through the CHEST Project (TRA2017-89139-C2-1-R). The respondent also wants to express his gratitude to Convergent Science for their kind support in the use of CONVERGE software for performing the CFD simulations. Parts of the work presented in this thesis have been supported in a collaborative framework with research partners at Argonne National Laboratory and their support is greatly acknowledged. / Pachano Prieto, LM. (2020). CFD modeling of combustion and soot production in Diesel sprays [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/142189

O objetivo desta tese é o estudo de formação de fuligem em chamas laminares de difusão. Para o modelo de formação de fuligem é escolhido um modelo semi-empírico de duas equações para prever a fração mássica de fuligem e o número de partículas de fuligem. O modelo descreve os processos de nucleação, de crescimento superficial e de oxidação das partículas. Para o modelo de radiação, a perda de calor por radiação térmica (gás e fuligem) é modelada considerando o modelo de gás cinza no limite de chama opticamente fina (OTA - Optically Thin Approximation). São avaliados diferentes modelos de cálculo das propriedades de transporte (detalhado e simplificado). Em relação à cinética química, tanto modelos detalhados quanto reduzidos são utilizados. No presente estudo, é explorada a técnica automática de redução conhecida como Flamelet Generated Manifold (FGM), sendo que esta técnica é capaz de resolver cinética química detalhada com tempos computacionais reduzidos. Para verificar o modelo de formação de fuligem foram realizados uma variedade de experimentos numéricos, desde chamas laminares unidimensionais adiabáticas de etileno em configuração tipo jatos opostos (counterflow) até chamas laminares bidimensionais com perda de calor de etileno em configuração tipo jato (coflow). Para testar a limitação do modelo os acoplamentos de massa e energia entre a fase sólida e a fase gasosa são investigados e quantificados para as chamas contra-corrente Os resultados mostraram que os termos de radiação da fase gasosa e sólida são os termos de maior importancia para as chamas estudas. Os termos de acoplamento adicionais (massa e propriedade termodinâmicas) são geralmente termos de efeitos de segunda ordem, mas a importância destes termos aumenta conforme a quantidade de fuligem aumenta. Como uma recomendação geral o acoplamento com todos os termos deve ser levado em conta somente quando a fração mássica de fuligem, YS, for igual ou superior a 0.008. Na sequência a formação de fuligem foi estudada em chamas bi-dimensionais de etileno em configuração jato laminar usando cinética química detalhada e explorando os efeitos de diferentes modelos de cálculo de propriedades de transporte. Foi encontrado novamente que os termos de radiação da fase gasosa e sólida são os termos de maior importância e uma primeira aproximação para resolver a chama bidimensional de jato laminar de etileno pode ser feita usando o modelo de transporte simplificado. Finalmente, o modelo de fuligem é implementado com a técnica de redução FGM e diferentes formas de armazenar as informações sobre o modelo de fuligem nas tabelas termoquímicas (manifold) são testadas A melhor opção testada neste trabalho é a de resolver todos os flamelets com as fases sólida e gasosa acopladas e armazenar as taxas de reação da fuligem por área de partícula no manifold. Nas simulações bidimensionais estas taxas são então recuperadas para resolver as equações adicionais de formação de fuligem. Os resultados mostraram uma boa concordância qualitativa entre as predições do FGM e da solução detalhada, mas a grande quantidade de fuligem no sistema ainda introduz alguns desafios para a obtenção de bons resultados quantitativos. Entretanto, este trabalho demonstrou o grande potencial do método FGM em predizer a formação de fuligem em chamas multidimensionais de difusão de etileno em tempos computacionais reduzidos. / The objective of this thesis is to study soot formation in laminar diffusion flames. For soot modeling, a semi-empirical two equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth and oxidation. For flame radiation, the radiant heat losses (gas and soot) is modelled by using the grey-gas approximation with Optically Thin Approximation (OTA). Different transport models (detailed or simplified) are evaluated. For the chemical kinetics, detailed and reduced approaches are employed. In the present work, the automatic reduction technique known as Flamelet Generated Manifold (FGM) is being explored. This reduction technique is able to deal with detailed kinetic mechanisms with reduced computational times. To assess the soot formation a variety of numerical experiments were done, from one-dimensional ethylene counterflow adiabatic flames to two-dimensional coflow ethylene flames with heat loss. In order to assess modeling limitations the mass and energy coupling between soot solid particles and gas-phase species are investigated and quantified for counterflow flames. It is found that the gas and soot radiation terms are of primary importance for flame simulations. The additional coupling terms (mass and thermodynamic properties) are generally a second order effect, but their importance increase as the soot amount increases As a general recommendation the full coupling should be taken into account only when the soot mass fraction, YS, is equal to or larger than 0.008. Then the simulation of soot is applied to two-dimensional ethylene co-flow flames with detailed chemical kinetics and explores the effect of different transport models on soot predictions. It is found that the gas and soot radiation terms are also of primary importance for flame simulations and that a first attempt to solve the two-dimensional ethylene co-flow flame can be done using a simplified transport model. Finally an implementation of the soot model with the FGM reduction technique is done and different forms for storing soot information in the manifold is explored. The best option tested in this work is to solve all flamelets with soot and gas-phase species in a coupled manner, and to store the soot rates in terms of specific surface area in the manifold. In the two-dimensional simulations, these soot rates are then retrieved to solve the additional equations for soot modeling. The results showed a good qualitative agreement between FGM solution and the detailed solution, but the high amount of soot in the system still imposes some challenges to obtain good quantitative results. Nevertheless, it was demonstrated the great potential of the method for predicting soot formation in multidimensional ethylene diffusion flames with reduced computational time.

O objetivo desta tese é o estudo de formação de fuligem em chamas laminares de difusão. Para o modelo de formação de fuligem é escolhido um modelo semi-empírico de duas equações para prever a fração mássica de fuligem e o número de partículas de fuligem. O modelo descreve os processos de nucleação, de crescimento superficial e de oxidação das partículas. Para o modelo de radiação, a perda de calor por radiação térmica (gás e fuligem) é modelada considerando o modelo de gás cinza no limite de chama opticamente fina (OTA - Optically Thin Approximation). São avaliados diferentes modelos de cálculo das propriedades de transporte (detalhado e simplificado). Em relação à cinética química, tanto modelos detalhados quanto reduzidos são utilizados. No presente estudo, é explorada a técnica automática de redução conhecida como Flamelet Generated Manifold (FGM), sendo que esta técnica é capaz de resolver cinética química detalhada com tempos computacionais reduzidos. Para verificar o modelo de formação de fuligem foram realizados uma variedade de experimentos numéricos, desde chamas laminares unidimensionais adiabáticas de etileno em configuração tipo jatos opostos (counterflow) até chamas laminares bidimensionais com perda de calor de etileno em configuração tipo jato (coflow). Para testar a limitação do modelo os acoplamentos de massa e energia entre a fase sólida e a fase gasosa são investigados e quantificados para as chamas contra-corrente Os resultados mostraram que os termos de radiação da fase gasosa e sólida são os termos de maior importancia para as chamas estudas. Os termos de acoplamento adicionais (massa e propriedade termodinâmicas) são geralmente termos de efeitos de segunda ordem, mas a importância destes termos aumenta conforme a quantidade de fuligem aumenta. Como uma recomendação geral o acoplamento com todos os termos deve ser levado em conta somente quando a fração mássica de fuligem, YS, for igual ou superior a 0.008. Na sequência a formação de fuligem foi estudada em chamas bi-dimensionais de etileno em configuração jato laminar usando cinética química detalhada e explorando os efeitos de diferentes modelos de cálculo de propriedades de transporte. Foi encontrado novamente que os termos de radiação da fase gasosa e sólida são os termos de maior importância e uma primeira aproximação para resolver a chama bidimensional de jato laminar de etileno pode ser feita usando o modelo de transporte simplificado. Finalmente, o modelo de fuligem é implementado com a técnica de redução FGM e diferentes formas de armazenar as informações sobre o modelo de fuligem nas tabelas termoquímicas (manifold) são testadas A melhor opção testada neste trabalho é a de resolver todos os flamelets com as fases sólida e gasosa acopladas e armazenar as taxas de reação da fuligem por área de partícula no manifold. Nas simulações bidimensionais estas taxas são então recuperadas para resolver as equações adicionais de formação de fuligem. Os resultados mostraram uma boa concordância qualitativa entre as predições do FGM e da solução detalhada, mas a grande quantidade de fuligem no sistema ainda introduz alguns desafios para a obtenção de bons resultados quantitativos. Entretanto, este trabalho demonstrou o grande potencial do método FGM em predizer a formação de fuligem em chamas multidimensionais de difusão de etileno em tempos computacionais reduzidos. / The objective of this thesis is to study soot formation in laminar diffusion flames. For soot modeling, a semi-empirical two equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth and oxidation. For flame radiation, the radiant heat losses (gas and soot) is modelled by using the grey-gas approximation with Optically Thin Approximation (OTA). Different transport models (detailed or simplified) are evaluated. For the chemical kinetics, detailed and reduced approaches are employed. In the present work, the automatic reduction technique known as Flamelet Generated Manifold (FGM) is being explored. This reduction technique is able to deal with detailed kinetic mechanisms with reduced computational times. To assess the soot formation a variety of numerical experiments were done, from one-dimensional ethylene counterflow adiabatic flames to two-dimensional coflow ethylene flames with heat loss. In order to assess modeling limitations the mass and energy coupling between soot solid particles and gas-phase species are investigated and quantified for counterflow flames. It is found that the gas and soot radiation terms are of primary importance for flame simulations. The additional coupling terms (mass and thermodynamic properties) are generally a second order effect, but their importance increase as the soot amount increases As a general recommendation the full coupling should be taken into account only when the soot mass fraction, YS, is equal to or larger than 0.008. Then the simulation of soot is applied to two-dimensional ethylene co-flow flames with detailed chemical kinetics and explores the effect of different transport models on soot predictions. It is found that the gas and soot radiation terms are also of primary importance for flame simulations and that a first attempt to solve the two-dimensional ethylene co-flow flame can be done using a simplified transport model. Finally an implementation of the soot model with the FGM reduction technique is done and different forms for storing soot information in the manifold is explored. The best option tested in this work is to solve all flamelets with soot and gas-phase species in a coupled manner, and to store the soot rates in terms of specific surface area in the manifold. In the two-dimensional simulations, these soot rates are then retrieved to solve the additional equations for soot modeling. The results showed a good qualitative agreement between FGM solution and the detailed solution, but the high amount of soot in the system still imposes some challenges to obtain good quantitative results. Nevertheless, it was demonstrated the great potential of the method for predicting soot formation in multidimensional ethylene diffusion flames with reduced computational time.

O objetivo desta tese é o estudo de formação de fuligem em chamas laminares de difusão. Para o modelo de formação de fuligem é escolhido um modelo semi-empírico de duas equações para prever a fração mássica de fuligem e o número de partículas de fuligem. O modelo descreve os processos de nucleação, de crescimento superficial e de oxidação das partículas. Para o modelo de radiação, a perda de calor por radiação térmica (gás e fuligem) é modelada considerando o modelo de gás cinza no limite de chama opticamente fina (OTA - Optically Thin Approximation). São avaliados diferentes modelos de cálculo das propriedades de transporte (detalhado e simplificado). Em relação à cinética química, tanto modelos detalhados quanto reduzidos são utilizados. No presente estudo, é explorada a técnica automática de redução conhecida como Flamelet Generated Manifold (FGM), sendo que esta técnica é capaz de resolver cinética química detalhada com tempos computacionais reduzidos. Para verificar o modelo de formação de fuligem foram realizados uma variedade de experimentos numéricos, desde chamas laminares unidimensionais adiabáticas de etileno em configuração tipo jatos opostos (counterflow) até chamas laminares bidimensionais com perda de calor de etileno em configuração tipo jato (coflow). Para testar a limitação do modelo os acoplamentos de massa e energia entre a fase sólida e a fase gasosa são investigados e quantificados para as chamas contra-corrente Os resultados mostraram que os termos de radiação da fase gasosa e sólida são os termos de maior importancia para as chamas estudas. Os termos de acoplamento adicionais (massa e propriedade termodinâmicas) são geralmente termos de efeitos de segunda ordem, mas a importância destes termos aumenta conforme a quantidade de fuligem aumenta. Como uma recomendação geral o acoplamento com todos os termos deve ser levado em conta somente quando a fração mássica de fuligem, YS, for igual ou superior a 0.008. Na sequência a formação de fuligem foi estudada em chamas bi-dimensionais de etileno em configuração jato laminar usando cinética química detalhada e explorando os efeitos de diferentes modelos de cálculo de propriedades de transporte. Foi encontrado novamente que os termos de radiação da fase gasosa e sólida são os termos de maior importância e uma primeira aproximação para resolver a chama bidimensional de jato laminar de etileno pode ser feita usando o modelo de transporte simplificado. Finalmente, o modelo de fuligem é implementado com a técnica de redução FGM e diferentes formas de armazenar as informações sobre o modelo de fuligem nas tabelas termoquímicas (manifold) são testadas A melhor opção testada neste trabalho é a de resolver todos os flamelets com as fases sólida e gasosa acopladas e armazenar as taxas de reação da fuligem por área de partícula no manifold. Nas simulações bidimensionais estas taxas são então recuperadas para resolver as equações adicionais de formação de fuligem. Os resultados mostraram uma boa concordância qualitativa entre as predições do FGM e da solução detalhada, mas a grande quantidade de fuligem no sistema ainda introduz alguns desafios para a obtenção de bons resultados quantitativos. Entretanto, este trabalho demonstrou o grande potencial do método FGM em predizer a formação de fuligem em chamas multidimensionais de difusão de etileno em tempos computacionais reduzidos. / The objective of this thesis is to study soot formation in laminar diffusion flames. For soot modeling, a semi-empirical two equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth and oxidation. For flame radiation, the radiant heat losses (gas and soot) is modelled by using the grey-gas approximation with Optically Thin Approximation (OTA). Different transport models (detailed or simplified) are evaluated. For the chemical kinetics, detailed and reduced approaches are employed. In the present work, the automatic reduction technique known as Flamelet Generated Manifold (FGM) is being explored. This reduction technique is able to deal with detailed kinetic mechanisms with reduced computational times. To assess the soot formation a variety of numerical experiments were done, from one-dimensional ethylene counterflow adiabatic flames to two-dimensional coflow ethylene flames with heat loss. In order to assess modeling limitations the mass and energy coupling between soot solid particles and gas-phase species are investigated and quantified for counterflow flames. It is found that the gas and soot radiation terms are of primary importance for flame simulations. The additional coupling terms (mass and thermodynamic properties) are generally a second order effect, but their importance increase as the soot amount increases As a general recommendation the full coupling should be taken into account only when the soot mass fraction, YS, is equal to or larger than 0.008. Then the simulation of soot is applied to two-dimensional ethylene co-flow flames with detailed chemical kinetics and explores the effect of different transport models on soot predictions. It is found that the gas and soot radiation terms are also of primary importance for flame simulations and that a first attempt to solve the two-dimensional ethylene co-flow flame can be done using a simplified transport model. Finally an implementation of the soot model with the FGM reduction technique is done and different forms for storing soot information in the manifold is explored. The best option tested in this work is to solve all flamelets with soot and gas-phase species in a coupled manner, and to store the soot rates in terms of specific surface area in the manifold. In the two-dimensional simulations, these soot rates are then retrieved to solve the additional equations for soot modeling. The results showed a good qualitative agreement between FGM solution and the detailed solution, but the high amount of soot in the system still imposes some challenges to obtain good quantitative results. Nevertheless, it was demonstrated the great potential of the method for predicting soot formation in multidimensional ethylene diffusion flames with reduced computational time.

The probability density function (PDF) approach is a powerful technique for large eddy simulation (LES) based modeling of turbulent reacting flows. In this approach, the joint-PDF of all reacting scalars is estimated by solving a PDF transport equation, thus providing detailed information about small-scale correlations between these quantities. The objective of this work is to further develop the LES/PDF approach for studying flame stabilization in supersonic combustors, and for soot modeling in turbulent flames. Supersonic combustors are characterized by strong shock-turbulence interactions which preclude the application of conventional Lagrangian stochastic methods for solving the PDF transport equation. A viable alternative is provided by quadrature based methods which are deterministic and Eulerian. In this work, it is first demonstrated that the numerical errors associated with LES require special care in the development of PDF solution algorithms. The direct quadrature method of moments (DQMOM) is one quadrature-based approach developed for supersonic combustion modeling. This approach is shown to generate inconsistent evolution of the scalar moments. Further, gradient-based source terms that appear in the DQMOM transport equations are severely underpredicted in LES leading to artificial mixing of fuel and oxidizer. To overcome these numerical issues, a new approach called semi-discrete quadrature method of moments (SeQMOM) is formulated. The performance of the new technique is compared with the DQMOM approach in canonical flow configurations as well as a three-dimensional supersonic cavity stabilized flame configuration. The SeQMOM approach is shown to predict subfilter statistics accurately compared to the DQMOM approach. For soot modeling in turbulent flows, an LES/PDF approach is integrated with detailed models for soot formation and growth. The PDF approach directly evolves the joint statistics of the gas-phase scalars and a set of moments of the soot number density function. This LES/PDF approach is then used to simulate a turbulent natural gas flame. A Lagrangian method formulated in cylindrical coordinates solves the high dimensional PDF transport equation and is coupled to an Eulerian LES solver. The LES/PDF simulations show that soot formation is highly intermittent and is always restricted to the fuel-rich region of the flow. The PDF of soot moments has a wide spread leading to a large subfilter variance. Further, the conditional statistics of soot moments conditioned on mixture fraction and reaction progress variable show strong correlation between the gas phase composition and soot moments. / text

[ES] La formación de hollín en los motores de inyección directa de gasolina se rige por complejas interacciones entre procesos físicos, que varían en función del modo de funcionamiento del motor. Con la aplicación de normativas cada vez más estrictas, junto con ensayos de homologación más rigurosos, el reto de cumplir los límites de emisión de partículas se ha vuelto excepcionalmente exigente. En consecuencia, los fabricantes han recurrido a métodos de post-tratamiento de partículas a la salida del motor, como los filtros de partículas de gasolina, aunque a costa del rendimiento del motor y del consumo de combustible. Además, existe una creciente preocupación por la calidad del aire ambiente debido a las emisiones de partículas de los vehículos, con estudios que indican potenciales riesgos cancerígenos para la salud humana. Las herramientas de modelización ofrecen la ventaja de disminuir la necesidad de costosos ensayos de calibración experimental destinados a encontrar la estrategia óptima para mitigar las emisiones de partículas dentro de unos límites aceptables. En esta Tesis, se desarrolla un marco de modelización exhaustivo que aborda tres fuentes principales de formación de hollín: la mezcla inadecuada que da lugar a zonas ricas, las películas de combustible en las puntas de los inyectores y las películas de combustible en las paredes de la cámara de combustión debido al impacto del chorro. Al incorporar submodelos físicos para abordar estas vías de formación de hollín, se consigue una descripción completa de la distribución del tamaño de las partículas, teniendo en cuenta que las distintas fuentes contribuyen de forma diferente a las emisiones totales de partículas. El número, la masa y el tamaño de las partículas se determinan mediante la utilización de un mecanismo de reacción detallado para resolver las reacciones químicas que se producen en las regiones ricas definidas por el modelo multizona que se integra con un solucionador estocástico de partículas. El enfoque se complementó con simulaciones CFD tridimensional no reactivo para validar las formulaciones de los diversos submodelos. A continuación, se utilizaron ensayos experimentales realizados en un banco de pruebas de motores para evaluar las predicciones del modelo y también sirvieron como herramienta de validación para determinados submodelos en los casos en los que surgían ambigüedades en las simulaciones CFD. Además, se utilizaron ensayos ópticos realizados en un motor distinto del considerado en este estudio para obtener más información sobre los fenómenos en el cilindro que conducen a la formación de hollín. La herramienta ha demostrado su capacidad para predecir con exactitud la distribución del tamaño de las partículas en distintas condiciones de funcionamiento, captando eficazmente los cambios en la configuración de los parámetros del motor y las condiciones termodinámicas. Se identificó que las principales fuentes que contribuían a las emisiones de partículas procedían de la disminución de la calidad de la mezcla, especialmente evidente en cargas más elevadas debido al enriquecimiento del combustible, y de la humectación de la punta del inyector, que se manifestaba en todas las condiciones. Por consiguiente, el modelo de mezcla, basado en los parámetros de mezcla pertinentes, resultó eficaz para generar distribuciones de dosados, mientras que el modelo de evaporación de la película del inyector calculó la masa de la película de combustible en la punta de forma satisfactoria. Al integrar estos aspectos con un enfoque de postprocesamiento exhaustivo, que incorpora la eficiencia del recuento de partículas y consideraciones sobre el volumen de escape, se logró un alto grado de concordancia entre las predicciones numéricas y los datos de partículas medidos. Como resultado, la herramienta puede aprovecharse para simular las emisiones de partículas en condiciones transitorias de conducción, facilitando así el desarrollo de motores en el futuro próximo. / [CA] La formació de sutge en els motors d'injecció directa de gasolina es regeix per complexes interaccions entre processos físics, que varien en funció de la manera de funcionament del motor. Amb l'aplicació de reglaments cada vegada més estrictes i normes d'assaig més rigoroses, el repte de complir els límits d'emissió de partícules s'ha tornat excepcionalment exigent. En conseqüència, els fabricants han recorregut a mètodes de tractament de partícules posteriors a l'eixida del motor, com els filtres de partícules de gasolina, encara que a costa del rendiment del motor i del consum de combustible. A més, existeix una creixent preocupació per la qualitat de l'aire ambient degut a les emissions de partícules dels vehicles, amb estudis que indiquen riscos cancerígens potencials per a la salut humana. Les eines de modelització ofereixen l'avantatge de disminuir la necessitat de costosos assajos de calibratge experimental destinats a trobar l'estratègia òptima per a mitigar les emissions de partícules dins d'uns límits acceptables. En aquesta Tesi, es desenvolupa un marc de modelització exhaustiu que aborda tres fonts principals de formació de sutge: la mescla inadequada que dona lloc a zones riques, les pel·lícules de combustible en les puntes dels injectors i les pel·lícules de combustible en les parets de la cambra de combustió a causa de l'impacte del doll. En incorporar submodels físics per a abordar aquestes vies de formació de sutge, s'aconsegueix una descripció completa de la distribució de la grandària de les partícules, tenint en compte que les diferents fonts contribueixen de manera diferent a les emissions totals de partícules. El número, la massa i la grandària de les partícules es determinen mitjançant la utilització d'un mecanisme de reacció detallat per a resoldre les reaccions químiques que es produeixen a les regions riques definides pel model multizona que s'integra amb un solucionador estocàstic de partícules. L'enfocament es va complementar amb simulacions CFD tridimensional no reactiu per a validar les formulacions dels diversos submodels associats a les diferents vies. A continuació, es van utilitzar assajos experimentals realitzats en un banc de proves de motors per a avaluar les prediccions del model i també van servir com a eina de validació per a determinats submodels en els casos en els quals sorgien ambigüitats en les simulacions CFD. A més, es van utilitzar assajos òptics realitzats en un motor diferent del considerat en aquest estudi per a obtenir més informació sobre els fenòmens en el cilindre que condueixen a la formació de sutge. L'eina de modelització ha demostrat la seua capacitat per a predir amb exactitud la distribució de la grandària de les partícules en diferents condicions de funcionament, captant eficaçment els canvis en la configuració dels paràmetres del motor i les condicions termodinàmiques. Es va identificar que les principals fonts que contribuïen a les emissions de partícules procedien de la disminució de la qualitat de la mescla, especialment evident en càrregues més elevades a causa de l'enriquiment del combustible, i de la humectació de la punta de l'injector, que es manifestava en totes les condicions. Per consegüent, el model de mescla, basat en els paràmetres de mescla pertinents, va resultar eficaç per a generar distribucions de dosatges, mentre que el model d'evaporació de la pel·lícula de l'injector va calcular la massa de la pel·lícula de combustible en la punta de manera satisfactòria. En integrar aquests aspectes amb un enfocament de postprocessament exhaustiu, que incorpora l'eficiència del recompte de partícules i consideracions sobre el volum de fuita, es va aconseguir un alt grau de concordança entre les prediccions numèriques i les dades de partícules mesurades. Com a resultat, l'eina pot aprofitar-se per a simular les emissions de partícules en condicions transitòries de conducció, facilitant així el desenvolupament de motors en el futur pròxim. / [EN] Soot formation in gasoline direct injection engines is governed by complex interactions among physical processes, which vary depending on the engine operating mode. With the implementation of increasingly stringent regulations and more rigorous testing standards, the challenge of meeting particle emission limits has become exceptionally demanding. Consequently, manufacturers have resorted to post-engine outlet particle treatment methods, such as gasoline particulate filters, albeit at the expense of engine performance and fuel consumption. Moreover, there is a growing concern regarding environmental air quality due to vehicle particle emissions, with studies indicating potential cancerous risks to human health. Modeling tools offer the benefit of diminishing the need for expensive experimental calibration campaigns aimed at finding the optimal strategy for mitigating particle emissions within acceptable limits. In this Thesis, a comprehensive modeling framework is developed that addresses three primary sources of soot formation: inadequate mixing resulting in rich pockets, fuel films on injector tips, and fuel films on combustion chamber walls due to spray impingement. By incorporating physical sub-models to address these pathways of soot formation, a comprehensive depiction of the particle size distribution becomes achievable, considering that the various sources contribute differently to overall particle emissions. Particle number, mass, and size are determined through the utilization of a detailed reaction mechanism to solve the chemical reactions occurring in rich regions defined by the multi-zone model that is integrated with a stochastic particle solver. The approach was complemented by non-reactive 3D CFD simulations to validate the formulations of the various sub-models associated with different pathways. Experimental measurements conducted on an engine test bench were then utilized to evaluate the model predictions and also served as a validation tool for certain sub-models in cases where ambiguity arose in the 3D CFD simulations. Furthermore, optical engine experiments conducted on a different engine from the one considered in this study were employed to gain further insights into the in-cylinder phenomena driving soot formation. The modeling tool has demonstrated its capability to accurately predict the particle size distribution across various operating conditions, effectively capturing changes in engine parameter settings and thermodynamic conditions. The primary sources contributing to particle emissions stemmed from diminished mixture quality, particularly evident at higher loads due to fuel enrichment, and injector tip wetting, which manifested across all conditions. Therefore, the mixing model, based on pertinent mixing parameters proved effective in generating equivalence ratio distributions, while the injector film evaporation model computed fuel film mass on the injector tip in a commendable way. By integrating these aspects with a thorough post-processing approach, incorporating particle counting efficiency and exhaust volume considerations, a high degree of agreement was achieved between numerical predictions and measured particle data. As a result, the tool can be leveraged to simulate particle emissions under transient driving conditions, thereby facilitating engine development in the foreseeable future. / Abboud, R. (2024). A Comprehensive Modeling Toolchain for Particle Emissions in GDI Engines [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/207344