Des bâtiments de grande hauteur sont construits avec un poids et un amortissement structurel de plus en plus faibles en lien avec l'évolution des techniques de construction et des matériaux. La connaissance des charges de vent dynamiques est un enjeu important pour la conception des grands bâtiments afin de garantir leur sécurité structurelle. L'objectif de cette thèse est d’évaluer la capacité de la simulation numérique des grandes échelles (LES) à prédire les charges de vent sur les structures et d’étudier l'influence des conditions d’entrée d’une simulation LES sur ces charges. Des expériences ont été menées à échelle réduite dans la soufflerie atmosphérique NSA du CSTB afin de documenter l’écoulement atmosphérique modélisé, de caractériser son interaction avec un bâtiment et les charges de vent statiques et dynamiques résultantes. Le sillage du bâtiment a été caractérisé grâce à des mesures PIV. Les efforts globaux et les pressions locales ont été mesurés par une balance et des prises de pression à haute fréquence. Ces expériences en soufflerie ont permis de développer un générateur de conditions amont (GCA) pour la simulation LES, visant à reproduire les principales caractéristiques de la turbulence dans la couche limite. La base de données constituée a également permis de qualifier les résultats des simulations LES réalisées avec le code OpenFOAM dans la configuration de l’expérience. L’utilisation du nouveau GCA et d’un générateur dégradé qui ne respecte pas toutes les caractéristiques de l'écoulement a permis de montrer la nécessité de bien reproduire les caractéristiques du vent incident pour accéder aux charges dynamiques sur le bâtiment. / High-rise buildings are built with increasingly low weight and structural damping in relation to the evolution of construction techniques and materials. The understanding of dynamic wind loads is an important issue for the design of high-rise buildings in order to guarantee their structural safety. The objective of the present work is to assess the ability of large eddy simulation (LES) to predict wind loads on structures and to investigate the influence of the inflow boundary conditions of a LES simulation on these loads. Experiments were carried out at a small scale in the NSA atmospheric wind tunnel of CSTB to document the modeled atmospheric boundary layer, to characterize its interaction with a building and the resulting static and dynamic wind loads. The wake flow around the building has been characterized by PIV measurements. Global and local wind loads were measured by a high frequency force balance and high frequency pressure taps. These wind tunnel experiments allowed for the development of an inflow turbulence generator for the LES simulation,which was aimed at reproducing the main characteristics of turbulence in the boundary layer. The database also made it possible to assess the quality of the results of the LES simulations carried out with the OpenFOAM code in the same configuration as the experiment. The use of both the new turbulence generator and a degraded one that does not account for all the characteristics of the flow has made it possible to show the necessity to reproduce the characteristics of the upstream wind flow in order to access the dynamic wind loads on the building.

En hiver, par condition anticyclonique hivernale, les régions montagneuses urbanisées font l'objet de conditions atmosphériques stables et découplées, qui conduisent à des épisodes de forte pollution particulaire. Ce travail de thèse traite de la caractérisation des facteurs de pollution de l'air aux particules par condition stable en vallée alpine encaissée, du point de vue de la dynamique atmosphérique. Ce travail repose sur des simulations numériques idéalisées et en conditions réelles à l'aide du modèle Weather Research and Forecasting (WRF), les particules étant modélisées par un champ de traceur passif.Dans une première partie, des simulations numériques à haute résolution ont été réalisées, en utilisant deux configurations différentes de vallées tridimensionnelles idéalisées ouvrant sur une plaine. La première configuration correspond à une vallée de largeur constante dans la direction de l'axe de la vallée ouvrant sur une plaine. Dans la seconde configuration, la vallée est composée de deux sections, la section aval, qui ouvre sur la plaine, étant plus étroite que la section amont. Cette configuration est appelée "pooling". Quelle que soit la configuration, la variation du profil vertical de température le long de l'axe de la vallée conduit à d'un gradient de pression horizontal générant un vent de vallée. Les configurations de type "pooling" sont associées à des concentrations de polluants plus élevées que dans le cas où la largeur de la vallée ne varie pas, car le vent de vallée est plus faible dans la configuration "pooling". L'impact de la largeur de la section aval sur la concentration de polluants dans la section amont est remarquable: la ventilation de la section amont peut être complètement bloquée pendant la majeure partie de la simulation (six heures) lorsque le rapport de largeur de la section amont à la section aval est égal à 10.La deuxième partie concerne une situation réelle en vallée alpine, durant un épisode fortement pollué de février 2015. La dynamique atmosphérique et la concentration de particules fines (PM) ont été modélisées dans une section de la vallée de l'Arve autour de la ville de Passy par le modèle WRF-Chem durant cet épisode. Pour cela, le cadastre d'émission préparé par l'agence de qualité de l'air de la région Auvergne Rhône-Alpes a été implémenté dans le modèle. La dynamique de la couche limite a été comparée à des profils verticaux de vent et température mesurés lors d'une campagne de terrain durant l'épisode pollué, et la concentration de PM comparée aux données de stations de qualité de l'air. La contribution des vallées tributaires à la pollution atmosphérique dans cette section de la vallée de l'Arve a d'abord été étudiée. Les résultats montrent que cette contribution est très faible par rapport à celle des sources locales; ceci implique que les émissions $in-situ$ sont principalement responsables de la forte concentration de particules polluantes enregistrée dans la vallée car le fond de vallée est découplé de l'atmosphère au-dessus de la vallée. Les processus de ventilation et leur influence sur la concentration de particules (PM) ont ensuite été étudiés. L'analyse montre qu'en raison de la ventilation limitée autour de la ville de Passy, la variabilité horaire de la concentration de PM dans la ville est contrôlée par les émissions. / During wintertime anticyclonic regimes, urbanized mountain areas often experience stable and decoupled atmospheric conditions, resulting in severe episodes of particulate air pollution. This study deals with the characterization of drivers of particulate air pollution in deep alpine valleys under such stable conditions from the point of view of atmospheric dynamics. The work has been carried out through the implementation of both idealized and real case numerical simulations using the Weather Research and Forecasting (WRF) model. Particulate air pollution has been modelled by implementing passive tracers in the simulations.In a first part, high-resolution numerical simulations have been performed using two different configurations of three-dimensional idealized valleys opening onto a plain. The first configuration corresponds to a valley of constant width in the along-valley direction, directly opening onto a plain. The second one consists in a valley of varying width in that direction, with an upstream section of larger width than the downstream section, which opens onto the plain. The latter configuration is referred to as a pooling case. The change in the vertical temperature profile along the valley axis in all configurations results in a horizontal pressure gradient that leads to the development of an along-valley flow. Configurations in which the upstream section of the valley opens onto a narrower one have shown to be prone to higher pollutant concentrations than in the case of a valley opening directly onto a plain, due to the weaker down-valley flow developing in the variable width configurations. The impact of the downstream valley section on the concentration of pollutants in the upstream valley section is remarkable: the ventilation of the upstream valley section can be completely blocked for most of the simulation (six hours) when the ratio of the upstream to downstream valley section width is equal to 10.The second part deals with a real case situation, focusing on a section of the alpine Arve River valley around the city of Passy. The numerical modelling of the atmospheric dynamics and particulate matter (PM) concentration in that section has been performed during a strongly polluted wintertime episode of February 2015. For this purpose the emission inventory prepared by the air quality agency of the région Auvergne Rhône-Alpes has been implemented in the WRF-Chem model. The boundary layer dynamics has been compared to wind and temperature profiles collected during a field campaign conducted during the episode, and PM concentration has been compared to data recorded by air quality stations. The contribution of the valley tributaries to air pollution within the section considered has been first studied. Results show that this contribution is very small compared to that from local sources; this implies that $in-situ$ emissions are primarily responsible for the high PM concentration recorded in the valley since the valley bottom is decoupled from the atmosphere above the valley. The ventilation characteristics of the valley and their influence on the particulate matter (PM) concentration have next been investigated. The analysis indicates that due to the limited ventilation around the city of Passy, the hourly variability of PM concentration therein is driven by that of the emissions.

En période hivernale, le golfe du Lion est sujet à des conditions de vents régionaux forts (Mistral et/ou Tramontane) qui transportent des masses d'air continentales froides et sèches au dessus de la mer. Ces événements, les Cold Air Outbreaks (CAO) , conduisent à d'intenses échanges air-mer et donc à un pompage de chaleur qui favorise la formation d'eaux denses et le déclenchement de la convection océanique profonde. La bonne représentation de ces échanges air-mer intenses dans les modèles de climat et de prévision numérique du temps reste à l'heure actuelle une problématique majeure. Elle est au cœur du projet ASICS-MED centré sur compréhension des mécanismes de formation d'eaux denses en Méditerranée et qui s'inscrit dans le cadre de la thématique " Échanges air-mer intenses " du programme HyMeX dédié à l'étude du cycle de l'eau en Méditerranée. Les processus qui s'opèrent au sein de la couche limite atmosphérique marine (CLAM) et de la couche de mélange océanique (CMO) interagissent entre eux à différentes échelles spatiales et temporelles. La compréhension de l'évolution globale de la CLAM mais également des mécanismes locaux nécessitent la prise en compte de l'ensemble des processus. L'étude présentée ici est consacrée à la structure moyenne et turbulente de la CLAM en conditions de vents forts. L'objectif est de déterminer comment l'organisation du champ turbulent est impactée lors d'épisodes de CAO et d'estimer les flux de surface associés à ces conditions de vents forts. La méthodologie adoptée est basée sur l'utilisation conjointe d'observations aéroportées collectées lors de la campagne de mesure SOP2 d'HyMeX et de simulations numériques. La campagne de mesure SOP2 d'HyMeX qui a eu lieu au cours de l'hiver 2013 dans le golfe du Lion a permis de documenter grâce à l'avion de recherche ATR42 la structure moyenne et turbulente de la CLAM lors de 11 épisodes de CAO. Une analyse spectrale s'appuyant sur un modèle analytique a été réalisé sur 181 paliers (i.e. segments de vol rectilignes et stabilisés en altitude). Les profils verticaux des échelles turbulentes caractéristiques ainsi que la forme du spectre de la vitesse verticale ont permis de mettre en évidence un allongement des structures énergétiques dans l'axe du vent moyen associé à l'organisation du champ turbulent sous la forme de rouleaux longitudinaux. Une description unidirectionnelle du champ turbulent tridimensionnel peut conduire à une représentativité limitée des structures cohérentes au sein des échantillons. Cependant, la connaissance des profils de flux sur toute l'épaisseur de la CLAM est nécessaire pour l'estimation des échanges air-mer. Une méthode de correction des flux turbulents calculés par eddy correlation a été appliqué afin de prendre en compte les erreurs systématique et aléatoire relatives à la mesure et au traitement de données. Cette correction a permis de déterminer les meilleures estimations possibles des flux extrapolés à la surface avec une marge d'incertitude pour les 11 épisodes de CAO documentés lors de la campagne SOP2 d'HyMeX. La comparaison de ces estimations aéroportées aux autres sources d'information dérivées de paramétrisations des flux a permis de mettre en évidence une sous-estimation systématique du flux de chaleur latente en conditions de vents forts. Une approche numérique a permis de compléter l'analyse de la structure moyenne et turbulente de la CLAM lors d'épisodes de CAO. / During winter, local strong winds (Mistral or Tramontana) occurred in the Gulf of Lion which bring cold and dry continental air over a warmer sea. Those events, the cold air outbreaks, can lead to intense air-sea interactions which favour dense water formation and deep oceanic convection. The representation of air-sea exchanges is a fundamental aspect of of climate modelling and numerical weather forecasting. The ASICS-MED project aims to identify fine-scale processes leading to dense water formation and is a part of the "Intense air-sea exchanges" topic of the HyMeX program devoted to hydrological cycle in the Mediterranean. The processes occurring within the marine atmospheric boundary layer (MABL) and the oceanic mixing layer (ML) interact with one another at different spatial and temporal scales. Understanding the overall evolution of the MABL but also the local mechanisms requires taking into account all the processes. The study presented here is devoted to the mean and turbulent structure of the MABL under strong wind conditions. The objective is to determine how the organization of the turbulent field is impacted during CAO events and to estimate the surface fluxes associated with these strong wind conditions. The methodology adopted is based on the joint use of airborne observations collected during the HyMeX-SOP2 field campaign and numerical simulations. The HyMeX-SOP2 field campaign took place during the winter of 2013 in the Gulf of Lion. The research aircraft ATR42 was operated to document the mean and turbulent structure of the MABL during 11 CAO events. A spectral analysis based on an analytic model was carried out on 181 legs (i.e. stacked straight and level runs stabilized in altitude). The vertical profiles of the turbulent characteristic scales as well as the shape of the vertical velocity spectrum revealed an elongation of the energy structures in the mean wind direction associated with the organization of the turbulent field into longitudinal rolls. A unidirectional sampling of the three-dimensional turbulent field may lead to a limited representativeness of the coherent structures within the samples. However, knowledge of kinematic fluxes profiles over the entire thickness of the CLAM is necessary to estimate air-sea exchanges. A correction method was applied to turbulent fluxes calculated by eddy correlation in order to take into account systematic and random errors related to measurement and data processing. This correction made it possible to determine the best possible estimates of the extrapolated surface fluxes with a margin of uncertainty for the 11 CAO events documented during the HyMeX-SOP2 field campaign. The comparison of these airborne estimates with the other sources of information derived from bulk parameterizations show a systematic underestimation of the latent heat flux under strong wind conditions. A numerical approach allowed to complete the analysis of the mean and turbulent structure of the MABL during CAO events. The numerical study, based on the non-hydrostatic Meso-NH model, focuses on an episode of strong Tramontana with winds greater than 25m/s documented during the HyMeX-SOP2 field campaign. In a first step, a one-dimensional framework made it possible to understand the forcing terms necessary to reproduce in a realistic way the development of the observed MABL. This reference configuration allowed, in a second time, a Large-Eddy Simulation of the CAO event. This simulation has been validated using airborne data and has allowed to deepen the description of the turbulent field as well as the evolution of the coherent structures oriented in the axis of the mean wind.

Understanding the interactions between ocean and atmosphere in regions of oceanographic fronts is of vital importance for the improvement of numerical models for weather and climate forecasting. In the South Atlantic Ocean (SAO) the meeting between the warm waters of the Brazil Current (BC) and the cold waters of the Malvinas (Falkland) Current (MC) in the region known as the Brazil-Malvinas Confluence (BMC), results in intense mesoscale oceanic activity and, for this reason, this region is considered one of the most energetic of the Global Ocean. The interactions resulting from the thermal contrast in regions oceanographic fronts of the OAS are investigated in this work through estimates of heat fluxes based on data collected in situ and by satellite. The results of this study show that the response to the thermal contrasts found in the ocean is in the form of heat fluxes and these fluxes are critical in modulating the atmospheric boundary layer (ABL). Estimation based on data collected in situ show that in the warm side (north) of the oceanographic front the fluxes are more intense (latent heat: 62 W/m² and sensible heat: 0.6 W/m²) than in the cold side (south) (latent heat: 5.8 W/m² and sensible heat: -13.8 W/m²). In the South Atlantic Current (SAC) along the 30° S parallel, heat fluxes are directly related to the meandering characteristic of the current. The data collected in situ, in addition to allow heat flux estimates at a better spatial resolution, were used to develop a new method for estimating the heat energy exchanged between the atmosphere and the ocean caused by the presence of mesoscale oceanic structures. This methodology consists in the comparison of a radiosonde profile taken over waters of the structure of interest and another taken over waters which do not belong to this structure. The methodology was used to estimate the heat energy transfer between the atmosphere and the ocean over the top of three structures sampled in the OAS. The estimation of the heat energy transferred by a warm eddy detached from the BC points to an energy in the latent (sensible) form of 1.6 1017 J (-2.8 1016 J) which corresponds to approximately 0.011 % of the total heat energy of the eddy transferred to the atmosphere during the field experiment and 0.78 % transferred during the supposed lifetime of the eddy (3 months). Along the CSA two oceanic structures were studied: (i) a cold meander that receives from the atmosphere energy in the latent (sensible) form of 1.4 106 J/m2 (5.4 105 J/m2), and (ii) warmer waters associated with a detached eddy from the Agulhas Current (AC) that transfer to the atmosphe heat energy of approximately 4 106 J/m2 an 5.7 106 J/m2 in the latent and sensible forms, respectively. The estimation of heat energy transfer on top of mesoscale oceanic structures clearly demonstrate the importance of these structures for the heat exchanges between the ocean and the atmosphere and must be taken into account in future works about this subject in the SAO. / A compreensão das interações entre oceano e atmosfera em regiões de frentes oceanográficas é de vital importância para o melhoramento de modelos numéricos de previsão do tempo e clima. No Oceano Atlântico Sul (OAS) o encontro entre as águas quentes da Corrente do Brasil (CB) com as águas frias da Corrente das Malvinas (CM), na região denominada Confluência Brasil-Malvinas (CBM), resulta em intensa atividade oceânica de mesoescala e, por esse motivo, essa região é considerada uma das mais energéticas do Oceano Global. As interações resultantes do contraste termal ao longo de regiões de frentes oceanográficas no OAS são investigadas neste trabalho através de estimativas de fluxos de calor baseadas em dados de satélite e dados coletados in situ. Os resultados do trabalho demonstram que a resposta aos contrastes termais encontrados no oceano se dá na forma de fluxos de calor e que esses fluxos são fundamentais na modulação da Camada Limite Atmosférica (CLA). As estimativas com base em dados coletados in situ demonstram que no lado quente (norte) da frente oceanográfica os fluxos são mais intensos (calor latente: 62 W/m² e calor sensível: 0,6 W/m²) que nos lado frio (sul) (calor latente: 5,8 W/m² e calor sensível: -13,8 W/m²). Na Corrente Sul Atlântica (CSA), ao longo do paralelo de 30° S, os fluxos de calor estão diretamente relacionados a característica meandrante da corrente. Os dados coletados in situ, além de possibilitarem estimativas de fluxo de calor com uma melhor resolução espacial, foram usados no desenvolvimento de uma nova metodologia para estimativa da energia calorífica trocada entre oceano e atmosfera em virtude da presença de estruturas oceânicas de mesoescala. Essa metodologia consiste na comparação entre um perfil de radiossonda tomado sobre águas da estrutura de interesse e outro tomado sobre águas que não pertencem a essa estrutura. A metodologia desenvolvida foi utilizada para determinar a transferência de energia calorífica entre oceano e atmosfera em três estruturas amostradas no OAS. A estimativa da energia calorífica transferida por um vórtice quente desprendido da CB aponta para uma energia na forma latente (sensível) de 1,6 1017 J (-2,8 1016 J) que corresponde a aproximadamente 0,011 % da energia calorífica total do vórtice transferida durante o experimento de campo e de 0,78 % da energia do vórtice transferidos durante o tempo suposto de vida do vórtice (3 meses). Ao longo da CSA, duas estruturas oceânicas foram estudadas: (i) um meandro frio que recebe da atmosfera uma energia na forma latente (sensível) de 1,4 106 J/m2 (5,4 105 J/m2) e (ii) águas mais quentes associadas a um vórtice desprendido da Corrente das Agulhas (CA) que transferem para a atmosfera uma energia calorífica de aproximadamente 4 106 J/m2 e 5,7 106 J/m2 nas formas latente e sensível, respectivamente. As estimativas da transferência de energia calorífica sobre estruturas oceânicas de mesoescala demonstram claramente a importância destas nas trocas de calor entre o oceano e a atmosfera e devem ser levadas em consideração em trabalhos futuros sobre o tema no OAS.

Doppler radar wind profilers are already widely used to measure atmospheric winds throughout the free troposphere and stratosphere. Several methods have been developed to quantify atmospheric turbulence with such radars, but to date they have remained largely un-tested; this thesis presents the first comprehensive validation of one such method. Conventional in-situ measurements of turbulence have been concentrated in the surface layer, with some aircraft and balloon platforms measuring at higher altitudes on a case study basis. Radars offer the opportunity to measure turbulence near continuously, and at a range of altitudes, to provide the first long term observations of atmospheric turbulence above the surface layer. Two radars were used in this study, a Mesosphere-Stratosphere-Troposphere (MST) radar, at Capel Dewi, West Wales, and the Facility for Ground Based Atmospheric Measurements (FGAM) mobile boundary layer profiler. In-situ measurements were made using aircraft and tethered-balloon borne turbulence probes. The spectral width method was chosen for detailed testing, which uses the width of a radar's Doppler spectrum as a measure of atmospheric velocity variance. Broader Doppler spectra indicate stronger turbulence. To obtain Gaussian Doppler spectra (a requirement of the spectral width method), combination of between five and seven consecutive spectra was required. Individual MST spectra were particularly non-Gaussian, because of the sparse nature of turbulence at its observation altitudes. The width of Gaussian fits to the Doppler spectrum were compared to those from the `raw' spectrum, to ensure that non-atmospheric signals were not measured. Corrections for non-turbulent broadening, such as beam broadening, and signal processing, were investigated. Shear broadening was found to be small, and the errors in its calculation large, so no corrections for wind shear were applied. Beam broadening was found to be the dominant broadening contribution, and also contributed the largest uncertainty to spectral widths. Corrected spectral widths were found to correlate with aircraft measurements for both radars. Observing spectral widths over time periods of 40 and 60 minutes for the boundary layer profiler and MST radar respectively, gave the best measure of turbulence intensity and variability. Median spectral widths gave the best average over that period, with two-sigma limits (where sigma is the standard deviation of spectral widths) giving the best representation of the variability in turbulence. Turbulent kinetic energies were derived from spectral widths; typical boundary layer values were 0.13 m 2.s (-2) with a two-sigma range of 0.04-0.25 m 2.s (-2), and peaked at 0.21 m 2.s (-2) with a two-sigma range of 0.08-0.61 m 2.s (-2). Turbulent kinetic energy dissipation rates were also calculated from spectral widths, requiring radiosonde measurements of atmospheric stability. Dissipation rates compared well width aircraft measurements, reaching peaks of 1x10 (-3) m 2.s (-3) within 200 m of the ground, and decreasing to 1-2x10 (-5) m 2.s (-3) near the boundary layer capping inversion. Typical boundary layer values were between 1-3x10 (-4) m 2.s (-3). Those values are in close agreement with dissipation rates from previous studies.

Title: Urban Ventilation Dependence on Geometric Configuration Author: RNDr. Ing. Libor Kukačka Department: Department of Atmospheric Physics Supervisor: prof. RNDr. Zbyněk Jaňour, DrSc., Academy of Sciences of the Czech Republic, Institute of Thermomechanics, v. v. i. Abstract: The main goal of the thesis is to investigate the impact of urban geometry on the urban ventilation using wind-tunnel modelling. To measure the pollutant transport, both advective and turbulent, within complex urban geometries with a high temporal resolution a special measurement method was developed. At first, the pollution of a simplified urban area was simulated by a ground-level point source and the ventilation of the intersection with respect to four wind directions was studied. Later, the pollution of other simplified and complex urban areas was simulated by a ground-level line source and the ventilation of three different street canyons with respect to wind direction perpendicular and oblique to their along-canyon axis was investigated. The clear impact of urban complexity and wind direction on street canyon ventilation is demonstrated at lateral and top openings of all investigated canyons and the intersection. Whilst the dominance of the pollutant advection is demonstrated at the eaves of pitched roofs, at the roof ridges...

In this thesis, we have evaluated a stochastic Lagrangian model for computing particle deposition rates with prospects to use for atmospheric deposition applications.  The model is one-dimensional and models the particle dynamics in the boundary layers near walls and obstacles by simulating the coherent turbulent structures and Brownian motion governing the wall-normal transport. The deposition model is used with a hybrid deterministic/stochastic particle dispersion model governing the dynamics in the turbulent bulk flow. We used a steady-state RANS k-ϵ turbulence model to simulate the turbulent fluid flow in a neutral atmospheric boundary layer (ABL) using the with inflow boundary conditions by Richards & Hoxey (1993). The turbulence model is solved with the SIMPLE algorithm using the OpenFOAM software. The mean-field characteristic of the turbulent flow in the computational domain is exported and used for the particle model. The particle model is a Lagrangian Langevin-type model, consisting of a system of stochastic differential equations. The particle model was solved using a weakly first order a-stable scheme. We evaluated the deposition model by computing the deposition rate for a range of particle sizes and compared our results with collected experimental wind tunnel data. The numerical experiment was done in a computational domain based on the ABL model by Hargreaves & Wright (2007), a rectangular domain with a logarithmic wind profile. We used a particle source near the inflow boundary with an instantaneously release at the initial time. Results showed disagreement with the experimental data and was only valid for medium sized particles. However, time restrictions led to the analysis being cut short and only a single simulation was conducted. A definite conclusion on the suitability of the method could not be made based solely on this single results. Some uncertainties were identified and discussed for further potential work on the evaluation of the method. However, one conclusion was drawn on the performance of the method. The computational cost was concluded to be too high with the first order particle scheme used and higher order schemes is required for any practical use of the method for atmospheric deposition applications.

Accurate forecasts of wind speed and power production are of great value for wind power producers. In Southwest Iceland, wind power installations are being planned by various entities. This study aims to create optimal wind speed and wind power production forecasts for wind power production in Southwest Iceland by applying statistical post-processing methods to a deterministic HARMONIE-AROME forecast at a single point in space. Three such methods were implemented for a 22 month-long set of forecast-observation samples in 1h resolution: Temporal Smoothing (TS), Observational Distributions on Discrete Intervals (ODDI - a relatively simple classification algorithm) and Quantile Regression Forest (QRF - a relatively complicated Machine Learning Algorithm). Wind power forecasts were derived directly from forecasts of wind speed using an idealized power curve. Four different metrics were given equal weight in the evaluation of the methods: Root Mean Square Error (RMSE), Miss Rate of the 95-percent forecast interval (MR95), Mean Median Forecast Interval Width (MMFIW - a metric to measure the forecast sharpness) and Continuous Ranked Probability Score (CRPS). Of the three methods, TS performed inadequately while ODDI and QRF performed significantly better, and similarly to each other. Both ODDI and QRF predict wind speed and power production slightly more accurately than deterministic AROME in terms of their Root Mean Square Error. In addition to an overall evaluation of all three methods, ODDI and QRF were evaluated conditionally. The results indicate that QRF performs significantly better  than ODDI at forecasting wind speed and wind power at wind speeds above 13 m/s. Else, no strong discrepancies were found between their conditional performance. The results of this study are limited by a relatively scarce data set and correspondingly short time series. The results indicate that applying statistical post-processing methods of varying complexity to deterministic wind speed forecasts is a viable approach to gaining a probabilistic insight into the wind power potential at a given location.

The effect of building geometry on the wind environment of cities is such that it can cause problems like wind danger, discomfort and poor ventilation of airborne pollutants. Computational fluid dynamics (CFD) can play a role in assessing changes in wind environment caused by building projects before realisation at little cost. However, the current state-of-the-art methods, RANS and LES, force a steep trade-off between accuracy and computational cost, and neither method is truly predictive. Time-resolved adaptive direct finite element simulation (DFS) is a method for CFD that is predictive and automatically optimises the mesh for a goal quantity, making it both efficient and accurate. In this thesis, DFS was implemented in FEniCS and used on basic validation cases to provide a proof of concept for the use of this method in the building aerodynamics, on resources freely available to anyone. The results show that the method is accurate to within 10% of the validation data with respect to the goal quantity. Visually, the expected flow features are clearly identifiable. DFS was successfully applied to a relatively complicated building geometry, with a total computation time of about 120 core-hours. We conclude that DFS has significant potential as a method for evaluating urban wind environments. Furthermore, because of its ease of use and lack of parameters, DFS can play an important role in helping architects, designers and students understand the effect of urban geometries on the wind environment. This report provides a basis for further research on DFS for building aerodynamics, as validation on more diverse urban geometries is still necessary. / Effekten av byggnaders form och geometri är så viktig att den kan ge problem för ventilation av t.ex. föroreningar, för energieffektivitet, och för vindfaror med t.ex. hög vindhastihet som kan vara farligt eller skapa obehag. Beräkningsströmningsdynamik (CFD) kan ha en roll i bedömningen av byggnadsprojekt i ett tidigt skede till liten kostnad. Dock är de etablerade och ledande metodikerna, RANS och LES, inte prediktiva och tvingar fram en kompromiss mellan beräkningskosnad och noggrannhet. Vår metodik “Time-resolved adaptive direct finite element simulation” (DFS) är en metod för CFD som är prediktiv och automatiskt optimerar beräkningsnätet (och därmed beräkningskostnaden) för en given målkvantitet, som ger både effektivitet och noggrannhet. I denna avhandling implementerades DFS i FEniCS och användes i grundläggande valideringsfall för att ge ett proof of conceptför användning av denna metod i byggnadsaerodynamik, på resurser som är fritt tillgängliga för alla. Resultaten visar att metoden är korrekt inom 10% av valideringsdata med avseende på målkvantiteten. Visuellt är de förväntade flödesfunktionerna tydligt identifierbara. DFS applicerades framgångsrikt på en relativt komplicerad byggnadsgeometri med en total beräkningstid på cirka 120 kärntimmar, vilket är en försumbar kostnad. Vi drar slutsatsen att DFS har en betydande potential som metod för utvärdering av stadsvindmiljöer. Dessutom, på grund av dess användarvänlighet och frihet från parametrar, kan DFS spela en viktig roll för att hjälpa arkitekter, designers och studenter att förstå effekterna av stadsgeometrier på vindmiljön. Denna rapport ger en grund för vidare forskning om DFS för aerodynamik, eftersom validering av mer olika stadsgeometrier fortfarande är nödvändig.

La physique de la couche limite atmosphérique en domaine océanique est principalement régie par les processus couplés liés au vent, à l’état de mer local, et à des effets de flottabilité. Leur compréhension reste néanmoins parcellaire et leurs descriptions théoriques et stochastiques sont pour le moins lacunaires, lorsqu’elles ne sont tout simplement pas mises à mal par les rares observations. Dans un contexte d’exploitation croissante de la ressource éolienne offshore, la mise en place de méthodes numériques visant à une description plus fine des propriétés turbulentes de cette couche limite sera une étape déterminante dans la réduction des coûts et l’optimisation des structures pour des rendements de récupération d’énergie améliorés. Ainsi, un outil numérique a été mis en place afin d’étudier le couplage entre un écoulement atmosphérique et l’état de mer. Un code Large-Eddy Simulation massivement parallèle pour la simulation des écoulements atmosphériques incompressibles développé par P. Sullivan au National Center for Atmospheric Research est couplé à un code spectral d’états de mer non-linéaires développé au Laboratoire de recherche en Hydrodynamique, Energétique et Environnement Atmosphérique. De nombreuses configurations de vents et d’états de mer sont modélisées. On montre que les lois semi empiriques souvent utilisées pour représenter la distribution verticale de la vitesse moyenne du vent sont une bonne approximation dans les situations où un petit état de mer est soumis à un fort vent. Néanmoins, dans le cas de houles très rapides se propageant dans des zones de faible vent, la création d’un jet de vent par la houle invalide ces lois semi-empiriques. / Modelling the dynamic coupling of ocean-atmosphere systems requires a fundamental and quantitative understanding of the mechanisms governing the windwave interactions: despite numerous studies, our current understanding remains quite incomplete and, in certain conditions, sparse field observations contradict the usual theoretical and stochastic models. Within the context of a growing exploitation of the offshore wind energy and the development of met ocean models, a fine description of this resource is a key issue. Field experiments and numerical modelling have revealed that atmospheric stability and wave effects, including the dynamic sea surface roughness, are two major factors affecting the wind field over oceans. A numerical tool has been implemented in order to study the coupling between an atmospheric flow and the seastate. A massively parallel large-eddy simulation developed by P. Sullivan at the National Center for Atmospheric Research is then coupled to a High-Order Spectral wave model developed at the Hydrodynamics,Energetics & Atmospheric Environment Laboratory in Ecole Centrale de Nantes. Numerous configurations of wind and sea states are investigated. It appears that, under strongly forced wind conditions above a small sea state, the semi-empirical laws referred to as standards in the international guidelines are a good approximation for the vertical profile of the mean wind speed. However, for light winds overlying fast-moving swell, the presence of a wave induced wind jet is observed, invalidating the use of such logarithmic laws.