In Diesel engines, the internal flow characteristics in the fuel injection nozzles, such as the turbulence level and distribution, the cavitation pattern and the velocity profile affect significantly the air-fuel mixture in the spray and subsequently the combustion process. Since the possibility to observe experimentally and measure the flow inside real size Diesel injectors is very limited, Computational Fluid Dynamics (CFD) calculations are generally used to obtain the relevant information. The work presented within this thesis is focused on the study of cavitation in real size automotive injectors by using a commercial CFD code. It is divided in three major phases, each corresponding to a different complementary objective. The first objective of the current work is to assess the ability of the cavitation model included in the CFD code to predict cavitating flow conditions. For this, the model is validated for an injector-like study case defined in the literature, and for which experimental data is available in different operating conditions, before and after the start of cavitation. Preliminary studies are performed to analyze the effects on the solution obtained of various numerical parameters of the cavitation model itself and of the solver, and to determine the adequate setup of the model. It may be concluded that overall the cavitation model is able to predict the onset and development of cavitation accurately. Indeed, there is satisfactory agreement between the experimental data of injection rate and choked flow conditions and the corresponding numerical solution.This study serves as the basis for the physical and numerical understanding of the problem. Next, using the model configuration obtained from the previous study, unsteady flow calculations are performed for real-size single and multi-hole sac type Diesel injectors, each one with two types of nozzles, tapered and cylindrical. The objective is to validate the model with real automotive cases and to ununderstand in what way some physical factors, such as geometry, operating conditions and needle position affect the inception of cavitation and its development in the nozzle holes. These calculations are made at full needle lift and for various values of injection pressure and back-pressure. The results obtained for injection rate, momentum flux and effective injection velocity at the exit of the nozzles are compared with available CMT-Motores Térmicos in-house experimental data. Also, the cavitation pattern inside the nozzle and its effect on the internal nozzle flow is analyzed. The model predicts with reasonable accuracy the effects of geometry and operating conditions. / Patouna, S. (2012). A CFD STUDY OF CAVITATION IN REAL SIZE DIESEL INJECTORS [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/14723 / Palancia

El proceso de atomización desde una vena o lámina líquida hasta multitud de gotas dispersas en un medio gaseoso ha sido un fenómeno de interés desde hace varias décadas, especialmente en el campo de los motores de combustión interna alternativos. Multitud de estudios experimentales han sido publicados al respecto, pues una buena mezcla de aire-combustible asegura una evaporación y combustión mucho más eficientes, aumentando la potencia del motor y reduciendo la cantidad de contaminantes emitidos. Con el auge de las técnicas computacionales, muchos modelos han sido desarrollados para estudiar este proceso de atomización y mezcla. Uno de los últimos modelos que han aparecido es el llamado ELSA (Eulerian-Lagrangian Spray Atomization), que utiliza un modelo Euleriano para la parte densa del chorro y cambia a un modelo Lagrangiano cuando la concentración de líquido es suficientemente pequeña, aprovechando de esta manera las ventajas de ambos. En el presente trabajo se ha desarrollado un modelo puramente Euleriano para estudiar la influencia de la geometría interna de la tobera de inyección en el proceso de atomización y mezcla. Se ha estudiado únicamente el proceso de inyección diésel. Este modelo permite resolver en un único dominio el flujo interno y el externo, evitando así las comunes simplificaciones y limitaciones de la interpolación entre ambos dominios resueltos por separado. Los resultados actuales son prometedores, el modelo predice con un error aceptable la penetración del chorro, el flujo másico y de cantidad de movimiento, los perfiles de velocidad y concentración, así como otros parámetros característicos del chorro. / Martí Gómez-Aldaraví, P. (2014). DEVELOPMENT OF A COMPUTATIONAL MODEL FOR A SIMULTANEOUS SIMULATION OF INTERNAL FLOW AND SPRAY BREAK-UP OF THE DIESEL INJECTION PROCESS [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/43719 / TESIS / Premios Extraordinarios de tesis doctorales

Plug nozzle, a passive altitude adaptive nozzle, for futuristic SSTO applications, exhibits greater efficiency as compared to conventional nozzles over a wide range of altitudes. The plug nozzle comprises of a primary nozzle and a contoured plug; an under–expanded jet exiting the primary nozzle is allowed to further expand over the plug surface for altitude adaptation. At design condition the flow expands correctly to the ambient conditions on the full length plug surface, while at off design conditions the flow adapts to the ambient conditions through wave interactions within the nozzle core jet. Based on thrust to weight considerations, the full length plug is truncated and this results in a base flow rich in flow physics. In addition, the base flow exhibits an interesting transitional behaviour from open wake to a closed wake because of the wave interactions within the nozzle core jet. The plug surface flow can further exhibit flow complexities because of wave interactions resulting from the shear layer emanating from the splitter plates, in case of clustered plug flows. Considering these flow complexities, the design of the plug nozzles and analysing the associated flows can be a challenge to the aerodynamic community. An attempt has been made in understanding this class of flows in this thesis. This objective has been accomplished using both experimental and computational tools. In the present work, both the linear and annular plug nozzle geometries have been analysed for a wide range of pressure ratios spanning from 5to 80. The linear and annular nozzles have been designed for similar flow conditions and their respective design pressure ratios are 60and 66. From the experimental and computational results, it has been shown that the computational solver performs well in predicting the wave interactions on the plug surface. In addition the limitations of the computational solver in predicting the plug base flows in general has been brought out. This limitation in itself need not be considered as a serious handicap in the design and analysis of plug nozzle flows; this is because the plug base contribution to the thrust is very minimal, as has been brought out in this thesis. Apart from this the high quality experimental data generated is also of immense value to the CFD community as this also serves as a valuable data base for CFD code validation. For analysis, the plug flow field has been categorized into three different regimes based on the primary nozzle lip expansion fan extent. The flow field is categorised based on the reflection of the primary nozzle lip expansion fan from plug surface, base region shear layer and symmetry line downstream of the base region recirculation bubble. This flow division is particularly helpful in understanding the base wake characteristics with increasing pressure ratio. The base lip pressure and the base pressure variation have been discussed with respect to the primary nozzle lip expansion fan extent. In the open wake regime (or for low pressure ratios) the wave interactions within the core jet flow impinge on the base region shear layer. Because of these interactions it is difficult to propose an empirical model for open wake base pressure. In the closed wake regime (for higher pressure ratios), the base region recirculation bubble is completely under the shower of primary nozzle lip expansion fan. Hence the base lip pressure and base pressure are frozen with respect to stagnation conditions. Based on these insights it was possible to propose empirical models for linear and annular closed wake base pressure. Along with these, a mathematical model defining a reference pressure ratio PR∗, beyond which the closed wake base pressure is expected to be more than the ambient pressure has also been proposed. This is expected to serve as a good design parameter. In case of linear plug flows, this also serves the purpose of base wake transition, for the cases considered in this thesis. The flow expansion process or the primary nozzle lip expansion fan extent was also useful in understanding the differences between the linear and annular plug nozzle flow fields. In a linear plug nozzle, the flow expands only in the streamwise direction while in an annular plug nozzle the flow expands both along the streamwise and azimuthal directions. The flow expands at a faster rate in case of annular nozzle as against linear nozzle. Hence differences are observed between the linear and annular nozzle on plug and base surfaces. On the annular plug surface more wave interactions are observed because of faster expansion. With regard to base characteristics, faster expansion in annular plug nozzle, with respect to linear nozzle, results in a lower base lip pressure, lower base pressure and higher wake transition pressure ratio. The realistic cluster plug configurations have also been considered for the present studies. The effects of clustering on the plug nozzle flow field have been brought out by considering two different linear cluster nozzles and one annular cluster nozzle. The differences in the flow field of a simple and cluster plug nozzle has been discussed. In case of simple plug nozzle wave interactions are observed only in the stream wise direction, while in case of cluster plug nozzle three dimensional wave interactions are observed because of the splitter plates. Along the splitter plate differential end conditions introduce a curved recompression shock on the plug surface. This recompression shock in turn induces a streamwise vortex and also a secondary shock. It has been observed that differences between the simple and cluster plug surface pressure field are because of three dimensional wave interactions. Regarding the base pressure, differences between the simple and cluster geometries were observed for shorter truncation plug lengths (20% length plug). While for longer plug lengths (more than 34% length) the effects of clustering were reduced on the base pressure. Regarding the transition pressure ratio, differences were observed between simple and clustered plug nozzles for all the plug lengths considered. In addition, the performance of the plug nozzles has been carried out. From the analysis it was found that the primary nozzle and plug surface are major contributors towards thrust. The base surface contributes only about 2– 3% of the thrust at design condition. Hence, from a design point of view, a computational solver can be a useful tool considering its efficacy on the plug surface and in the primary nozzle.

Plug nozzle, a passive altitude adaptive nozzle, for futuristic SSTO applications, exhibits greater efficiency as compared to conventional nozzles over a wide range of altitudes. The plug nozzle comprises of a primary nozzle and a contoured plug; an under–expanded jet exiting the primary nozzle is allowed to further expand over the plug surface for altitude adaptation. At design condition the flow expands correctly to the ambient conditions on the full length plug surface, while at off design conditions the flow adapts to the ambient conditions through wave interactions within the nozzle core jet. Based on thrust to weight considerations, the full length plug is truncated and this results in a base flow rich in flow physics. In addition, the base flow exhibits an interesting transitional behaviour from open wake to a closed wake because of the wave interactions within the nozzle core jet. The plug surface flow can further exhibit flow complexities because of wave interactions resulting from the shear layer emanating from the splitter plates, in case of clustered plug flows. Considering these flow complexities, the design of the plug nozzles and analysing the associated flows can be a challenge to the aerodynamic community. An attempt has been made in understanding this class of flows in this thesis. This objective has been accomplished using both experimental and computational tools. In the present work, both the linear and annular plug nozzle geometries have been analysed for a wide range of pressure ratios spanning from 5to 80. The linear and annular nozzles have been designed for similar flow conditions and their respective design pressure ratios are 60and 66. From the experimental and computational results, it has been shown that the computational solver performs well in predicting the wave interactions on the plug surface. In addition the limitations of the computational solver in predicting the plug base flows in general has been brought out. This limitation in itself need not be considered as a serious handicap in the design and analysis of plug nozzle flows; this is because the plug base contribution to the thrust is very minimal, as has been brought out in this thesis. Apart from this the high quality experimental data generated is also of immense value to the CFD community as this also serves as a valuable data base for CFD code validation. For analysis, the plug flow field has been categorized into three different regimes based on the primary nozzle lip expansion fan extent. The flow field is categorised based on the reflection of the primary nozzle lip expansion fan from plug surface, base region shear layer and symmetry line downstream of the base region recirculation bubble. This flow division is particularly helpful in understanding the base wake characteristics with increasing pressure ratio. The base lip pressure and the base pressure variation have been discussed with respect to the primary nozzle lip expansion fan extent. In the open wake regime (or for low pressure ratios) the wave interactions within the core jet flow impinge on the base region shear layer. Because of these interactions it is difficult to propose an empirical model for open wake base pressure. In the closed wake regime (for higher pressure ratios), the base region recirculation bubble is completely under the shower of primary nozzle lip expansion fan. Hence the base lip pressure and base pressure are frozen with respect to stagnation conditions. Based on these insights it was possible to propose empirical models for linear and annular closed wake base pressure. Along with these, a mathematical model defining a reference pressure ratio PR∗, beyond which the closed wake base pressure is expected to be more than the ambient pressure has also been proposed. This is expected to serve as a good design parameter. In case of linear plug flows, this also serves the purpose of base wake transition, for the cases considered in this thesis. The flow expansion process or the primary nozzle lip expansion fan extent was also useful in understanding the differences between the linear and annular plug nozzle flow fields. In a linear plug nozzle, the flow expands only in the streamwise direction while in an annular plug nozzle the flow expands both along the streamwise and azimuthal directions. The flow expands at a faster rate in case of annular nozzle as against linear nozzle. Hence differences are observed between the linear and annular nozzle on plug and base surfaces. On the annular plug surface more wave interactions are observed because of faster expansion. With regard to base characteristics, faster expansion in annular plug nozzle, with respect to linear nozzle, results in a lower base lip pressure, lower base pressure and higher wake transition pressure ratio. The realistic cluster plug configurations have also been considered for the present studies. The effects of clustering on the plug nozzle flow field have been brought out by considering two different linear cluster nozzles and one annular cluster nozzle. The differences in the flow field of a simple and cluster plug nozzle has been discussed. In case of simple plug nozzle wave interactions are observed only in the stream wise direction, while in case of cluster plug nozzle three dimensional wave interactions are observed because of the splitter plates. Along the splitter plate differential end conditions introduce a curved recompression shock on the plug surface. This recompression shock in turn induces a streamwise vortex and also a secondary shock. It has been observed that differences between the simple and cluster plug surface pressure field are because of three dimensional wave interactions. Regarding the base pressure, differences between the simple and cluster geometries were observed for shorter truncation plug lengths (20% length plug). While for longer plug lengths (more than 34% length) the effects of clustering were reduced on the base pressure. Regarding the transition pressure ratio, differences were observed between simple and clustered plug nozzles for all the plug lengths considered. In addition, the performance of the plug nozzles has been carried out. From the analysis it was found that the primary nozzle and plug surface are major contributors towards thrust. The base surface contributes only about 2– 3% of the thrust at design condition. Hence, from a design point of view, a computational solver can be a useful tool considering its efficacy on the plug surface and in the primary nozzle.

This thesis studies the influence of internal nozzle flow characteristics over a large spectrum of experimental conditions and diagnostics. Experiments were carried out for two nozzle geometries---cylindrical and conical single hole nozzles---and three different fuels. Two of the fuels are pure components---n-heptane and n-dodecane---while the third fuel consists of a three-component surrogate to better represent the physical and chemical properties of diesel fuel. Measurements include a complete hydraulic characterization consisting of instantaneous injection rate and spray momentum flux measurements; a high-speed visualization of isothermal liquid spray; a high-speed visualization of the evaporative inert spray, imaging liquid and vapor phases simultaneously and finally, a high-speed visualization of the high temperature reactive spray, imaging vapor phase and OH* chemiluminescence for each injection event. All high-temperature diagnostics were performed in a continuous flow test chamber that allows an accurate control on a wide range of thermodynamic conditions (up to 1000 K and 15 MPa). The experimental findings from this work, and the large database obtained (available for download at: http://www.cmt.upv.es/DD01.aspx), could be used to validate CFD models that could help the community understand the fundamental driving mechanisms behind these observations. / En esta tesis se estudia la influencia del flujo interno sobre un amplio espectro de condiciones y diagnósticos experimentales. Se realizaron experimentos para dos geometrías de tobera---toberas cilíndrica y cónica de un único orificio---y tres combustibles. Dos de los combustibles son puros---n-heptano y n-dodecano--- mientras el tercero es un combustible sustituto que consiste en una mezcla de tres componentes que busca representar mejor las propiedades físicas y químicas del diesel. Las medidas incluyen una caracterización hidráulica completa, compuesta por tasa de inyección y cantidad de movimiento instantáneas; una visualización de alta velocidad del chorro líquido isotermo; una visualización de alta velocidad del chorro inerte evaporativo, con captura simultánea de las fases líquida y vapor y, finalmente, una visualización del chorro reactivo a alta temperatura, con captura de la fase vapor y la quimioluminiscencia del radical OH* para cada evento de inyección. Todos los diagnósticos en condiciones de alta temperatura fueron realizados en una maqueta de alta presión y temperatura de flujo constante que permite controlar con precisión un rango amplio de condiciones termodinámicas (hasta 1000 K y 15 MPa). Los resultados experimentales y la gran base de datos obtenida en este trabajo (disponible en: http://www.cmt.upv.es/DD01.aspx), podrían ser utilizados para validar modelos CFD detallados que podrían ayudar a la comunidad científica a entender mejor los mecanismos fundamentales que producen los resultados observados. / Aquesta tesi estudia la influència del flux intern sobre un gran espectre de condicions i diagnòstics experimentals. Es van realitzar experiments per a dos geometries de tovera---toveres ci¿líndrica i cónica amb un únic orifici---i tres combustibles. Dos dels combustibles són purs---n-heptà i n-dodecà--- mentre el tercer combustible consisteix en una mescla de tres components que formen un combustible substitut que busca representar millor les propietats físiques i químiques del dièsel. Les mesures inclouen una caracterització hidràulica completa, composta per taxa d'injecció i quantitat de moviment instantanis; visualització d'alta velocitat del doll líquid isoterme; visualització d'alta velocitat del doll inert evaporatiu, capturant simultàniament les fases líquid i vapor i, finalment, una visualització del doll reactiu a alta temperatura, capturant la fase vapor i la quimioluminiscència del radical OH per a cada esdeveniment d'injecció. Tots els diagnòstics en condicions d'alta temperatura van ser realitzats en una insta¿lació d'alta pressió i temperatura amb flux constant que permet controlar amb precisió un ampli rang de condicions termodinàmiques (fins a 1000 K i 15 MPa). Els resultats experimentals i la gran base de dades obtinguda en aquest treball (disponible a la web en: http://www.cmt.upv.es/dd01.aspx), podrien ser utilitzats per tal de validar models CFD detallats que podrien ajudar a la comunitat científica a entendre millor els mecanismes fonamentals que produeixen aquestes observacions. / Viera Sotillo, JP. (2017). Experimental study of the effect of nozzle geometry on the performance of direct-injection diesel sprays for three different fuels [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/81857 / TESIS

Many industrial applications acquire necessary thermal energy through the combustion process. The basic element of each combustion appliance is a burner and one~part~of~it~is a~nozzle system that supplies fuel to a combustion chamber. The geometry of the fuel nozzle significantly affects the intensity of mixing the fuel with the combustion air and thus the stability of the combustion. The main subject of~this diploma thesis is~determination of~velocity coefficients for nozzles with different geometries. The knowledge of~correct values of~velocity coefficients is a key parameter for the design of~the burner and~its subsequent operation. For the calculation of~velocity coefficients, the exit nozzle velocities were used. For chosen nozzles, a~theoretical exit nozzle velocities were calculated. They were compared with the actual exit nozzle velocities, which were measured experimentally using the hot-wire anemometry, and with velocities, which were calculated using the CFD simulation method.

Convergent-Divergent (C-D) nozzles are used in rocket engines to produce thrust as a reaction to the acceleration of hot combustion chamber gases in the opposite direction. To maximize the engine performance at high altitudes, large area ratio, bell-shaped or contoured nozzles are used. At lower altitudes, the exit pressure of these nozzles is lower than the ambient pressure. During this over-expanded condition, the nozzle-internal flow adapts to the ambient pressure through an oblique shock. But the boundary layer inside the divergent portion of the nozzle is unable to withstand the pressure rise associated with the shock, and consequently flow separation is induced. Numerical simulation of separated flows in rocket nozzles is challenging because the existing turbulence models are unable to correctly predict shock-induced flow separation. The present thesis addresses this problem. Axisymmetric, steady-state, Reynolds-Averaged Navier-Stokes (RANS) simulations of a conical nozzle and three sub-scale contoured nozzles were carried out to numerically predict flow separation in over-expanded rocket nozzles at different nozzle pressure ratios (NPR). The conical nozzle is the JPL 45◦-15◦ and the contoured nozzles are the VAC-S1, the DLR-PAR and the VAC-S6-short. The commercial CFD code ANSYS FLUENT 13 was first validated for simulation of separated cold gas flows in the VAC-S1 nozzle. Some modeling issues in the numerical simulations of flow separation in rocket nozzles were determined. It is recognized that compressibility correction, nozzle-lip thickness and upstream-extension of the external domain are the sources of uncertainty, besides turbulence modeling. In high-speed turbulent flows, compressibility is known to affect dissipation rate of turbulence kinetic energy. As a consequence, a reduction in the spreading rate of supersonic mixing layers occurs. Whereas, the standard turbulence models are developed and calibrated for incompressible flows and hence, do not account for this effect. ANSYS FLUENT uses the compressibility correction proposed by Wilcox [1] which modifies the turbulence dissipation terms based on turbulent Mach number. This, as shown in this thesis, may not be appropriate to the prediction of flow separation in rocket nozzles. Simulation results of the standard SST model, with and without the compressibility correction, are compared with the experimental data at NPR=22 for the DLR-PAR nozzle. Compressibility correction is found to cause under-prediction of separation location and hence its use in the prediction of flow separation is not recommended. In the literature, computational domains for the simulation of DLR subscale nozzles have thick nozzle-lips whereas for the VAC subscale nozzles they have no nozzle-lip. Effect of nozzle-lip thickness on flow separation is studied in the DLR-PAR nozzle by varying its nozzle-lip thickness. It is found that nozzle-lip thickness significantly influences both separation location and post-separation pressure recovery by means of the recirculation bubbles formed at the nozzle-lip. Usually, experimental values of free stream turbulence are unknown. So conventionally, to minimize solution dependence on the boundary conditions specified for the ambient flow, the computational domain external to the nozzle is extended in the upstream direction. Its effect on flow separation is studied in the DLR-PAR nozzle through simulations conducted with and without this domain extension. No considerable effect on separation location and pressure recovery is found. The two eddy-viscosity based turbulence models, Spalart-Allmaras (SA) model and Shear Stress Transport (SST) model, are well known to predict separation location better than other eddy-viscosity models, but with moderate success. Their performances, in terms of predicting separation location and post-separation wall pressure distribution, were compared with each other and evaluated against experimental data for the conical and two contoured nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the range of NPRs and to the type of nozzle. Depending on NPR, the SST model either under-predicts or over-predicts Free Shock Separation (FSS). Moreover, it also fails to capture Restricted Shock Separation (RSS). With compressibility correction, it under-predicts separation at all NPRs to a greater extent. Even though RSS is captured by using compressibility correction, the transition from FSS to RSS is over-predicted [2]. Early efforts by few researchers to improve predictions of nozzle flow separation by realizability corrections to turbulence models have not been successful, especially in terms of capturing both the separation types. Therefore, causes of turbulence modeling failure in predicting nozzle flow separation correctly were further investigated. It is learnt that limiting of the shear stress inside boundary layer, due to Bradshaw’s assumption, and over-prediction of jet spreading rate are the causes of SST model’s failure in predicting nozzle flow separation correctly. Based on this physical reasoning, values of the a 1 parameter and the two diffusion coefficients σk,2 and σω,2 were empirically modified to match the predicted wall pressure distributions with experimental data of the DLR-PAR and the VAC-S6-short nozzles. The results confirm that accurate prediction of flow separation in rocket nozzles indeed depends on the correct prediction of spreading rate of the supersonic separation-jet. It is demonstrated that accurate RANS simulation of flow separation in rocket nozzles over a wide range of NPRs is feasible by modified values of the diffusion coefficients in turbulence model.

Vegetable oil methyl esters obtained by transesterification of vegetable oils are considered to be suitable alternative fuels for diesel engines. However, higher viscosity, surface tension and boiling temperatures of biodiesels may adversely affect spray characteristics as compared to those of diesel. Thus, spray characteristics of Jatropha Methyl Ester (JME) are studied by comparing them to those of diesel in a high-pressure chamber with optical access to simulate the actual in-cylinder conditions. Also, the effect of inner-nozzle cavitation on JME and diesel sprays is studied by utilizing two nozzles, one with sharp entry-radius and the other with larger entry-radius. Finally, spray characteristics of surrogate fuels such as n-dodecane and n-hexadecane are also studied. The first part of the work concerning precise measurements of inner-nozzle geometry revealed that one of the nozzles has a hole diameter of 190-µm and entry-radius of around 70-µm, while the other has a hole diameter of 208-µm and entry-radius of around 10-µm. Injection rate-shape and coefficient of discharge for JME and diesel flow through the two nozzles were then measured. It was observed that while the coefficients of discharge (Cd) are almost identical for JME and diesel, the nozzle with entry radius of 10-µm exhibited around 20% lower Cd than that of the entry-radius of 70-µm. This observation coupled with insight from complementary CFD simulations of inner-nozzle flow showed that the lower Cd of the nozzle with entry-radius of 10-µm could be attributed to inner-nozzle cavitation. The second part of the work involved measurement of non-evaporating spray characteristics including spray-tip penetration, spray-cone angle and droplet size measurement under realistic operating conditions using techniques such as Shadowgraphy and Particle/Droplet Imaging Analysis (PDIA). The non-evaporating spray of the fuels are studied by injecting them using a common-rail fuel injection system into the high-pressure chamber maintained at room temperature. Experimental results show that JME is associated with a slightly faster spray-tip penetration and narrow spray-cone angle indicating inferior spray atomization which is confirmed by around 5% larger droplet sizes. Slower spray-tip penetration, wider spray-cone angle and around 5% smaller droplet sizes are observed for the spray from the cavitating nozzle. Thus, the inner nozzle cavitation is observed to improve the atomization of diesel and JME sprays. The differences in spray characteristics of JME and diesel reduce as the injection pressure increases. The spray-tip penetrations of both surrogates are observed to almost match that of diesel. The third part of the work involved measurements of evaporating spray liquid length, vapour penetration and spread angle for JME, diesel and surrogates at conditions of 50 bar chamber pressure and 900 K temperature. It is observed that the JME exhibits around 16% longer liquid length than that of diesel. The liquid length of n-dodecane is significantly lower than that of diesel and liquid length of n-hexadecane is around 20% higher than that of n-dodecane mimicking the trend of JME and diesel. The liquid length of n-hexadecane is very close to that of diesel at all the three test conditions. Interestingly, the vapour penetration and spread angle for all the fuels is observed to be almost identical. As the cold spray and evaporating spray characteristics of n-hexadecane match well with those of diesel, n-hexadecane can be chosen as a pure component surrogate for diesel. Finally, an analytical model for predicting the spray vapour penetration is assessed with the experimentally-observed trends of penetration and spray spread angle. The model indicated that the effect of fuel density variation is compensated by the corresponding variation in injection velocity for a given injection pressure to result in a similar vapour penetration. Overall, the present work, in addition to studying the effect of fuel physical properties and cavitation on sprays, has generated a comprehensive experimental database on non-evaporating and evaporating sprays of biodiesel, diesel, and pure component surrogates, which would aid significantly in validation of CFD simulations.

[ES] La creciente preocupación por los problemas medioambientales, la disponibilidad de combustibles fósiles unido a la gran demanda de vehículos, han llevado a los gobiernos a regular las emisiones emitidas a la atmósfera. Existen propuestas de adoptar fuentes de energía renovables. Sin embargo, la sustitución de los combustibles derivados del petróleo no será fácil, rápida o rentable, y el transporte propulsado por motores de combustión interna (ICE) seguirá destacando en los próximos años. La eficiencia de la combustión y el rendimiento del motor están influenciados por el complejo proceso de inyección. La inyección directa de gasolina (GDI) aumenta el ahorro de combustible y cumple los requisitos de emisiones contaminantes, aunque queda potencial por descubrir. Por ello, ha sido objeto de estudio en los últimos años y, en consecuencia, de la presente Tesis. Este trabajo tiene como motivación mejorar el entendimiento en el campo del GDI. La compleja naturaleza transitoria del proceso de inyección hace que el estudio experimental sea un desafío. La Mecánica de Fluidos Computacional (CFD) surge como una potente alternativa a los experimentos y ha sido adoptada para esta investigación. Bajo este contexto, el objetivo de la presente Tesis es desarrollar una metodología predictiva para la caracterización hidráulica del inyector, capaz de ser aplicada a las actuales y futuras generaciones de inyectores GDI, independientemente de las características del inyector y del software de estudio. Una vez validada, el objetivo posterior es utilizar los resultados para analizar el comportamiento del chorro. Este enfoque busca seguir los pasos de la comunidad científica sustituyendo la práctica experimental. La validación de la metodología se lleva a cabo mediante su aplicación en dos inyectores GDI solenoides multi-orificio diferentes. Además, se han utilizado dos códigos CFD comerciales: CONVERGE y StarCCM+. La metodología predictiva se centra en el estudio del flujo interno y el campo cercano para caracterizar hidráulicamente el inyector. El problema a tratar se define como un sistema multifásico en un marco Euleriano y considerando un único fluido. El tratamiento del flujo multifásico se realiza mediante el enfoque Volume-of-Fluid (VOF). Además, se emplea el Homogeneous Relaxation Model (HRM) para considerar el intercambio de masa entre las fases líquida y vapor debido a cavitación y flash boiling. La turbulencia se ha tratado a partir de los enfoques Reynolds-Averaged Navier-Stokes (RANS) y Large Eddy Simulations (LES). Por otro lado, en cuanto al estudio del flujo externo, se ha adoptado el Discrete Droplet Model (DDM). La atomización y el chorro están influenciados por la geometría de la tobera, por lo que la estrategia de acoplamiento del flujo interno y externo complementa los análisis. Se han adoptado enfoques de acoplamiento unidireccional y mapeado, utilizando como parámetros de entrada los datos de flujo interno de la validada metodología. Esta Tesis aporta una nueva y valiosa metodología predictiva con una elevada precisión a la hora de caracterizar el proceso de inyección en comparativa con datos experimentales. Por otro lado, es directamente trasferible a distintos códigos de cálculo así como aplicable a inyectores con características dispares sin perjudicar las exigencias del modelo. La correcta caracterización del flujo interno ha permitido emplear los datos obtenidos para analizar el comportamiento del chorro eliminando la necesidad de usar datos experimentales. Los resultados obtenidos capturan el comportamiento macroscópico del chorro con una precisión comparable a los experimentos. Aunque todavía hay muchos retos que afrontar, la presente Tesis supone un gran avance en el campo del GDI. El remarcable progreso se debe al desarrollo y uso de una metodología totalmente predictiva, que permite prescindir de la mayoría de los experimentos para contribuir a una mayor y más amplia visión de la física del proceso de inyección. / [CA] La creixent preocupació pels problemes ambientals, la limitada disponibilitat de combustibles fòssils, acompanyat a la gran demanda de vehicles, ha portat el govern a regular els nivells d'emissions emesos a l'atmosfera. Existeixen propostes d'adoptar fonts d'energia renovables. Tanmateix, la substitució dels combustibles líquids derivats del petroli no es durà a terme de forma fàcil, ràpida o rentable, i el transport propulsat per motors de combustió interna (ICE) continuarà destacant en els pròxims anys. L'eficiència de la combustió i el rendiment del motor són fortament influenciats pel complex procés d'injecció. La injecció directa de gasolina (GDI) augmenta l'estalvi de combustible i complix amb els requisits d'emissions, encara que queda molt potencial per descobrir. Per això, aquest ha sigut objecte d'investigació en els últims anys i, com a conseqüència, d'aquesta Tesi. Aquest treball té com a motivació millorar l'enteniment en el camp del GDI. La complexa natura transitòria de la injecció fa que l'estudi experimental siga força complex. La Mecànica de Fluids Computacional (CFD) sorgeix com una potent alternativa als experiments, i ha sigut adoptada per aquesta investigació. Baix aquest mateix context, es proposa com a objectiu principal d'aquesta Tesi el desenvolupament d'una metodologia predictiva per a la caracterització hidràulica de l'injector, capaç de ser aplicada a les actuals i futures generacions d'injectors GDI (independentment de les característiques de l'injector i del software d'estudi). Una vegada validada, el posterior objectiu és analitzar el comportament de l'esprai. Aquest enfocament busca seguir els passos de la comunitat científica substituint la pràctica experimental. La validació de la metodologia ha sigut duta a terme mitjançant la seva aplicació en dos injectors GDI solenoides multi-orifici. A més, s'han utilitzat dos software CFD comercials: CONVERGE i StarCCM+. La metodologia predictiva se centra en l'estudi del flux intern i el camp proper per tal de caracteritzar hidràulicament l'injector. El problema a tractar es defineix en base a un sistema multi-fàsic en un marc Eulerià i considerant un únic fluid. El tractament del fluid multi-fàsic es realitza mitjançant l'aproximació Volume-of-Fluid (VOF). A més, s'utilitza el Homogeneous Relaxation Model (HRM) per tal de considerar l'intercambi de massa entre les fases líquida i vapor degut als fenòmens de cavitació i flash boiling. La turbulència s'ha tractac a través dels enfocaments Reynolds-Averaged Navier-Stokes (RANS) i Large Eddy Simulations (LES). Pel que fa a l'estudi del fluix extern, s'ha adoptat el Discrete Droplet Model (DDM). Sent conscients que el comportament l'atomització i l'esprai estan influenciats per la geometria de la tovera, l'estratègia d'acoblament del flux intern i extern complementa les anàlisis. S'han adoptat els enfocaments d'acoblament unidireccional i mapejat, utilitzant com a paràmetres d'entrada les dades del flux intern obtingudes amb la validada metodologia. Aquesta Tesi aporta una nova i valuosa metodologia predictiva amb una elevada precisió a l'hora de caracteritzar el procés d'injecció en comparativa amb dades experimentals. És directament transferible a diversos codis de càlcul així com aplicable a injectors amb característiques dispars sense perjudicar les exigències del model. La correcta caracterització del flux intern ha permès utilitzar les dades obtingudes per tal d'analitzar el comportament de l'esprai, eliminant la necessitat d'emprar dades experimentals. Els resultats obtinguts d'aquest estudi capturen el comportament macroscòpic de l'esprai amb una precisió comparable als experiments. Encara que queden molts reptes per afrontar, aquesta Tesi aporta un important avanç al camp del GDI. La ruptura prové del desenvolupament i ús d'una metodologia completament predictiva, que substitueix els experiments requerits i així contribueix a una millor i més ampla visió de la física del procés d'injecció. / [EN] Concerns about climate change, availability of fuel resources and the high demand for vehicles, have led governments to regulate the level of pollution emitted by engines into the atmosphere. There is a strong desire to adopt renewable and sustainable energy sources. However, the substitution of liquid fuels derived from petroleum will not emerge easily, quickly or economically, and Internal Combustion Engines (ICE) will continue to excel for the next few years. Combustion efficiency and engine performance are strongly influenced by the complex fuel injection process. Gasoline Direct Injection (GDI) strategies increase fuel economy and meet emission requirements, although many challenges remain, which has therefore been one of the main research objectives in recent years and of this Thesis. The present research aims to provide a better understanding in the field of GDI. The transient and complex nature of the injection process makes the experimental study of GDI quite challenging. Therefore, Computational Fluid Dynamics (CFD) emerges as a powerful alternative adopted for this research. In this context, the main objective of the present Thesis is to develop a predictive methodology capable of being applied to current and future generations of GDI injectors, regardless of the injector features and the software employed, for the hydraulic characterization of the injector. Once validated, the subsequent goal is to employ the obtained results to analyze the behavior of the spray downstream of the injector. The approach attempts to follow the footsteps of the research community to avoid experimental practice. The predictive methodology has been validated through its application to two multi-hole solenoid GDI injectors with different features. In addition, the mentioned methodology has been evaluated using diverse commercial software: CONVERGE and StarCCM+. The methodology focuses on the study of the internal and near-field flow to hydraulically characterize the injector. So the problem to be addressed is a multi-phase system, performed in an Eulerian framework, modeled through a single-fluid approach. The multi-phase flow is treated by means of the Volume-of-Fluid (VOF) approach. Homogeneous Relaxation Model (HRM) is employed to consider the mass exchange between liquid and vapor fuel phases, due to cavitation and flash boiling. The turbulence treatment has been performed from both Reynolds-Averaged Navier-Stokes (RANS) and Large Eddy Simulations (LES) approaches. Regarding the external flow study, the Discrete Droplet Model (DDM) has been adopted. In addition, being aware that atomization and spray behavior is greatly influenced by the nozzle geometry, the coupling strategy of the internal and external flow complements the analyses. One-way coupling and mapping approaches have been adopted, using as input parameters the internal flow data obtained from the already validated methodology. Accordingly, this Thesis provides a new and valuable predictive methodology, which has demonstrated a high accuracy in characterizing the flow behavior during the injection process through comparison with experimental data. It has also proven to be directly transferable to different CFD software and applicable to injectors with dissimilar characteristics without compromising the requirements of the model. The correct internal flow characterization has made it possible to employ the obtained data to analyze the spray patterns, which eliminates the need to consider experimental data. The outcomes of this study macroscopically capture the jet behavior with an accuracy comparable to experiments under different operating conditions. Although there are still many challenges to face, the present Thesis brings a breakthrough in the field of GDI. The quantum leap arises from the development and use of a fully predictive methodology, allowing to avoid most experiments to contribute to a greater and broader vision of the injection process physics. / María Martínez García has been founded through a grant from the Government of Generalitat Valenciana with reference ACIF/2018/118 and financial support from the European Union. These same institutions, Government of Generalitat Valenciana and the European Union, supported through a grant for pre-doctoral stays out of the Comunitat Valenciana with reference BEFPI/2020/057 the research carried out during the stay at Aerothermochemistry and Combustion Systems Laboratory, Swiss Federal Institute of Technology, ETH Zurich, Switzerland. Special gratitude from the author to both institutions, Government of Generalitat Valenciana and the European Union, for making this dream possible / Martínez García, M. (2022). Computational Study of the Injection Process in Gasoline Direct Injection (GDI) Engines [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/185180 / TESIS