Understanding the physics of mixing for two fluids is a complicated problem and has always been an interesting phenomenon to study. Surface discharge is the oldest, least expensive and simplest way of discharging industrial or domestic wastewater into rivers and estuaries. Because of the lower degree of dilution in surface discharges, critical conditions are more likely to occur. Having a better understanding of the mixing phenomenon in these cases will help to predict the environmental effects more accurately. In this study, surface discharges of jets into waterbodies with or without crossflows were investigated numerically and experimentally. Three-dimensional (3-D) Computational Fluid Dynamics (CFD) models were developed for studying the surface discharge of jets into water bodies using different turbulence models. Reynolds stress turbulence models and spatially filtered Large Eddy Simulation (LES) were used in the numerical models. The effects of inclusion of free surface water in the CFD models on the performance of the numerical model results were investigated. Numerical model results were compared with the experimental data in the literature as well as the experimental works performed in this study. Experimental works for buoyant and non-buoyant surface discharge of jets into crossflow and stagnant water were conducted in this study. A new setup was designed and built in the Civil Engineering Hydraulics Laboratory at the University of Ottawa to perform the desired experiments. Stereoscopic Particle Image Velocimetry (Stereo-PIV) was used to measure the instantaneous spatial and temporal 3-D velocity distribution on several planes of measurement downstream of the jet with the frequency of 40 Hz. Averaged 3-D velocity distribution was extracted on different planes of measurement to show the transformation of the velocity vectors from a “jet-like” to a “plume-like” flow regime. Averaged 3-D velocity distribution and streamlines illustrated the flow transformation of the surface jets. Experimental results detected the formation and evolution of vortices in the surface jet’s flow structure over the measurement zone. Additional turbulent flow characteristics such as the turbulent kinetic energy (k), turbulent kinetic energy dissipation rate (ϵ), and turbulent eddy viscosity (υt) were calculated using the measured time history of the 3-D velocity field.

L'Océan Austral est une zone clef pour la circulation océanique tant à cause de l'intensité du courant circumpolaire antarctique qu'en tant que région de formation des masses d'eaux abyssales de l'océan global. Pour modéliser l'océan et prévoir les changements climatiques futurs, il est important de comprendre les processus de mélange diapycnal qui lient ces eaux abyssales aux couches supérieures. Dans l'Océan Austral, des courants profonds et intenses s'écoulent sur une topographie accidentée, ce qui génère des ondes internes de relief très énergétiques. Actuellement, la dissipation de l'énergie induite par ces ondes de relief est la candidate principale pour expliquer les forts taux de mélange observés à ces latitudes. L'objet du présent travail de thèse est de comprendre comment les ondes internes de relief sont dissipées et affectent la circulation et le mélange diapycnal dans l'océan abyssal. Nous examinons l'impact de ces ondes sur le mélange profond au moyen d'une combinaison d'expertise de terrain, de simulations non hydrostatiques bi-dimensionnelles et de calculs théoriques. Sur la gamme de paramètres étudiés, nous montrons, en présence des ondes de relief, une intensification du taux de dissipation d'énergie cinétique turbulente sur une profondeur de 1000 m au-dessus de la topographie, atteignant typiquement ~20 mW/m2. Nous montrons également comment les ondes participent à des interactions triadiques impliquant des oscillations inertielles qui sont amplifiées par intéractions résonantes contrôlées par les ondes de relief. Finalement, nous préparons de futures études tri-dimensionnelles en concevant un cadre numérique et en décrivant des outils théoriques adaptés à ce problème. Nos résultats préliminaires en trois dimensions montrent qui le confinement méridien de la topographie réduit significativement l'émissions d'ondes internes de relief. / The Southern Ocean plays a key role in global ocean circulation by connecting the major ocean basins with the intense Antarctic Circumpolar Current and as a formation region for abyssal water masses of the global ocean. Understanding the diapycnal mixing processes that link these abyssal waters to the overlying layers is essential both for ocean modelling and for predicting future climate change. In the Southern Ocean, deep reaching currents impinge on rough topography and create highly energetic internal lee waves. The dissipation of the energy of these internal lee waves is the main candidate for explaining the high mixing rates between waters of different densities observed at these latitudes. The purpose of this study is to understand the fate of the internal lee wave energy and how it affects the circulation and diapycnal mixing in the abyssal ocean. We first study the impact of internal lee waves on deep mixing with the combination of field expertise, two-dimensional non hydrostatic numerical simulations and theoretical developments. Over the range of parameters studied, an enhanced bottom turbulent kinetic energy dissipation is observed in the bottom 1000 m, typically reaching $sim$ 20 mW.m$^{-2}$. We further show that internal lee waves undergo non-dissipative wave-wave interactions that can be rationalized as resonant triad interactions between the bottom emitted internal lee waves, inertial oscillations and linear combinations of these two waves. We then build a three-dimensional model configuration and specific diagnostic methods that pave the way for future investigations in three dimensions. Preliminary results with the three-dimensional numerical configuration show that the meridional confinement of the topography notably reduces the emission of internal lee waves.

Depuis une dizaine d’années, l’évolution de la puissance des ordinateurs a permis de développer l’utilisation, dans les études d’ingénierie, des simulations 3D CFD (Computational Fluid Dynamics) pour l’étude de l’atmosphère à petite échelle, en particulier pour la dispersion de polluants sur des sites industriels et urbains complexes. Compte tenu de la complexité des domaines à étudier et des ressources de calcul généralement disponibles, ces études sont la plupart du temps réalisées à l’aide des modèles RANS (Reynolds Averaged Navier-Stokes), et particulièrement avec le modèle de fermeture k – e. Différents travaux de validation de l’approche RANS k – e ont mis en évidence quelques limitations à reproduire la dynamique de l’écoulement et de la dispersion dans des configurations géométriques complexes. Le travail de recherche réalisé dans le cadre de cette thèse a pour objectif une caractérisation expérimentale fine de l’écoulement et de la dispersion turbulente autour d’un obstacle bidimensionnel placé dans une couche limite de surface, afin d’évaluer la validité des modèles RANS en vue de leur application pour l’étude de la dispersion atmosphérique.Dans un premier temps, nous avons utilisé des techniques d’anémométrie à fil chaud, d’anémométrie laser Doppler et d’anémométrie par image de particules, pour déterminer le champ de vitesse dans une couche limite de surface rugueuse et autour d’un obstacle bidimensionnel de section carrée. Une attention particulière a été portée sur l’analyse des termes de l’équation évolutive de l’énergie cinétique turbulente (ECT) et sur la détermination de la viscosité turbulente vt. Différentes approches ont également été utilisées pour estimer le taux de dissipation e de l’énergie cinétique turbulente. Nous avons mis en évidence que ces différentes approches fournissent des résultats comparables dans le cas de la couche limite, tandis que seule la technique estimant e comme le résidu de l’ECT est applicable dans le sillage de l’obstacle. De plus, nos mesures ont permis d’évaluer les paramétrisations du modèle k – e et de montrer que la valeur du coefficient Cμ = 0.09 ne semble pas adaptée dans le cas de la couche limite, conduisant à une surestimation de vt, alors que cette valeur apparait plus adaptée dans le cas de l’obstacle. Une étude de sensibilité, portant la détermination de la constante σk du modèle k – e, indique une contribution non négligeable des termes de corrélation entre la vitesse et la pression dans le sillage de l’obstacle.Dans un deuxième temps, nous avons étudié la dispersion d’un scalaire passif, en mesurant les différents moments statistiques de la concentration, au moyen d’un détecteur à ionisation de flamme. Nous avons également déterminé les flux turbulents de masse, par un couplage entre les mesures de vitesse et de concentration, en prenant soin de contrôler les influences réciproques des deux techniques de mesure. Ces mesures nous ont permis de tester la validité de différents modèles de fermeture de l’équation d’advection-diffusion pour estimer les flux dans le sens vertical et dans le sens longitudinal. Nous avons également pu déterminer expérimentalement le coefficient de diffusivité turbulente Dt, nous permettant d’évaluer un nombre de Schmidt turbulent Sct, afin de mettre en évidence que la valeur Sct = 0.7 est adaptée à la majorité des zones étudiées, excepté dans la zone de recirculation induite par l’obstacle. Enfin, nous nous sommes intéressés aux différents termes de l’équation de la variance de la concentration et plus particulièrement à son taux de dissipation. À nouveau, les mesures nous ont permis de tester un modèle de fermeture disponible dans la littérature et de montrer la bonne cohérence entre le modèle et l’expérience. / In the last decades, there has been an increasing use of Computational Fluid Dynamics (CFD)simulations to evaluate the impact of atmospheric pollutants dispersion in within industrial and urban sites. Given the high geometrical complexity of these sites, these simulations are mainly performed adopting a Reynolds Averaged Navier-Stokes (RANS) approach and using k−e closure models. As is well known from previous studies, RANS k−e simulations are affected by some limitations that prevent them correctly reproducing the dynamics of the flow and the pollutant dispersion in complex geometrical configurations. The aim of the PhD is to provide a detailed experimental characterization of the flow and the turbulent dispersion around an idealized two-dimensional obstacle placed within a boundary layer flow. This is subsequently used to analyse the reliability of RANS closure models as predictive tools for the atmospheric dispersion of airborne pollutants. Initially we focus on the flow dynamics of a boundary layer flow developing over a rough wall and in the wake of a 2D obstacle. The velocity field is investigated experimentally by means of different measurement techniques, namely Hot Wire Anemometry (HWA), Laser Doppler Anemometry (LDA) and Stereo-Particle Imagery Velocimetry (PIV). A particular attention was devoted to the estimate of the turbulent viscosity nt as well as of the terms composing the turbulent kinetic energy budget (TKE), including its rate of dissipation e which was determined adopting different approaches. These measurements allowed us to analyse the accuracy of the parameterizations included in a standard k−e closure model. Our analysis show that a value of the coefficient Cμ = 0.09 leads to significant overestimation of nt in a boundary layer flow. Conversely, adopting Cμ = 0.09 provides a good agreement between modeled and direct estimates of nt in the wake of the obstacle. As a second step, we studied the dispersion of a passive scalar emitted by a ground level line source. To that purpose we measured the first four order moments of the concentration probability density function by mean of a flame ionization detector (FID). Furthermore, the coupling of the FID system with the LDA or HWA system allowed us to directly measure the turbulent mass transfer, i.e. the correlation between velocity and concentration fluctuations. The combination of these two techniques was carefully analyzed, in order to quantify eventual mutual disturbances of one measurement technique on the other. The measurements of the velocity/concentration correlations allowed us to determine experimentally the turbulent diffusivity Dt and the turbulent Schmidt number Sct , and therefore to test the accuracy of different closure models of the advection-distribution equation. Our results show that the value of the turbulent Schmidt number is approximately equal to 0.7 in most of the domain, except in the recirculation zone on the wake of the obstacle. Experimental data provide also a complete description of the spatial distribution of the concentration variance, and of the term composing its budget (with a focus on its dissipation). As for the velocity field, we test the reliability of different closure model proposed in the literature of the turbulent mass transfer terms, enlightening the shortcomings of simple gradient-law closer models.