Losses associated with turbulent flows dissipate a significant amount of generated energy. Such losses originate from the drag force, which is often described as the sum of the pressure drag and the friction drag. This thesis sets out to explore the hypothesis that passive wall motions driven by fluid mechanical forces are able to reduce the friction drag in fully developed turbulent boundary layers. Firstly, the streamwise and spanwise opposition controls proposed by Choi et al. (1994, Journal of Fluid Mechanics) are revisited to identify beneficial wall motions. Near-wall streamwise or spanwise velocity fluctuations are measured along a detection plane parallel to the wall (sensing). For streamwise control, the wall velocities are set to be equivalent to the measured streamwise velocity fluctuations, whereas for spanwise control, the wall velocities are set to have the same magnitude but opposite direction as the measured spanwise velocity fluctuations (actuation). Direct numerical simulations of canonical turbulent channel flows are carried out at low (Reτ ≈ 180) and intermediate (Reτ ≈ 1000) Reynolds numbers to quantify the effect of the distance between the wall and the detection plane. The investigation reveals the primary differences between the mechanisms underlying the two active in-plane controls. The modified flow features and turbulence statistics show that the streamwise control amplifies the most energetic streamwise velocity fluctuations and damps the near-wall vorticity fluctuations. In comparison, the spanwise control induces near-wall vorticity in order to counteract the quasi-streamwise vortices of the near-wall cycle and suppress turbulence production. Although, the working principles of the active controls are fundamentally different, both achieve drag reduction by mitigating momentum transfer between the velocity components. Secondly, two theoretical passive compliant wall models are proposed, the aim being to sustain beneficial wall motions identified by active flow control simulations. In the proposed models, streamwise or spanwise in-plane wall motions are governed by an array of independent one-degree-of-freedom damped harmonic oscillators. Unidirectional wall motions are driven by local streamwise or spanwise wall shear stresses. A weak coupling scheme is implemented to investigate the interaction between the compliant surface models and the turbulent flow in the channel by means of direct numerical simulations. A linear analytical solution of the coupled differential equation system is derived for laminar pulsatile channel flows allowing verification and validation of the numerical model. The obtained analytical solution is utilised to map the parameter space of the passive controls and estimate the effect of the wall motions. It is shown that depending on the control parameters, the proposed compliant walls decrease or increase the vorticity fluctuations at the wall similarly to the active controls. This is confirmed by direct numerical simulations. On the one hand, when the control parameters are chosen appropriately, the passive streamwise control damps the near-wall vorticity fluctuations and sustains the same drag reduction mechanism as the active streamwise control. This leads to modest, 3.7% and 2.3% drag reductions at low and intermediate Reynolds numbers. On the other hand, the spanwise passive control is not capable of increasing the near-wall vorticity fluctuations as dictated by the active spanwise control. For this reason, passive spanwise wall motions can increase the friction drag by more than 50%. The results emphasise the necessity of anisotropy for a practical compliant wall design. The present work demonstrates for the first time that passive wall motions can decrease friction drag in fully turbulent wall-bounded flows. The thesis sheds light on the working principle of an active streamwise control, and proposes a passive streamwise control exploiting the same drag reduction mechanism. An analytical model is developed to give a ready prediction of the statistical behaviour of passive in-plane wall motions. Whereas streamwise passive wall motions are found beneficial when the control parameters are chosen appropriately, solely spanwise passive wall motions lead to a drag penalty.

Des mélanges de gaz et de particules sont présents dans divers environnements géophysiques. De tels mélanges chauds sont générés par des éruptions volcaniques explosives et comprennent des écoulements de conduit, des jets et des panaches, ainsi que des courants de densité pyroclastiques. La concentration de particules dans ces mélanges volcaniques peut varier fortement, allant de concentrations élevées (>50 % en volume) dans les écoulements denses fluidisés à des concentrations très faibles dans les suspensions diluées dans lesquelles les particules sont mises en suspension par la phase gazeuse turbulente. Une limite de concentration inférieure à ~% en volume dans les suspensions diluées a été suggérée par des études récentes, car des concentrations plus élevées nécessiteraient une énergie cinétique turbulente excessive. L'objectif principal de cette thèse est d'étudier expérimentalement le comportement d'un écoulement d'air turbulent dans un cylindre avec des concentrations de particules croissantes, pour différents nombres de Reynolds et en utilisant différents types de particules. Les nombres de Reynolds des mélanges gaz-particules dans les expériences atteignaient ~106. Une première série d'expériences a été menée avec des billes de verre de différentes tailles allant de 75-80 μm jusqu'à 2 mm, pour un total de huit tailles de particules. Au-dessus d'un seuil de concentration moyenne de 0.5-3 % en volume, qui augmentait avec le nombre de Reynolds, le comportement de l'écoulement a montré une transition d'une suspension homogène de particules (sous la concentration maximale) vers une séparation en une partie basale dense et une partie supérieure diluée contenant la concentration maximale des particules. Ce seuil de concentration a été détecté à l'aide de mesures de pression et d'une méthode impliquant une sphère dont la densité était légèrement inférieure à la densité apparente des particules et qui pouvait donc flotter au-dessus de la partie basale dense, si celle-ci était présente. Des vidéos à haute vitesse ont révélé que l'apparition de la concentration maximale de particules coïncidait avec l'émergence d’amas de particules dans la partie turbulente diluée. Dans une deuxième partie de la thèse, les expériences ont été répétées pour cinq gammes de tailles de particules de céramique et elles ont révélé le même comportement général que pour les billes de verre. Pour les deux types de particules, une concentration maximale a pu être détectée pour presque toutes les tailles de particules et a montré une augmentation avec le nombre de Reynolds à la puissance 1/5 (billes de verre) ou 0.4 (billes de céramique). Compte tenu du nombre de Reynolds des particules, la concentration maximale des particules augmente ensuite jusqu'à la puissance de 1/6 pour les particules de céramique et de verre. Ces résultats ouvrent de nouvelles perspectives sur la structure des mélanges gaz-particules volcaniques et ils fournissent également des contraintes pour les données d'entrée et de sortie des simulations numériques et pour les observations géophysiques. / Mixtures consisting of gas and particles can be found in various geophysical environments. Hot mixtures are generated by explosive volcanic eruptions and include conduit flows, jets and buoyant plumes, and pyroclastic density currents. The particle concentration within these volcanic mixtures can vary highly, from high concentrations (>50 vol. %) in dense fluidized flows to very low concentrations in dilute suspensions in which the particles are suspended by the turbulent gas phase. A concentration limit of less than ~1 vol. % in dilute suspensions was suggested by recent studies, as higher concentrations would require excessive turbulent kinetic energy. The main objective of this thesis was to investigate experimentally the behavior of a turbulent air flow in a pipe with increasing particle concentrations, for different Reynolds numbers and using different types of particles. The Reynolds numbers of the gas-particle mixtures in the experiments were up to ~106. A first set of experiments was conducted with glass beads of varying sizes from 75-80 μm up to 2 mm, for eight particle size ranges in total. Above a bulk concentration threshold of ~0.5-3 vol. %, which increased with the Reynolds number, the flow behavior changed from a homogeneous suspension of particles (below the maximum concentration) to a separation into a dense basal part and an upper dilute part carrying the maximum concentration of particles. This concentration threshold was detected with pressure measurements and a method that involved a ball of a slightly lower density than the bulk density of the particles, which could thus float over the dense basal part, if present. High-speed videos revealed that the occurrence of the maximum particle concentration coincided with the emergence of particle clusters in the dilute turbulent part. In a second part of the thesis, the experiments were repeated for five ceramic particle size ranges and they yielded the same general behavior as for the glass beads. For both types of particles, a maximum concentration could be detected for almost all particle size ranges and showed an increase with the mixture Reynolds number to the power 1/5 (glass beads) or 0.4 (ceramic beads). Considering the particle Reynolds number the maximum particle concentration then increase to the power 1/6 for both ceramic and glass particles. These results give new insights about the structure of volcanic gas-particle mixtures and they also provide constraints for input and output data of numerical simulations and for geophysical observations.

Among the different types of flows encountered in practical applications, the physics of turbulent separated flows and turbulence induced secondary motion are not fully understood despite the large amount of previous experimental and numerical work. The objectives of this work are to study theoretically and computationally the conditions at the separation and reattachment point, the numerical simulation of turbulence induced secondary motion in non-circular ducts, and to provide a comprehensive test of different RANS models of these types of flow. In a theoretical study of flow separation, a Lagrangian approach was first used to derive an Eulerian criterion, which associates separation and reattachment points to a critical point in the eigenvalues of the Cauchy-Green tensor. A turbulent separated boundary layer under the influence of an adverse pressure gradient was simulated using DNS and LES techniques. A bootstrapping method was used to obtain high fidelity results at a relatively high Reynolds number with which the performance of some of the most commonly used eddy-viscosity turbulence models was evaluated. The DNS and LES results were used to assess the consistency of the different terms in the k−e , ζ −f , k −ω and Spalart-Allmaras models. Different wall-modelling techniques were employed for the calculation of separated boundary layers. The exact values of the modelled terms were calculated using the reference DNS and LES dataset. These results were used for both a priori and a posteriori tests. It was determined that the eddy-viscosity assumption works well, and that anisotropic effects are not significant in separated boundary layer. For the secondary flow calculation in non-circular ducts, direct numerical simulations of turbulent flow in square and skewed ducts were carried out to determine the effect of the duct (rhombus) included angle on both the mean and turbulence energy budgets. Two skewed ducts, with included angles of 30 and 60 degrees, were simulated. The capability of different turbulence models to predict the secondary velocity field was investigated. Results obtained from a non-linear stress-strain constitutive relation was found to be fairly accurate for the flows at the range of Reynolds number considered in this study. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2010-11-26 13:52:18.361

A new hybrid approach to model high Reynolds number wall-bounded turbulent flows is developed based on coupling the two-level simulation (TLS) approach in the inner region with conventional large eddy simulation (LES) away from the wall. This new approach is significantly different from previous near-wall approaches for LES. In this hybrid TLS-LES approach, a very fine small-scale (SS) mesh is embedded inside the coarse LES mesh in the near-wall region. The SS equations capture fine-scale temporal and spatial variations in all three cartesian directions for all three velocity components near the wall. The TLS-LES equations are derived based on defining a new scale separation operator. The TLS-LES equations in the transition region are obtained by blending the TLS large-scale and LES equations. A new incompressible parallel flow solver is developed that accurately and reliably predicts turbulent flows using TLS-LES. The solver uses a primitive variable formulation based on an artificial compressibility approach and a dual time stepping method. The advective terms are discretized using fourth-order energy conservative finite differences. The SS equations are also integrated in parallel, which reduces the overall cost of the TLS-LES approach. The TLS-LES approach is validated and investigated for canonical channel flows, channel flow with adverse pressure gradient and asymmetric plane diffuser flow. The results suggest that the TLS-LES approach yields very reasonable predictions of most of the crucial flow features in spite of using relatively coarse grids.

Cette thèse s’inscrit dans le contexte de la simulation et du contrôle des vibrations de structures montées sur ressort qui peuvent apparaître sous l’effet de l’interaction avec l’écoulement de sillage instationnaire. Le contrôle de ce phénomène, appelé vibrations induites par vortex (VIV), est un enjeu critique dans l’optimisation de nombreux systèmes. Une méthode de frontières immergées (IBM) a été intégrée dans l’algorithme PISO du code OpenFOAM, dédié à la simulation d’écoulements fluides incompressibles. La méthode IBM permet une représentation précise de corps fixes ou en mouvement, tout en conservant des maillages structurés conduisant à des algorithmes plus précis et efficaces en termes de performances numériques. Pour calculer la divergence de l’équation de quantité de mouvement dans une boucle PISO et l’interpolation des flux, un calcul hybride orignal a été proposé avec une résolution analytique utilisant l’équation de la fonction noyau des quantités impliquant le terme force de l’IBM (quantités singulières). La méthode a été étendu au formalisme d’écoulements en régimes turbulents. Une loi de paroi a été intégrée permettant de modéliser la couche limite à grand nombre de Reynolds. Le travail de validation a été réalisé au regard des données expérimentales et numériques disponibles dans la littérature pour l’étude d’écoulements autour de cylindres et de sphères, sur une large gamme de nombres de Reynolds. Avec l’objectif de développer des lois de contrôle optimal pour le VIV, basées sur les mécanismes d’instabilité linéaire du système couplé dans le cadre de la théorie du contrôle, un solveur adjoint a été développé et validé. / This thesis is related to the simulation and the control of the vortex induced vibrations phenomenon (VIV), which can result from the fluid structure interactions between an unsteady wake and the body, when the shedding frequency in the wake is close to the natural frequency of the body. The control of VIV is a critical issue when optimizing many systems. An Immersed Boundaries Method (IBM) was implemented into the PISO algorithm as a new library of OpenFOAM, in order to perform reliable simulations of incompressible flows around bluff bodies.To compute the divergence of the momentum equation and the interpolation of the fluxes, an hybrid calculation with an analytical resolution of the quantities involving the force term (singular quantities) has been proposed. The mesh convergence of several errors was shown by means of a manufactured solution, allowing to analyze both the errors irelated to the discretization and to the IBM. The new algorithm was subsequently extended to the RANS and DDES formalism proposed in OpenFOAM for the simulation of turbulent flows. A wall law was integrated into theIBM method to model the boundary layers that develop around the bodies at large Reynolds numbers. Various 2D and 3D well-documented test cases of academic flows around fixed or moving solid bodies (cylinderand sphere) have been simulated and carefully validated against existing data from the literature in a large range of Reynolds numbers. With the objective of developing optimal control laws for VIV, based on the linear instability mechanisms of the coupled system within the framework of the control theory, a new adjoint solver was also developed and validated in OpenFOAM.

Cette thèse vise à démêler les principaux mécanismes impliqués dans les écoulements transitoires et turbulents. L’idée centrale est d'utiliser une technique d’optimisation non linéaire pour étudier l’origine et le rôle des structures cohérentes habituellement observées dans ces écoulements. Cette méthode a été utilisée dans trois contextes différents. Tout d’abord, un écoulement laminaire linéairement stable a été considéré et l'optimisation a été utilisée pour calculer les perturbations les plus amplifiées parmi toutes les perturbations capables de déclencher une transition vers la turbulence. Une fois que la turbulence est bien établie, une optimisation non linéaire entièrement 3D maximisant l'énergie cinétique turbulente est utilisée pour étudier les structures cohérentes qui peuplent l’écoulement turbulent et les mécanismes responsables de la croissance et de l’échange d’énergie (optimale) sont étudiés. Ensuite, une approche de type système dynamique est appliquée aux équations du mouvement. La géométrie de l’espace des phases est étudiée en utilisant la théorie de la croissance transitoire pour évaluer l’importance des variétés stable et instable dans la dynamique. Dans le même cadre, un algorithme de minimisation non linéaire est utilisé pour calculer les connexions hétérocliniques parmi les solutions invariantes des équations de Navier-Stokes. / This thesis aims at unraveling the main mechanisms involved in transitional and turbulent flows. The central idea is that of using a nonlinear optimization technique to investigate the origin and role of coherent structures usually observed in these flows. This method has been used in three different contexts. First, a linearly stable laminar flow has been considered and the optimization has been used to compute the most amplified perturbations among all disturbances able to trigger transition to turbulence. Once turbulence is well established, a fully 3D nonlinear optimization maximizing the turbulent kinetic energy is used to study coherent structures populating turbulent shear flow as well as investigate the mechanisms responsible for the energy (optimally) growth and exchange. Then, a dynamical system approach is applied to fluid flow equations. The geometry of the state space is investigated by using transient growth theory to reveal the importance of the stable and unstable manifold. In the same framework, a nonlinear minimization algorithm is used to compute heteroclinic connections among invariant solutions of the Navier-Stokes equations.

This work uses data from the three components of wind and temperature, collected in an experimental site in the Amazon rainforest for 10 months during GOAmazon Project. A total of 10 levels of sensors were deployed on a micrometeorological tower. Focusing on the nocturnal boundary layer, an analysis of the temporal scales of the motion, using the multiresolution decomposition, has shown that the contributions from horizontal, nonturbulent fluctuations with long temporal scales, can be as significant as purely turbulent fluctuation. On, weakly stable nights the dominant temporal scales of the flow are those associated with, downward (sweeps) and upward events (ejections), which occur with fullydeveloped turbulence, , having dominant time scales between 10 and 100 s. Through the analysis of two-point correlations, it was possible to show that horizontal events with long time scales propagate from the top to within the canopy, being detected at different times, and. The vertical component correlations are larger at the upper canopy, not showing any time delay. The occurrence of positive sensible heat flux near the surface, with times scales larger than 100 s, has been identified in the study of overall averages. It is hypothesized that on very stable nights, non-turbulent modes associated with longer time scales, and referred as "submeso" have great impact on the horizontal components, becoming an important cause of the flow near the forest floor. In these situations, the most relevant time scales are longer than 300 a and dominate the almost entire vertical profile. For these cases, the correlations of turbulent variables decay rapidly a, being between 0.2 and 0.3 for the horizontal components while not exceeding 0.1 for the vertical component does. This reinforces the hypothesis that, for this scenario, the most correlated events between the top of the canopy and its interior are horizontal in nature. / O presente estudo utiliza dados das três componentes do vento e da temperatura, coletados em um sítio experimental na Floresta Amazônica, durante 10 meses, através do Projeto GOAmazon, no qual 10 níveis de sensores foram dispostos ao longo de uma torre micrometeorológica. Com ênfase na camada limite noturna, uma análise das escalas temporais do movimento, usando o método de decomposição em multiresolução, mostrou que as contribuições devido aos movimentos horizontais, de escalas mais longas, podem ser tão significativas quanto as puramente turbulentas. Em noites que a camada limite é fracamente estável, movimentos descendentes (varreduras) e ascendentes (ejeções), associados à turbulência bem desenvolvida, são os responsáveis pela intensidade turbulenta existente, com suas escalas temporais mais significativas entre 10 e 100 s. Através da análise das correlações entre pontos foi possível mostrar que estes eventos propagam-se desde a copa até o interior do dossel, sendo detectados em diferentes instantes de tempo, mais intensamente nas componentes horizontais, conforme alcançam os níveis mais profundos e que estas decaem verticalmente, sendo bem correlacionadas até aproximadamente 0,8 h. Já para a componente vertical as correlações são altas em todos os níveis do perfil e não apresentam atraso. Além disso, a ocorrência de fluxos de calor sensível positivo próximo a superfície, em escalas temporais maiores que 100 s, foi identificada no estudo das médias gerais. Em noites de condições de estabilidade alta, modos não turbulentos, associados a escalas temporais mais longas, chamados submeso têm grande impacto nas componentes horizontais do movimento e tornam-se os principais causadores dos fluxos. Nestas situações, as escalas temporais mais relevantes são maiores que 300 s e dominam praticamente todo o perfil vertical. Para estes casos as correlações das variáveis turbulentas decaem rapidamente e apresentam atraso quase nulo, entretanto as componentes horizontais apresentam correlações entre 0,2 e 0,3 nos níveis mais baixos do dossel, enquanto a componente vertical não ultrapassa 0,1. Isto reforça a hipótese de que, para este cenário, os eventos correlacionados entre a copa e o interior do dossel são horizontais.

Les écoulements turbulents de gaz denses, qui sont d’un grand intérêt pour un large éventail d'applications, sont le siège de phénomènes physiques encore peu connus et difficiles à étudier par des approches expérimentale. Dans ce travail, nous étudions pour la première fois l’influence des effets de gaz denses sur la structure de la turbulence compressible à l’aide de simulations numériques. Le fluide considéré est le PP11, un fluorocarbure lourd, dont le comportement thermodynamique a été représenté à l’aide de différentes lois d’état, afin de quantifier la sensibilité des solutions aux choix de modélisation. Nous avons considéré d’abord la décroissance d’une turbulence homogène isotrope compressible. Les fluctuations de température sont négligeables, alors que celles de la vitesse du son sont importantes à cause de leur forte dépendance de la densité. Le comportement particulier de la vitesse du son modifie de manière significative la structure de la turbulence, conduisant à la formation de shocklets de détente. L’analyse de la contribution des différentes structures à la dissipation d’énergie et à la génération d’enstrophie montre que, pour un gaz dense, les régions de forte dilatation jouent un rôle similaire à celles de forte compression, contrairement aux gaz parfaits, dans lesquels le comportement est fortement dissymétrique. Ensuite, nous avons mené des simulations numériques pour une configuration de canal plan en régime supersonique, pour plusieurs valeurs des nombres de Mach et de Reynolds. Les résultats confirment la validité de l’hypothèse de Morkovin. L’introduction d’une loi d’échelle semi-locale prenant en compte le variations de densité et viscosité, permet de comparer les profils des grandeurs turbulentes (contraintes de Reynolds, anisotropie, budgets d’énergie) avec ces observés en gaz parfait. Les variables thermodynamiques, quant à elles, présentent une évolution très différente pour un gaz parfait et pour un gaz dense, la chaleur spécifique élevée de ce dernier conduisant à un découplage des effets dynamiques et thermiques et à un comportement proche de celui d’un fluide incompressible avec des propriétés variables. / Dense gas turbulent flows, of great interest for a wide range of engineering applications, exhibit physical phenomena that are still poorly understood and difficult to reproduce experimentally. In this work, we study for the first time the influence of dense gas effects on the structure of compressible turbulence by means of numerical simulations. The fluid considered is PP11, a heavy fluorocarbon, whose thermodynamic behavior is described by means of different equations of state to quantify the sensitivity of solutions to modelling choices. First, we considered the decay of compressible homogeneous isotropic turbulence. Temperature fluctuations are found to be negligible, whereas those of the speed of sound are large because of the strong dependence on density. The peculiar behavior of the speed of sound significantly modifies the structure of the turbulence, leading to the occurrence of expansion shocklets. The analysis of the contribution of the different structures to energy dissipation and enstrophy generation shows that, for a dense gas, high expansion regions play a role similar to high compression ones, unlike perfect gases, in which the observed behaviour is highly asymmetric. Then, we carried out numerical simulations of a supersonic turbulent channel flow for several values of Mach and Reynolds numbers. The results confirm the validity of the Morkovin' hypothesis. The introduction of a semi-local scaling, taking into account density and viscosity variations across the channel, allow to compare the wall-normal profiles of turbulent quantities (Reynolds stresses, anisotropy, energy budgets) with those observed in ideal gases. Nevertheless, the thermodynamic variables exhibit a different evolution between perfect and dense gases, since the high specific heats of the latter lead to a decoupling of dynamic and thermal effects, and to a behavior close to that of variable property incompressible fluids.

Cette thèse vise à développer et évaluer la méthode de Galerkin discontinue (DG) pour la simulationdes grandes échelles (LES) des écoulements turbulents. L’approche DG présente un nombre d’avantages intéressants pour la LES : ordre élevé, stencil compact, prise en compte des maillages non structurés et expression de la solution numérique dans une base de polynômes permettant l’utilisation de modèles de turbulence multi-échelle. Parmi ce type de modèles, nous nous sommes intéressés ici à la méthode Variational Multiscale (VMS) qui consiste à séparer les échelles résolues dans la base de polynômes pour restreindre l’influence du modèle à une gamme réduite d’échelles. Les modèles considérés ont été paramétrés en prenant en compte les fonctions de transfert spécifiques aux discrétisations DG. La précision de la méthode pour la représentation de phénomènes turbulents variés a été évaluée à travers la réalisation de DNS de configurations académiques. Enfin, l’approche VMS/DGa été éprouvée sur des configurations simples à haut nombre de Reynolds. Il apparaît que cette méthodologie permet la représentation précise des phénomènes turbulents pour un coût réduit en terme de degrés de liberté. / This work focuses on the development of the Discontinuous Galerkin (DG) method for the large-eddy simulation (LES) of turbulents flows. The DG method shows some interesting properties for LES : high-order of accuracy, compact stencil, unstructured meshes and amodal polynomial basis which can be used to implement multiscale turbulence models. We consider in this work the Variational Multiscale approach (VMS), which consists in splitting the resolved scales into two components using the modal basis in order to restrict the action of the model to a given range of small scales. The models have been tuned using the transfer functions of the DG hp-discretizations. The accuracy of the DG method for the representation of turbulent phenomena has been assessed through DNS of free and wall-bounded canonical flows. Finally, the VMS/DG approach has been assessed for simple configurations at high Reynolds numbers. We have shown that this particular approach allows for an accurate representation of turbulent flows for coarse discretizations.

La ventilation joue un rôle important dans le confort thermique des occupants d'un bâtiment en climat chaud, en contribuant au rafraîchissement de l'air qui les entoure. Qu'elle soit mécanique ou naturelle, la ventilation doit être maîtrisée pour ne pas gêner l'occupant et respecter des normes ou réglementations en vigueur. Ces gênes sont liées à des vitesses d'air ou à une intensité turbulente trop élevée. Les concepteurs doivent alors faire appel à l'outil numérique pour une prédiction fine des écoulements d'air. La simulation de configurations à l'échelle du bâtiment peut se faire par une approche moyennée des équations de Navier-Stokes en complément d'un modèle de turbulence. Ces simulations sont utilisées par les chercheurs comme des outils de dimensionnement, ou encore, d'optimisation des composants de ventilation. De plus, la forme des bouches de ventilation peut contribuer passivement à l'optimisation de certains phénomènes aérauliques. L'amélioration de ces formes peut ainsi se faire par l'utilisation de méthodes d'optimisation de forme. L'optimisation topologique par ajout de matière permet de trouver des formes pour optimiser des fonctionnelles objectifs définies sur le fluide ou sur ses frontières. C'est sur cette méthode que ces travaux de thèse se concentrent pour proposer un outil de contrôle des écoulements d'air dans le bâtiment par la recherche de formes optimales de bouches de ventilation. Ces travaux de thèse proposent une contribution à l'optimisation topologique d'écoulements turbulents dans le bâtiment. Dans un premier temps, la méthode par ajout de matière est appliquée pour minimiser les pertes de charge dans une conduite d'aération en forme de Té. Le modèle adjoint développé est soumis à l'hypothèse de turbulence gelée. Dans un second temps le modèle adjoint complet est proposé pour le modèle de turbulence standard k-epsilon pour la réduction des pertes de charge d'une part et de l'intensité turbulente d'autre part. Enfin, ces outils sont appliqués à l'optimisation de forme de bouches de ventilation. Les résultats montrent ainsi un bon potentiel de l'optimisation topologique par ajout de matière pour l'orientation des écoulements d'air mais ne garantissent pas la maîtrise des vitesses d'air dans la pièce. De plus, la minimisation de l'intensité turbulente grâce à l'approche complète développée a contribué à la réduction du taux d'insatisfaction lié à une intensité turbulente trop élevée dans la pièce. / Ventilation plays a key role in thermal comfort of building's occupants in hot climates by refreshing air surrounding them. Mechanical or natural ventilation must be controlled for two reasons: do not disturb the bulding's occupants and comply with the regulations in force. Discomfort is linked to too high air velocities or turbulent intensity. Designers can use the numerical tools for a finer prediction of airflow. The simulation of configurations at the building scale can be done using averaged Navier-Stokes equations approach in addition to a turbulence model. These simulations are used by researchers as sizing tools or for the optimization of ventilation components. In addition, the shape of the ventilation nozzle can passively contributes to the optimization of some aeraulics phenomena. The improvement of these ventilation components can be achieved by the use of shape optimization methods. Topological optimization by addition of material makes it possible for the optimization of cost functions defined on the fluid or on its boundaries. The main objective of this manuscript is to propose a tool to control airflows in building by the search for optimal shape of ventilation nozzle. This work proposes a contribution to the topological optimization of turbulent flows in buildings. In a first step, topological optimization by adding material is applied to minimize pressure losses in a T-shaped pipe. The developed model is subjected to the hypothesis of the frozen turbulence. In a second step, the complete adjoint model is proposed for the standard turbulence model k-epsilon for the minimization of the total pressure losses on the one hand and the turbulent intensity on the other hand. Finally, these tools are applied to the shape optimization of ventilation nozzle. The results of topological optimization by adding virtual material show good potential for the orientation of the airflows but does not guarantee the control of the air velocities in the room. Moreover, the minimization of turbulent intensity through the complete approach contributed to the reduction of the dissatisfaction rate due to excessive turbulent intensity in the room.