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.

Flow-induced vibration caused by fully developed pipe flow has been recognized, but not fully investigated under turbulent conditions. This thesis focuses on the development of a numerical Fluid-Structure Interaction (FSI) model that will help define the relationship between pipe wall vibration and the physical characteristics of turbulent flow. Commercial FSI software packages are based on Reynolds Averaged Navier-Stokes (RANS) fluid models, which do not compute the instantaneous fluctuations in turbulent flow. This thesis presents an FSI approach based on Large Eddy Simulation (LES) flow models, which do compute the instantaneous fluctuations in turbulent flow. The results based on the LES models indicate that these fluctuations contribute to the pipe vibration. It is shown that there is a near quadratic relationship between the standard deviation of the pressure field on the pipe wall and the flow rate. It is also shown that a strong relationship between pipe vibration and flow rate exists. This research has a direct impact on the geothermal, nuclear, and other fluid transport industries.

High speed planing hulls are currently widely used for example in recreational and emergency vessel applications. However, very little CFD research has been done for planing vessels, especially for those with stepped hulls. A validated CFD method for planing stepped hulls could be a valuable improvement for the design phase of such hulls. In this thesis, a CFD method for stepped hulls, with a primary focus on two-step hulls, is developed using STAR-CCM+. As a secondary objective, porpoising instability of two-step hulls is investigated. The simulations are divided into two parts: In the first part a method is developed and validated with existing experimental and numerical data for a simple model scale planing hull with one step. In the second part the method is applied for two two-step hulls provided with Hydrolift AS. A maximum two degrees of freedom, trim and heave, are used, as well as RANS based k-w SST turbulence model and Volume of Fluid (VOF) as a free surface model. The results for the one-step hull mostly corresponded well with the validation data. For the two-step hulls, validation data did not exists and they were first simulated with a fixed trim and sinkage and compered between each other. In the simulations with free trim and heave both hulls experienced unstable porpoising behavior.

Early in the ship design process, naval architects must often evaluate and compare multiple hull forms for a specific set of requirements. Analytical tools are useful for quick comparisons, but they usually specialize in a specific hull type and are therefore not adequate for comparing dissimilar hull types. Scale model hydrodynamic testing is the traditional evaluation method, and is applicable to most hull forms. Scale model tests are usually performed on the largest model possible in order to achieve the most accurate performance predictions. However, such testing is very resource intensive, and is therefore not a cost effective method of evaluating multiple hull forms. This thesis explores the testing of small scale models. It is hypothesized that although the data acquired by these tests will not be accurate enough for performance predictions, they will be accurate enough to rank the performance of the multiple hull forms being evaluated.

L'objectif de la présente étude est d'analyser la modélisation de la turbulence de simulations en moyenne de Reynolds (RANS) dans le cadre d'écoulements de type turbomachines, en utilisant des simulations aux grandes échelles (SGE) comme référence. L'étude porte sur deux cas test: un décollement de coin dans une grille d'aubes rectiligne, et un écoulement de jeu pour un aubage isolé dans un jet. Deux lois de comportement, la loi de comportement de Boussinesq et la loi de comportement quadratique (quadratic constitutive relation ou QCR), sont analysées, avec deux versions du modèle de turbulence k-omega de Wilcox. Les lois de comportement étudiées reposent sur deux hypothèses: une hypothèse d'alignement entre le tenseur de Reynolds et un tenseur construit à partir de l'écoulement moyen, et une hypothèse sur la viscosité turbulente. L'hypothèse d'alignement est étudiée à partir de la SGE, pour laquelle les deux tenseurs sont indépendamment connus, en utilisant un indicateur construit sur le produit scalaire des tenseurs. Les résultats sont présentés sous forme d'une fonction de répartition de la valeur de l'indicateur pour le domaine complet, puis pour trois sous-domaines d'intérêt: l'entrée, une région où l'écoulement interagit fortement avec les parois, et une région où l'écoulement est fortement tourbillonnaire. L'hypothèse d'alignement n'est que rarement valide pour la loi de comportement de Boussinesq. Pour la QCR, les résultats sont meilleurs en entrée, comparé à la loi de Boussinesq. Il ne sont cependant pas meilleurs pour les régions où l'écoulement est plus tourbillonnaire. Une amélioration de la loi de comportement est nécessaire pour pouvoir faire progresser la modélisation turbulente en RANS. En revanche, l'utilisation de l'énergie cinétique turbulente et du taux de dissipation spécifique semble correcte pour estimer la valeur de la viscosité turbulente. L'analyse de la modélisation de l'équation d'énergie cinétique turbulente (ECT) est réalisée au travers d'une comparaison terme à terme avec l'équation d'ECT résolue par la SGE. Les résultats SGE présentent une turbulence qui n'est pas à l'équilibre : la production et la dissipation ne sont pas superposées, et le terme de transport est important. Pour le RANS, la turbulence est à l'équilibre : la production et la dissipation sont superposées, et le terme de transport est de faible intensité. Un modèle de turbulence qui prend en compte le déséquilibre est nécessaire pour améliorer ce point. En dernier lieu, une nouvelle formulation hybride RANS/SGE est proposée, fondée sur la distance à la paroi en unités de paroi. La formulation est validée dans un canal bi-périodique et un premier essai est réalisé sur le cas de décollement de coin, mais d'autres analyses sont nécessaires avant que cette formulation ne soit fonctionnelle. / The present study aims at analysing turbulence modelling in Reynolds-averaged Navier-Stokes (RANS) simulations, in the context of turbomachinery flows, using large-eddy simulations (LES) as references. Two test cases are considered: a corner separation (CS) flow in a linear compressor cascade, and a tip-leakage (TL) flow of a single blade in a jet. Two constitutive relations, the Boussinesq constitutive relation and the quadratic constitutive relation (QCR), are investigated, with two versions of Wilcox's $k-\omega$ turbulence model. The studied constitutive relations rely on two hypotheses: an alignment hypothesis between the Reynolds stress tensor and a mean flow tensor, and an hypothesis on the turbulent viscosity. The alignment hypothesis is investigated using LES, where both the tensors are known independently, with an indicator built on the inner product of the tensors. The results are presented as probability density functions of the indicator value for the entire domain first, and then for three specific areas of interest: the inlet area, similar to a boundary-layer flow, an area of strong interaction between the flow and the walls (CS: passage area, TL: tip clearance) and an area of highly vortical flow (CS: separation wake, TL: tip-leakage vortex). The alignment hypothesis is rarely verified in any area for the Boussinesq constitutive relation. For the QCR, the results are improved for the inlet areas compared to the Boussinesq constitutive relation, but no significant improvement is found in the highly vortical regions. An improvement of the constitutive relation is needed in order to improve the RANS turbulence modelling. In contrast, the use of the turbulent kinetic energy and the specific dissipation rate appears quite correct to estimate the turbulent viscosity. The modelling of the RANS turbulent kinetic energy (TKE) budget equation is investigated through a term to term comparison with the resolved LES TKE budget equation. The LES presents a turbulence that is not at equilibrium, with the production and the dissipation not superimposed, and an important amount of transport. This differs from the RANS models, at equilibrium: the production and the dissipation are superimposed, with a small amount of transport. The development of a non-equilibrium turbulence model for RANS simulations could improve this aspect of turbulence modelling. Finally, a new hybrid RANS-LES formulation, based on the wall distance in wall units, is also proposed. It is validated on a bi-periodical channel flow, and a first attempt is made on the corner separation case, but further investigations are still needed for the model to be fully operational.

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.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada. / To study and forecast the environmental impacts that a hydroelectric power plant reservoir may produce, a numerical simulation that takes into account the effects of stratification and scalar mixing (such as temperature, salinity or water-soluble substance) is required. This work proposes a methodology for the study of the environmental fluid flow phenomena, mainly for flows in which the knowledge of the interaction between stratification and mixing can give important notions of the phenomena that occur. For this, a numerical simulation tool for 3D environmental flow is developed. A tetrahedral mesh generator of the reservoir based on the terrain topology and an algebraic turbulence model based on the Richardson number are the main tools developed. The main difficulty in tetrahedral mesh generation of a reservoir is nonuniform distribution of the points related to the huge ratio between the horizontal and vertical scales of the reservoir. In this type of point distributions, conventional tetrahedron mesh generation algorithm may become unstable. For this reason, a unstructured tetrahedral mesh generator is developed and the methodology used to obtain conforming elements is described. Triangular surface mesh generation using the Delaunay triangulation and the construction of the tetrahedra from the triangular surface mesh are the main steps to the mesh generator. The hydrodynamic simulation of reservoirs with a turbulence model provides a useful tool that is computationally viable for engineering purposes. Furthermore, the turbulence model based on the Richardson number takes into account the effects of interaction between turbulence and stratification. The algebraic model is the simplest among the various models of turbulence, but provides realistic results with the fitting of a small amount of parameters. Eddy-Viscosity/Diffusivity models for stratified turbulent flows models are incorporated. Using the Finite Element Method (FEM) approximation the Reynolds-averaged Navier-Stokes (RANS) and mean scalar transport equations are approximated. The convective terms are discretized employing the Semi-Lagrangian method, and the spatial discretization is based on the Galerkin method. The computational results are compared with the results available in the literature. Finally, the simulation of the flow in a branch of a reservoir is presented.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada. / To study and forecast the environmental impacts that a hydroelectric power plant reservoir may produce, a numerical simulation that takes into account the effects of stratification and scalar mixing (such as temperature, salinity or water-soluble substance) is required. This work proposes a methodology for the study of the environmental fluid flow phenomena, mainly for flows in which the knowledge of the interaction between stratification and mixing can give important notions of the phenomena that occur. For this, a numerical simulation tool for 3D environmental flow is developed. A tetrahedral mesh generator of the reservoir based on the terrain topology and an algebraic turbulence model based on the Richardson number are the main tools developed. The main difficulty in tetrahedral mesh generation of a reservoir is nonuniform distribution of the points related to the huge ratio between the horizontal and vertical scales of the reservoir. In this type of point distributions, conventional tetrahedron mesh generation algorithm may become unstable. For this reason, a unstructured tetrahedral mesh generator is developed and the methodology used to obtain conforming elements is described. Triangular surface mesh generation using the Delaunay triangulation and the construction of the tetrahedra from the triangular surface mesh are the main steps to the mesh generator. The hydrodynamic simulation of reservoirs with a turbulence model provides a useful tool that is computationally viable for engineering purposes. Furthermore, the turbulence model based on the Richardson number takes into account the effects of interaction between turbulence and stratification. The algebraic model is the simplest among the various models of turbulence, but provides realistic results with the fitting of a small amount of parameters. Eddy-Viscosity/Diffusivity models for stratified turbulent flows models are incorporated. Using the Finite Element Method (FEM) approximation the Reynolds-averaged Navier-Stokes (RANS) and mean scalar transport equations are approximated. The convective terms are discretized employing the Semi-Lagrangian method, and the spatial discretization is based on the Galerkin method. The computational results are compared with the results available in the literature. Finally, the simulation of the flow in a branch of a reservoir is presented.