Afin de mieux comprendre les mécanismes liés à l’érosion du sol sous l’effet du vent, le transport de particules solides dans un écoulement de couche limite turbulente à l’échelle d’une soufflerie est étudié à l’aide de simulations numériques. La présence d’une ou plusieurs collines Gaussiennes au sol permet d’étudier les effets de la topographie sur le transport, le dépôt et la réémission de particules solides. L’écoulement du fluide porteur est résolu par la Simulation des Grandes Échelles (SGE). Des modèles de paroi pour la vitesse du fluide sont implémentés afin de mieux représenter l’écoulement proche d’une colline. Le mouvement des particules est pris en compte par un suivi Lagrangien. Des modèles d’envol et de rebond sont développés et utilisés pour prendre en compte l’émission et l’impact au sol des particules. Dans la première partie, l’écoulement au-dessus de collines transversales est simulé et validé par des comparaisons avec différentes expériences. Selon Oke [1988], l’écoulement dans la canopée urbaine peut être schématiquement caractérisé par différents régimes en fonction du positionnement relatif des obstacles. Ce concept est appliqué au cas des dunes, assimilées à des collines dans notre étude. L’accent est mis sur la zone de recirculation (ZR) formée derrière ces collines. Les variations de la ZR sont examinées en fonction de différents paramètres dont la configuration des collines et le nombre de Reynolds. De plus, une étude portant sur la sous couche rugueuse est effectuée de façon à déterminer l’effet de la rugosité due à la couche de particules solides au sol. La seconde partie du travail porte sur la simulation des particules au-dessus des collines. L’objectif est l’amélioration des modélisations concernant l’envol, le rebond et le couplage entre le fluide et les particules. Un premier travail de validation est réalisé en utilisant le modèle complet de transport des particules solides. En particulier, l’évolution du flux d’émission des particules, estimé par le modèle d’envol, en fonction du nombre de Shields, donne des résultats comparables aux modèles classiques de saltation et aux expériences de la littérature. Au-dessus des collines, le transport des particules solides est étudié par des profils de concentration et de vitesse moyenne. Pour analyser les résultats, deux cartographies sont réalisées. La première donne l’intensité des événements locaux et instantanés qui seraient à l’origine de l’évacuation des particules piégées au sein de la ZR. La seconde montre la distribution des particules déposées au sol. Ces résultats permettent d’identifier des zones sujettes à l’érosion et à l’accumulation autour des collines. Enfin, les flux des particules piégées et déposées à l’intérieur de la ZR sont quantifiés et comparés aux flux des particules émises en amont. Ces flux, bien que faibles par rapport au flux entrant, contribueraient aux migrations des dunes et à l’avancée des déserts. / The transport of solid particles inside a laboratory-scale turbulent boundary-layer is studied by numerical simulations, to obtain a better understanding of the mechanisms associated with wind erosion of soil. The presence of one or several Gaussian hills allows a study of the topographic effects on the transport, deposition and re-emission of solid particles. The carrier fluid motion is resolved in a Large Eddy Simulation (LES). Wall models are implemented to better account for the effects of turbulent flow near the terrain. Particle trajectories are calculated using a Lagrangian tracking. Take-off and rebound models are developed in order to take into account particle emissions and impacts at the wall. In the first part, the flow over transversal Gaussian hills is simulated and validated by comparison with different experiments. According to Oke [1988], the flow inside an urban canopy can be schematically characterised into different flow regimes depending on the relative localisation of the obstacles at the ground. This concept is applied to the case of sand dunes, assimilated to 2D hills in this study. The focus is on the recirculation zone (RZ) on the lee side, which has the characteristic of increasing the residence time and the interaction fluid/particle in general, particle trapping and deposition in particular. The variations of RZ with different hill geometries and Reynolds numbers are examined. A study on the roughness sublayer is conducted in order to determine the roughness effects due to the layer of solid particles on the wall. The second part of the work is devoted to the simulation of solid particle transport over the Gaussian hills. The objective is to improve the modelling of particle take-off, rebound and the two-way coupling between the fluid and the particle. A first work of validation is conducted by using the complete model of solid particle transport developed in this thesis. In particular, the evolution of particle emission flux predicted by the take-off model is in accordance with classical saltation models and experiments from the literature. Over the Gaussian hills, analysis of particle transport is conducted using concentration and mean velocity fields. Two mappings are realised. The first indicates the intensity of the local and instantaneous flow structures that arguably regulate the re-entrainment of particles trapped inside the RZ. The second shows the accumulation of particles on the wall. These results highlight zones prone to wind erosion and particle deposition around the hills. Last but not least, the fluxes of particle trapping and deposition inside the RZ are quantified and compared to the incoming flux from upstream. These fluxes, albeit relatively weak in comparison to the incoming one, contribute potentially to dune migrations and desertification.

Cette thèse s'inscrit dans un projet de modélisation numérique de la dispersion atmosphérique de polluants à travers le code de mécanique des fluides numérique Code_Saturne. L'objectif est de pouvoir simuler la dispersion atmosphérique de polluants en environnement complexe, c'est-à-dire autour de centrales, sites industriels ou en milieu urbain. Dans ce contexte, nous nous concentrons sur la modélisation de la dispersion des polluants à micro-échelle, c'est-à-dire pour des distances de l'ordre de quelques mètres à quelques kilomètres et correspondant à des échelles de temps de l'ordre de quelques dizaines de secondes à quelques dizaines de minutes : on parle de modélisation en champ proche. L’approche suivie dans ces travaux de recherche suit une formulation hybride eulérienne/lagrangienne, où les champs dynamiques moyens relatifs au fluide porteur (pression, vitesse, température, turbulence) sont calculés via une approche eulérienne et sont ensuite fournis au solveur lagrangien. Ce type de formulation est couramment utilisé dans la littérature atmosphérique pour son efficacité numérique. Le modèle lagrangien stochastique considéré dans nos travaux est le Simplified Langevin Model (SLM), développé par Pope (1985,2000). Ce modèle appartient aux méthodes communément appelées méthodes PDF (Probability Density Function), et, à notre connaissance, n'a pas été exploité auparavant dans le contexte de la dispersion atmosphérique. Premièrement, nous montrons que le SLM respecte le critère dit de mélange homogène (Thomson, 1987). Ce critère, essentiel pour juger de la bonne qualité d'un modèle lagrangien stochastique, correspond au fait que si des particules sont initialement uniformément réparties dans un fluide incompressible, alors elles doivent le rester. Nous vérifions le bon respect du critère de mélange homogène pour trois cas de turbulence inhomogène représentatifs d'une large gamme d'applications pratiques : une couche de mélange, un canal plan infini, ainsi qu'un cas de type atmosphérique mettant en jeu un obstacle au sein d'une couche limite neutre. Nous montrons que le bon respect du critère de mélange homogène réside simplement en la bonne introduction du terme de gradient de pression en tant que terme de dérive moyen dans le modèle de Langevin (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Nous discutons parallèlement de l'importance de la consistance entre champs eulériens et lagrangiens dans le cadre de telles formulations hybrides eulériennes/lagrangiennes. Ensuite, nous validons le modèle dans le cas d'un rejet de polluant ponctuel et continu, en conditions de vent uniforme et turbulence homogène. Dans ces conditions, nous disposons en effet d'une solution analytique nous permettant une vérification précise. Nous observons que dans ce cas, le modèle lagrangien discrimine bien les deux différents régimes de diffusion de champ proche et champ lointain, ce qui n'est pas le cas d'un modèle eulérien à viscosité turbulente (Bahlali et al., 2018b).Enfin, nous travaillons sur la validation du modèle sur plusieurs campagnes expérimentales en atmosphère réelle, en tenant compte de la stratification thermique de l'atmosphère et de la présence de bâtiments. Le premier programme expérimental considéré dans nos travaux concerne le site du SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), dans la banlieue sud de Paris, et met en jeu une stratification stable de la couche limite atmosphérique. La seconde campagne étudiée est l'expérience MUST (Mock Urban Setting Test). Réalisée aux Etats-Unis, dans le désert de l'Utah, cette expérience a pour but de représenter une ville idéalisée, au travers d'un ensemble de lignées de conteneurs. Deux rejets ont été simulés et analysés, respectivement en conditions d'atmosphère neutre et stable (Bahlali et al., 2018a) / This Ph.D. thesis is part of a project that aims at modeling pollutant atmospheric dispersion with the Computational Fluid Dynamics code Code_Saturne. The objective is to simulate atmospheric dispersion of pollutants in a complex environment, that is to say around power plants, industrial sites or in urban areas. In this context, the focus is on modeling the dispersion at micro-scale, that is for distances of the order of a few meters to a few kilometers and corresponding to time scales of the order of a few tens of seconds to a few tens of minutes: this is also called the near field area. The approach followed in this thesis follows a hybrid Eulerian/Lagrangian formulation, where the mean dynamical fields relative to the carrier fluid (pressure, velocity, temperature, turbulence) are calculated through an Eulerian approach and are then provided to the Lagrangian solver. This type of formulation is commonly used in the atmospheric literature for its numerical efficiency. The Lagrangian stochastic model considered in our work is the Simplified Langevin Model (SLM), developed by Pope (1985,2000). This model belongs to the methods commonly referred to as PDF (Probability Density Function) methods, and, to our knowledge, has not been used before in the context of atmospheric dispersion. First, we show that the SLM meets the so-called well-mixed criterion (Thomson, 1987). This criterion, essential for any Lagrangian stochastic model to be regarded as acceptable, corresponds to the fact that if particles are initially uniformly distributed in an incompressible fluid, then they must remain so. We check the good respect of the well-mixed criterion for three cases of inhomogeneous turbulence representative of a wide range of practical applications: a mixing layer, an infinite plane channel, and an atmospheric-like case involving an obstacle within a neutral boundary layer. We show that the good respect of the well-mixed criterion lies simply in the good introduction of the pressure gradient term as the mean drift term in the Langevin model (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Also, we discuss the importance of consistency between Eulerian and Lagrangian fields in the framework of such Eulerian/Lagrangian hybrid formulations. Then, we validate the model in the case of continuous point source pollutant dispersion, under uniform wind and homogeneous turbulence. In these conditions, there is an analytical solution allowing a precise verification. We observe that in this case, the Lagrangian model discriminates well the two different near- and far-field diffusion regimes, which is not the case for an Eulerian model based on the eddy-viscosity hypothesis (Bahlali et al., 2018b).Finally, we work on the validation of the model on several experimental campaigns in real atmosphere, taking into account atmospheric thermal stratification and the presence of buildings. The first experimental program considered in our work has been conducted on the `SIRTA' site (Site Instrumental de Recherche par Télédétection Atmosphérique), in the southern suburb of Paris, and involves a stably stratified surface layer. The second campaign studied is the MUST (Mock Urban Setting Test) experiment. Conducted in the United States, in Utah's desert, this experiment aims at representing an idealized city, through several ranges of containers. Two cases are simulated and analyzed, respectively corresponding to neutral and stable atmospheric stratifications (Bahlali et al., 2018a)

O formato do perfil de velocidades do vento varia de acordo com as características locais da superfície terrestre e de rugosidade do terreno, parâmetros que definem o perfil da Camada-Limite Atmosférica (CLA). As características do escoamento do ar atmosférico sobre e ao redor de acidentes geográficos, tais como morros e colinas, são de grande interesse para aplicações relacionadas à Engenharia de Turbinas e Parques Eólicos. No topo de morros, ocorre a aceleração do vento, fenômeno que pode representar um fator decisivo para a instalação de aerogeradores. Este trabalho dedica-se ao estudo do comportamento da CLA como função da inclinação e rugosidade superficial da elevação, fazendo uso da Dinâmica de Fluidos Computacional (CFD) para construir perfis de velocidade do vento e de intensidade de turbulência. O problema de fechamento das Equações Médias de Reynolds (RANS) é contornado com o uso do modelo de turbulência k-ω SST; os resultados numéricos obtidos são comparados com dados experimentais medidos em túnel de vento sobre modelos em escala dos morros. São testados oito modelos de morros com declives que variam de 25° a 64° para dois tipos de categorias de terreno, em 2D e 3D, e são aplicados dois códigos analíticos para representar o perfil de velocidades de entrada. Resultados numéricos para os perfis de velocidade apresentam diferença inferior a 4% em relação aos respectivos dados obtidos experimentalmente. Os perfis de intensidade de turbulência apresentam diferença máxima na casa dos 7% em comparação aos dados experimentais, o que é explicado pelo fato de que não é possível inserir o perfil de entrada de intensidade de turbulência nas simulações numéricas. Em alternativa, foi usado um valor constante resultado da média dos valores dos perfis usados no túnel de vento. Os modelos de morro em 3D apresentam maior concordância nos resultados de velocidade que os modelos em 2D e que ademais quanto maior é a inclinação do morro maior é a concordância com as medições experimentais. / The shape of the wind velocity profile changes according to local features of terrain shape and roughness, which are parameters responsible for defining the Atmospheric Boundary Layer (ABL) profile. Air flow characteristics over and around landforms, such as hills, are of considerable importance for applications related to Wind Farm and Turbine Engineering. The air flow is accelerated on top of hills, which can represent a decisive factor for Wind Turbine placement choices. The present work focuses on the study of ABL behavior as a function of slope and surface roughness of hill-shaped landforms, using the Computational Fluid Dynamics (CFD) to build wind velocity and turbulent intensity profiles. Reynolds-Averaged Navier-Stokes (RANS) equations are closed using the SST k-ω turbulence model; numerical results are compared to experimental data measured in wind tunnel over scale models of the hills under consideration. Eight hill models with slopes varying from 25° to 64° were tested for two types of terrain categories in 2D and 3D, and two analytical codes are used to represent the inlet velocity profiles. Numerical results for the velocity profiles show differences under 4% when compared to their respective experimental data. Turbulent intensity profiles show maximum differences around 7% when compared to experimental data, this can be explained by not being possible to insert inlet turbulent intensity profiles in the simulations. Alternatively, constant values based on the averages of the turbulent intensity at the wind tunnel inlet were used. The 3D models present greater concordance in the speed results than the 2D models and that in addition the greater the slope of the hill, the greater the agreement with the experimental measurements.

A Lagrangian particle model was applied to estimating emission strength and air pollutant concentrations specifically for the short-range dispersion of an air pollutant in the atmospheric boundary layer. The model performance was evaluated with experimental data. The model was then used as the platform of parametric uncertainty analysis, in which effects of uncertainties in five parameters (Monin-Obukhov length, friction velocity, roughness height, mixing height, and the universal constant of the random component) of the model on mean ground-level concentrations were examined under slightly and moderately stable conditions. The analysis was performed under a probabilistic framework using Monte Carlo simulations with Latin hypercube sampling and linear regression modeling. In addition, four studies related to the Lagrangian particle modeling was included. They are an alternative technique of formulating joint probability density functions of velocity for atmospheric turbulence based on the Koehler-Symanowski technique, analysis of local increments in a multidimensional single-particle Lagrangian particle model using the algebra of Ito integrals and the Wagner-Platen formula, analogy between the diffusion limit of Lagrangian particle models and the classical theory of turbulent diffusion, and evaluation of some proposed forms of the Lagrangian velocity autocorrelation of turbulence.

<p>Accurate modeling of land-atmosphere interaction would help us understand the persistent weather conditions and further contribute to the skill of seasonal climate prediction. In this study, seasonal variations in radiation and precipitation forcing are included in a stochastic soil water balance model to explore the seasonal evolution of soil moisture probabilistic structure. The theoretical results show soil moisture tends to exhibit bimodal behavior only in summer when there are strong positive feedback from soil moisture to subsequent rainfall. Besides the statistical analysis of soil moisture – rainfall feedback, simplified mixed-layer models, coupled with soil-plant-atmosphere continuum, are also used to study heat flux partitioning, cloud initiation, and strength of moist convection. Approximate analytical solutions to the mixed-layer model are derived by applying Penman-Monteith approach, which help explain the roles of equilibrium evaporation and vapor pressure deficit in controlling the diurnal evolution of boundary layer. Results from mixed-layer model also define four regimes for possible convection in terms of cloud/no-cloud formation and low/high convection intensity. Finally, cloud-topped mixed-layer model is developed to simulate the boundary-layer dynamics after the cloud formation, when the evaporative and radiative cooling other than surface heat flux may significantly contribute to the growth of the boundary layer.</p> / Dissertation

No description available.

La dynamique de la couche limite atmosphérique d'une vallée alpine est influencée par le relief environnant et par l’écoulement de grande échelle qui la surmonte. La paramétrisation de cette circulation atmosphérique requiert donc de caractériser finement ces effets. C'est l’objectif de ce travail de thèse : comprendre l’influence du relief environnant une vallée sur les bilans de masse et de chaleur au travers d’une section de cette vallée, par conditions stables et sèches et lorsque le vent synoptique est faible mais non négligeable. Le travail s’appuie sur des simulations numériques.Plusieurs vallées idéalisées ont tout d’abord été considérées: une vallée infiniment longue (bidimensionnelle) et une vallée tridimensionnelle, qualifiée de supérieure, ouvrant soit sur une plaine (cas “vallée-plaine”), soit sur une autre vallée, qualifiée d’inférieure. Cette seconde vallée est soit plus large (cas “drainage”) ou plus étroite (cas “quasi-stagnation”).Dans les vallées tridimensionnelles, deux régimes principaux ont été identifiés, quelle que soit le cas considéré : un régime transitoire, avant que le vent de vallée (descendant) ne se développe, puis un régime quasi-stationnaire, quand le vent de vallée est complètement développé. La présence d’une vallée inférieure réduit la variation de température le long de la vallée, de sorte que le vent de vallée y est plus faible que dans le cas vallée-plaine. En conséquence, la durée du régime transitoire augmente et est maximum pour le cas quasi-stagnation. Lorsque la vallée inférieure est très étroite, la variation de température peut même changer de signe, conduisant à un vent de vallée montant, de la vallée inférieure vers la vallée supérieure. Durant ce régime transitoire, le taux de refroidissement moyenné sur le volume de la vallée est maximum, sa valeur dépendant du cas considéré. En conclusion, les cas drainage et quasi-stagnation conduisent à une couche limite dans la vallée supérieure plus froide et plus profonde que dans le cas vallée-plaine.Dans le régime quasi-stationnaire, le taux de refroidissement moyenné sur le volume de la vallée est plus faible que dans le régime transitoire et varie peu en fonction du cas considéré. En effet, lorsque la vallée inférieure devient plus étroite, le réchauffement lié aux effets advectifs diminue car la vitesse du vent de vallée diminue, de sorte que la contribution (refroidissante) du flux de chaleur sensible diminue également. La conservation de la masse dans la couche limite de la vallée supérieure est assurée par un équilibre entre la convergence des vents de pente au sommet de la couche limite (alimenté par un courant de retour au-dessus (et en sens inverse) du vent de vallée descendant) et la divergence du vent de vallée, les effets de subsidence loin des parois de la vallée jouant un rôle négligeable.Le cas réaliste de la vallée de l’Arve autour de Passy durant une période d’observation intensive de la campagne de mesures PASSY-2015 a permis de caractériser l’impact des vallées environnant Passy sur les bilans de masse et de chaleur dans la vallée. Une couche d’air froid persistante se forme en fond de vallée, suite à l’advection d’air chaud associée au passage d’une crête anticyclonique au-dessus de l’Europe. Les écoulements le long des vallées tributaires présentent une grande variabilité durant la phase persistante de l’épisode, dépendant de la variabilité de l’écoulement à grande échelle, et ont un impact majeur sur l’intensité de la couche d’air froid et la hauteur de l’inversion qui la surmonte. La forte stratification près du sol conduit à leur décollement au-dessus du fond de vallée, les empêchant d'y pénétrer. L’évolution de l’écoulement à grande échelle durant l’épisode a un profond impact sur la dynamique proche du fond de vallée. Durant la nuit en effet, la canalisation de cet écoulement réduit la variation de température le long de la vallée contrôlant le vent de vallée, favorisant la stagnation de l’air. / Alpine valleys are rarely closed systems, implying that the atmospheric boundary layer of a particular valley is influenced by the surrounding terrain and large-scale flows. A detailed characterisation and quantification of these effects is required in order to design appropriate parameterisation schemes for complex terrains. The focus of this work is to improve the understanding of the effects of surrounding terrain (plains, valleys or tributaries) on the heat and mass budgets of the stable boundary layer of a valley, under dry and weak large-scale wind conditions. Numerical simulations using idealised and real frameworks are performed to meet this goal. Several idealised terrains (configurations) were considered: an infinitely long valley (i.e. two-dimensional), and upstream valleys opening either on a plain (valley-plain), on a wider valley (draining) or on a narrower valley (pooling). In three-dimensional valleys, two main regimes can be identified for all configurations: a transient regime, before the down-valley flow develops, followed by a quasi-steady regime, when the down-valley flow is fully developed. The presence of a downstream valley reduces the along-valley temperature difference, therefore leading to weaker down-valley flows. As a result, the duration of the transient regime increases compared to the respective valley-plain configuration. Its duration is longest for pooling configuration. For strong pooling the along-valley temperature difference can reverse, forcing up-valley flows from the narrower towards the wider valley. In this regime, the volume-averaged cooling rate is found maximum and its magnitude dependent on the configuration considered. Therefore pooling and draining induce colder and deeper boundary layers than the respective valley-plain configurations. In the quasi-steady regime the cooling rate is smaller than in the transient regime, and almost independent of the configuration considered. Indeed, as the pooling character is more pronounced, the warming contribution from advection to the heat budget decreases because of weaker down-valley flows, and so does the cooling contribution from the surface sensible heat flux. The mass budget of the valley boundary layer was found to be controlled by a balance between the convergence of downslope flows at the boundary layer top and the divergence of down-valley flows along the valley axis, with negligible contributions of subsidence far from the slopes. The mass budget highlighted the importance of the return current above the down-valley flow, which may contribute significantly to the inflow of air at the top of the boundary layer. A case-study of a persistent cold-air pool event which occurred in February 2015 in the Arve River Valley during the intensive observation period 1 of the PASSY-2015 field campaign, allowed to quantify the effects of neighbouring valleys on the heat and mass budgets of a real valley atmosphere. The cold-air pool persisted because of warm air advection at the valley top, associated with the passage of an upper-level ridge over Europe. The contributions from each tributary valley to the mass and heat budgets of the valley atmosphere were found to vary from day to day within the persistent stage of the cold-air pool, depending on the large-scale flow. Tributary flows had significant impact on the height of the inversion layer and the strength of the cold-air pool, transporting a significant amount of mass within the valley atmosphere throughout the night. The strong stratification of the near-surface atmosphere prevented the tributary flows from penetrating down to the valley floor. The evolution of the large-scale flow during the episode had a profound impact on the near-surface circulation of the valley. The channelling of the large-scale flow at night, can lead to the decrease of the horizontal temperature difference driving the near-surface down-valley flow, favouring the stagnation of the air close to the ground.

A estrutura do vento varia de acordo com as características do terreno e com a rugosidade da superfície terrestre, desacelerando proporcionalmente sua intensidade de acordo com a proximidade do solo, o que determina a constituição da camada limite atmosférica (C.L.A.). As características do escoamento sobre e no entorno de características topográficas obtusas tais como morros são de grande interesse em muitas aplicações, especialmente aquelas ligadas à engenharia de vento. Esta pesquisa foca a investigação da C.L.A. sobre terrenos complexos, analisando a estrutura do escoamento turbulento, a separação e a recomposição do mesmo. Para tanto, dois métodos são empregados e comparados na presente investigação para identificar a influência da orografia complexa sobre o escoamento do vento: padrões ou códigos de carga de vento e análises experimentais em túnel de vento. Nove modelos experimentais de morros isolados, sendo quatro simétricos bidimensionais, quatro simétricos tridimensionais e um assimétrico, todos com a consideração de dois tipos de terreno, categoria I – plano e categorias III-IV – medianamente rugosas são analisados. A partir de uma simulação experimental da camada limite em túnel de vento, é possível parametrizar os efeitos do vento sobre terrenos complexos (MILLER, 1995). A definição de variáveis, tais como: o perfil de velocidade do vento, a intensidade de turbulência, os efeitos topográficos na velocidade do vento constituem-se em elementos fundamentais para cálculos estruturais de edificações situadas no entorno. Foram desenvolvidas simulações experimentais no túnel de vento de camada limite Prof. Joaquim Blessmann, da Universidade Federal do Rio Grande do Sul. Na superfície do modelo assimétrico, na radial principal do mesmo, foram localizados nove perfis de medição, contendo cada um, vinte alturas de investigação. Os demais modelos tiveram como foco de análise o cume dos mesmos, também com perfis definidos nas mesmas alturas. As medições da velocidade do vento e da intensidade da turbulência foram procedidas por um sistema de anemometria de fio quente. Os dados obtidos em túnel de vento foram confrontados com expressões empíricas calculadas para os mesmos pontos segundo cinco códigos ou padrões de carga de vento, pontuando também as correlações entre os mesmos e entre estes com a norma brasileira NBR 6123 (1988). Os modelos, códigos ou padrões analisados foram: Jackson e Hunt (1975) e Davenport, Surry e Lemelin (1988), models e as normas a seguir nominadas:Norma Brasileira: cargas de vento em Edificações - NBR 6123 (1988); European Standard:Eurocode1: Basis of Design and Actions on Structures, CEN TC 250: 2002; Australian/New Zealand Standard: Minimum Design Loads on Structures, AS/NZS 1170.2: 2002; Architectural Institute of Japan AIJ: 2004; American Society of Civil Engineering Standard, ASCE 7-95 (ASCE 7-95), Minimum Design Loads for Buildings and Other Structures;National Building Code of Canada, 2005, (NRCC 2005). Além do estudo comparativo acerca da resposta estimada pelos códigos nominados, esta pesquisa constitui-se em um banco de dados de medições em pontos localizados nos morros nominados, em túnel de vento. O comparativo entre os padrões mostra a inexistência de uma harmonização entre os mesmos para a consideração dos parâmetros a serem empregues para o cálculo de cargas de vento. Estas diferenças na definição dos parâmetros básicos para o carregamento de vento em estruturas determinam grandes dificuldades na unificação de formatos recomendados na previsão das cargas de vento. Comparativamente aos dados experimentais, os padrões, em geral, mostram-se conservadores para os dois tipos de morros analisados, simétrico e assimétrico, para os dois tipos de terreno, tanto considerando-se análises em 2D ou 3D. / The structure of the wind varies with the characteristics of the terrain and roughness land surface, slowing its intensity proportionally according to the surrounding terrain, which determines the onset of boundary layer (ABL). The characteristics of the flow over and around topographic features such as hills are of great interest in many applications, especially those related to wind engineering. From an experimental simulation of the boundary layer wind tunnel, it is possible to parameterize the effects of wind over complex terrain (MILLER, 1995). The definition of variables, such as the profile of wind speed, the turbulence intensity, the topographic effects on wind speed are key elements in structural calculations for buildings situated around the area. This research focuses on the investigation of the ABL complex terrain conditions, analyzing the structure of turbulent flow and characterization of separation and reattachment of the flow. Experimental simulations were developed in the wind tunnel of the atmospheric boundary layer Prof. Joaquim Blessmann, Federal University of Rio Grande do Sul in nine models of hills, four symmetrical two-dimensional, four- symmetrical threedimensional and one asymmetrical, all considering two types of terrain, category I - plan and Category III-IV- moderately rough. The surface of the asymmetric model was measured in nine profiles ploted on the main radial of the hill, with twenty heights each, and the other models were examined at top of the hill. The measurements were performed with a system of hot wire anemometry to measure the wind velocity and intensity of turbulence. The identification of the data obtained in the wind tunnel were confronted with empirical expressions for the same points, in order to establish the correlations between patterns and among these with NBR 6123 (1988). Two models and five codes of wind loads are analyzed: Jackson and Hunt (1975) and e Davenport, Surry e Lemelin (1988) Models and Brazilian Association of Technical Standards: Wind Load on Buildings, NBR 6123 (1988); European Standard: Eurocode1: Basis of Design and Actions on Structures, CEN TC 250: 2002; Australian/New Zealand Standard: Minimum Design Loads on Structures, AS/NZS 1170.2: 2002; Architectural Institute of Japan AIJ: 2004; American Society of Civil Engineering Standard, ASCE 7-95 (ASCE 7-95), Minimum Design Loads for Buildings and Other Structures; National Building Code of Canada, 2005, (NRCC 2005) codes. This study focuses not only the comparison of the response estimated by international codes nominees, but also a data bank of wind tunnel data to validate this tool based on empirical expressions. The comparison of the patterns shows a lack of consideration for international harmonization of the parameters to be employed for the calculations of wind loads. These differences in defining the basic parameters for the wind loading on structures determines difficulties to unify the formats recommended in the prediction of wind loads. Compared to the experimental data, the patterns will generally show up conservative for both types of mounts analyzed, symmetrical and asymmetrical, for both types of terrain, both considering 2D or 3D.

O formato do perfil de velocidades do vento varia de acordo com as características locais da superfície terrestre e de rugosidade do terreno, parâmetros que definem o perfil da Camada-Limite Atmosférica (CLA). As características do escoamento do ar atmosférico sobre e ao redor de acidentes geográficos, tais como morros e colinas, são de grande interesse para aplicações relacionadas à Engenharia de Turbinas e Parques Eólicos. No topo de morros, ocorre a aceleração do vento, fenômeno que pode representar um fator decisivo para a instalação de aerogeradores. Este trabalho dedica-se ao estudo do comportamento da CLA como função da inclinação e rugosidade superficial da elevação, fazendo uso da Dinâmica de Fluidos Computacional (CFD) para construir perfis de velocidade do vento e de intensidade de turbulência. O problema de fechamento das Equações Médias de Reynolds (RANS) é contornado com o uso do modelo de turbulência k-ω SST; os resultados numéricos obtidos são comparados com dados experimentais medidos em túnel de vento sobre modelos em escala dos morros. São testados oito modelos de morros com declives que variam de 25° a 64° para dois tipos de categorias de terreno, em 2D e 3D, e são aplicados dois códigos analíticos para representar o perfil de velocidades de entrada. Resultados numéricos para os perfis de velocidade apresentam diferença inferior a 4% em relação aos respectivos dados obtidos experimentalmente. Os perfis de intensidade de turbulência apresentam diferença máxima na casa dos 7% em comparação aos dados experimentais, o que é explicado pelo fato de que não é possível inserir o perfil de entrada de intensidade de turbulência nas simulações numéricas. Em alternativa, foi usado um valor constante resultado da média dos valores dos perfis usados no túnel de vento. Os modelos de morro em 3D apresentam maior concordância nos resultados de velocidade que os modelos em 2D e que ademais quanto maior é a inclinação do morro maior é a concordância com as medições experimentais. / The shape of the wind velocity profile changes according to local features of terrain shape and roughness, which are parameters responsible for defining the Atmospheric Boundary Layer (ABL) profile. Air flow characteristics over and around landforms, such as hills, are of considerable importance for applications related to Wind Farm and Turbine Engineering. The air flow is accelerated on top of hills, which can represent a decisive factor for Wind Turbine placement choices. The present work focuses on the study of ABL behavior as a function of slope and surface roughness of hill-shaped landforms, using the Computational Fluid Dynamics (CFD) to build wind velocity and turbulent intensity profiles. Reynolds-Averaged Navier-Stokes (RANS) equations are closed using the SST k-ω turbulence model; numerical results are compared to experimental data measured in wind tunnel over scale models of the hills under consideration. Eight hill models with slopes varying from 25° to 64° were tested for two types of terrain categories in 2D and 3D, and two analytical codes are used to represent the inlet velocity profiles. Numerical results for the velocity profiles show differences under 4% when compared to their respective experimental data. Turbulent intensity profiles show maximum differences around 7% when compared to experimental data, this can be explained by not being possible to insert inlet turbulent intensity profiles in the simulations. Alternatively, constant values based on the averages of the turbulent intensity at the wind tunnel inlet were used. The 3D models present greater concordance in the speed results than the 2D models and that in addition the greater the slope of the hill, the greater the agreement with the experimental measurements.

abstract: This thesis focuses on an improved understanding of the dynamics at different length scales of wind farms in an atmospheric boundary layer (ABL) using a series of visualization studies and Fourier, wavelet based spectral analysis using high fidelity large eddy simulation (LES). For this purpose, a robust LES based neutral ABL model at very high Reynolds number has been developed using a high order spectral element method which has been validated against the previous literature. This ABL methodology has been used as a building block to drive large wind turbine arrays or wind farms residing inside the boundary layer as documented in the subsequent work. Studies conducted in the thesis involving massive periodic wind farms with neutral ABL have indicated towards the presence of large scale coherent structures that contribute to the power generated by the wind turbines via downdraft mechanisms which are also responsible for the modulation of near wall dynamics. This key idea about the modulation of large scales have seen a lot of promise in the application of flow past vertically staggered wind farms with turbines at different scales. Eventually, studies involving wind farms have been progressively evolved in a framework of inflow-outflow where the turbulent inflow is being fed from the precursor ABL using a spectral interpolation technique. This methodology has been used to enhance the understanding related to the multiscale physics of wind farm ABL interaction, where phenomenon like the growth of the inner layer, and wake impingement effects in the subsequent rows of wind turbines are important owing to the streamwise heterogeneity of the flow. Finally, the presence of realistic geophysical effects in the turbulent inflow have been investigated that influence the flow past the wind turbine arrays. Some of the geophysical effects that have been considered include the presence of the Coriolis forces as well as the temporal variation of mean wind magnitude and direction that might occur due to mesoscale dynamics. This study has been compared against field experimental results which provides an important step towards understanding the capability of the mean data driven LES methodology in predicting realistic flow structures. / Dissertation/Thesis / Turbulent flow past a large periodic wind farm / Doctoral Dissertation Mechanical Engineering 2018