L'objectif de cette thèse est d'étudier numériquement l'aéroacoustique à faibles nombres de Mach(M inf 0.3) pour un écoulement de couche limite turbulente épaisse affleurant une cavité, sur la base de simulations numériques à grandes échelles (LES). Un profil de vitesse en loi puissance et pour une couche limite d'équilibre ont servi comme conditions en entrée du domaine de calcul. La couche limite d'équilibre, sans et avec gradient de pression adverse, a été résolue par une approche asymptotique basée sur une formulation déficitaire avec un nouveau modèle de longueur de mélange. Ce dernier a été validé pour améliorer les comparaisons avec les expériences et les simulations numériques directes. Des simulations LES ont permis de regarder l'influence de l'épaisseur de la couche limite turbulente amont sur le mode d'oscillation d'une cavité L/D=4. Un accord satisfaisant avec les expériences d'Haigermoser et l'émergence du mode de cisaillement a été obtenu pour la vitesse amont de 5.8m/s. Le mode était de type sillage pour les deux autres cas tests (20et40m/s). Finalement, une simulation 3D a montré que le mode de sillage est un artefact du calcul 2D. En utilisant l'analogie de Lighthill-Curle et les champs de pression in stationnaire issus de la simulation, nous avons déterminé les niveaux de pression sonore dans le champ proche et lointain. Conformément aux expériences d'Haigermoser, une faible directivité vers l'amont est trouvée. Le mode de sillage influence très fortement les niveaux de pression acoustique. / The objective of this thesis is to study numerically the aeroacoustics of low Mach number (M inf 0.3) fow with thick turbulent boundary layer past a cavity based on Large Eddy Simulation (LES). Velocity profiles from power law and equilibrium turbulent boundary layer were imposed as inlet conditions on the computational domain. The equilibrium turbulent boundary layer profles (zero and adverse pressure gradient) have been generated using a symptotic approach with an improved mixing length model. A good agreement is observed between the computed boundary layer profiles and the profiles obtained from experiments and direct numerical simulations. LES results present the infuence of the thickness of the incoming turbulent boundary layers on the mode of oscillation in the shallow cavity of L/D=4. An agreement with the experiments of Haigermoser and the shear mode have been found for the upstream velocity 5.8m/s. Wake mode was observed for the other two test cases at 20 and 40m/s. A 3D cavity simulation is performed to show that the wake mode observed in the 2D calculations is an artifact. The hydrodynamic pressure feld obtained from the 2D simulation is used as an input to the acoustic analogy (Lighthill-Curle's analogy), to compute the acoustic pressure feld at the near and far feld of the cavities. Conforming the experiments of Haigermoser, a weak directivity of sound propagation was observed. Shear mode infuences the sound pressure levels strongly.

Améliorer la représentation de la couche limite stable constitue un des grands challenges de la prévision numérique du temps et du climat. Sa représentation est clé pour la prévision du brouillard, du gel des surfaces, des inversions de température, du jet de basse couche et des épisodes de pollution. De plus, à l'échelle climatique, la hausse de la température moyenne globale de l'air en surface impacte davantage les régions polaires : améliorer la représentation de la couche limite stable est un enjeu important pour réduire les incertitudes autour des projections climatiques. Depuis une quinzaine d'années, les exercices d'intercomparaison de modèles GABLS ont montré que le mélange turbulent dans la couche limite stable est généralement surestimé par les modèles de prévision du temps. En effet, de nombreux modèles intensifient artificiellement l'activité de leur schéma de turbulence afin d'éviter une décroissance inévitable du mélange lorsque la stabilité dépasse un seuil critique en terme de nombre de Richardson gradient. Ce problème numérique et théorique n'est pas en accord avec de nombreuses observations et simulations à haute résolution qui montrent une activité turbulente séparée en deux régimes : un régime faiblement stable dans lequel l'atmosphère est turbulente de manière continue et intense, et un régime très stable dans lequel la turbulence est très intermittente, anisotrope et faible en intensité. Ces travaux de thèse s'articulent autour de deux parties dont l'objectif principal est d'améliorer la paramétrisation de la turbulence dans le modèle atmosphérique de recherche Méso-NH développé conjointement par Météo-France et le Laboratoire d'Aérologie, et dans le modèle opérationnel AROME. Cette étude utilise une méthodologie communément employée dans le développement de paramétrisations qui consiste à comparer des simulations à très haute résolution qui résolvent les structures turbulentes les plus énergétiques (LES) à des simulations uni-colonnes d'un modèle méso-échelle. Plusieurs simulations 3D couvrant différents régimes de stabilité de l'atmosphère sont réalisées avec Méso-NH. Les limites du modèle LES en stratification stable sont documentées. Une première partie répond à la problématique de la surestimation du mélange dans le régime faiblement stable. Une expression originale pour la longueur de mélange est formulée. La longueur de mélange est un paramètre clé pour les schémas de turbulence associés à une équation pronostique pour l'énergie cinétique turbulente. Cette longueur de mélange non-locale combine un terme de cisaillement vertical du vent horizontal à une formulation existante qui repose sur la flottabilité. Le nouveau schéma est évalué dans des simulations 1D par rapport aux LES d'une part ; et dans le modèle opérationnel AROME par rapport aux observations de l'ensemble du réseau opérationnel de Météo-France d'autre part. Une deuxième partie apporte des éléments d'évaluation d'un schéma combinant deux équations pronostiques pour les énergies cinétiques et potentielles turbulentes. En condition stable, le flux de chaleur négatif contribue à la production d'énergie potentielle turbulente. L'interaction entre les deux équations d'évolution permet, via une meilleure prise en compte de l'anisotropie et d'un terme à contre gradient dans le flux de chaleur, de limiter la destruction de l'énergie turbulente dans les modèles. Dans les cas simulés, cette nouvelle formulation ne montre pas un meilleur comportement par rapport à un schéma à une équation pour l'énergie cinétique turbulente car le mécanisme d'auto-préservation n'est pas dominant par rapport au terme de dissipation. Il conviendra d'améliorer la paramétrisation du terme de dissipation dans le régime très stable. / The modeling of the stable atmospheric boundary layer is one of the current challenge faced by weather and climate models. The stable boundary layer is a key for the prediction of fog, surface frost, temperature inversion, low-level jet and pollution peaks. Furthermore, polar regions, where stable boundary layer predominates, are one of the region with the largest temperature rise : the stable boundary layer modeling is crucial for the reduction of the spread of climate predictions. Since more than 15 years, the GABLS models intercomparison exercices have shown that turbulent mixing in the stable boundary layer is overestimated by numerical weather prediction models. Numerous models artificially strengthen the activity of their turbulence scheme to avoid a laminarization of the flow at a critical value of the gradient Richardson number. The existence of this threshold is only a theoretical and a numerical issues. Numerous observations and high-resolution numerical simulations do not support this concept and show two different regimes : the weakly stable boundary layer that is continuously and strongly turbulent; and the very stable boundary layer globally intermittent with a highly anisotropic and very weak turbulence. This thesis aims at improving the turbulence scheme within the atmospheric research model Méso-NH developped by Météo-France and the Laboratoire d'Aérologie, and the operational weather forecast model AROME. We use a traditional methodology based on the comparison of high-resolution simulations that dynamically resolve the most energetic turbulent eddies (Large-Eddy Simulations) to single-column simulations. Several LES covering the weakly and the very stable boundary layer were performed with Méso-NH. The limits of applicability of LES in stratified conditions are documented. The first part of the study deals with the overmixing in the weakly stable boundary layer. We propose a new diagnostic formulation for the mixing length which is a key parameter for turbulence schemes based on a prognostic equation for the turbulent kinetic energy. The new formulation adds a local vertical wind shear term to a non-local buoyancy-based mixing length currently used in Méso-NH and in the French operational model AROME. The new scheme is evaluated first in single-column simulations with Méso-NH and compared to LES, and then in the AROME model with respect to observations collected from the operational network of Météo-France. The second part presents a theoretical and numerical evaluation of a turbulence scheme based on two prognostic equations for the turbulent kinetic and potentiel energies. In stratified conditions, the heat flux contributes to the production of turbulent potential energy. The laminarization of the flow is then limited by a reduction of the destruction of the turbulent kinetic energy by a better representation of the anisotropy and a counter-gradient term in the heat flux. On the simulated cases, this new formulation behaves similarly than the scheme with one equation for the turbulent kinetic energy because the self-preservation mechanism is not dominant compared to the dissipation term. Further research should improve the turbulent kinetic energy dissipation closure in the very stable regime.

Dans ce travail, nous avons abordé le phénomène de la Convection Thermique Turbulente sous des angles nouveaux et innovants. <br /><br />Le premier système que nous abordons expérimentalement est un Canal Vertical dit Infini au sein duquel règne un gradient de température moyen constant. Le flux de masse dans ce canal est nul. L'image qui se dégage de nos mesures est celle d'un écoulement majoritairement inertiel, où les coefficients dissipatifs (la viscosité en l'occurrence) n'interviennent que pour fixer une longueur de cohérence L. Cette longueur est celle sur laquelle les panaches thermiques peuvent êtres considérés comme en « chute libre ». Le transport horizontal (d'impulsion et de chaleur) est entièrement dû aux fluctuations. La « longueur de mélange » associée est petite devant la largeur du canal. Par contre, le transport de chaleur vertical est dû à des structures cohérentes, les panaches. <br />Ces panaches, nous les retrouvons dans une étude Lagrangienne de l'écoulement au centre d'une cellule de Rayleigh Bénard. La sonde est une sphère de 2 cm de diamètre qui a la même densité que le fluide que nous utilisons, équipée de thermomètres et d'un émetteur radio. Elle est transportée par les panaches, ce qui nous permet une étude statistique de ceux-ci.

De nouveaux modèles d'atmosphère sont présentés, incluant les profils de raie d'hélium neutre améliorés de Beauchamp (1995) et le formalisme de probabilité d'occupation pour ce même atome. Ces modèles sont utilisés pour calculer une grille de spectres synthétiques correspondant à des atmosphères riches en hélium et contenant des traces d'hydrogène. Cette grille est utilisée pour déterminer les paramètres atmosphériques principaux des étoiles de notre échantillon, soient la température effective, la gravité de surface et l'abondance d'hydrogène. Notre échantillon contient des spectres visibles de haut rapport signal-sur-bruit pour 102 naines blanches riches en hélium, dont 29 ont été observés au cours de ce projet, ce qui en fait le plus grand échantillon de spectres de qualité de naines blanches riches en hélium. Des spectres synthétiques ont été calculés en utilisant différentes valeurs du paramètre α de la théorie de la longueur de mélange dans le but de calibrer empiriquement la valeur de ce paramètre pour les DB. Afin d'améliorer la précision sur les paramètres atmosphériques de quelques étoiles, nous avons utilisé des spectres couvrant la raie Hα pour mieux déterminer l'abondance d'hydrogène. Finalement, nous avons calculé la distribution de masse de notre échantillon et la fonction de luminosité des DB. La distribution de masse montre une coupure à 0.5 fois la masse solaire qui est prédite par les modèles d'évolution stellaire et dévoile une masse moyenne significativement plus élevée pour les étoiles de type DBA. La masse moyenne de l'ensemble des DB et DBA est très proche de celle des DA. La fonction de luminosité nous permet de calculer que le rapport du nombre de DB sur le nombre de DA vaut environ 25%. / New model atmospheres are presented, including improved neutral helium lines from Beauchamp (1995) and the occupation probability formalism for that atom. These models are used to compute a grid of synthetic spectra for helium rich atmospheres with different hydrogen abundances. This grid is used to determine the principal atmospheric parameters of the stars in our sample, e.g. effective temperature, surface gravity and hydrogen abundance. There are 102 high quality spectra of helium-rich white dwarfs in our sample, making it the largest sample of this kind. 29 of these spectra were observed for this project. Synthetic spectra using different values of the α parameter from the mixing length theory have been calculated in order to determine the correct value of this parameter for DB model atmospheres. Finally, we have computed the mass distribution of our sample and the DB luminosity funtion. The mass distribution shows a clear cutoff at 0.5 solar masses which is predicted by stellar evolution theory and gives a significantly higher mean mass for the DBA stars of the sample. However, the global mean mass of our sample is very close to that of DA stars. With our luminosity function, we found a number ratio of DB stars over DA stars of about 25%.

De nouveaux modèles d'atmosphère sont présentés, incluant les profils de raie d'hélium neutre améliorés de Beauchamp (1995) et le formalisme de probabilité d'occupation pour ce même atome. Ces modèles sont utilisés pour calculer une grille de spectres synthétiques correspondant à des atmosphères riches en hélium et contenant des traces d'hydrogène. Cette grille est utilisée pour déterminer les paramètres atmosphériques principaux des étoiles de notre échantillon, soient la température effective, la gravité de surface et l'abondance d'hydrogène. Notre échantillon contient des spectres visibles de haut rapport signal-sur-bruit pour 102 naines blanches riches en hélium, dont 29 ont été observés au cours de ce projet, ce qui en fait le plus grand échantillon de spectres de qualité de naines blanches riches en hélium. Des spectres synthétiques ont été calculés en utilisant différentes valeurs du paramètre α de la théorie de la longueur de mélange dans le but de calibrer empiriquement la valeur de ce paramètre pour les DB. Afin d'améliorer la précision sur les paramètres atmosphériques de quelques étoiles, nous avons utilisé des spectres couvrant la raie Hα pour mieux déterminer l'abondance d'hydrogène. Finalement, nous avons calculé la distribution de masse de notre échantillon et la fonction de luminosité des DB. La distribution de masse montre une coupure à 0.5 fois la masse solaire qui est prédite par les modèles d'évolution stellaire et dévoile une masse moyenne significativement plus élevée pour les étoiles de type DBA. La masse moyenne de l'ensemble des DB et DBA est très proche de celle des DA. La fonction de luminosité nous permet de calculer que le rapport du nombre de DB sur le nombre de DA vaut environ 25%. / New model atmospheres are presented, including improved neutral helium lines from Beauchamp (1995) and the occupation probability formalism for that atom. These models are used to compute a grid of synthetic spectra for helium rich atmospheres with different hydrogen abundances. This grid is used to determine the principal atmospheric parameters of the stars in our sample, e.g. effective temperature, surface gravity and hydrogen abundance. There are 102 high quality spectra of helium-rich white dwarfs in our sample, making it the largest sample of this kind. 29 of these spectra were observed for this project. Synthetic spectra using different values of the α parameter from the mixing length theory have been calculated in order to determine the correct value of this parameter for DB model atmospheres. Finally, we have computed the mass distribution of our sample and the DB luminosity funtion. The mass distribution shows a clear cutoff at 0.5 solar masses which is predicted by stellar evolution theory and gives a significantly higher mean mass for the DBA stars of the sample. However, the global mean mass of our sample is very close to that of DA stars. With our luminosity function, we found a number ratio of DB stars over DA stars of about 25%.