La politique européenne recommande aux industriels de l'aéronautique de concevoir des aéronefs capables de proposer une réduction de 50% du bruit perçu, de 50% de la consommation et de 80% des émissions de gaz NOx d'ici 2020. Les ingénieurs sont donc amenés à envisager des solutions innovantes basées sur le contrôle des écoulements, en particulier celui des décollements puisqu'ils ont généralement des effets délétères sur les performances aérodynamiques. Les travaux présentés dans ce mémoire sont dédiés au contrôle en boucle fermée du décollement d'une rampe arrondie par un jet synthétique réglable en fréquence et en amplitude. Le mémoire s'articule autour de deux axes principaux. Le premier axe a permis d'étudier la réponse du décollement à une perturbation périodique par simulations numériques de type LES afin de poursuivre la compréhension de la phy- sique impliquée dans le contrôle et de définir des objectifs de contrôle. Les fréquences et amplitudes optimales au regard de différents critères ont ainsi été identifiées et analysées. Le second axe a été consacré au contrôle en boucle fermée de la réponse quasi-statique du système par l'application de stratégies de régulations PID, extremums-seekings optimisés et logique floue. Leur conception a nécessité l'identification d'un modèle boîte-noire non-linéaire NARX. Les algorithmes non-linéaires, en particulier celui par logique floue, se sont montrés particulièrement performants. La démarche a été validée grâce à une application en temps réel d'une régulation extremum-seeking sur simulation RANS instationnaire.

The primary requirements for a modern industrial gas turbine consist of a continuous trend of an increasing efficiency combined with very low emissions in a robust, cost-effective manner. To fulfil these tasks a high turbine inlet temperature together with advanced dry low NOX combustion chambers are employed. These dry low NOX combustion chambers generate a rather flat temperature profile compared to previous generation gas turbines, which have a rather parabolic temperature profile before the nozzle guide vane. This means that the nozzle guide vane endwall heat load for modern gas turbines is much higher compared to previous generation gas turbines. Therefore the prediction of the nozzle guide vane flow field and endwall heat transfer is crucial for the engineering task of the design layout of the vane endwall cooling system. The present study is directed towards establishing new in-depth aerodynamic and endwall heat transfer knowledge for an advanced nozzle guide vane of a modern industrial gas turbine. To reach this objective the physical processes and effects which cause the different flow fields and the endwall heat transfer pattern in a baseline configuration, a combustion chamber variant, a heat shield variant without and with additional cooling air and a cavity variant without and with additional cooling air have been investigated. The variants, which differ from the simplified baseline configuration, apply design elements which are commonly used in real modern gas turbines. This research area is crucial for the nozzle guide vane endwall heat transfer, especially for the advanced design of the nozzle guide vane of a modern industrial gas turbine and has so far hardly been investigated in the open literature. For the experimental aerodynamic and endwall heat transfer research of the baseline configuration of the advanced nozzle guide vane geometry a new low pressure, low temperature test facility has been developed, designed and constructed, since no experimental heat transfer data exist in the open literature for this type of vane configuration. The new test rig consists of a linear cascade with the baseline configuration of the advanced nozzle guide vane geometry with four upscaled airfoils and three flow passages. For the aerodynamic tests the two middle airfoils and the hub and the tip endwall are instrumented with pressure taps to monitor the Mach number distribution. For the heat transfer tests the temperature distribution on the hub endwall is measured via thermography. The analysis of these measurements, including comparisons to research in the open literature shows that the new test rig generates accurate and reproducible results which give confidence that it is a reliable tool for the experimental aerodynamic and heat transfer research on the advanced nozzle guide vane of a modern industrial gas turbine. Previous own research work together with the numerical analysis performed in another part of the project as well as conclusions from a detailed literature study lead to the conclusion that advanced Navier-Stokes CFD tools with the v2-f turbulence model are most suitable for the calculation of the flow field and the endwall heat transfer of turbine vanes and blades. Therefore this numerical tool, validated against different vane and blade geometries and for different flow conditions, has been chosen for the numerical aerodynamic and endwall heat transfer research of the advanced nozzle guide vane of a modern industrial gas turbine. The evaluation of the numerical and experimental investigations of the baseline configuration of the advanced design of a nozzle guide vane shows the flow field of an advanced mid-loaded airfoil design with the features to reduce total airfoil losses. For the hub endwall of the baseline configuration of the advanced design of a nozzle guide vane the flow characteristics and heat transfer features of the classical vane endwall secondary flow model can be detected with a very weak intensity and geometric extension compared to the studies of less advanced vane geometries in the open literature. A detailed analysis of the numerical simulations and the experimental data showed very good qualitative and quantitative agreement for the three-dimensional flow field and the endwall heat transfer. These findings, together with the evaluations obtained from the open literature, lead to the conclusions that selected CFD software Fluent together with the applied v2-f turbulence model exhibits a high level of general applicability and is not tuned to a special vane or blade geometry. Therefore the CFD code Fluent with the v2-f turbulence model has been selected for the research of the influence of the several geometric variants of the baseline configuration on the flow field and the hub endwall heat transfer of the advanced nozzle guide vane of a modern industrial gas turbine. Most of the vane endwall heat transfer research in the open literature has been carried out only for baseline configurations of the flow path between combustion chamber and nozzle guide vane. Such a simplified geometry consists of a long, planar undisturbed approach length upstream of the nozzle guide vane. The design of real modern industrial gas turbines however requires often significant variations from this baseline configuration consisting of air-cooled heat shields and purged cavities between the combustion chamber and the nozzle guide vane. A detailed evaluation of the flow field and the endwall heat transfer shows major differences between the baseline and the heat shield configuration. The heat shield in front of the airfoil of the nozzle guide vane influences the secondary flow field and the endwall heat transfer pattern strongly. Additional cooling air, released under the heat shield has a distinctive influence as well. Also the cavity between the combustion chamber and the nozzle guide vane affects the secondary flow field and the endwall heat transfer pattern. Here the influence of additional cavity cooling air is more decisive. The results of the detailed studies of the geometric variants are applied to formulate guidelines for an optimized design of the flow path between the combustion chamber and the nozzle guide vane and the nozzle guide vane endwall cooling configuration of next-generation industrial gas turbines. / QC 20100917

Cette thèse étudie essentiellement des phénomènes touchant aux domaines aéronautique et spatial. Elle traite par simulations numériques des écoulements dans les tuyères de moteur-fusée, les écoulements d'arrière-corps d'avions, et les écoulements supersoniques en présence d'interactions onde de choc/couche limite. La motivation principale de ce travail a été de comprendre les différents facteurs qui gouvernent ces écoulements. Ce travail, soutenu par le groupe SNECMA, est consacré plus particulièrement à l'étude et à la modélisation des phénomènes en proche paroi en utilisant le code industriel N3S-Natur. L'objectif visé est d'améliorer les chaînes de conception et de développement industriels (RANS) et de développer de nouvelles méthodes de calcul (LES/DES) pour une meilleure maîtrise des systèmes énergétiques. Ainsi, les modèles mis au point dans cette étude ont été d'abord validés sur une variété d'écoulement simple (couche limite compressible, jet supersonique,...) avant d'être appliqués aux écoulements complexes (tuyères avec film de refroidissement, arrière-corps 3D, interaction visqueuse et décollement tridimensionnel). Une partie de l‘étude a été consacrée aux phénomènes transitoires d'amorçage rapide des tuyères propulsives. En ce qui concerne les aspects instationnaires, il a été montré, à travers l'étude d'une interaction onde de choc/couche limite, que la DES (Detached Eddy Simulation), tout comme la LES (Large Eddy Simulation), constitue un outil pertinent, permettant une prédiction fine des caractéristiques moyennes et fluctuantes des écoulements supersoniques décollés. Les résultats obtenus dans ce travail confirment la portée et l'intérêt scientifiques des études en aérodynamique supersonique.

Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests.

This thesis deals with the near-field of confluent jets, which can be of interest in many engineering applications such as design of a ventilation supply device. The physical effect of interaction between multiple closely spaced jets is studied using experimental and numerical methods. The primary aim of this study is to explore a better understanding of flow and turbulence behavior of multiple interacting jets. The main goal is to gain an insight into the confluence of jets occurring in the near-field of multiple interacting jets. The array of multiple interacting jets is studied when they are placed on a flat and a curved surface. To obtain the boundary conditions at the nozzle exits of the confluent jets on a curved surface, the results of numerical prediction of a cylindrical air supply device using two turbulence models (realizable 𝑘 − 𝜖 and Reynolds stress model) are validated with hot-wire anemometry (HWA) near different nozzles discharge in the array. A single round jet is then studied to find the appropriate turbulence models for the prediction of the three-dimensional flow field and to gain an understanding of the effect of the boundary conditions predicted at the nozzle inlet. In comparison with HWA measurements, the turbulence models with low Reynolds correction (𝑘 − 𝜖 and shear stress transport [SST] 𝑘 − 𝜔) give reasonable flow predictions for the single round jet with the prescribed inlet boundary conditions, while the transition models (𝑘 − 𝑘l − 𝜔𝜔 and transition SST 𝑘 − 𝜔) are unable to predict the flow in the turbulent region. The results of numerical prediction (low Reynolds SST 𝑘 − 𝜔 model) using the prescribed inlet boundary conditions agree well with the HWA measurement in the nearfield of confluent jets on a curved surface, except in the merging region. Instantaneous velocity measurements are performed by laser Doppler anemometry (LDA) and particle image velocimetry (PIV) in two different configurations, a single row of parallel coplanar jets and an inline array of jets on a flat surface. The results of LDA and PIV are compared, which exhibit good agreement except near the nozzle exits. The streamwise velocity profile of the jets in the initial region shows a saddle back shape with attenuated turbulence in the core region and two off-centered narrow peaks. When confluent jets issue from an array of closely spaced nozzles, they may converge, merge, and combine after a certain distance downstream of the nozzle edge. The deflection plays a salient role for the multiple interacting jets (except in the single row configuration), where all the jets are converged towards the center of the array. The jet position, such as central, side and corner jets, significantly influences the development features of the jets, such as velocity decay and lateral displacement. The flow field of confluent jets exhibits asymmetrical distributions of Reynolds stresses around the axis of the jets and highly anisotropic turbulence. The velocity decays slower in the combined regio  of confluent jets than a single jet. Using the response surface methodology, the correlations between characteristic points (merging and combined points) and the statistically significant terms of the three design factors (inlet velocity, spacing between the nozzles and diameter of the nozzles) are determined for the single row of coplanar parallel jets. The computational parametric study of the single row configuration shows that spacing has the greatest impact on the near-field characteristics.

La contribution majeure de cette thèse consiste en un formalisme et une méthodologie permettant de réaliser une analyse de stabilité globale des écoulements turbulents. La dynamique de ces écoulements est modélisée à l'aide des équations moyennées RANS, on s'intéresse ainsi à l'évolution des grandes échelles turbulentes. Un formalisme global est adopté permettant d'analyser des écoulements complexes. Une approche de type discrète est proposée, où les équations sont d'abord discrétisées puis linéarisées par différences finies. Cette approche permet d'adopter une stratégie générique vis à vis du système d'équations utilisées, comme le choix d'un modèle turbulent, et évite une linéarisation analytique fastidieuse des équations. Par ailleurs, cette méthode permet également l'utilisation systématique d'un code de simulation numérique afin de réaliser une étude de stabilité linéaire. Enfin, on démontre que l'analyse de la sensibilité à des perturbations stationnaires peut être réalisée grâce à ce formalisme et ce pour des écoulements laminaires et turbulents. Cette analyse détermine les zones où un contrôle stationnaire permettrait de réduire les instationnarités observées, facilitant la conception de stratégies efficaces de contrôle en boucle ouverte. La méthode est testée en premier lieu sur deux écoulements laminaires, où l'on reproduit les résultats obtenus par de précédentes études sur la dynamique d'oscillateur du sillage d'un cylindre bidimensionnel ainsi que sur la dynamique d'amplificateur de bruit d'une couche limite. La robustesse et la validité de notre méthode sont ensuite analysées sur un cas d'écoulement compressible et turbulent dans une cavité profonde. La précision des gradients de sensibilité est vérifiée, et la physique de l'écoulement, modes instables, propriétés acoustiques, impact de la modélisation de la turbulence, est détaillée. Afin de mieux appréhender la portabilité ainsi que la valeur ajoutée de notre méthode, on présentera ensuite trois cas d'études réalisées à l'aide de nos outils. On s'intéressera en premier lieu au phénomène de buffet sur un profil bidimensionnel, puis on présentera des résultats obtenus sur la caractérisation comme amplificateur de bruit d'un cas d'interaction de choc-couche limite, enfin une analyse du screech dans les jets sous détendus sera proposée. Enfin, on présente en dernier lieu une étude de la dynamique turbulente du sillage derrière un cylindre en forme de D.

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.

This thesis is concerned with the study of parallel mixing of two dissimilar gases under compressible conditions in the confined environment. A number of numerical studies are reported in the literature for the compressible mixing of two streams of gases where (1) both the streams are of similar gases at the same temperatures, (2) both the streams are at different temperatures with similar gases, and (3) dissimilar gases are with nearly equal temperatures. The combination of dissimilar gases at large temperature difference, mixing under compressible conditions, as in the case of scramjet propulsion, has not been adequately addressed numerically. Also many of the earlier studies have used two dimensional numerical simulation and showed good match with the experimental results on mixing layers that are inherently three dimensional in nature. In the present study, both two-dimensional (2-d) and three dimensional (3-d) studies are reported and in particular the effect of side wall on the three dimensionality of the flow field is analyzed, and the reasons of the good match of two dimensional simulations with experimental results have been discussed. Both two dimensional and three dimensional model free simulations have been conducted for a flow configuration on which experimental results are available. In this flow configuration, the mixing duct has a rectangular cross section with height to width ratio of 0.5. In the upper part of the duct hydrogen gas at a temperature of 103 K is injected through a single manifold of two Ludweig tubes and in the lower part of the duct nitrogen gas at a temperature of 2436 K is supplied through an expansion tube, both the gases are at Mach numbers of 3.1 and 4.0 respectively. Measurements in the experiment are limited to wall pressures and heat flux. The choice of this experimental condition gives an opportunity to study the effect of large temperature difference on the mixing of two dissimilar gases with large molecular weights under compressible conditions. Both two dimensional and three dimensional model free simulations are carried out using higher order numerical scheme (4th order spatial and 2nd order temporal) to understand the structure and evolution of supersonic confined mixing layer of similar and dissimilar gases. Two dimensional simulations are carried out by both SPARK (finite difference method) and OpenFOAM (finite volume method based open source software that was specially picked out and put together), while 3D model free simulations are carried out by OpenFOAM. A fine grid structure with higher grid resolution near the walls and shear layer is chosen. The effect of forcing of fluctuations on the inlet velocity shows no appreciable change in the fully developed turbulent region of the flow. The flow variables are averaged after the attainment of statistical steady state established through monitoring the concentration of inert species introduced in the initial guess. The effect of side wall on the flow structure on the mixing layer is studied by comparing the simulation results with and without side wall. Two dimensional simulations show a good match for the growth rate of shear layer and experimental wall pressures. Three dimensional simulations without side wall shows 14% higher growth rate of shear layer than that of two dimensional simulations. The wall pressures predicted by these three dimensional simulations are also lower than that predicted using two dimensional simulations (6%) and experimental (9%) results in the downstream direction of the mixing duct. Three dimensionality of the flow is thought of as a cause for these differences. Simulations with the presence of side wall show that there is no remarkable difference of three dimensionality of the flow in terms of the variables and turbulence statistics compared to the case without side walls. However, the growth rate of shear layer and wall surface pressures matches well with that predicted using two dimensional simulations. It has been argued that this good match in shear layer growth rate occurs due to formation of oblique disturbances in presence of side walls that are considered responsible for the decrease in growth rate in 3-d mixing layers. The wall pressure match is argued to be good because of hindrance from side wall in the distribution of momentum in third direction results in higher wall pressure. The effect of dissimilar gases at large temperature difference on the growth rate reduction in compressible conditions is studied. Taking experimental conditions as baseline case, simulations are carried out for a range of convective Mach numbers. Simulations are also carried out for the same range of convective Mach numbers considering the mixing of similar gases at the same temperature. The normalized growth rates with incompressible counterpart for both the cases show that the dissimilar gas combination with large temperature difference shows higher growth rate. This result confirms earlier stability analysis that predicts increased growth rate for such cases. The growth rate reduction of a compressible mixing layer is argued to occur due to reduced pressure strain term in the Reynolds stress equation. This reduction also requires the pressure and density fluctuation correlation to be very near to unity. This holds good for a mixing layer formed between two similar gases at same temperature. For dissimilar gases at different temperatures this assumption does not hold well, and pressure-density correlation coefficient shows departure from unity. Further analysis of temperature density correlation factor, and temperature fluctuations shows that the changes in density occur predominantly due to temperature effects, than due to pressure effects. The mechanism of density variations is found to be different for similar and dissimilar gases, while for similar gases the density variations are due to pressure variations. For dissimilar gases density variation is also affected by temperature variations in addition to pressure variations. It has been observed that the traditional k-ε turbulence model within the RANS (Reynolds Averaged Navier Stokes) framework fails to capture the growth rate reduction for compressible shear layers. The performance of k-ε turbulence model is tested for the mixing of dissimilar gases at large temperature difference. For the experimental test case the shear layer growth rate and wall pressures show good match with other model free simulations. Simulations are further carried out for a range of convective Mach numbers keeping the mixing gases and their temperatures same. It has been observed that a drop in the growth rate is well predicted by RANS simulations. Further, the compressibility option has been removed and it has been observed that for the density and temperature difference, even for incompressible case, the drop in growth rate exists. This behaviour shows that the decrease in growth rate is mainly due to the interaction of temperature and species mass fraction on density. Also it can be inferred that RANS with k-ε turbulence model is able to capture the compressible shear layer growth rate for dissimilar gases at high temperature difference. The mixing of heat and species is governed by the values of turbulent Prandtl and Schmidt numbers respectively. These numbers have been observed to vary for different flow conditions, while affecting the flow field considerable in the form of temperature and species distribution. Model free simulations are carried out on an incompressible convective Mach number mixing layer, and the results are compared with that of a compressible mixing layer to study the effect of compressibility on the values of turbulent Prandtl / Schmidt numbers. It has been observed that both turbulent Prandtl and Schmidt numbers show an almost constant value in the mixing layer region for incompressible case. While, for a compressible case, both turbulent Prandtl and Schmidt numbers show a continuous variation within the mixing layer. However, the turbulent Lewis number is observed to be near unity for both incompressible and compressible cases. The thesis is composed of 8 chapters. An introduction of the subject with critical and relevant literature survey is presented in chapter 1. Chapter 2 describes the mathematical formulation and assumptions along with solution methodology needed for the simulations. Chapter 3 deals with the two and three dimensional model free simulations of the non reacting mixing layer. The effect of the presence of side wall is studied in chapter 4. Chapter 5 deals with the effect of compressibility on the mixing of two dissimilar gases at largely different temperatures. The performance of k-ε turbulence model is checked for dissimilar gases in Chapter 6. Chapter 7 is concerned with the effect of compressibility on turbulent Prandtl and Schmidt numbers. Finally concluding remarks are presented in chapter 8. The main aim of this thesis is the exploration of parallel mixing of dissimilar gases under compressible conditions for both two and three dimensional cases. The outcome of the thesis is (a) a finding that the presence of sidewall in a mixing duct does not make flow field two dimensional, instead it causes the formation of oblique disturbances and the shear layer growth rate is reduced, (b) that it has been shown that the growth rates of dissimilar gases are affected far more by large temperature difference than by compressibility as in case of similar gases, (c) that the growth rates of compressible shear layers formed between dissimilar gases are better predicted using k-εturbulence model and (d) that for compressible mixing conditions the turbulent Prandtl and Schmidt numbers vary continuously in the mixing layer region necessitating the use of some kind of model instead of assuming constant values.

This thesis focuses on ingestion of foreign objects into standard turboprop engine GE H80 situated in aircraft Let L-410 Turbolet. Aim of this study is to create methodology of numerical simulation of particle movement inside the engine, which could be used during design process of Inertial Particle Separator device. Thesis consists of backward-facing step benchmark study which validates used methodology. Second part describes flow field calculation and numerical setup. The last part is dedicated to particle tracking analysis. Simulated trajectories are visually investigated, and coordinates of particle impacts at 1st rotor of a compressor are correlated to position of real observed damage.

A key component of research in the aerospace industry constitutes hypersonic flights (M>5) which includes the design of commercial high-speed aircrafts and development of rockets. Computational analysis becomes more important due to the difficulty in performing experiments and reliability of its results at these harsh operating conditions. There is an increasing demand from the industry for the accurate prediction of wall-shear and heat transfer with a low computational cost. Direct Numerical Simulations (DNS) create the standard for accuracy, but its practical usage is difficult and limited because of its high cost of computation. The usage of Reynold's Averaged Navier Stokes (RANS) simulations provide an affordable gateway for industry to capitalize its lower computational time for practical applications. However, the presence of existing RANS turbulence closure models and associated wall functions result in poor prediction of wall fluxes and inaccurate solutions in comparison with high fidelity DNS data. In recent years, machine learning emerged as a new approach for physical modeling. This thesis explores the potential of employing Machine Learning (ML) to improve the predictions of wall fluxes for hypersonic turbulent boundary layers. Fine-grid RANS simulations are used as training data to construct a suitable machine learning model to improve the solutions and predictions of wall quantities for coarser meshes. This strategy eliminates the usage of wall models and extends the range of applicability of grid sizes without a significant drop in accuracy of solutions. Random forest methodology coupled with a bagged aggregation algorithm helps in modeling a correction factor for the velocity gradient at the first grid points. The training data set for the ML model extracted from fine-grid RANS, includes neighbor cell information to address the memory effect of turbulence, and an optimal set of parameters to model the gradient correction factor. The successful demonstration of accurate predictions of wall-shear for coarse grids using this methodology, provides the confidence to build machine learning models to use DNS or high-fidelity modeling results as training data for reduced-order turbulence model development. This paves the way to integrate machine learning with RANS to produce accurate solutions with significantly lesser computational costs for hypersonic boundary layer problems.