Comparison of turbulence models for the computational fluid dynamics simulation of wind turbine wakes in the atmospheric boundary layer

An elliptic computational fluid dynamics wake model based on the actuator disk concept is used to simulate a wind turbine, approximated by a disk upon which a distribution of forces, defined as axial momentum sources, is applied on an incoming non-uniform shear flow. The rotor is supposed to be uniformly loaded with the exerted forces estimated as a function of the incident wind speed, thrust coefficient and rotor diameter. The model is assessed in terms of wind speed deficit and added turbulence intensity for different turbulence models and is validated from experimental measurements of the Sexbierum wind turbine experiment.

Atmospheric boundary layer development over a narrow coastal plain during onshore flow

A study of the structure of the daytime atmospheric boundary layer during onshore flow over a narrow coastal plain is presented. The main emphasis of the study is on the nature and causes of heating and cooling observed in the boundary layer temperature profiles. Measurements included vertical temperature profiles above at least two sites derived from radiosondes and aircraft, as well as surface estimates of radiative and sensible heat fluxes. Surface meteorological and pilot balloon data were also available, providing further evidence of short-term changes in atmospheric boundary layer structure. The Manawatu case was representative of autumnal anticyclonic conditions with weak pressure gradients, and illustrated typical diurnal development of a convective boundary layer over a coastal plain bordered by mountain ranges, with a transition from a stable nocturnal situation to a well-mixed profile in the afternoon. The profiles show surface input of heat propagating upwards through the boundary layer during the day, as well as entrainment of heat at the top associated with shear induced turbulence and/or penetrative convection. Applying a one-dimensional model, estimates of boundary layer heat budget components were obtained for four time periods during the day. Later periods were affected by cumulus cloud development at the top of the boundary layer, resulting in significant changes in individual components. Input of sensible heat from the surface decreased, while the addition of heat to the boundary layer from both cloud condensation and advection increased.

Characteristics of ABL during active and weak phases of monsoon

Understanding of the Atmospheric Boundary Layer (ABL) is imperative in the arena of the monsoon field. Here, the features of the ABL are studied employing Conserved Variable Analysis (CVA) using equivalent potential temperature and humidity. In addition, virtual potential temperature and wind are used during active and weak phases of monsoon. The analysis is carried out utilising the radiosonde observations during the monsoon months for two stations situated in the west coast of India. All these parameters show considerable variations during active and weak monsoon phases in both the stations. The core speed and core height vary with these epochs. The core speed is found to be more than 38 knots in the active monsoon phase around 1.2 km over Trivandrum and around 2 km over Mangalore. But during weak monsoon phase the core wind speed is decreased and core height is elevated over both stations. The wind direction shows an additional along shore component during weak monsoon period. The Convective Boundary Layer (CBL) height shows increase during weak monsoon phase over both stations due to less cloudiness and subsequent insolation. The CBL height during the southwest monsoon is more over Mangalore and is attributed by the orographic lifting in the windward side of the Western Ghats while the influence of the Ghats is less over Trivandrum.

Characteristics of oscillation in surface layer and above over Cochin

The oscillations in the Atmospheric Boundary Layer (ABL) are important because the transport mechanism from the surface to the upper atmosphere is governed by the ABL characteristics. The study was carried out using wind and temperature data observed at surface, 925 hPa and 850 hPa levels over Cochin and the different frequencies embedded in the boundary layer parameters are identified by employing wavelet technique. Surface boundary layer characteristics over the monsoon region are closely linked to the upper layer monsoon features. In this perception it is important to study the various oscillations in the surface boundary layer and the layer above. It is found that the wind and temperature at different levels show oscillations in Quasi Biweekly Mode (QBM) and Intra Seasonal Oscillation (ISO) bands as observed in a typical monsoon system. Amplitude of the oscillation varies with height. The amplitude of the QBM periodicity is more in the surface levels but in the upper levels the amplitude of the ISO periodicity is more than that of the QBM. From this, it is obvious that the controlling mechanism of QBM band is surface parameters such as surface friction and that for ISO band is associated with the active-break cycles of monsoon system

Variations of surface boundary layer parameters associated with Cyclone Gonu over the Arabian Sea using QuikSCAT data

This study attempted to quantify the variations of the surface marine atmospheric boundary layer (MABL) parameters associated with the tropical Cyclone Gonu formed over the Arabian Sea during 30 May–7 June 2007 (just after the monsoon onset). These characteristics were evaluated in terms of surface wind, drag coefficient, wind stress, horizontal divergence, and frictional velocity using 0.5◦ × 0.5◦ resolution Quick Scatterometer (QuikSCAT) wind products. The variation of these different surface boundary layer parameters was studied for three defined cyclone life stages: prior to the formation, during, and after the cyclone passage. Drastic variations of the MABL parameters during the passage of the cyclone were observed. The wind strength increased from 12 to 22 m s−1 in association with different stages of Gonu. Frictional velocity increased from a value of 0.1–0.6 m s−1 during the formative stage of the system to a high value of 0.3–1.4 m s−1 during the mature stage. Drag coefficient varied from 1.5 × 10−3 to 2.5 × 10−3 during the occurrence of Gonu. Wind stress values varied from 0.4 to 1.1 N m−2. Wind stress curl values varied from 10 × 10−7 to 45 × 10−7 N m−3. Generally, convergent winds prevailed with the numerical value of divergence varying from 0 to –4 × 10−5 s−1. Maximum variations of the wind parameters were found in the wall cloud region of the cyclone. The parameters returned to normally observed values in 1–3 days after the cyclone passage

Characteristic Study of the Boundary Layer Parameters over the Arabian Sea and the Bay of Bengal Using the QuikSCAT Dataset

The marine atmospheric boundary layer (MABL) plays a vital role in the transport of momentum and heat from the surface of the ocean into the atmosphere. A detailed study on the MABL characteristics was carried out using high-resolution surface-wind data as measured by the QuikSCAT (Quick scatterometer) satellite. Spatial variations in the surface wind, frictional velocity, roughness parameter and drag coe±cient for the di®erent seasons were studied. The surface wind was strong during the southwest monsoon season due to the modulation induced by the Low Level Jetstream. The drag coe±cient was larger during this season, due to the strong winds and was lower during the winter months. The spatial variations in the frictional velocity over the seas was small during the post-monsoon season (»0.2 m s¡1). The maximum spatial variation in the frictional velocity was found over the south Arabian Sea (0.3 to 0.5 m s¡1) during the southwest monsoon period, followed by the pre-monsoon over the Bay of Bengal (0.1 to 0.25 m s¡1). The mean wind-stress curl during the winter was positive over the equatorial region, with a maximum value of 1.5£10¡7 N m¡3, but on either side of the equatorial belt, a negative wind-stress curl dominated. The area average of the frictional velocity and drag coe±cient over the Arabian Sea and Bay of Bengal were also studied. The values of frictional velocity shows a variability that is similar to the intraseasonal oscillation (ISO) and this was con¯rmed via wavelet analysis. In the case of the drag coe±cient, the prominent oscillations were ISO and quasi-biweekly mode (QBM). The interrelationship between the drag coe±cient and the frictional velocity with wind speed in both the Arabian Sea and the Bay of Bengal was also studied.

A case study of boundary layer ventilation by convection and coastal processes

It is often assumed that ventilation of the atmospheric boundary layer is weak in the absence of fronts, but is this always true? In this paper we investigate the processes responsible for ventilation of the atmospheric boundary layer during a nonfrontal day that occurred on 9 May 2005 using the UK Met Office Unified Model. Pollution sources are represented by the constant emission of a passive tracer everywhere over land. The ventilation processes observed include shallow convection, turbulent mixing followed by large-scale ascent, a sea breeze circulation and coastal outflow. Vertical distributions of tracer are validated qualitatively with AMPEP (Aircraft Measurement of chemical Processing Export fluxes of Pollutants over the UK) CO aircraft measurements and are shown to agree impressively well. Budget calculations of tracers are performed in order to determine the relative importance of these ventilation processes. Coastal outflow and the sea breeze circulation were found to ventilate 26% of the boundary layer tracer by sunset of which 2% was above 2 km. A combination of coastal outflow, the sea breeze circulation, turbulent mixing and large-scale ascent ventilated 46% of the boundary layer tracer, of which 10% was above 2 km. Finally, coastal outflow, the sea breeze circulation, turbulent mixing, large-scale ascent and shallow convection together ventilated 52% of the tracer into the free troposphere, of which 26% was above 2 km. Hence this study shows that significant ventilation of the boundary layer can occur in the absence of fronts (and thus during high-pressure events). Turbulent mixing and convection processes can double the amount of pollution ventilated from the boundary layer.

Synoptic controls on boundary-layer characteristics

We report the characteristics of the three-dimensional, time evolving, atmospheric boundary layer that develops beneath an idealised, dry, baroclinic weather system. The boundary-layer structure is forced by thermal advection associated with the weather system. Large positive heat fluxes behind the cold front drive a vigorous convective boundary layer, whereas moderate negative heat fluxes in the warm sector between the cold and warm fronts generate shallow, stably stratified or neutral boundary layers. The forcing of the boundary-layer structure is quantified by forming an Eulerian mass budget integrated over the depth of the boundary layer. The mass budget indicates that tropospheric air is entrained into the boundary layer both in the vicinity of the high-pressure centre, and behind the cold front. It is then transported horizontally within the boundary layer and converges towards the cyclone's warm sector, whence it is ventilated out into the troposphere. This cycling of air is likely to be important for the ventilation of pollution out of the boundary layer, and for the transformation of the properties of large-scale air masses.

Boundary-layer friction in midlatitude cyclone

Results from an idealized three-dimensional baroclinic life-cycle model are interpreted in a potential vorticity (PV) framework to identify the physical mechanisms by which frictional processes acting in the atmospheric boundary layer modify and reduce the baroclinic development of a midlatitude storm. Considering a life cycle where the only non-conservative process acting is boundary-layer friction, the rate of change of depth-averaged PV within the boundary layer is governed by frictional generation of PV and the flux of PV into the free troposphere. Frictional generation of PV has two contributions: Ekman generation, which is directly analogous to the well-known Ekman-pumping mechanism for barotropic vortices, and baroclinic generation, which depends on the turning of the wind in the boundary layer and low-level horizontal temperature gradients. It is usually assumed, at least implicitly, that an Ekman process of negative PV generation is the mechanism whereby friction reduces the strength and growth rates of baroclinic systems. Although there is evidence for this mechanism, it is shown that baroclinic generation of PV dominates, producing positive PV anomalies downstream of the low centre, close to developing warm and cold fronts. These PV anomalies are advected by the large-scale warm conveyor belt flow upwards and polewards, fluxed into the troposphere near the warm front, and then advected westwards relative to the system. The result is a thin band of positive PV in the lower troposphere above the surface low centre. This PV is shown to be associated with a positive static stability anomaly, which Rossby edge wave theory suggests reduces the strength of the coupling between the upper- and lower-level PV anomalies, thereby reducing the rate of baroclinic development. This mechanism, which is a result of the baroclinic dynamics in the frontal regions, is in marked contrast with simple barotropic spin-down ideas. Finally we note the implications of these frictionally generated PV anomalies for cyclone forecasting.

Advances in understanding mineral dust and boundary layer processes over the Sahara from Fennec aircraft observations

The Fennec climate program aims to improve understanding of the Saharan climate system through a synergy of observations and modelling. We present a description of the Fennec airborne observations during 2011 and 2012 over the remote Sahara (Mauritania and Mali) and the advances in the understanding of mineral dust and boundary layer processes they have provided. Aircraft instrumentation aboard the UK FAAM BAe146 and French SAFIRE Falcon 20 is described, with specific focus on instrumentation specially developed and relevant to Saharan meteorology and dust. Flight locations, aims and associated meteorology are described. Examples and applications of aircraft measurements from the Fennec flights are presented, highlighting new scientific results delivered using a synergy of different instruments and aircraft. These include: (1) the first airborne measurement of dust particles sized up to 300 microns and associated dust fluxes in the Saharan atmospheric boundary layer (SABL), (2) dust uplift from the breakdown of the nocturnal low-level jet before becoming visible in SEVIRI satellite imagery, (3) vertical profiles of the unique vertical structure of turbulent fluxes in the SABL, (4) in-situ observations of processes in SABL clouds showing dust acting as CCN and IN at −15 °C, (5) dual-aircraft observations of the SABL dynamics, thermodynamics and composition in the Saharan heat low region (SHL), (6) airborne observations of a dust storm associated with a cold-pool (haboob) issued from deep convection over the Atlas, (7) the first airborne chemical composition measurements of dust in the SHL region with differing composition, sources (determined using Lagrangian backward trajectory calculations) and absorption properties between 2011 and 2012, (8) coincident ozone and dust surface area measurements suggest coarser particles provide a route for ozone depletion, (9) discrepancies between airborne coarse mode size distributions and AERONET sunphotometer retrievals under light dust loadings. These results provide insights into boundary layer and dust processes in the SHL region – a region of substantial global climatic importance.

Multiyear measurements of the oceanic and atmospheric boundary layers at the Brazil-Malvinas confluence region

This study analyzes and discusses data taken from oceanic and atmospheric measurements performed simultaneously at the Brazil-Malvinas Confluence (BMC) region in the southwestern Atlantic Ocean. This area is one of the most dynamical frontal regions of the world ocean. Data were collected during four research cruises in the region once a year in consecutive years between 2004 and 2007. Very few studies have addressed the importance of studying the air-sea coupling at the BMC region. Lateral temperature gradients at the study region were as high as 0.3 degrees C km(-1) at the surface and subsurface. In the oceanic boundary layer, the vertical temperature gradient reached 0.08 degrees C m(-1) at 500 m depth. Our results show that the marine atmospheric boundary layer (MABL) at the BMC region is modulated by the strong sea surface temperature (SST) gradients present at the sea surface. The mean MABL structure is thicker over the warmside of the BMC where Brazil Current (BC) waters predominate. The opposite occurs over the coldside of the confluence where waters from the Malvinas (Falkland) Current (MC) are found. The warmside of the confluence presented systematically higher MABL top height compared to the coldside. This type of modulation at the synoptic scale is consistent to what happens in other frontal regions of the world ocean, where the MABL adjusts itself to modifications along the SST gradients. Over warm waters at the BMC region, the MABL static instability and turbulence were increased while winds at the lower portion of the MABL were strong. Over the coldside of the BC/MC front an opposite behavior is found: the MABL is thinner and more stable. Our results suggest that the sea-level pressure (SLP) was also modulated locally, together with static stability vertical mixing mechanism, by the surface condition during all cruises. SST gradients at the BMC region modulate the synoptic atmospheric pressure gradient. Postfrontal and prefrontal conditions produce opposite thermal advections in the MABL that lead to different pressure intensification patterns across the confluence.

Detection and Modeling of Boundary Layer Height

This thesis tackles the problem of the automated detection of the atmospheric boundary layer (BL) height, h, from aerosol lidar/ceilometer observations. A new method, the Bayesian Selective Method (BSM), is presented. It implements a Bayesian statistical inference procedure which combines in an statistically optimal way different sources of information. Firstly atmospheric stratification boundaries are located from discontinuities in the ceilometer back-scattered signal. The BSM then identifies the discontinuity edge that has the highest probability to effectively mark the BL height. Information from the contemporaneus physical boundary layer model simulations and a climatological dataset of BL height evolution are combined in the assimilation framework to assist this choice. The BSM algorithm has been tested for four months of continuous ceilometer measurements collected during the BASE:ALFA project and is shown to realistically diagnose the BL depth evolution in many different weather conditions. Then the BASE:ALFA dataset is used to investigate the boundary layer structure in stable conditions. Functions from the Obukhov similarity theory are used as regression curves to fit observed velocity and temperature profiles in the lower half of the stable boundary layer. Surface fluxes of heat and momentum are best-fitting parameters in this exercise and are compared with what measured by a sonic anemometer. The comparison shows remarkable discrepancies, more evident in cases for which the bulk Richardson number turns out to be quite large. This analysis supports earlier results, that surface turbulent fluxes are not the appropriate scaling parameters for profiles of mean quantities in very stable conditions. One of the practical consequences is that boundary layer height diagnostic formulations which mainly rely on surface fluxes are in disagreement to what obtained by inspecting co-located radiosounding profiles.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

A oscilação Madden - Julian na Amazônia Oriental: variáveis superfíciais

Tendo como foco as múltiplas escalas de tempo que atuam na Amazônia, este trabalho foi desenvolvido com o objetivo de investigar a possível influencia da Oscilação Madden – Julian (OMJ) em elementos turbulentos da CLP. A OMJ foi identificada a partir de 30 anos de dados de reanálise de radiação de onda longa (ROL) e componente zonal do vento (u). As grandezas turbulentas foram estudadas a partir da variância, covariância e coeficiente de correlação de um conjunto de dados de resposta rápida coletado na torre micrometeorológica de Caxiuanã (PA), e tratados com a Transformada em Ondeletas (TO) para se obter a contribuição de cada escala para estes momentos estatísticos. A análise dos 30 anos de dados de ROL e u mostrou que a ocorrência da OMJ está ligada com o fenômeno do El Niño/Oscilação Sul (ENOS), bem como influência do ENOS no tempo da região amazônica pode estar associado a presença ou não da OMJ. Foi observado que anos de El Niño tendem a desfavorecer a ocorrência da OMJ e anos de La Niña tendem a favorecer o desenvolvimento da oscilação. Caso uma OMJ se desenvolva durante um episodio de El Niño, a oscilação pode influenciar a temperatura, a velocidade do vento e a precipitação de forma diferente ao do El Niño. A análise por fase da OMJ mostrou que, em Belém, há diferença significativa na temperatura máxima e na precipitação entre cada fase, porém, a temperatura mínima e o módulo do vento apresentaram pouca diferença. Os fluxos cinemáticos turbulentos analisados, por escala, em três horários distintos, foram mais diferentes durante o período diurno, principalmente w’T’ e w’q’. A diferença entre fase ativa e fase inativa foi reduzindo com passar do dia, durante o período de transição dia – noite, poucas escalas tiveram diferença significativa, e durante a noite, nenhuma escala teve nível de confiança acima ou igual a 95%. Estes resultados indicam que a convecção diurna é o mecanismo responsável por esta diferença e como a OMJ atua como uma grande célula convectiva, a convecção local é amplificada, explicando a grande diferença observada entre as fases durante o período diurno.