L’interaction onde de choc/couche limite en écoulement transsonique autour de profils aérodynamiques est étudiée numériquement utilisant différentes classes de modélisation de la turbulence. Les approches utilisées sont celles de modèles URANS et de méthodes hybrides RANS-LES. L’emploi d’une correction de compressibilité pour les fermetures à une équation est aussi évalué. Premièrement, la séparation intermittente induite par le choc sur un profil supercritique en conditions d’incidence proches de l’angle critique d’apparition du tremblement est analysée. Suite à des simulations URANS, la modélisation statistique la mieux adaptée est étudiée et utilisée dans l’approche DDES (Delayed Detached-Eddy Simulation). L’étude de la topologie de l’écoulement, des pressions pariétales et champs de vitesse statistiques montrent que les principales caractéristiques de l’oscillation auto-entretenue du choc sont capturées par les simulations. De plus, la DDES prédit des fluctuations secondaires de l’écoulement qui n’apparaissent pas en URANS. L’étude de l’interface instationnaire RANS-LES montre que la DDES évite le MSD (modeled stress depletion) pour les phases de l’écoulement attaché ou séparé. Le problème de la ‘zone grise’ et de son influence sur les résultats est considéré. Les conclusions de l’étude sur le profil supercritique est ensuite appliquées à l’étude numérique d’un profil transsonique laminaire. Dans ce contexte, l’effet de la position de la transition de la couche limite sur les caractéristiques de deux régimes d’interaction choc/couche limite sélectionnés est étudié. En conditions de tremblement, les simulations montrent une forte influence du point de transition sur l’amplitude du mouvement du choc et sur l’instationnarité globale de l’écoulement. / Shock wave/boundary layer interactions arising in the transonic flow over airfoils are studied numerically using different levels of turbulence modeling. The simulations employ standard URANS models suitable for aerodynamics and hybrid RANS-LES methods. The use of a compressibility correction for one-equation closures is also considered. First, the intermittent shock-induced separation occurring over a supercritical airfoil at an angle of attack close to the buffet onset boundary is investigated. After a set of URANS computations, a scale-resolving simulation is performed using the best statistical approach in the context of a Delayed Detached-Eddy Simulation (DDES). The analysis of the flow topology and of the statistical wall-pressure distributions and velocity fields show that the main features of the self-sustained shock-wave oscillation are predicted by the simulations. The DDES also captures secondary flow fluctuations which are not predicted by URANS. An examination of the unsteady RANS-LES interface shows that the DDES successfully prevents modeled-stress depletion whether the flow is attached or separated. The gray area issue and its impact on the results are also addressed. The conclusions from the supercritical airfoil simulations are then applied to the numerical study of a laminar transonic profile. Following a preliminary characterization of the airfoil aerodynamics, the effect of the boundary layer transition location on the properties of two selected shock wave/boundary layer interaction regimes is assessed. In transonic buffet conditions, the simulations indicate a strong dependence of the shock-wave motion amplitude and of the global flow unsteadiness on the tripping location.

Cette thèse vise à une meilleure compréhension de la dynamique du voilier et à la validation des outils numériques de prédiction de performances et d’optimisation par l'étude expérimentale in situ du problème aéro-élastique d'un gréement. Une instrumentation est développée sur un voilier de 8m de type J80 pour la mesure dynamique des efforts dans le gréement, de la forme des voiles en navigation, du vent et des attitudes du bateau. Un effort particulier est apporté à la mesure des caractéristiques géométriques et mécaniques des éléments du gréement, la calibration des capteurs et au système d'acquisition des données. Les principaux résultats montrent que le voilier instrumenté est un outil adapté pour les mesures instationnaires et soulignent l'amplitude de variation d'effort rencontrée en mer (20 à 50% de l'effort moyen dans une houle modérée). En outre, les variations du signal d'effort sont déphasées avec l'angle d'assiette, créant un phénomène d'hystérésis. Le comportement dynamique d'un voilier en mouvement diffère ainsi de l'approche quasi-statique. Les simulations numériques proviennent du code ARAVANTI, couplage implicite d’un code structure éléments finis ARA et d’un code fluide parfait, limitant son domaine de validité aux allures de près Les résultats de simulation sont très proches des cas stationnaires et concordent bien avec les mesures en instationnaire dans une houle de face. L'expérimentation numérique d'un gréement soumis à des oscillations harmoniques en tangage souligne l'importance de l'approche Interaction Fluide Structure (IFS) et montre que l’énergie échangée par le système avec la houle est reliée à la fréquence réduite et l'amplitude du mouvement. Certaines informations n'étant pas disponibles sur le voilier instrumenté, une expérience contrôlée en laboratoire est développée. Elle consiste en un carré de tissu tenu par deux lattes en oscillation forcée. Les mesures sur cette « voile oscillante » permettent d'étudier les phénomènes IFS avec décollement et sont utilisées pour la validation du couplage ARA-ISIS entre un code fluide Navier-Stokes (RANS) et le même code structure. / This work presents a full scale experimental study on the aero-elastic wind/sails/rig interaction in real navigation condition with the aim to give a reliable database of unsteady measurement. This database is used for the investigation of the dynamic behavior and loads in the rigging and for an experimental validation of an unsteady Fluid Structure Interaction (FSI) model. An inboard instrumentation system has been developed on a 8 meter yacht (J80 class) to simultaneously and dynamically measure the navigation parameters, yacht's motion, sails flying shape, wind and loads in the rigging. A special effort is made on mechanical and geometrical characteristics measurement, sensors calibration and data acquisition system synchronization. Results show that the instrumented boat is a reliable tool to measure the unsteady phenomena in navigation. Dynamic measurements at sea underline the load variation encountered, which represent 20 to 50% of the mean value in a moderate sea state. Oscillations of loads exhibit phase shift with the trim angle, reason for an hysteresis phenomenon, which shows that the dynamic behavior of a sail plan subject to yacht motion clearly deviates from the quasi-steady theory. Simulations are made with ARAVANTI, an implicit coupling of a Finite Element Method structural model ARA and an inviscid fluid model which restricts the simulation domain to upwind conditions. The simulation results compare very well with the experimental data for steady sailing conditions and show a good agreement in unsteady conditions (head swell). Numerical investigation of a sail plan submitted to harmonic pitching motion underlines the importance of FSI modeling and shows that the energy exchanged by the system with the swell increases with the motion reduced frequency and amplitude. Some information is not accessible on the instrumented boat and requires developing a controlled test case in laboratory. The experiment consists of a spinnaker fabric square mounted on two carbon battens moved in forced oscillation. This test case is used to study FSI phenomena with a separated flow and gives experimental results for the validation of the coupling ARA-ISIS of a RANS fluid model with the same structure model.

Cette thèse est principalement dédiée à la simulation des écoulements massivement décollés dans le domaine spatial. Nous avons restreint notre étude aux écoulements d'arrière-corps, pour lesquels ces décollements sont imposés par des changements brutaux de la géométrie. Dans le domaine spatial, le caractère fortement compressible des écoulements rencontrés impose l'utilisation de schémas numériques robustes. D'un autre coté, la simulation fine de la turbulence impose des schémas d'ordre élevé et peu dissipatifs. Ces deux spécifications, apparemment contradictoires, doivent pourtant coexister au sein d'une même simulation. Les modèles de turbulence ainsi que les schémas de discrétisation sont indissociables et leur couplage doit impérativement être considéré. Les schémas numériques doivent garder leur précision formelle dans des géométries complexes et des maillages très irréguliers imposés par le contexte industriel. Cette étude analyse le schéma de discrétisation utilisé dans le code de calcul FLUSEPA développé par Airbus Defence & Space. Ce schéma est robuste et précis pour des écoulements avec chocs et il présente une faible sensibilité au maillage (l'ordre 3 étant conservé même sur des maillages fortement perturbés). Malheureusement, le schéma possède une trop faible résolvabilité liée à un niveau de dissipation trop élevé pour envisager des simulations hybrides RANS/LES. Pour pallier à cet inconvénient, nous nous sommes penchés vers une solution basée sur un recentrage conditionnel et local : dans les zones dominées par des structures tourbillonnaires, une fonction analytique assure un recentrage local lorsque la stabilité numérique le permet. Cette condition de stabilité assure le couplage entre le schéma et le modèle. De cette manière, les viscosités laminaire et tourbillonnaire sont les seules à jouer un rôle dans les régions dominées par la vorticité et servent aussi à stabiliser le schéma numérique. Cette étude présente de plus une comparaison qualitative et quantitative de plusieurs modèles hybrides RANS/LES, à égalité de maillage et de schéma utilisés Pour cela, un certain nombre d'améliorations (notamment de leur capacité à résoudre les instabilités de Kelvin-Helmohlotz sans retard), proposées dans la littérature ou bien introduites dans cette thèse, sont prises en compte. Les applications numériques étudiées concernent des géométries allant de la marche descendante au lanceur spatial complet à échelle réduite. / This thesis is mainly dedicated to the simulation of massively separated flows in the space domain. We restricted our study to afterbody flows, where the separation is imposed by abrupt geometry changes. In the space domain, highly compressible flows require the use of robust numerical schemes. On the other hand, the simulation of turbulence imposes high-order low dissipative numerical schemes. These two specifications, apparently contradictory, must coexist within the same simulation. The coupling between turbulence models and discretization schemes is of the utmost importance and must be considered. Numerical schemes should keep their formal accuracy on complex geometries and on very irregular meshes imposed by the industrial context. In this research, we analyze the discretization scheme implemented in the FLUSEPA solver, developed by Airbus Defence & Space. Such a scheme is robust and accurate for flows with shocks and exhibits a low sensitivity to the grid (the third order of accuracy being ensured, even on highly irregular grids). Unfortunately, the scheme possesses a too low resolvability related to a too high numerical dissipation for RANS/LES simulations. To circumvent this problem, we considered a conditional and local re-centering strategy: in regions dominated by vortical structures, an analytic function provides local re-centering when a numerical stability condition is satisfied. This stability condition ensures the coupling between the numerical scheme and the model. In this way, only the turbulent and the laminar viscosities play a role in regions dominated by vorticity, and also allow to stabilize the numerical scheme. This study provides also a qualitative and quantitative assessment of several hybrid RANS/LES models, using the same grids and discretization scheme. For this purpose some recent improvements (improving their ability to trigger the Kelvin-Helmohlotz instabilities without delay), proposed in the litterature or suggested in this work, are taken into account. Numerical applications include geometrical configurations ranging from a backward facing step to realistic launcher configurations.

Mieux comprendre les spécificités de la dynamique des couches limites à grand nombre de Reynolds malgré les contraintes métrologiques et son coût de simulation numérique est crucial. A titre d'exemple, cette dynamique peut déterminer plus de la moitié de la traînée d'un avion en croisière. Décrire la turbulence pariétale peut guider la résolution numérique d'une partie des fluctuations à un coût maîtrisé par des stratégies WMLES (simulation des grandes échelles avec modèle de paroi). Les présentes analyses physiques de couches limites turbulentes incompressibles à gradient de pression nul et à grand nombre de Reynolds s'appuient sur des simulations numériques avancées. Après validation d'une base de données, le frottement moyen pariétal est décomposé selon l'identité FIK (Fukagata et al. (2002)), dont l'application malgré le développement spatial est discutée. Une analyse spectrale montre que les grandes échelles (\lambda_x > \delta) contribuent à environ la moitié du frottement vers Re_\theta = 10^4. Les limitations de l'identité FIK motivent la dérivation d'une décomposition physique de la génération du frottement dont le comportement asymptotique est alors relié à la production d'énergie cinétique turbulente dans la zone logarithmique. Pour mieux reconstruire les spectres spatiaux, une nouvelle méthode d'estimation de la vitesse de convection turbulente en fonction de la longueur d'onde des fluctuations, adaptée au développement spatial et à des signaux temporels de durée finie, est dérivée, interprétée et évaluée à Re_\theta = 13000. Certaines des conclusions éclairent des modifications d'une stratégie WMLES, le mode III de la méthode ZDES. / Better understanding the specificities of the dynamics of high-Reynolds number boundary layers despite metrological constraints and its numerical simulation cost is crucial. For instance, this dynamics can determine more than half of the drag of a cruising aircraft. Describing wall turbulence can guide the numerical resolution of some of the fluctuations at a limited cost by WMLES strategies (wall-modelled large eddy simulation). The present physical analyses of zero-pressure gradient incompressible turbulent boundary layers at high Reynolds number rely on advanced numerical simulations. After validating a database, mean skin friction is decomposed by means of the FIK identity (Fukagata et al. (2002)), whose application despite the spatial growth is discussed. A spectral analysis shows that the large scales (\lambda_x > \delta) contribute approximately half of the friction near Re_\theta = 10^4. The limitations of the FIK identity motivate the derivation of a physical decomposition of the generation of friction whose asymptotic behaviour is then related to turbulent kinetic energy production in the logarithmic layer. In order to better reconstruct spatial spectra, a new method to estimate the turbulent convection velocity as a function of the wavelength of the fluctuations, adapted to spatial growth and to temporal signals of finite duration, is derived, interpreted, and assessed at Re_\theta = 13000. Some of the conclusions enlighten modifications to a WMLES strategy, mode III of the ZDES method.

La modélisation à l’aide de l’approche RANS (Reynolds Averaged Navier Stokes) a été menée en trois dimensions, en considérant les régimes permanent et non permanent, dans le but de simuler les écoulements au niveau des jonctions, bifurcations et déversoirs d’orage. A travers ces trois exemples d’étude, plusieurs stratégies de modélisation (différentes combinaison lois de paroi/modèles de turbulence, différents algorithmes de couplage pression/vitesse, différents schémas de discrétisation spatiale, différentes conditions aux limites, différents types et tailles de maille, etc.) ont été testées et évaluées à l’aide de plusieurs indicateurs de performance (de type RMS –Root Mean Square) en s’appuyant sur les données de vitesses (vitesses moyennes in situ et champs de vitesses obtenus par PIV en laboratoire), hauteurs d’eau, débits (répartition de débits en bifurcation en laboratoire ou débits déversés in situ). Les résultats obtenus sont transposables aux autres cas de singularités et ouvrages spéciaux rencontrés en réseau d’assainissement et montrent que : i) les résultats des simulations 3D sont sensibles à la rugosité, aux conditions limites de hauteur et de vitesse ; ii) dans les trois cas d’étude, les schémas de discrétisation du second ordre et l’algorithme de couplage pression/vitesse PISO sont appropriés ; ii) la loi de paroi scalable couplée aux modèles de turbulence de type k-ε pour le cas des jonctions (avec un débit latéral inférieur ou égal au débit principal) ou des déversoirs semblent convenir, tandis que le modèle de turbulence RSM associé à la loi de paroi enhanced ou scalable permet de mieux représenter les écoulements à travers les bifurcations ou au niveau des jonctions lorsque le débit latéral est dominant. Sur la base de ces résultats, une méthodologie de modélisation plus générale définie en six étapes et fondée sur le guide proposée par Jakeman et al. (2006) a été mise au point. La méthodologie ainsi définie a été utilisée pour i) améliorer l’instrumentation du déversoir OTHU (Observatoire de Terrain en Hydrologie Urbaine) situé à Ecully, à partir de la simulation de sa courbe de fonctionnement et en fournissant les incertitudes sur les débits déversés obtenus ; ii) simuler l’implémentation des capteurs débitmétriques à l’aval d’une jonction. Elle a permis de concevoir et de dimensionner un nouveau dispositif de maîtrise des flux d’eau et de polluants déversés (technologie DSM – dispositif de surveillance et de maîtrise des flux déversés au milieu naturel). Ce dispositif a fait l’objet d’un dépôt de brevet international. Enfin, la mise en œuvre de cette méthodologie a été à l’origine de la découverte d’une nouvelle structure d’écoulement dans la branche latérale d’une bifurcation à 90°. L’analyse des résultats des simulations des écoulements mettant en évidence cette nouvelle structure a montré qu’il était possible de prédire l’apparition de cette dernière à partir du rapport d’aspect et du nombre de Froude. / The understanding of sewer flows behaviour is a key component for better urban drainage monitoring and management. However, these flows are conveyed across singularities (such as bends, drops, deviations, etc.) and special structures (combined sewer overflows –CSOs–, channels junction, dividing flow structures, etc.). These singularities and specific structures exhibit complex geometries, leading to open channel turbulent, three-dimensional and multiphase (pollutants and storm and sewer waters) flows. Using three-dimensional CFD (Computational Fluid Dynamics) platform allows a better understanding of mechanisms of contaminants transport through these structures and singularities, leading to a better sewer monitoring. In this thesis, 3D-RANS (Reynolds Averaged Navier Stokes) modelling approach under steady-state conditions is used in order to study flows within CSOs, junctions and bifurcations. Through these three structures, several modelling strategies (wall law/turbulence model combination, velocity/pressure coupling algorithm, spatial discretisation schemes, boundary conditions, computational mesh –shape and size of cells–, etc.) are tested and evaluated thanks to performance indicators (such as RMS –Root Mean Square– indicators) based on velocities (in situ mean velocities and PIV velocity fields obtained in laboratory), water depths and discharge (discharge repartition for bifurcation in laboratory or in situ overflow discharge, for CSOs). Results deriving from these tests are transposable to other singularities or special structures encountered in sewer network and suggest that: i) simulated CFD results are sensitive to the roughness coefficient; ii) for the three studied structures, second-order discretisation schemes and SIMPLE or PISO velocity/pressure coupling algorithm are appropriate; iii) scalable wall function associated to the group of k-ε turbulence model for junction flows (with a lateral inflow lower or equal to the main inflow) or for CSOs is appropriate, whereas RSM turbulence model associated to enhanced wall function allows a better representation of bifurcation flows or junctions flows when the lateral inflow is greater than the main inflow. Based on these results and on Jakeman et al. (2006) guideline, a six steps-methodology focused on the using of RANS approach modelling has been proposed. This six steps-methodology is used in order to i) enhance the monitoring of an OTHU (Observatoire de Terrain en Hydrologie Urbaine) CSO located at Ecully accounting for uncertainties on overflow discharge values ;ii) simulate the performance of flowmeters downstream of a junction, defining the best location for these sensors. This methodology is used to design the new overflow discharge measurement device. This device is an international patent. Finally, the application of the methodology led to point out a new flow structure, occurring in the downstream lateral branch of a 90° bifurcation.

L’écoulement de jeu dans les compresseurs axiaux est étudié à l’aide de la Zonal Detached Eddy Simulation (ZDES). L’objectif consiste à évaluer la capacité de méthodes hybrides URANS/LES à simuler l’écoulement de jeu d’un compresseur axial réaliste afin de mieux comprendre la physique de cet écoulement, notamment son comportement au vannage ainsi que l’effet de sillages venant du stator amont sur le rotor aval. Après avoir choisi la méthode hybride ZDES, un banc d’essai numérique est défini afin de simuler le premier rotor du compresseur de recherche CREATE. Ce banc a la particularité de pouvoir prendre en compte les effets instationnaires venant de la roue directrice d’entrée (RDE), notamment son sillage ainsi que les tourbillons générés en pied et en tête. Basé sur des critères de maillage ZDES, il est utilisé pour évaluer cette méthode comparativement aux méthodes classiques RANS et URANS. La ZDES est validée par étape jusqu’à une analyse spectrale de l’écoulement de jeu se basant sur des données expérimentales. Elle s’est révélée capable de capturer plus précisément l’intensité et la position des phénomènes instationnaires rencontrés en tête du rotor, notamment le tourbillon de jeu. Les densités spectrales de puissance analysées montrent que cela est dû en partie à une meilleure prise en compte du transfert d’énergie des grandes vers les petites structures de l’écoulement avant leur dissipation. De plus, l’écart entre les approches s’accentue lorsque le tourbillon de jeu traverse le choc en tête. Proche pompage, les effets d’interaction entre le choc, le tourbillon de jeu, la couche limite carter et le tourbillon venant de la tête de la RDE sont amplifiés. Le décollement de la couche limite carter s’accentue et une inversion locale de l’écoulement est observée. De plus, le tourbillon de jeu s’élargit et est dévié vers la pale adjacente, ce qui intensifie le phénomène de double écoulement de jeu. L’interaction du tourbillon venant de la tête de la RDE avec le choc et le tourbillon de jeu du rotor est ensuite étudiée au point de dessin. Un battement du tourbillon de jeu est rencontré lors de l’interaction de ce tourbillon avec le tourbillon de tête de la RDE, ce qui diminue le double écoulement de jeu. / The tip-leakage flow in axial compressors is studied with the Zonal Detached Eddy Simulation (ZDES). This study aims at evaluating the capability of hybrid URANS/LES methods to simulate the tip-leakage flow within a realistic axial compressor in order to better understand the involved physics, especially the behaviour of the flow near surge and the effects of stator wakes on the downstream rotor. Once the ZDES method is chosen, a numerical test bench is defined to simulate the first rotor of the research compressor CREATE. This bench takes into account the unsteady effects of the Inlet Guide Vane (IGV), such as its wake as well as vortices generated at the IGV hub and tip. It is based upon ZDES meshing criteria and is used to evaluate this method compared to classic RANS and URANS approaches. A method validation is carried out up to a spectral analysis compared to experimental data. The ZDES is capable to capture more accurately the intensity and position of the unsteady phenomena encountered in the tip region, especially the tip-leakage vortex. The power spectral densities highlight that this partly originates from a better capture of the energy transfer from large to small structures until their dissipation. The discrepancy between the methods is accentuated as the tip-leakage vortex crosses the shock. Near the surge line, the interactions between the shock, the tip-leakage vortex, the boundary layer developing on the shroud and the vortex generated by the IGV tip are amplified. The boundary layer on the shroud separates earlier and a local flow inversion occurs. Besides, the tip-leakage vortex widens and is deflected toward the adjacent blade. This strengthens the double leakage. At the design operating point, the interaction of the IGV tip vortex with the shock and the rotor tip vortex is studied. A vortex flutter is observed as the IGV tip vortex arrives on the rotor blade and stretches the rotor tip vortex. This phenomenon decreases the double leakage.

Ce travail de thèse concerne l'étude d'extracteurs d'air hybrides éoliens. Il se décompose en des investigations amont et appliquées. Dans la partie amont, des simulations fines ont été effectuées en conduite carrée avec et sans rotation, pour des nombres de Reynolds turbulents de l'ordre de 600, afin d'analyser l'impact de la rotation sur la turbulence. Elles ont montré que la rotation rompt la symétrie de l'écoulement. La partie appliquée est dédiée à la conception d'une nouvelle géométrie d'extracteur d'air. Cette géométrie a été proposée à partir de l'analyse de simulations RANS. Ses performances ont été confirmées par des mesures expérimentales sur banc d'essais. Les tests en soufflerie d'un système de captage d'énergie éolienne, conçu pour l'extracteur, ont mis en évidence son adéquation au régime de fonctionnement de l'extracteur. Les essais expérimentaux de l'extracteur complet, montrent que le système de captage apporte une part significative de l'énergie. Des essais en soufflerie ont permis d'observer le comportement global de l'extracteur. / This PhD work concerns the study of hybrid air extractors. It is composed of upstream and applied investigations. In the upstream part, fine simulations are realized in square duct flow with and without rotation to analyse the impact of rotation on turbulence. It is found that rotation removes symmetry property of the flow with turbulent Reynolds number of 600. The applied part is dedicated to the conception of a new air extractor geometry. This geometry is proposed from the analyse of RANS simulations. Its performances are confirmed by experimental measurements on test rig. Wind tunnel tests of a wind power capturing system, designed for the extractor, show a good adequation to the operating regime of the extractor. Experimental investigations on the complete air extractor, show the wind power capturing system brings a significant part of the energy. Wind tunnel tests allow to observe the complete air extractor behaviour.

La piscine fait partie des établissements publics les plus fréquentés dans notre société. En effet, il ne s’agit pas uniquement d’un lieu de pratique d'activités physiques, mais également un espace de détente, de jeu, d’éducation et de lien familial. Il est de toute évidence essentiel, de fournir un environnement intérieur confortable et sain pour ses occupants. Cependant, en raison de sa dimension, son besoin excessif en énergie et la complexité des phénomènes physiques évoluant à l’intérieur, il est difficile de parvenir à un équilibre optimum entre : qualité de l’air intérieur, confort thermique des occupants et efficacité énergique du bâtiment. Il faut pour cela, parvenir à une description des mécanismes qui façonnent la structure de l’écoulement de l’air par une analyse profonde de ces phénomènes qui sont à l'origine des transferts de chaleur et de masse mis en jeu à l’intérieur. Ainsi, l’objectif visé de cette thèse est de présenter une étude numérique thermo aéraulique, par CFD en régime stationnaire et transitoire, qui permet d’évaluer le comportement dynamique, thermique et thermodynamique des différents phénomènes physiques qui évoluent à l’intérieur de la piscine intérieure semi-olympique de l’université Bishop’s (Sherbrooke, Canada) afin d’améliorer la qualité de l’air intérieur et le confort thermique ainsi que son rendement énergétique. Les simulations sont réalisées avec le logiciel libre OpenFOAM en utilisant une approche RANS. Une étude thermo-aéraulique par CFD a d’abord été réalisée sur une cavité rectangulaire avec plancher chauffé, afin d’appréhender les simulations thermo aérauliques. Cela a abouti à la détermination de la meilleure configuration d’aération pour une qualité de l’air et un confort thermique optimum. Plusieurs simulations CFD du flux d'air tridimensionnel avec transfert de chaleur et de masse ont été aussi effectuées ultérieurement pour la piscine, afin d’évaluer les effets des conditions climatiques extérieures et ceux des nageurs sur l'atmosphère intérieure. En adoptant plusieurs modèles de turbulence de type RANS, la comparaison des résultats obtenus avec les données expérimentales de référence a permis de valider le code OpenFOAM. Les données expérimentales ont été recueillies au préalable au sein de la piscine de l’Université Bishop’s à l’aide d’un dispositif conçu et adapté aux conditions internes propre à la piscine et qui est équipé de plusieurs capteurs pour la mesure de : température, humidité relative et vitesse. Enfin, une étude thermo-aéraulique de la piscine en régime turbulent transitoire pour une durée de 24 heures pour les jours typiques d'été et d'hiver a été réalisée afin de prédire l’évolution de la distribution des paramètres tels que la vitesse, la température et l'humidité relative. Une analyse statistique a permis de montrer que les conditions climatiques extérieures n'ont pas d'effet sur l'environnement interne de celle-ci. D’ailleurs, sa très bonne isolation thermique démontrée par un calcul détaillé des pertes thermiques à travers son enveloppe confirme ce constat. D’autre part, l’évaluation de la qualité de l'air intérieur et le confort thermique des occupants a révélé que ces derniers sont inacceptables. Suite auxquels, un ajustement des paramètres de conditionnement de l’air a été apporté pour fin d’amélioration. / Abstract : The swimming pool is one of the most popular public establishments in our society and is not just a place for physical activities but also a space for relaxation, play, education and family ties. It is therefore important to ensure a healthy and comfortable indoor environment for the occupants. However, given the size, energy requirement and complexity of the physical phenomena that take place within such space, it is difficult to achieve an optimum balance between interior air quality, thermal comfort of occupants and energy efficiency of the building. This requires a description of the mechanisms, which determine the structure of the airflow by a profound analysis of these phenomena, which are the origin of the heat and mass transfers involved inside such spaces. The objective of this thesis is to present a numerical thermo-ventilation study using CFD (computational fluid dynamic) in stationary and transient regime that allows to evaluate the dynamic, thermal and thermodynamic behaviors of the various phenomena that take place inside the semi-Olympic closed swimming pool at Bishop's University (Sherbrooke, Qc, Canada). The aim is to improve the indoor air quality and thermal comfort of occupants as well as its energy efficiency. The simulations are carried out using OpenFOAM (Open Field Operation and Manipulation) using a Reynolds-Averaged Navier-Stokes (RANS) approach. To do this, a CFD thermo-ventilation study was first carried out on a rectangular cavity with heated floor in order to understand the thermo-ventilation simulations. This has led to the determination of the best ventilation configuration for optimum air quality and thermal comfort. Several CFD simulations of the three-dimensional airflow with heat and mass transfer were also carried out later for the indoor swimming pool to evaluate the effects of outdoor climatic conditions and swimmers on the indoor atmosphere of the pool. By adopting several RANS turbulence models, the comparison of the results obtained with the experimental data allowed to validate the OpenFOAM code. The experimental data were collected in the pool at Bishop's University using a device designed and adapted to the pool’s internal conditions. The devise is equipped with several sensors to measure temperature, relative humidity and velocity. Finally, a thermo-ventilation study of the swimming pool in transient turbulent regime for a duration of 24 hours for typical days of summer and winter was conducted in order to predict the distribution of the various parameters such as velocity, temperature and relative humidity. A statistical analysis showed that the external climatic conditions have no effect on the internal environment of the swimming pool. Moreover, its good thermal insulation demonstrated by a detailed calculation of the thermal losses through building envelope confirms this observation. On the other hand, the evaluation of the indoor air quality and the thermal comfort of occupants revealed that the conditions inside the pool are unacceptable. After which, an adjustment of the air conditioning parameters was made for improvements.

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

Ce travail est une contribution à la simulation numérique d'écoulements turbulents dans un but d'application industrielle. Nous nous intéressons dans un premier temps à une nouvelle combinaison "VMS-LES/procédure dynamique" pour la simulation d'écoulements autour de cylindres circulaire et carré. L'approche VMS-LES mise en œuvre a pour originalité d'utiliser une procédure de moyennage sur des volumes finis agglomérés dans le but de séparer les échelles, l'approche dynamique étant celle introduite par Germano en LES. Une approche hybride RANS/VMS-LES est ensuite évaluée sur le problème du cylindre circulaire à des nombres de Reynolds élevés. Cette approche introduit un paramètre d'hybridation qui selon la résolution locale de grille privilège le modèle RANS ou celui VMS-LES. Enfin, les performances d'un algorithme de Schwarz deux-niveau, qui utilise les méthodes de déflation et de balancing, sont examinées d'un point de vue efficacité et scalabilité dans le cas de simulations VMS-LES. / This work is a contribution to the numerical simulation of turbulent flows with the aim of industrial application. At first, we focus on a new combination "VMS-LES/dynamic procedure" for the simulation of flows around circular and square cylinders. The VMS-LES approach adopted in this work is original in using an averaging procedure over agglomerated finite volumes in order to separate the scales, the dynamic approach being the one introduced by Germano in LES. A RANS/VMS-LES hybrid approach is then evaluated on the circular cylinder test case at high Reynolds numbers. This approach introduces a hybridization parameter which privileges the RANS model or the VMS-LES model according to the grid resolution. Finally, the performance of a two-level Schwarz algorithm, which uses the deflation and balancing methods, are examined in terms of efficiency and scalability in the context of VMS-LES simulations.