As complexidades das configurações, da geometria e da topologia dos escoamentos encontrados em sistemas tais como câmaras de combustão de turbo-máquinas requerem métodos numéricos eficientes para modelar os escoamentos turbulentos reativos existentes. Este trabalho de modelagem numérica está relacionado com o desenvolvimento de ferramentas numéricas eficientes necessárias ao estudo dos escoamentos turbulentos reativos nessas aplicações. Na primeira parte desta tese, após relembrar os elementos de base que permitem a descrição de um escoamento turbulento reativo prée-misturado, são apresentadas as diversas etapas que presidem a elaboração do programa computacional. Inicialmente, introduz-se a técnica de compressibilidade artificial, que consiste em modificar as equações governantes acrescentando termos não-estacionários no pseudo-tempo, destinados a introduzir uma velocidade do som finita no domínio de cálculo. Em seguida é exposta a discretização por volumes finitos das equações governantes em malhas não-estruturadas bidimensionais híbridas, que contêm elementos triangulares e/ou quadriláteros. A adição de termos de dissipação artificial bem com o uso de um método do tipo Runge-Kutta para integrar os termos temporais no pseudo-tempo das equações governantes, que completam o método, são comentados em seguir. O programa de cálculo, desenvolvido na linguagem C++, foi validado considerando-se o problema de uma chama plana turbulenta, propagando-se livremente em uma turbulência "congelada". Este fenômeno foi modelado com um termo de produção química médio do tipo flamelet, no qual um corte representando extinção foi introduzido com a finalidade de selecionar uma só velocidade possível de propagação. Os resultados obtidos para uma chama lenta (0,5 m/s) ou rápida (10 m/s) foram comparados com aqueles obtidos na literatura. A precisão dos resultados é bastante satisfatória e os gradientes elevados de pressão estática, característicos do termo fonte empregado, são capturados corretamente, tanto no que se refere a sua intensidade quanto à sua localização, dentro da frente de chama média. Uma análise de sensibilidade e precisão da malha foi igualmente conduzida, assim como um estudo das características de convergência em função do valor do coeficiente de compressibilidade artificial. Tendo sido validado, o programa de cálculo foi, em seguida, utilizado com a finalidade de estudar esteiras inertes ou reativas encontradas à jusante de um obstáculo de seção transversal triangular, colocado no meio de um escoamento turbulento em um canal. Um modelo de turbulência do tipo k-e, acoplado a dois tipos de leis de parede foi utilizado. O desprendimento turbilhonar, típico dos escoamentos de esteiras inertes considerados, foi previsto corretamente com, em particular, uma freqüência de Strouhal de desprendimento com boa concordância com os resultados experimentais disponíveis. As esteiras reativas foram calculadas com um modelo de flamelet, onde os transportes turbulentos da variável de avanço da combustão foram calculados por aproximação do tipo gradiente. Três formas do termo de produção química médio, existentes na literatura, foram utilizados. Aparentemente, a estrutura da esteira reativa prevista é fortemente dependente da formulação adotada para o termo fonte médio. Embora os resultados obtidos sejam comparados àqueles apresentados por outro autores, em uma configuração de escoamento similar, a comparação com resultados experimentais mostra que é necessário uma melhoria na modelagem física de forma a prever este tipo de escoamento reativo de maneira satisfatória.

Ce travail de thèse est consacré à l'étude numérique et théorique de la modulation de la turbulence par des particules. Cette étude s'appuie sur des résultats issus de simulations de type Euler/Lagrange qui résolvent directement les équations instantanées de la phase gazeuse et effectuent un suivi de trajectoires des particules. La configuration étudiée représente une nappe de particules injectées à haute vitesse dans une turbulence homogène isotrope décroissante. Le mouvement des particules est supposé uniquement gouverné par la force de traînée visqueuse. Le chargement en particules est suffisamment important pour que les particules influent sur la phase gazeuse (couplage inverse) mais suffisamment faible pour pouvoir négliger les collisions interparticulaires. Une analyse des équations de transport des principales grandeurs moyennes de l'écoulement est menée pour déterminer les effets directs et indirects des particules sur la turbulence fluide. L'étude des transferts d'énergie entre phases montre que la présence des particules tend à détruire la turbulence gazeuse au centre de la nappe et à l'augmenter à la périphérie. Ce dernier effet est causé par la forte corrélation entre la distribution de particules et la vitesse instantanée du gaz. Le modèle k-epsilon est ensuite étudié et la validité de ses hypothèses de fermeture en écoulement diphasique est éprouvée à l'aide de tests a priori. une nouvelle formulation de type viscosité turbulente, fonction des paramètres diphasiques, est utilisée pour modéliser le tenseur de Reynolds du gaz. Une équation de Langevin diphasique est également testée pour modéliser les équations de vitesse de dérive et de covariance des fluctuations de vitesse fluide-particules.

For turbulent flow calculations, some of the well-known turbulence models in the literature are applied on a previously developed Navier-Stokes solver designed to handle laminar flows. A finite volume formulation, which is cell-based for inviscid terms and cell-vertex for viscous terms, is used for numerical discretization of the Navier-Stokes equations in conservative form. This formulation is combined with one-step, explicit time marching Lax-Wendroff numerical scheme that is second order accurate in space. To minimize non-physical oscillations resulting from the numerical scheme, second and fourth order artificial smoothing terms are added. To increase the convergence rate of the solver, local time stepping technique is applied. Before applying turbulence models, Navier-Stokes solver is tested for a case of subsonic, laminar flow over a flat plate. The results are in close agreement with Blasius similarity solutions. To calculate turbulent flows, Boussinesq eddy-viscosity approach is utilized. The eddy viscosity (also called turbulent viscosity), which arises as a consequence of this approach, is calculated using Cebeci-Smith, Michel et. al., Baldwin-Lomax, Chien&rsquo / s k-epsilon and Wilcox&rsquo / s k-omega turbulence models. To evaluate the performances of these turbulence models and to compare them with each other, the solver has been tested for a case of subsonic, laminar - transition fixed - turbulent flow over a flat plate. The results are verified by analytical solutions and empirical correlations.

El cultivo protegido bajo invernadero tiene altos costes energéticos derivados de la ventilación forzada. En cambio, la ventilación natural puede ser una solución barata que reduzca el consumo energético. No obstante, un diseño que permita la ventilación natural del invernadero supone un reto debido a la complejidad del fenómeno físico. Frente a otros métodos de diseño, como los métodos experimentales o analíticos, en los últimos tiempos la dinámica de fluidos computacional (CFD por sus siglas en inglés) se ha convertido en la herramienta más utilizada para estudiar este tipo de fenómenos, gracias a su relativo bajo coste y a la rapidez en la obtención de resultados. No obstante, los modelos CFD deben refrendarse mediante validaciones realizadas a través de datos experimentales. Un análisis bibliográfico detallado del uso de la CFD aplicada a invernaderos muestra, en general, que los trabajos carecen de suficientes datos experimentales, seguramente debido al alto coste de los sensores para adquirirlos y la dificultad que conlleva el trabajo en condiciones de campo. Además, se observa que las simulaciones CFD en invernaderos no se realizan con un procedimiento sistematizado. Por ello, la presente tesis, por un lado describe un dispositivo y unos métodos sencillos y baratos para obtener datos atmosféricos, y por otro, propone una visión crítica sobre la investigación realizada hasta el momento, con el fin de sistematizar la manera de generar modelos CFD aplicados a la ventilación natural de invernaderos. Finalmente, la tesis se complementa con un ejemplo sobre un caso práctico. Para ello, en primer lugar, se realizó una revisión bibliográfica de las diferentes guías de buenas prácticas en diferentes campos de la tecnología principalmente en edificación, para sistematizar y adaptar las recomendaciones para generar modelos CFD en invernaderos. En segundo lugar, se desarrolló un sistema de adquisición de datos sencillo, consistente en una red de sensores, que permite medir simultáneamente la velocidad y dirección de viento en 20 puntos. Este sistema de sensores fue calibrado y probado en campo satisfactoriamente obteniendo una precisión similar a los anemómetros comerciales con un precio 30 veces superior. En tercer lugar, se generaron 24 modelos del flujo de aire alrededor del invernadero, resultantes de la combinación de cuatro modelos de turbulencia (k- ¿; RNG k-¿; k-¿ y RSM); dos esquemas de discretización (primer orden y segundo orden) y tres velocidades de viento exterior (3; 3,5 y 4 m/s) con el fin de analizar sus diferencias y demostrar sus ventajas e inconvenientes. Por este motivo, en cuarto lugar, se comprobó su capacidad de ajustarse a los datos de campo, validando los modelos con un análisis de regresión lineal sobre los datos experimentales. Con este estudió se reveló que los modelos SST k-¿ y el RSM (segundo orden) son los que mejor representan el flujo de ventilación y se demostró que el modelo k-¿ estándar (primer orden), el más utilizado en la bibliografía, no sólo ofrece resultados diferentes al resto de modelos, sino que su rendimiento es pobre para predecir el flujo de ventilación. / Granell Ruiz, R. (2014). Análisis del flujo ambiental y propuesta metodológica para simulaciones CFD aplicadas a la ventilación natural de invernaderos [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/37194 / TESIS

Purpose of this thesis is to check influence of turbulent model used for simulation of flow close to primary elementi inserted into piping. The goal is to check if results computed by these models are equal and how precise is their prediction.

In this thesis, we have evaluated a stochastic Lagrangian model for computing particle deposition rates with prospects to use for atmospheric deposition applications.  The model is one-dimensional and models the particle dynamics in the boundary layers near walls and obstacles by simulating the coherent turbulent structures and Brownian motion governing the wall-normal transport. The deposition model is used with a hybrid deterministic/stochastic particle dispersion model governing the dynamics in the turbulent bulk flow. We used a steady-state RANS k-ϵ turbulence model to simulate the turbulent fluid flow in a neutral atmospheric boundary layer (ABL) using the with inflow boundary conditions by Richards & Hoxey (1993). The turbulence model is solved with the SIMPLE algorithm using the OpenFOAM software. The mean-field characteristic of the turbulent flow in the computational domain is exported and used for the particle model. The particle model is a Lagrangian Langevin-type model, consisting of a system of stochastic differential equations. The particle model was solved using a weakly first order a-stable scheme. We evaluated the deposition model by computing the deposition rate for a range of particle sizes and compared our results with collected experimental wind tunnel data. The numerical experiment was done in a computational domain based on the ABL model by Hargreaves & Wright (2007), a rectangular domain with a logarithmic wind profile. We used a particle source near the inflow boundary with an instantaneously release at the initial time. Results showed disagreement with the experimental data and was only valid for medium sized particles. However, time restrictions led to the analysis being cut short and only a single simulation was conducted. A definite conclusion on the suitability of the method could not be made based solely on this single results. Some uncertainties were identified and discussed for further potential work on the evaluation of the method. However, one conclusion was drawn on the performance of the method. The computational cost was concluded to be too high with the first order particle scheme used and higher order schemes is required for any practical use of the method for atmospheric deposition applications.

<p>This is a thesis work proposed by Sweco System in order to carry out a study related to the heating system of a circular fuel oil storage tank or cistern. The study tank is a 23m diameter and 18m height with a storage capacity of around 7500m3 of Eo5 heavy fuel oil. The content ought to be at a minimum storage temperature of 50ºC so that the fuel oil is fluid enough and operation labors can be adequately performed. In fact, these types of heavy fuel oils have fairly high viscosities at lower temperatures and the heating and pumping system can be compromised at temperatures below the pour point. For this purpose a heating system is installed to maintain the fluid warm. So far the system was operated by an oil burner but there are plans to its replacement by a District Heating-heat exchanger combo. Thereby, tank heating needs, flow and thermal patterns and heat transfer within it are principally studied.</p><p> </p><p>Tank boundaries are studied and their thermal resistances are calculated in order to dimension heat supply capacity. The study implies Finite Elements (Comsol Multiphysics) and Finite Volume (Fluent) analysis to work out some stationary heat transfer by conduction cases on some parts and thermal bridges present on these boundaries. Afterwards both cooling and heating processes of the fuel oil are studied using several strategies: basic models and Computational Fluid Dynamics (CFD). CFD work with Fluent is focused on optimizing inlet and outlet topologies. Understanding the cooling process is sought as well; Fluent CFD transient models are simulated in this way as well. Additionally the effect of filling levels is taken into account leading to a multiphase (fuel oil and air) flow cases where especially heating coupling of both phases is analyzed.</p><p> </p><p>Results show that maximum heat supply needs are around 80kW when the tank temperature is around 60ºC and 70kW when it is around 50ºC. Expectedly the main characteristic of the flow turns out to be the buoyancy driven convective pattern. K-ε turbulence viscous models are applied to both heating and cooling processes showing thermal stratification, especially at the bottom of the tank. Hotter fluid above follows very complex flow patterns. During the heating processes models used predict fairly well mixed and homogenous temperature distribution regardless small stratification at the bottom of the tank. In this way no concrete inlet-outlet configuration shows clear advantages over the rest. Due to the insulation of the tank, low thermal conductivity of the fluid and vast amount of mass present in the tank, the cooling process is slow (fluid average temperature drops around 5.7 ºC from 60ºC in 15 days when the tank is full and ambient temperature is considered to be at -20ºC) and lies somewhere in the middle between the solid rigid and perfect mixture cooling processes. However, due to stratification some parts of the fluid reach minimum admissible temperatures much faster than average temperature does. On the other hand, as expected, air phase acts as an additional thermal resistance; anyhow the cooling process is still faster for lower filling levels than the full one.</p> / <p>El presente proyecto fue propuesto por Sweco Systems para llevar a cabo un estudio relacionado con el sistema de calefacción de una cisterna o tanque de almacenamiento de fuel oil circular. Dicho tanque tiene 23 m de diámetro  y 18 m de altura con una capacidad de almacenamiento de alrededor de 7500 m³ de Eo5 fuel oil pesado. El contenido mantenerse a una temperatura mínima de 50 ºC de manera que el fuel oil es suficientemente fluido para que las labores de operación puedan ser ejecutadas adecuadamente. De hecho, estos tipos de fuel oil pesado tienen altas viscosidades a bajas temperaturas y, por tanto, tanto los sistemas de calefacción y como el de bombeo pueden verse comprometidosr a temperaturas por debajo del pour point. Con este fin un sistema de calefacción es instalado para mantener el fluido suficientemente caliente. Hasta el momento, el sistema era operado por un quemador de fuel, sin embargo, hay planes que éste sea sustituido por un combo intercambiador de calor-District Heating. Por lo tanto, principalmente son estudiadas las necesidades de calefacción así como los flujos térmicos y fluidos.</p><p>Se estudian las fronteras del tanque, y sus respectivas resistencias térmicas son calculadas con el fin de dimensionar la capacidad necesaria de suministro de calor. El estudio implica Elementos Finitos (Comsol Multiphysics) y Volúmenes Finitos (Fluent) para elaborar análisis estacionarios de transferencia de calor por conducción en algunos casos. Existen puentes térmicos en las paredes y su importancia es también anallizada. Posteriormente se estudian tanto los procesos de calentamiento y enfriamiento del fuel oil utilizando diversas estrategias: modelos básicos y Dinámica de Fluidos Computacional (CFD). El trabajo con CFD se centra en la optimización de topologías de entradas y salidas del sistema. También es solicitado entender el proceso de enfriamiento; En este sentido, se simulan modelos CFD transitorios de Fluent. Además, el efecto de los niveles de llenado se tiene en cuenta dando lugar a estudios de flujo multifase (fuel oil y aire), haciendo hincapié en el análisis de acoplamiento de transferencia de calor entre las dos fases.</p><p>Los resultados muestran que las necesidades de calefacción máximas son de alrededor de 80kW cuando la temperatura del tanque es de alrededor de 60 º C y 70kW cuando está alrededor de 50 ºC. Como era de esperar, la principal característica de este tipo de flujos es la convección natural resultante de las fuerzas de flotabilidad. Se aplican modelos turbulentos k-ε a los procesos de calentamiento y enfriamiento, mostrando estratificación térmica, sobre todo en la parte inferior de la cisterna. El líquido más caliente que se sitúa encima muestra complejos patrones de flujo. Durante los procesos de calentamiento, los modelos utilizados predicen un buen mezclado y distribución homogénea de la temperatura independientemente de esta pequeña estratificación en la parte inferior de la cisterna. De esta manera, ninguna concreta configuración de entradas-salidas simuladas muestra claras ventajas sobre el resto. Debido al aislamiento de la cisterna, la baja conductividad térmica del fluido y la gran cantidad de masa presente en el tanque el proceso de enfriamiento es lento (la temperatura media del fluido desciende 5.7 º C desde 60 º C en 15 días cuando el tanque está lleno y la temperatura ambiente es de -20 º C) y se encuentra en algún lugar en medio de los procesos de enfriamiento del sólido rígido y perfecta mezcla. Sin embargo, debido a la estratificación, algunas partes el líquido alcanzan la temperatura mínima admisible mucho más rápido que la media de temperatura. Por otra parte, como se esperaba, la fase de aire actúa como una resistencia térmica adicional, de todos modos, el proceso de enfriamiento es aún más rápido para niveles de llenado más bajos que el lleno.</p>

This is a thesis work proposed by Sweco System in order to carry out a study related to the heating system of a circular fuel oil storage tank or cistern. The study tank is a 23m diameter and 18m height with a storage capacity of around 7500m3 of Eo5 heavy fuel oil. The content ought to be at a minimum storage temperature of 50ºC so that the fuel oil is fluid enough and operation labors can be adequately performed. In fact, these types of heavy fuel oils have fairly high viscosities at lower temperatures and the heating and pumping system can be compromised at temperatures below the pour point. For this purpose a heating system is installed to maintain the fluid warm. So far the system was operated by an oil burner but there are plans to its replacement by a District Heating-heat exchanger combo. Thereby, tank heating needs, flow and thermal patterns and heat transfer within it are principally studied.   Tank boundaries are studied and their thermal resistances are calculated in order to dimension heat supply capacity. The study implies Finite Elements (Comsol Multiphysics) and Finite Volume (Fluent) analysis to work out some stationary heat transfer by conduction cases on some parts and thermal bridges present on these boundaries. Afterwards both cooling and heating processes of the fuel oil are studied using several strategies: basic models and Computational Fluid Dynamics (CFD). CFD work with Fluent is focused on optimizing inlet and outlet topologies. Understanding the cooling process is sought as well; Fluent CFD transient models are simulated in this way as well. Additionally the effect of filling levels is taken into account leading to a multiphase (fuel oil and air) flow cases where especially heating coupling of both phases is analyzed.   Results show that maximum heat supply needs are around 80kW when the tank temperature is around 60ºC and 70kW when it is around 50ºC. Expectedly the main characteristic of the flow turns out to be the buoyancy driven convective pattern. K-ε turbulence viscous models are applied to both heating and cooling processes showing thermal stratification, especially at the bottom of the tank. Hotter fluid above follows very complex flow patterns. During the heating processes models used predict fairly well mixed and homogenous temperature distribution regardless small stratification at the bottom of the tank. In this way no concrete inlet-outlet configuration shows clear advantages over the rest. Due to the insulation of the tank, low thermal conductivity of the fluid and vast amount of mass present in the tank, the cooling process is slow (fluid average temperature drops around 5.7 ºC from 60ºC in 15 days when the tank is full and ambient temperature is considered to be at -20ºC) and lies somewhere in the middle between the solid rigid and perfect mixture cooling processes. However, due to stratification some parts of the fluid reach minimum admissible temperatures much faster than average temperature does. On the other hand, as expected, air phase acts as an additional thermal resistance; anyhow the cooling process is still faster for lower filling levels than the full one. / El presente proyecto fue propuesto por Sweco Systems para llevar a cabo un estudio relacionado con el sistema de calefacción de una cisterna o tanque de almacenamiento de fuel oil circular. Dicho tanque tiene 23 m de diámetro  y 18 m de altura con una capacidad de almacenamiento de alrededor de 7500 m3 de Eo5 fuel oil pesado. El contenido mantenerse a una temperatura mínima de 50 ºC de manera que el fuel oil es suficientemente fluido para que las labores de operación puedan ser ejecutadas adecuadamente. De hecho, estos tipos de fuel oil pesado tienen altas viscosidades a bajas temperaturas y, por tanto, tanto los sistemas de calefacción y como el de bombeo pueden verse comprometidosr a temperaturas por debajo del pour point. Con este fin un sistema de calefacción es instalado para mantener el fluido suficientemente caliente. Hasta el momento, el sistema era operado por un quemador de fuel, sin embargo, hay planes que éste sea sustituido por un combo intercambiador de calor-District Heating. Por lo tanto, principalmente son estudiadas las necesidades de calefacción así como los flujos térmicos y fluidos. Se estudian las fronteras del tanque, y sus respectivas resistencias térmicas son calculadas con el fin de dimensionar la capacidad necesaria de suministro de calor. El estudio implica Elementos Finitos (Comsol Multiphysics) y Volúmenes Finitos (Fluent) para elaborar análisis estacionarios de transferencia de calor por conducción en algunos casos. Existen puentes térmicos en las paredes y su importancia es también anallizada. Posteriormente se estudian tanto los procesos de calentamiento y enfriamiento del fuel oil utilizando diversas estrategias: modelos básicos y Dinámica de Fluidos Computacional (CFD). El trabajo con CFD se centra en la optimización de topologías de entradas y salidas del sistema. También es solicitado entender el proceso de enfriamiento; En este sentido, se simulan modelos CFD transitorios de Fluent. Además, el efecto de los niveles de llenado se tiene en cuenta dando lugar a estudios de flujo multifase (fuel oil y aire), haciendo hincapié en el análisis de acoplamiento de transferencia de calor entre las dos fases. Los resultados muestran que las necesidades de calefacción máximas son de alrededor de 80kW cuando la temperatura del tanque es de alrededor de 60 º C y 70kW cuando está alrededor de 50 ºC. Como era de esperar, la principal característica de este tipo de flujos es la convección natural resultante de las fuerzas de flotabilidad. Se aplican modelos turbulentos k-ε a los procesos de calentamiento y enfriamiento, mostrando estratificación térmica, sobre todo en la parte inferior de la cisterna. El líquido más caliente que se sitúa encima muestra complejos patrones de flujo. Durante los procesos de calentamiento, los modelos utilizados predicen un buen mezclado y distribución homogénea de la temperatura independientemente de esta pequeña estratificación en la parte inferior de la cisterna. De esta manera, ninguna concreta configuración de entradas-salidas simuladas muestra claras ventajas sobre el resto. Debido al aislamiento de la cisterna, la baja conductividad térmica del fluido y la gran cantidad de masa presente en el tanque el proceso de enfriamiento es lento (la temperatura media del fluido desciende 5.7 º C desde 60 º C en 15 días cuando el tanque está lleno y la temperatura ambiente es de -20 º C) y se encuentra en algún lugar en medio de los procesos de enfriamiento del sólido rígido y perfecta mezcla. Sin embargo, debido a la estratificación, algunas partes el líquido alcanzan la temperatura mínima admisible mucho más rápido que la media de temperatura. Por otra parte, como se esperaba, la fase de aire actúa como una resistencia térmica adicional, de todos modos, el proceso de enfriamiento es aún más rápido para niveles de llenado más bajos que el lleno.

Spark ignited engines are still important for conventional as well as for hybrid power trains and are thus objective to optimization. Today a lot of functionalities arise from software solutions, which have to be calibrated. Modern engine technologies provide an extensive variability considering their valve train, fuel injection and load control. Thus, calibration efforts are really high and shall be reduced by introduction of virtual methods. In this work a physical 0D combustion model is set up, which can cope with a new generation of spark ignition engines. Therefore, at first cylinder thermodynamics are modeled and validated in the whole engine map with the help of a real-time capable approach. Afterwards an up to date turbulence model is introduced, which is based on a quasi-dimensional k-epsilon-approach and can cope with turbulence production from large scale shearing. A simplified model for ignition delay is implemented which emphasizes the transfer from laminar to turbulent flame propagation after ignition. The modeling is completed with the calculation of overall heat release rates in a 0D entrainment approach with the help of turbulent flame velocities. After validation of all sub-models, the 0D combustion prediction is used in combination with a 1D gas exchange analysis to virtually calibrate the modern engine torque structure and the ECU function for exhaust gas temperature with extensive simulations. / Moderne Ottomotoren spielen heute sowohl in konventionellen als auch hybriden Fahrzeugantrieben eine große Rolle. Aktuelle Konzepte sind hochvariabel bezüglich Ventilsteuerung, Kraftstoffeinspritzung und Laststeuerung und ihre Optimierungspotentiale erwachsen zumeist aus neuen Softwarefunktionen. Deren Applikation ist zeit- und kostenintensiv und soll durch virtuelle Methoden unterstützt werden. In der vorliegenden Arbeit wird ein physikalisches 0D Verbrennungsmodell für Ottomotoren aufgebaut und bis zur praktischen Anwendung geführt. Dafür wurde zuerst die Thermodynamik echtzeitfähig modelliert und im gesamten Motorenkennfeld abgeglichen. Der Aufbau eines neuen Turbulenzmodells auf Basis der quasidimensionalen k-epsilon-Gleichung ermöglicht anschließend, die veränderlichen Einflüsse globaler Ladungsbewegung auf die Turbulenz abzubilden. Für den Brennverzug wurde ein vereinfachtes Modell abgeleitet, welches den Übergang von laminarer zu turbulenter Flammenausbreitung nach der Zündung in den Vordergrund stellt. Der restliche Brennverlauf wird durch die physikalische Ermittlung der turbulenten Brenngeschwindigkeit in einem 0D Entrainment-Ansatz dargestellt. Nach Validierung aller Teilmodelle erfolgt die virtuelle Bedatung der Momentenstruktur und der Abgastemperaturfunktion für das Motorsteuergerät.

Spark ignited engines are still important for conventional as well as for hybrid power trains and are thus objective to optimization. Today a lot of functionalities arise from software solutions, which have to be calibrated. Modern engine technologies provide an extensive variability considering their valve train, fuel injection and load control. Thus, calibration efforts are really high and shall be reduced by introduction of virtual methods. In this work a physical 0D combustion model is set up, which can cope with a new generation of spark ignition engines. Therefore, at first cylinder thermodynamics are modeled and validated in the whole engine map with the help of a real-time capable approach. Afterwards an up to date turbulence model is introduced, which is based on a quasi-dimensional k-epsilon-approach and can cope with turbulence production from large scale shearing. A simplified model for ignition delay is implemented which emphasizes the transfer from laminar to turbulent flame propagation after ignition. The modeling is completed with the calculation of overall heat release rates in a 0D entrainment approach with the help of turbulent flame velocities. After validation of all sub-models, the 0D combustion prediction is used in combination with a 1D gas exchange analysis to virtually calibrate the modern engine torque structure and the ECU function for exhaust gas temperature with extensive simulations.:Contents 1 Introduction. 2 Thermodynamic modeling with real-time capability. 3 Quasi-dimensional modeling of turbulence and global charge motion. 4 Physical modeling of ignition delay. 5 Combustion modeling based on a 0D entrainment approach. 6 Virtual engine calibration with a quasi-dimensional combustion model. 7 Summary and outlook. / Moderne Ottomotoren spielen heute sowohl in konventionellen als auch hybriden Fahrzeugantrieben eine große Rolle. Aktuelle Konzepte sind hochvariabel bezüglich Ventilsteuerung, Kraftstoffeinspritzung und Laststeuerung und ihre Optimierungspotentiale erwachsen zumeist aus neuen Softwarefunktionen. Deren Applikation ist zeit- und kostenintensiv und soll durch virtuelle Methoden unterstützt werden. In der vorliegenden Arbeit wird ein physikalisches 0D Verbrennungsmodell für Ottomotoren aufgebaut und bis zur praktischen Anwendung geführt. Dafür wurde zuerst die Thermodynamik echtzeitfähig modelliert und im gesamten Motorenkennfeld abgeglichen. Der Aufbau eines neuen Turbulenzmodells auf Basis der quasidimensionalen k-epsilon-Gleichung ermöglicht anschließend, die veränderlichen Einflüsse globaler Ladungsbewegung auf die Turbulenz abzubilden. Für den Brennverzug wurde ein vereinfachtes Modell abgeleitet, welches den Übergang von laminarer zu turbulenter Flammenausbreitung nach der Zündung in den Vordergrund stellt. Der restliche Brennverlauf wird durch die physikalische Ermittlung der turbulenten Brenngeschwindigkeit in einem 0D Entrainment-Ansatz dargestellt. Nach Validierung aller Teilmodelle erfolgt die virtuelle Bedatung der Momentenstruktur und der Abgastemperaturfunktion für das Motorsteuergerät.:Contents 1 Introduction. 2 Thermodynamic modeling with real-time capability. 3 Quasi-dimensional modeling of turbulence and global charge motion. 4 Physical modeling of ignition delay. 5 Combustion modeling based on a 0D entrainment approach. 6 Virtual engine calibration with a quasi-dimensional combustion model. 7 Summary and outlook.