L'objectif de cette étude est d'analyser l'écoulement traversant le récepteur solaire haute température de la centrale PEGASE. Cet écoulement est turbulent et soumis à un très fort gradient de température. Nous utilisons la simulation des grandes échelles thermiques appliquée aux équations bas-Mach pour étudier l'écoulement turbulent dans une géométrie de canal plan bi-périodique avec températures imposées aux parois. Pour cette étude, nous avons réalisé des simulations pour deux nombres de Reynolds turbulents différents (Re_tau m = 180 et Re _taum = 395) et pour quatre rapports de température (T2=T1 = 1 ; 1; 04 ; 2 et 5). Après avoir validé notre modèle, nous avons réalisé une étude sur la modélisation sous-maille thermique qui conclut sur la nécessité d'utiliser un modèle sous-maille thermique dynamique pour un écoulement ayant une forte intensité turbulente soumis à un très fort gradient de température. Nous analysons l'influence du gradient de température sur les différents profils de vitesse et de température ainsi que sur les spectres d'énergie. Augmenter le gradient de température crée une dissymétrie de tous les profils. Une relaminarisation du côté chaud est visible. On remarque une nouvelle répartition des fluctuations à travers le canal ainsi que de l'énergie en fonction de la taille des échelles. Une création de fluctuations de vitesse et de température a également été mise évidence du côté froid et du côté chaud du canal. Ces effets sont dus, à la fois à une variation de la viscosité cinématique (effet direct de la température) et à une interaction entre les champs turbulents dynamique et thermique.

La demande croissante pour alléger les moteurs d'avions et diminuer les émissions polluantes de la propulsion aéronautique réclame à rendre plus compact le système de compression des moteurs, qui représente environ 40%-50% de la masse totale. Or, à taux de compression global égal, la réduction du nombre d'étage exige d'un point de vue aérodynamique une augmentation de la charge des aubes de compresseur par étage. La charge d'aube est aujourd'hui limitée car elle induit différents mécanismes de pertes tridimensionnelles très pénalisant. L'un des plus importants est le décollement de coin qui se forme à la jonction entre l'extrados de l'aube et le moyeu ou le carter. Bien que des travaux existent sur les mécanismes et paramètres intervenant dans le décollement de coin, il est encore difficile de proposer une méthode de contrôle efficace. Cela est principalement dû à deux raisons : (i) le manque de compréhension fine des mécanismes physiques, (ii) l'utilisation pour la conception de modèles de turbulence classiques de type RANS qui ne sont pas capables de prédire précisément le décollement de coin, car ils ne peuvent pas décrire correctement les mécanismes de transport turbulent. Des simulations de type RANS et LES sont présentées dans cette thèse sur une configuration de grille d'aubes de compresseur, et comparées avec les données expérimentales obtenues au LMFA (issues de travaux séparés). L'approche RANS surestime globalement le décollement de coin. Une amélioration significative est obtenue par la méthode LES, en particulier pour le coefficient de pression statique sur l'aube et les pertes de pression totale. Ces résultats montrent que la zone de décollement de coin, qui est la source principale des pertes, génère des tourbillons de grande échelle associés à de forts niveaux d'énergie. Les histogrammes bimodaux de la vitesse tangentielle qui ont été observés expérimentalement semblent confirmés par les résultats LES. En ce qui concerne les amplitudes des fluctuations de vitesse tangentielle, les résultats exprimentaux et ceux de la LES mettent en évidence deux pics sur certains profils perpendiculaires aux parois. Enfin, grâce à l'approche LES, les bilans de l'énergie cinétique turbulente sont calculés et analysés. Ils décrivent l'équilibre entre les termes de production, de dissipation et de transport. Une des perspectives de cette analyse est d'aider à améliorer la modélisation de la turbulence en approche RANS.

À l'allumage d'un lanceur spatial équipé de moteurs à propergol solide (MPS), une onde de souffle (ODS), caractérisée par un pic de grande amplitude et une large dépression, est générée pendant le transitoire de montée en pression. Puis le jet supersonique chaud rayonne un bruit de jet (BDJ) intense riche en composantes hautes fréquences. Cette thèse a été consacrée à la mise au point de méthodologies de simulation de ces phénomènes à l'aide de la plateforme de calcul CEDRE, en se comparant à des mesures expérimentales réalisées sur un petit MPS. Les phénomènes sont reproduits à l'aide de simulations aux grandes échelles (SGE). Les choix numériques (maillages, schémas) sont effectués à partir de cas-tests de validation. La méthodologie retenue pour l'ODS repose sur une approche directe simulant la génération et le transport de l'onde jusqu'aux capteurs expérimentaux. Ce calcul, qui estime aussi le transfert radiatif du jet, permet l'interprétation physique des phénomènes générateurs de l'ODS (effet piston du jet se développant). La prise en compte de la post combustion dans les premiers instants se révèle être un facteur de premier ordre concernant l'amplitude de l'ODS et son interaction avec le jet. Une approche hybride est adoptée pour le BDJ, chaînant une SGE du champ proche à la résolution des équations de Ffowcs Williams & Hawkings (FWH). Le bon accord du champ aérodynamique avec les mesures de la littérature incite à appliquer cette méthode au calcul de l'ODS. Finalement, celle-ci met en évidence le caractère fortement non-linéaire de l'ODS qui ne peut donc être rayonnée par FWH, et une interaction modèrée avec le jet, plus conforme aux observations expérimentales.

La thèse porte sur la simulation de la thermohydraulique et des transferts thermiques dans un faisceau d'aiguilles d'assemblage combustible de réacteur à neutrons rapides à caloporteur sodium.Des premiers calculs ont été réalisés par une approche moyennée de type RANS à l'aide du code industriel STAR-CCM+. De cette modélisation, il ressort une meilleure compréhension des transferts de chaleur opérés entre les aiguilles et le sodium. Les principales grandeurs macroscopiques de l'écoulement sont en accord avec les corrélations. Cependant, afin d'obtenir une description détaillée des fluctuations de température au niveau des fils espaceur, une approche plus détaillée de type LES et DNS est apparue indispensable. Pour la partie LES, le code TRIO_U a été utilisé. Concernant la partie DNS, un code de recherche a été utilisé. Ces approches requièrent des temps de calculs considérables qui ont nécessité des géométries représentatives mais simplifiées.L'approche DNS permet d'étudier l'écoulement à bas nombre de Prandtl, qui induit un comportement très différent du champ thermique relativement au champ hydraulique. Le calcul LES de l'assemblage montre que la présence du fil espaceur génère l'apparition de points chauds locaux (~20°C) en aval de celui-ci par rapport à l'écoulement sodium, au niveau de son contact avec l'aiguille. Les fluctuations de température au niveau des fils espaceur sont faibles (~1°C-2°C). En régime nominal, l'analyse spectrale montre l'absence de grande amplitude d'oscillations de température à basse fréquence (2-10 Hz); les conséquences sur la tenue mécanique des structures devront être analysées. / The thesis focuses on the numerical simulation of sodium flow in wire wrapped sub-assembly of Sodium-cooled Fast Reactor (SFR).First calculations were carried out by a time averaging approach called RANS (Reynolds- Averaged Navier-Stokes equations) using industrial code STAR-CCM+. This study gives a clear understanding of heat transfer between the fuel pin and sodium. The main variables of the macroscopic flow are in agreement with correlations used hitherto. However, to obtain a detailed description of temperature fluctuations around the spacer wire, more accurate approaches like LES (Large Eddy Simulation) and DNS (Direct Numerical Simulation) are clearly needed. For LES approach, the code TRIO_U was used and for the DNS approach, a research code was used. These approaches require a considerable long calculation time which leads to the need of representative but simplified geometry.The DNS approach enables us to study the thermal hydraulics of sodium that has very low Prandtl number inducing a very different behavior of thermal field in comparison to the hydraulic field. The LES approach is used to study the local region of sub-assembly. This study shows that spacer wire generates the local hot spots (~20°C) on the wake side of spacer wire with respect to the sodium flow at the region of contact with the fuel pin. Temperature fluctuations around the spacer wire are low (~1-2°C). Under nominal operation, the spectral analysis shows the absence of any dominant peak for temperature oscillations at low frequency (2-10Hz). The obtained spectra of temperature oscillations can be used as an input for further mechanical studies to determine its impact on the solid structures.

Un mascaret est une vague remontant contre le courant d’un fleuve lorsque la marée se propage dans un estuaire. À son passage, le mascaret induit une forte turbulence et un fort mélange dont les effets sur la vie de l’estuaire sont encore mal quantifiés. Ici, le phénomène est étudié expérimentalement et numériquement en utilisant un modèle d’onde positive se propageant contre un courant permanent.L’étude en laboratoire a permis de mesurer les variations de la surface libre, de la vitesse de l’écoulement ainsi que des échelles de turbulence. Lors de son passage, des fluctuations importantes de la surface libre et de la vitesse de l’écoulement sont observées, ainsi que des variations des échelles de turbulences. Des structures turbulentes semblent se former près du fond sous le front de l’onde et montent dans la colonne d’eau après le passage du front.La simulation numérique fut réalisée à partir de données expérimentales d’onde positive ondulée sur fond lisse. Une validation des méthodes numériques a été réalisée pour différente configuration. Les résultats des simulations d’onde positives donnent une cartographie détaillée de l’écoulement dans tout le canal. De plus, la simulation a permis d’identifier une inversion de la vitesse près des parois lors du passage des crêtes des ondes générant dans certaines configurations des structures turbulentes. / Tidal bores are surge waves propagating upstream rivers as the tide rushes into estuaries. They induce large turbulences and mixing of the river and estuary flow of which effects remain scarcely studied. Herein, tidal bores are investigated experimentally and numerically with an idealised model of positive surges propagating upstream an initially steady flow. The experimental work estimated flow changes and typical turbulent length scale evolution induced by undular bores with and without breaking roller. The bore passage was associated with large free surface and flow velocity fluctuations, together with some variations of the integral turbulent scales. Coherent turbulent structures appeared in the wake of leading wave near the bed and moved upward into the water column during the bore propagation. The numerical simulations were based on previous experimental work on undular bores. Some test cases were realised to verify the accuracy of the numerical methods. The results gave access to the detailed flow evolution during the bore propagation. Large velocity reversals were observed close to the no-slip boundaries. In some configurations, coherent turbulent structures appeared against the walls in the wake of the bore front.

Le bruit de cavité est un phénomène très fréquent dans le domaine des transports aériens.Il survient notamment lors de l’approche à l’atterrissage, où des interactions entre la cellule de l’aéronef et l’écoulement sont à l’origine de fortes émissions tonales. Il devient dès lors une source de pollution acoustique non-négligeable pour les populations résidant à proximité de zones aéroportuaires. Les études numériques et expérimentales décrites jusqu’à présent dans la littérature abordent essentiellement le cas des cavités rectangulaires. Pourtant, les cavités rencontrées en pratique dans l’industrie aéronautique impliquent des géométries souvent plus complexes. Lorsque ces cavités sont soumises à une excitation de nature aérodynamique, leur spécificité géométrique conduit le plus souvent à des réponses acoustiques assez éloignées des estimations issues de modèles académiques construits sur l’observation de cavités rectangulaires. Quelques travaux seulement abordent le cas des cavités cylindriques.Ce travail est consacré à l’étude aéroacoustique des cavités cylindriques, à l’initiative d’Airbus. Il s’inscrit dans le cadre du projet AEROCAV soutenu par la Fondation de Recherche pour l’Aéronautique & l’espace (FRAE). Son objectif est de déterminer les mécanismes impliqués dans les émissions acoustiques intenses et tonales pour les configurations étudiées.Une première partie présente les résultats expérimentaux issus des campagnes de mesures menées dans la soufflerie anéchoïque du Centre Acoustique du LMFA et de l’école Centrale de Lyon. Un modèle semi-empirique, reposant sur l’hypothèse d’une résonance acoustique pilotée par les instabilités présentes dans la couche de cisaillement à l’ouverture de la cavité,est construit à partir du modèle d’Elder (1978). Le modèle permet d’estimer les fréquences susceptibles de dominer l’acoustique rayonnée en champ lointain à partir de la donnée du champ moyen de vitesse longitudinale, que l’on mesure dans le plan de l’écoulement par Vélocimétrie par Imagerie des Particules (PIV).Une seconde partie est destinée au calcul direct du bruit rayonné par un écoulement laminaire ou turbulent affleurant une cavité cylindrique de référence. Il consiste à calculer le champ acoustique directement à partir de la résolution des équations tridimensionnelles de la mécanique des fluides. Le solver Alesia est présenté dans une version modifiée et adaptée à la mise en oeuvre d’une approche multidomaine d’ordre élevé faisant intervenir plusieurs maillages se recouvrant. Des techniques d’interpolation sont spécifiquement développées en vue d’assurer une communication bidirectionnelle entre les différents maillages, malgré des contraintes géométriques fortes. Un modèle d’excitation de l’écoulement est aussi développé afin de disposer de fluctuations dans l’écoulement incident, pour le cas turbulent. Ces deux points font l’originalité des calculs réalisés.Les simulations, menées sur une cavité de rapport d’aspect géométrique égal à 1 et soumise à un écoulement incident à Mach 0.2, montrent que le rayonnement acoustique peut être fidèlement reproduit numériquement. La couche de cisaillement est caractérisée par la présence de deux larges structures tourbillonnaires s’amplifiant lors de leur convection. Leur présence s’accompagne de fortes fluctuations de vitesse à l’origine d’un débit aérodynamique de fluide à l’ouverture qui excite la cavité acoustiquement. Une résonance forcée s’établit dans celle-ci, excitant la couche de mélange au voisinage du point de séparation. Ce couplage auto-entretenu est à l’origine du rayonnement acoustique intense et fortement tonal de la cavité. Il s’établit à une fréquence proche de la fréquence prédite par le modèle semi-empirique développé. / Cavity noise is a very frequent phenomenon in air transport. It occurs in particular during landing approaches, where airframe–flow interactions are responsible for strong tonal emissions. Accordingly, it turns to be a non negligible source of acoustic pollution for populations living near airport areas. Numerical and experimental studies reported in the literature tackle essentially the case of rectangular cavities. Nevertheless, cavities may often exhibit more complex shapes in practice. When subject to aerodynamic excitations, and because of their geometrical specificity, these cavities may have acoustic responses which can be rather far from estimations resulting from academic models designed for rectangular cavities. Only asmall number of studies tackle the case of cylindrical cavities.The present work requested by Airbus is dedicated to the study of aeroacoustics in cylindrical cavities. This work was been supported by the Fondation de Recherche pour l’Aéronautique& l’Espace (FRAE) under contract reference AEROCAV. It aims at discerning the mechanisms responsible for strong and tonal acoustic emissions for the studied configurations.Experimental data resulting from measurements performed in the anechoic wind-tunnel of the Centre Acoustique at ´Ecole Centrale de Lyon are presented in a first part. A semi-empirical model based on the hypothesis of a shear-layer driven acoustic resonance is constructed from the Elder model (1978). The model provides an estimation of the frequences which are likely to be predominant in the far-field acoustics, given the mean streamwise velocity field, currently measured in the flow plane by Particle Image Velocimetry (PIV).A second part deals with the direct computation of the noise radiated by a laminar or turbulent grazing flow over a standard cylindrical cavity. The method consists in the calculationof the acoustic field directly from the resolution of the tridimensional Navier–Stokes equations. The Alesia solver is presented in a modified form, adapted to the implementationof a high-order chimera method involving several overlapping grids. Interpolation techniques have been specifically developed to achieve a bidirectional communication between the meshes in spite of strong geometrical constraints. A flow excitation model has also been constructed in order to obtain fluctuations into the incoming flow in the turbulent case. These two last points make the present computations original. The simulations, which are performed on a cavity of geometric ratio taken as 1 and subject to a grazing flow of Mach 0.2, reveal that it is possible to retrieve the radiated noise numerically with high fidelity. They indicate the presence of two large amplifying vortices in the shearlayer. These vortices go with strong velocity fluctuations giving rise to an inflow of fluid at the cavity mouth which excites the cavity acoustically. A forced acoustic resonance occurs into the cavity, then destabilises the shear layer near the separation point. This self-sustained coupling is responsible for strong tonal radiations from the cavity. The frequency of the radiated noise is close to the one predicted by the semi-empirical model.

Lors d’un lancement spatial, le bruit des jets supersoniques chauds, générés par les moteurs-fusées au décollage et en interaction avec le pas de tir, est dommageable pour le lanceur et en particulier sa charge utile. Par conséquent, les acteurs du spatial cherchent à renforcer leur compréhension et leur maîtrise de cette ambiance acoustique, entre autres grâce à des méthodes et outils numériques. Toutefois, ils ne disposent pas d’une approche numérique globale capable de prendre en compte simultanément la génération fidèle du bruit, la propagation acoustique non-linéaire, les effets d’installation complexes et les géométries réalistes, pourtant inhérents aux applications spatiales. Dans cette optique, cette étude consiste à mettre en place et valider une méthodologie de simulation numérique par couplage fort Navier-Stokes − Euler, puis à l’appliquer à des cas réalistes de bruit de jet supersonique. L’objectif est d’affiner les capacités de prévision et de contribuer à la compréhension des mécanismes de génération de bruit dans de tels jets. Le solveur Navier-Stokes repose sur une méthode LES sur maillage non-structuré et le solveur acoustique sur une méthode de Galerkine discontinue d’ordre élevé sur maillage non-structuré. La méthodologie est tout d’abord évaluée sur des cas académiques visant à valider la simulation par couplage fort. Après des calculs préliminaires, la méthodologie est appliquée à la simulation du bruit d’un jet libre supersonique à Mach 3.1. Une méthode de déclenchement géométrique de la turbulence est implémentée sous la forme d’une marche à la paroi de la tuyère. La simulation aboutit à des estimations du bruit très proches des mesures réalisées au banc MARTEL et met en évidence des effets non-linéaires significatifs ainsi qu’un mécanisme singulier de rayonnement des ondes de Mach. Dans une démarche de progression vers des cas toujours plus réalistes, l’ensemble de l’approche numérique est finalement adaptée avec succès à la simulation du bruit d’un jet en présence d’un carneau. À terme, elle pourra être étendue à des configurations multi-jets réactifs, avec injection d’eau, voire à l’échelle 1. / During a space launch, the noise from hot supersonic jets, generated by rocket engines at liftoff and interacting with the launch pad, is harmful to the launcher and in particular its payload. Consequently, space actors are seeking to strengthen their understanding and control of this acoustic environment through numerical methods and tools, among the others. However, they do not dispose of a comprehensive numerical strategy that can simultaneously take into account accurate noise generation, nonlinear acoustic propagation, complex installation effects and realistic geometries, which are inherent to space applications. For this purpose, the present study consists in setting up and validating a numerical simulation methodology using a Navier-Stokes − Euler two-way coupling approach, then applying it to realistic cases of supersonic jet noise in order to improve prediction capabilities and contribute to the understanding of the noise generation mechanisms in such jets. The Navier-Stokes solver is based on an LES method on unstructured mesh and the acoustic solver on a high-order discontinuous Galerkin method on unstructured mesh. The methodology is first assessed on academic cases to validate the use of the two-way coupling. After preliminary computations, the methodology is applied to the simulation of the noise from a supersonic free jet at Mach 3.1. A geometric turbulence tripping method is implemented via a step at the nozzle wall. The computation leads to noise predictions very close to the experimental measurements performed at the MARTEL test bench and highlights significant nonlinear effects as well as a quite particular Mach waves radiation mechanism. Targeting even more realistic cases, the entire numerical approach is finally successfully adapted to the simulation of the noise from a supersonic jet configuration including a flame trench. In the future, it may be extended to configurations with clustered reactive jets, water injection devices or even at full scale.

Ces 50 dernières années, la majorité des paramètres de conception des moteurs cryotechniques ont été ajustés en l'absence d'une compréhension détaillée de la dynamique de flamme, en raison des limites des diagnostiques expérimentaux et des capacités de calcul. L'objectif de cette thèse est de réaliser des simulations numériques instationnaires d'écoulements réactifs transcritiques de haute fidélité, pour permettre une meilleure compréhension de la dynamique de flamme dans les moteurs cryotechniques et finalement guider leur amélioration. Dans un premier temps, la thermodynamique gaz-réel et son impact sur les schémas numériques sont présentés. Comme la Simulation aux Grandes Echelles (SGE) comporte des équations filtrées, les effets de filtrages induits par la thermodynamique gaz-réel sont ensuite mis en évidence dans une configuration transcritique type et un opérateur de diffusion artificiel, spécifique au gaz réel, est proposé pour lisser les gradients transcritiques en SGE. Dans un deuxième temps, une étude fondamentale du mélange turbulent et de la combustion dans la zone proche-injecteur des moteurs cryotechniques est menée grâce à la Simulation Numérique Directe (SND). Dans le cas non-réactif, les lâchers tourbillonnaires dans le sillage de la lèvre de l'injecteur jouent un rôle majeur dans le mélange turbulent et provoquent la formation de structures en peigne déjà observées expérimentalement dans des conditions similaires. Dans le cas réactif, la flamme reste attachée à la lèvre de l'injecteur, sans extinction locale, et les structures en peigne disparaissent. La structure de flamme est analysée et différents modes de combustion sont identifiés. Enfin, une étude de flamme-jet transcritique H2/O2, accrochée à un injecteur coaxial avec et sans retrait interne, est menée. Les résultats numériques sont d'abord validés par des données expérimentales pour l'injecteur sans retrait. Ensuite, la configuration avec retrait est comparée à la solution de référence sans retrait et à des données experimentales pour observer les effets de ce paramètre de conception sur l'efficacité de combustion.

Cette thèse présente la Simulation des Grandes Echelles (LES) de l’injection, de la pulvérisation et de la cavitation dans un injecteur pour les applications liées aux moteurs à combustion interne. Pour la modélisation de l’atomisation, on utilise le modèle ELSA (Eulerian Lagrangian Spray Atomization). Le modèle résout la fraction volumique du combustible liquide ainsi que la densité de surface d’interface liquide-gaz pour décrire le processus complet d’atomisation. Dans cette thèse, l’écoulement à l’intérieur de l’injecteur est également pris en compte pour une étude ultérieure de l’atomisation. L’étude présente l’application du modèle ELSA à un injecteur Diesel typique, à la fois dans le contexte de RANS et de LES.Le modèle est validé à l’aide de données expérimentales disponibles dans Engine Combustion Network (ECN). Le modèle ELSA, qui est normalement conçu pour les interfaces diffuses (non résolues), lorsque l’emplacement exact de l’interface liquide-gaz n’est pas pris en compte, est étendu pour fonctionner avec une formulation de type Volume of Fluid (VOF) de flux à deux phases, où l’interface est explicitement résolu. Le couplage est réalisé à l’aide de critères IRQ (Interface Resolution Quality), qui prennent en compte à la fois la courbure de l’interface et la quantité modélisée de la surface de l’interface. Le modèle ELSA est développé en premier lieu en considérant les deux phases comme incompressibles. L’extension à la phase compressible est également brièvement étudiée dans cette thèse. Il en résulte une formulation ELSA compressible qui prend en compte la densité variable de chaque phase. En collaboration avec l’Imperial College de Londres, la formulation de la fonction de densité de probabilité (PDF) avec les champs stochastiques est également explorée afin d’étudier l’atomisation. Dans les systèmes d’injection de carburant modernes, la pression locale à l’intérieur de l’injecteur tombe souvent en dessous de la pression de saturation en vapeur du carburant, ce qui entraîne une cavitation. La cavitation affecte le flux externe et la formulation du spray. Ainsi, une procédure est nécessaire pour étudier le changement de phase ainsi que la formulation du jet en utilisant une configuration numérique unique et cohérente. Une méthode qui couple le changement de phase à l’intérieur de l’injecteur à la pulvérisation externe du jet est développée dans cette thèse. Ceci est réalisé en utilisant le volume de formulation de fluide où l’interface est considérée entre le liquide et le gaz; le gaz est composé à la fois de vapeur et d’airambiant non condensable. / This thesis presents Large Eddy Simulation (LES) of fuel injection, atomization and cavitation inside the fuel injector for applications related to internal combustion engines. For atomization modeling, Eulerian Lagrangian Spray Atomization (ELSA) model is used. The model solves for volume fraction of liquid fuel as well as liquid-gas interface surface density to describe the complete atomization process. In this thesis, flow inside the injector is also considered for subsequent study of atomization. The study presents the application of ELSA model to a typical diesel injector, both in the context of RANS and LES. The model is validated with the help of experimental data available from Engine Combustion Network (ECN). The ELSA model which is normally designed for diffused (unresolved) interfaces, where the exact location of the liquid-gas interface is not considered, is extended to work with Volume of Fluid (VOF) type formulation of two phase flow, where interface is explicitly resolved. The coupling is achieved with the help of Interface Resolution Quality (IRQ) criteria, that takes into account both the interface curvature and modeled amount of interface surface. ELSA model is developed first considering both phases as incompressible, the extension to compressible phase is also briefly studied in this thesis, resulting in compressible ELSA formulation that takes into account varying density in each phase. In collaboration with Imperial College London, the Probability Density Function (PDF) formulation with Stochastic Fields is also explored to study atomization. In modern fuel injection systems, quite oftenthe local pressure inside the injector falls below the vapor saturation pressure of the fuel, resulting in cavitation. Cavitation effects the external flow and spray formulation. Thus, a procedure is required to study the phase change as well as jet formulation using a single and consistent numerical setup. A method is developed in this thesis that couples the phase change inside the injector to the external jet atomization. This is achieved using the volume of fluid formulation where the interface is considered between liquid and gas; gas consists of both the vapor and non condensible ambient air.