In this thesis, the ribbed ducts of the internal cooling passage in turbine blading are investigated to demonstrate the effects of high speed rotation. Rotation coupled with high temperature operating conditions alters the mean flow, turbulence, and heat transfer augmentation due to Coriolis and centrifugal buoyancy forces that arises from density stratification in the domain. Gas turbine engines operate in particle laden environments (sand, volcanic ash), and particulate matter ingested by the engine can make their way into the blade internal cooling passages over thousands of operating hours. These particulates can deposit on the walls of these cooling passages and degrade performance of the turbine blade. Large-Eddy Simulations (LES) with temperature dependent properties is used for turbulent flow and heat transfer in the ribbed cooling passages and Lagrangian tracking is used to calculate the particle trajectories together with a wall deposition model. The conditions used are Re=100,000, Rotation number, Ro = 0.0 and 0.2, and centrifugal Buoyancy parameters of Bo=0, 0.5, and 1.0. First, the independent effects of Coriolis and centrifugal buoyancy forces are investigated, with a focus on the additional augmentation obtained in heat transfer with the addition of centrifugal buoyancy. Coriolis forces are known to augment heat transfer at the trailing wall and attenuate the same at the leading wall. Phenomenological arguments stated that centrifugal buoyancy augments the effects of Coriolis forces in outward flow in the first pass while opposing the effect of Coriolis forces during inward flow in the second pass. In this study, it was found that in the first pass, centrifugal buoyancy had a greater effect in augmenting heat transfer at the trailing wall than in attenuating heat transfer at the leading wall. On the contrary, it aided heat transfer in the second half of the first pass at the leading wall by energizing the flow near the wall. Also, contrary to phenomenological arguments, inclusion of centrifugal buoyancy augmented heat transfer over Coriolis forces alone on both the leading and trailing walls of the second pass. Sand ingestion is then investigated, by injecting 200,000 particles in the size range of 0.5-175μm with 65% of the particles below 10 μm. Three duct wall temperatures are investigated, 950, 1000 and 1050 °C with an inlet temperature of flow and particles at 527 °C . The impingement, deposition levels, and impact characteristics are recorded as the particles move through the domain. It was found that the Coriolis force greatly increases deposition. This was made prevalent in the first pass, as 84% of the deposits in the domain occurred in the first pass for the rotating case, whereas only 27% of deposits occurred in the first pass for the stationary case with the majority of deposits occurring in the bend region. This was due to an increased interaction with the trailing wall in the rotating case whereas particles in the stationary case were allowed to remain in the mean flow and gain momentum, making rebounding from a wall during collision more likely than deposition. In contrast, the variation of wall temperatures caused little to no change in deposition levels. This was concluded to be a result of the high Reynolds number used in the flow. At high Reynolds numbers, the particles have a short residence times in the internal cooling circuit not allowing the flow and particles to heat up to the wall temperature. Overall, 87% of the injected particles deposited in the rotating duct whereas 58% deposited in the stationary duct. / Master of Science

In this dissertation, various VAT parameters are investigated to determine the effect of the overall efficiency of the turbine at a high Reynolds number. To increase the efficiency of the vertical axis turbines, 2D CFD simulations are completed in an effort to better understand the physics behind the operation of these turbines. Specifically, the effect of advance ratio, solidity, and wake interactions were investigated. Simulations were completed in OpenFOAM using the k-ω SST turbulence model at a nominal Reynolds number of 500,000 using a NACA 0015 airfoil. To simulate the motion of the turbine, Arbitrary Mesh Interfacing (AMI) was used. For all of the parameters tested, it was found that the geometric effective angle of attack seen by the turbine blades had a significant impact on the power extracted from the flow. The range of effective angles of attack was found to decrease as the advance ratio increased. In spite of this, a severe loss in the power coefficient occurred at an advance ratio of 2.5 during which the blade experienced dynamic stall. This effect was also seen when the number of turbine blades was changed to four, at a solidity of 1.08. This negative impact on performance was found to be due to the increase in the drag component of the tangential force when dynamic stall occurs. Results indicate that wake interactions between subsequent blades have a large impact on performance especially when the wake interaction alters the flow direction sufficiently to create conditions for dynamic stall. To improve the performance of the VAT in the presence of dynamic stall, calculations were completed of a static twisted blade profile using GenIDLEST and OpenFOAM. There was found to be no improvement in the lift coefficient when comparing the twisted blade profile with a 2D blade at the same median angle of attack as the twisted blade. To further see the effects of the twisted blade, an effective VAT pitching motion was given to the blade and again compared to a 2D blade with the same motion. In this case there was significant improvement seen in the performance of the twisted blade. / Master of Science

L’application de cette thèse est le contrôle actif de la combustion dans les brûleurs industriels à combustible liquide. Il s’agit d’explorer les possibilités de contrôle d’un spray par des jets gazeux auxiliaires. Deux familles d’actionneurs utilisant ce procédé ont été testées sur un atomiseur coaxial assisté par air. Le premier dispositif est appelé (Dev). Composé d’un unique jet actionneur, il vise à dévier le spray. La seconde configuration, appelée (Sw), est équipée de 4 jets auxiliaires tangents au spray afin de lui conférer un effet de swirl et d’en augmenter le taux d’expansion. Les mesures de granulométrie par PDA et les visualisations du spray par strioscopie démontrent un effet important du contrôle sur l’atomisation et la forme du spray. On observe en outre une déviation pouvant atteindre 30°avec l’actionneur (Dev) et une augmentation du taux d’expansion de 80% dans le cas (Sw). Des simulations du banc expérimental ont de plus été menées avec le code AVBP. L’écoulement de gaz est calculé par simulation aux grandes échelles (SGE ou LES en Anglais). L’approche lagrangienne est utilisée pour simuler la phase dispersée. Une attention particulière a été portée aux conditions d’injection du gaz et des gouttes dans le calcul. Ceci a abouti au développement d’une nouvelle condition limite caractéristique non réfléchissante (VFCBC) destinée à l’injection d’écoulements turbulents en LES compressible. Les résultats de LES présentent un bon accord avec les mesures expérimentales. Les effets du contrôle sur la dynamique des gouttes et sur la topologie du spray (forme, déviation, expansion) sont correctement décrits. / The present work focuses on active control of two-phase combustion in industrial burners. The generic method explored in this thesis consists in controlling the injected fuel spray with transverse air jets. Two families of these jet actuators are tested on a coaxial airblast atomizer. The first system (Dev) is used to modify the trajectory of the spray, while the second one (Sw) introduces swirl into the spray to modify its spreading rate and mixing with the surrounding air. Experimental characterisations of the controlled flow with Schlieren visualisations and Phase Doppler Anemometry (PDA) show that actuators induce important effect on the spray. The deviation angle reaches 30° for the actuator (Dev) and the expansion rate increases of 80 % in the swirl case (Sw). Simulations of the experiment are then performed with the CFD code AVBP. The gas flow is computed with Large Eddy Simulation (LES). A Lagrangian formulation is used to simulate droplets trajectories. A particular attention is given to the injection of the gas flow and the droplets in the calculations. Therefore, a new non-reflecting characteristic boundary condition (VFCBC) has been derived to inject turbulent flows in compressible LES. A good agreement is observed between simulation and experiment. Control effects on the spray topology ( features, deviation, spread rate) and on the droplets velocities and diameters are correctly described by the Lagrangian LES.

Les écoulements décollés sont fortement instationnaires, l'objectif de cette thèse a été de localiser et d'identifier les phénomènes à la source de ces instationnarités et de comprendre les processus physiques permettant le transfert de l'information de ces zones sources au reste de l'écoulement. Pour ce faire une analyse des résultats issus de simulations numériques a été réalisée. En étudiant la corrélation et la cohérence entre les positions de choc et les fluctuations de pression, l'interaction a pu être séparée en plusieurs parties distinctes. A l'aide de la théorie des caractéristiques définissant les directions et les cinématiques de propagation de l'information, les liens spatio-temporels entre ces différentes régions ont pu être déterminés. Les résultats de ces études couplés avec ceux issus des expériences ont montré clairement que les phénomènes se produisant à l'intérieur de la zone de recirculation existant en aval du choc de décollement gouvernent la dynamique de la totalité de l'interaction, aussi bien à basse fréquence qu'à moyenne fréquence. Ainsi les mouvements de choc apparaissent comme le miroir des phénomènes se produisant à l'intérieur de la zone décollée. Une représentation équivalente en fluide non visqueux permettant une description du comportement instationnaire de l'interaction a aussi été proposée. / Separated flow are often strongly unsteady; the aim of this thesis is to localize and identify the sources of the unsteadiness and to understand the physical phenomena governing the information transfer from these source zones to the rest of the flow. To do this, data used for this analysis have been obtained from numerical simulations (LES). Both cross-correlation and coherence between shock motion and pressure fluctuations have shown that the interaction can be split in several distinct zones. The theory of characteristics is used to define the information paths and the propagation velocities, so that the space-time links between these regions have been determined. Both numerical and experimental studies have clearly shown that phenomena present within the recirculation buble govern the whole of the interaction, at low and intermediate frequencies. Indeed the shock motions appears as the mirror of phenomena present in the separated zone. An inviscid equivalent scenario has been proposed to represent the interaction.

The spectral-element method (SEM) is used to study wall-bounded turbulent flowsin moderately complex geometries. The first part of the thesis is devoted to simulations of canonical flow cases, such as temporal K-type transitionand turbulent channel flow, to investigate general resolution requirements and computational efficiency of the numerical code nek5000. Large-eddy simulation (LES) is further performed of a plane asymmetric diffuser flow with an opening angle of 8.5 degrees, featuring turbulent flow separation. Good agreement with numerical studies of Herbst (2007) is obtained, and it is concluded that the use of a high-order method is advantageous for flows featuring pressure-induced separation. Moreover, it is shown, both a priori on simpler model problems and a posteriori using the full Navier--Stokes equations, that the numerical instability associated with SEM at high Reynolds numbers is cured either by employing over-integration (dealiasing) or a filter-based stabilisation, thus rendering simulations of moderate to high Reynolds number flows possible. The second part of the thesis is devoted to the first direct numerical simulation (DNS) of a truly three-dimensional, turbulent and separated diffuser flow at Re = 10 000 (based on bulk velocity and inflow-duct height), experimentally investigated by Cherry et al. (2008). The massively parallel capabilities of the spectral-element method are exploited by running the simulations on up to 32 768 processors. Very good agreement with experimental mean flow data is obtained and it is thus shown that well-resolved simulations of complex turbulent flows with high accuracy are possible at realistic Reynolds numberseven in complicated geometries. An explanation for the discovered asymmetry of the mean separated flow is provided and itis demonstrated that a large-scale quasi-periodic motion is present in the diffuser. In addition, a new diagnostic measure, based on the maximum vorticity stretching component in every spatial point, is designed and tested in a number of turbulent and transitional flows. Finally, Koopman mode decomposition is performed of a minimal channel flow and compared to classical proper orthogonal decomposition (POD). / QC 20111206

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

Multiphase flows are a common phenomenon. Rains, sediment transport in rivers, snow and dust storms, mud slides and avalanches are examples of multiphase flows occurring in nature. Blood flow is an example of multiphase flow in the human body, which is of vital importance for survival. Multiphase flows occur widely in industrial applications from hydrocarbon extrac-tion to fuel combustion in engines, from spray painting to spray drying, evaporators, pumps and pneumatic conveying. Predicting multiphase flows is of vital importance to understand natural phenomenon and to design and improve industrial processes. Separated flows and dispersed flows are two types of multiphase flows, which occur together in many industrial applications. Physical features of these two classes are different and the transition from one to another involves complex flow physics. Experimental studies of multiphase flows are not easy, as most real world phenomenon cannot be scaled down to laboratory models. Even for those phenomenon that can be demonstrated at lab-oratory scale, rescaling to real world applications requires mathematical models. There are many challenges in experimental measurements of multiphase flows as well. Measurement techniques well suited for single phase flows have constraints when measuring multiphase phenomenon. Un-certainty in experimental measurements poses considerable difficulties in validating numerical models developed for predicting these flows. Owing to the computational effort required, direct simulation of multiphase flows, even for small scale real world applications is out of present scope. Numerical methods have been developed for dealing with each class of flow separately, that in-volves use of models for phenomenon that is computationally demanding. Reynolds Averaged Navier-Stokes (RANS) methods for predicting multiphase flows place strong requirements on turbulence models, as information about fluctuating quantities in the field, that have significant effects on dispersed phase, is not available. Large Eddy Simulation (LES) gives better predictions than RANS as the instantaneous field data is available and large scale unsteadiness that effects the dispersed phase can be captured. Recent LES studies of multiphase flows showed that the sub-grid-scale (SGS) model used for the continuous phase has an effect on the evolution of the dispersed phase. In this work, LES of multiphase flows is performed using Explicit Filtering Large Eddy Sim-ulation method. In this method, spatial derivatives are computed using higher order compact schemes that have spectral-like resolution. SGS modeling is provided by the use of a filter with smoothly falling transfer function. This method is mathematically consistent and converges to a DNS as the grid is refined. It has been successfully applied to combustion and aero-acoustics and this work is the first application of the method to multiphase flows. Study of dispersed multiphase flows was carried out in this work. Modeling of the dispersed phase is kept simple since the in-tention was to evaluate the capability of explicit filtering LES method in predicting multiphase flows. Continuous phase is solved using a compressible formulation with explicit filtering method. Spatial derivatives are computed using fourth and sixth order compact schemes that use derivative splitting method proposed by Hixon & Turkel (2000a) and second order Runge-Kutta (RK2) time stepping. The grid is stretched as needed. Non-reflecting boundary conditions due to Poinsot & Lele (1992) are used to avoid acoustic reflections from boundaries. Buffer zones (Bogey & Bailly (2002)) are employed at outflow and lateral boundaries to damp vortical structures. The code developed for continuous phase is evaluated by studying round jets at Re =36,000 and comparing with experimental measurements of Hussein et al. (1994) and Panchapakesan & Lumley (1993). Simulations showed excellent agreement with experimental results. Rate of decay of axial velocity and the evolution of turbulence intensities on the centerline matched very well with measurements. Radial profiles of mean and fluctuating components of velocities exhibit self-similarity. A set of studies were then performed using this code to assess the effect of numerical scheme, grid refinement & stretching and simulation times on the predictions. Results from these simulations showed good agreements with experiments and established the code for use in multiphase flows under various simulation conditions. To assess the prediction of multiphase flows using this LES method, an evaporating spray ex-periment by Chen et al. (2006) was simulated. The experiment uses a nebuliser for generating a finely atomized spray of acetone, which avoids complex breakdown phenomenon associated with air blast atomizers and provides well defined boundary conditions for model evaluation. The neb-uliser sits upstream in a pipe carrying air and droplets travel along with air for a distance of 10 diameters before exiting into a wind tunnel with co-flowing air. Droplet breakdown, if any, takes place inside the pipe and the spray is finely atomized by the time it reaches pipe exit. One of the experimental cases at Re =31,600, with a mass loading of 1.1% and a jet velocity of 56 m/s is simulated. Particle size has a χsquared distribution with a Sauter mean diameter of 18µm. In the self-similar region, decay of centerline velocity and turbulence intensities matched well with ex-perimental results. Continuous phase exhibits self-similar behavior. A series of simulations were then performed to match the initial region of the spray by altering the inflow conditions in the sim-ulation. Simulation that matched the breakdown location of the experiment revealed the presence of a relaxation zone with a higher initial spreading rate, followed by a lower asymptotic spreading rate. Studies were performed to understand the effect of various phenomenon like evaporation and droplet size on this behavior. A study of breakdown region of particle-laden jets was performed to understand the presence of relaxation zone post breakdown. Flow conditions were similar to evaporating spray experiment except that particles do not evaporate, mass loading is 2% and jet Reynolds number Re =2000. A series of grid refinements were performed and on the largest grid, gird spacing Δy =7.5η, where ηis an estimate of the Kolmogorov length scale based on flow conditions. Decay of axial velocity on the centerline showed variations with grid refinement, tending to the experimentally measured value as the grid is refined. Variation of turbulence intensities along the centerline revealed a jump in axial velocity fluctuations at the breakdown location, while radial and azimuthal velocities showed a smooth increase to their asymptotic value. This jump was resolved on grid refinement and on fine grids axial velocity fluctuations followed the other two quantities closely in their rise to asymptotic state. Comparison of these quantities with a jet without particles revealed that the flow features are same for a jet with and without particles, and at the mass loading studied, particles have negligible effect on jet breakdown. Another study performed at a higher Reynolds number of Re =11,000, under similar flow conditions showed similar behavior. To assess the ability of predicting dispersed phase, simulations of particle-laden flows at low Stokes number were performed and compared against an experiment by Lau & Nathan (2014). The experiment studies variation of velocity and particle concentration along the centerline, and half widths of a jet velocity and concentration. Particles are injected into a pipe along with air, and the two phase flow is fully developed by the time it exits the pipe into a wind tunnel along with a co-flow. Particles are mono-disperse with a density of 1200 kg/m3. Mass loading is 40% so that particles have a significant effect on the continuous phase. Two cases at particle Stokes number of 1.4, one with Re =10,000, bulk velocity of 12 m/s and particle diameter of 20µm and another with Re =22,500, bulk velocity of 36 m/s and particle diameter of 10µm were simulated. Simulations of both the cases showed good match with experimental measurements of centerline decay for the continuous phase. For the dispersed case, simulations with larger particles showed good match with experimental results, while smaller particles showed differences. This was understood to be the effect of lateral migration which is prominent in case of smaller particles, the models for which have not been used in the present simulation study.

Turbulence is an open and challenging problem for mathematical approaches, physical modeling and numerical simulations. Numerical solutions contribute significantly to the understand of the nature and effects of turbulence. The focus of this thesis is the development of appropriate numerical methods for the computer simulation of turbulent flows. Many of the existing approaches to turbulence utilize analogies from kinetic theory. Degond & Lemou (J. Math. Fluid Mech., 4, 257-284, 2002) derived a k-✏ type turbulence model completely from kinetic theoretic framework. In the first part of this thesis, a numerical method is developed for the computer simulation based on this model. The Boltzmann equation used in the model has an isotropic, relaxation collision operator. The relaxation time in the collision operator depends on the microscopic turbulent energy, making it difficult to construct an efficient numerical scheme. In order to achieve the desired numerical efficiency, an appropriate change of frame is applied. This introduces a stiff relaxation source term in the equations and the concept of asymptotic preserving schemes is then applied to tackle the stiffness. Some simple numerical tests are introduced to validate the new scheme. In the second part of this thesis, alternative approaches are sought for more efficient numerical techniques. The Lattice Boltzmann Relaxation Scheme (LBRS) is a novel method developed recently by Rohan Deshmukh and S.V. Raghuram Rao for simulating compressible flows. Two different approaches for the construction of implicit sub grid scale -like models as Implicit Large Eddy Simulation (ILES) methods, based on LBRS, are proposed and are tested for Burgers turbulence, or Burgulence. The test cases are solved over a largely varying Reynolds number, demonstrating the efficiency of this new ILES-LBRS approach. In the third part of the thesis, as an approach towards the extension of ILES-LBRS to incompressible flows, an artificial compressibility model of LBRS is proposed. The modified framework, LBRS-ACM is then tested for standard viscous incompressible flow test cases.

Ce mémoire de thèse traite de l’analyse, l’intensification et l’optimisation du transfert thermique convectif dans les échangeurs à plaques et ailettes à pas décalés utilisés notamment dans le domaine automobile comme refroidisseurs d’air de suralimentation. Deux approches complémentaires sont abordées dans cette étude : des simulations numériques visant à l’analyse fine locale des caractéristiques de l’écoulement et des mécanismes de transfert, et une modélisation de type nodale permettant une caractérisation globale des performances thermo-aérauliques. Sur la gamme de Reynolds considérée différentes modélisations de la turbulence sont mises en œuvre et comparées. Ainsi des simulations aux grandes échelles (LES) permettent de qualifier des simulations de type RANS classiquement utilisées jusque-là : de fortes différences tant au niveau structuration de l’écoulement qu’au niveau performances globales sont ainsi mises en évidence selon le régime d’écoulement considéré. La mise au point d’un modèle nodal est ensuite abordée dans le but de mener des optimisations de géométries d’échangeurs non-conventionnels à pas décalés. Les différents scénarii d’optimisation considérés montrent l’intérêt de cette approche autorisant l’évaluation d’un nombre élevé de configurations géométriques. Dans une dernière partie une nouvelle géométrie innovante permettant de générer des tourbillons longitudinaux sur ce type d’ailettes est proposée et étudiée. / This thesis deals with the analysis, intensification and optimization of convective heat transfer in offset strip fins (OSF) heat exchangers used, for example, in the automotive field as water-cooled charge air coolers. Two complementary approaches are carried out in this study: CFD simulations to perform local fine analysis of the flow characteristics and transfer mechanisms, and a nodal type modeling allowing calculation of global aerothermal performance. Over the range of Reynolds numbers considered, different turbulence modeling approaches are implemented and compared: Large Eddy Simulations (LES) and RANS simulations which are usually used. The qualification of the RANS models shows that strong differences, both in the flow structure and at the overall performance evaluation level, may beobserved, depending on the flow regime considered. Then the development of a nodal model is presented. It aims at carrying out rapid optimization of geometries of unconventional OSF heat exchangers. The various optimization scenarios considered show the interest of this approach allowing the evaluation of a large number of geometric configurations. In a last part, an innovative new geometry that generates longitudinal vortices on this type of fins is proposed and studied.

Development of modern safety standards for hydrogen storage infrastructure requires fundamental insight into the physics of buoyant gas dispersion into ambient air. Also, from a practical engineering stand-point, flow patterns and dispersion of gas originating from orifices in the side wall of circular pipe or storage tank need to be studied. In this thesis, novel configurations were considered to investigate the evolution of turbulent jets issuing from realistic pipeline geometries. First, the effect of jet densities and Reynolds numbers on vertical jets were investigated, as they emerged from the side wall of a circular pipe, through a round orifice. The resulting jet flow was thus issued through a curved surface from a source whose original velocity components were nearly perpendicular to the direction of the ensuing jets. Particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) techniques were employed simultaneously to provide instantaneous and time-averaged flow fields of velocity and concentration. The realistic flow arrangement resulted in an asymmetric flow pattern and a significant deflection from the vertical axis of jets. The deflection was influenced by buoyancy, where heavier gases deflected more than lighter gases. These realistic jets experienced faster velocity decay, and asymmetric jet spreading compared to round jets due to significant turbulent mixing in their near field. In addition to that, horizontal multi-component jets issuing from a round orifice on the side wall of a circular tube were also investigated experimentally by the means of simultaneous velocity and concentration measurements. A range of Reynolds numbers and gas densities were considered to study the effects of buoyancy and asymmetry on the resulting flow structure. The realistic pipeline jets were always exhibited an asymmetry structure and found to deflect about the jet's streamwise axis in the near field. In the far field, the buoyancy dominated much closer to the orifice than expected in the axisymmetric round jet due to the realistic leak geometry along with the pipeline orientation considered in this study. In general, significant differences were found between the centreline trajectory, spreading rate, and velocity decay of conventional horizontal round axisymmetric jets issuing through flat plates and the pipeline leak-representative jets considered in the present study. Finally, the dispersion of turbulent multi-component jets issuing from high-aspect-ratio slots on the side wall of a circular tube were studies experimentally by employing simultaneous PIV and PLIF techniques. Two transversal & longitudinal oblong geometries in respect to the longitudinal axes of the tube , and with an aspect ratio of 10 were considered in this study. Both horizontal and vertical orientations along with broad range of Reynolds numbers and gas densities were considered to investigate the effects of buoyancy and asymmetry on the resulting flow structure. The ensuing jets were found to deflect along the jet streamwise axis, once more, due to the realistic pipeline leak-representative configuration. It was also found that increases in aspect ratio of these realistic jets caused a reduction in the angle of deflection, jet centreline decay rates and the width growth on both velocity and scalar fields compared to their round jets counterparts, most notably in the far field. These findings indicate that conventional jets (those that are issuing through flat surfaces) assumptions are inadequate to predict gas concentration, entrainment rates and, consequently, the extent of the flammability envelope of realistic gas leaks. Thus, extreme caution is required when using conventional jet assumptions to describe the physics of a buoyant jet emitted from realistic geometries. / Graduate