Validating CFD Predictions of Flow over an Escarpment Using Ground-Based and Airborne Measurement Devices

: Micrometeorological observations from a tower, an eddy-covariance (EC) station and an unmanned aircraft system (UAS) at the WINSENT test-site are used to validate a computational ﬂuid dynamics (CFD) model, driven by a mesoscale model. The observation site is characterised by a forested escarpment in a complex terrain. A two-day measurement campaign with a ﬂow almost perpendicular to the escarpment is analysed. The ﬁrst day is dominated by high wind speeds, while, on the second one, calm wind conditions are present. Despite some minor differences, the ﬂow structure, analysed in terms of horizontal wind speeds, wind direction and inclination angles shows similarities for both days. A real-time strategy is used for the CFD validation with the UAS measurement, where the model follows spatially and temporally the aircraft. This strategy has proved to be successful. Stability indices such as the potential temperature and the bulk Richardson number are calculated to diagnose atmospheric boundary layer (ABL) characteristics up to the highest ﬂight level. The calculated bulk Richardson values indicate a dynamically unstable region behind the escarpment and near the ground for both days. At higher altitudes, the ABL is returning to a near neutral state. The same characteristics are found in the model but only for the ﬁrst day. The second day, where shear instabilities are more dominant, is not well simulated. UAS proves its great value for sensing the ﬂow over complex terrains at high altitudes and we demonstrate the usefulness of UAS for validating and improving models.


Introduction
To achieve the objective of making the European Union climate-neutral by 2050, it will be necessary to maximize the deployment of renewable energy in the years to come [1]. The clean energy transition involves encouraging high levels of renewable energy penetrations. Harvesting power from wind power is one of the fastest-growing renewable energy methods. To meet the climate-neutral target, more wind energy capacity will be installed in the coming years. Independent of the exact number of new installations, most of them will be onshore [2]. A wind energy project includes site investigation, resource assessment, environmental impact assessment, technical and administrative procedures with local communities and, finally, construction work. The site investigation and resource assessment phases need accurate and reliable models for predicting the wind flow. During the last two decades, linear flow models, mainly based on the Jackson-Hunt theory [3], were extensively used in the wind industry for wind resource assessment and micro-siting. These models are, however, inaccurate for wind speed prediction in very complex terrain, which led to the use of computational fluid dynamics (CFD), capable of modelling the wind flow in details around specific orographic features. The spatial and temporal variability of the atmospheric boundary layer (ABL) makes the prediction of the wind flow by models one of the main challenges, especially in complex terrain. The commercial CFD software for wind resource assessment uses logarithmic wind profile and tunes parameters from the Monin-Obukhov similarity theory [4] to calculate annual energy prediction for site assessment. The logarithmic wind profiles are not site specific and don't have micro meteorological information of wind at the particular site. This is a major error source for the calculation of annual energy production. Hence, the variability of the atmosphere should be considered at the mesoscale. Only weather forecast models, driven by global models, can capture the dynamic processes of the atmosphere [5]. Today's mesoscale models are used increasingly for wind-energy applications with efforts in improving surface winds in complex terrain [6,7]. However, it has been observed that numerical noise is appearing with increasing terrain complexity [8]. A remedy is to use the data from mesoscale models of a given region as a realistic boundary condition for a more detailed microscale modelling. That way, the wind characteristics can be determined at a high resolution in the area of interest. This approach is referred in the literature as mesoscale-microscale modelling and have been successfully used in different domains. For instance, the improved predictability of pollutants such as NO x and NO 2 in urban areas by coupling a mesoscale model with a CFD model has been demonstrated in recent studies [9,10]. For wind assessment studies, the performance of coupled simulations is more and more recognized and has been extended. Even in the industry, it showed an improvement of wind power forecasting [11] .
All methods and models should be validated, ideally with in-situ and/or remote sensing measurements. The need for field campaigns to provide data for model validation in complex terrain has led to the funding of several projects over the last few years. Examples are the Kassel Experiment in Germany [12], the Alaiz Experiment in Spain [13], the double-ridge Perdigao in Portugal [14] or the Wind Forecast Improvement Project (WFIP2) in the U.S [15], all aiming at improving the understanding of physical processes and wind forecasts in complex terrain. Traditional in-situ platforms for atmospheric measurements are towers installed at a certain location. The vertical information is provided by instruments installed at several discrete levels. They can provide continuous recording of the investigated variables over long time periods. To follow the growth of wind turbine sizes, the meteorological tower should ideally become taller and have the same height as the wind turbines, making their installation technically complicated and more expensive. These limitations have led to the use of lidar (light detection and ranging) technologies for wind measurements. However, in complex environments, the spatial heterogeneity and transient features introduced by the terrain can lead to difficulties in interpreting the line of sight wind speed and reconstructing the wind field [16]. For instance, the central goal of the International Energy Agency Wind Task 32 community is developing frameworks for the deployment and use of wind lidars in complex flow. One promising approach is the use of two or more lidars to reconstruct the flow field in a complex terrain. Limitations of this approach are mainly the high costs and the uncertainties [17]. Other emerging techniques for wind-energy applications are flying platforms. Small unmanned aircraft systems (UAS) have become frequently used platforms for ABL research. They have the advantage of being versatile, even in highly complex terrain and the necessary infrastructure is small compared to masts and active remote sensing techniques. They can fly in the close vicinity of a wind turbine, even at very low altitudes, where the main vertical and horizontal inhomogeneities are occurring. The capability of UAS for turbulence characterization and providing high spatial-temporal sampling frequencies has been demonstrated in research works over the last several years. The use of these systems contributed to a better characterisation and understanding in aerosol formations [18,19], storm developments [20], atmospheric research [21][22][23], or for wind energy purposes [24,25].
In the research project WINSENT (Wind Science and Engineering Test Site in Complex Terrain), a field test-site is being built on the Swabian Alb in Southern Germany [26]. It will be the first wind energy test-site of this size in complex terrain. The first phase of the WINSENT project investigates the local wind flow at the test-site using different measurement equipments and simulation models. During the second stage, two wind turbines with a nominal output of around 750 kW will be installed next to the escarpment. The effect of the topography on the wind-turbines performances, wakes, and material fatigue will be studied both experimentally and numerically [27,28]. The on-site instruments consist of permanently installed towers equipped with ultrasonic anemometers and eddy-covariance (EC) micro-meteorological stations. Additionally, the multipurpose airborne sensor carrier MASC [29] and several lidars [30] carry out measurements over and downstream of the escarpment during different measurement campaigns. All of the equipment is going to provide new datasets to evaluate and validate the numerical models. In the framework of the WINSENT project, a numerical study is conducted with models, with a different range of scales such as the Weather Research and Forecasting (WRF) model [31], the Open Source Field Operation and Manipulation (OpenFOAM) software [32], and the solver FlOWer [33] .
In the present study, a numerical and experimental investigation of the flow at the WINSENT test-site is presented. Two flight measurement campaigns, which took place over two days in September 2018, are analysed. The first day is dominated by high wind speeds while the conditions on the second day are calmer. Measurements by the MASC, along with a meteorological tower and an EC station are used for the evaluation of our model. A real-time strategy is adopted for the comparison of the numerical results with the MASC measurements. The performance of our methodology is presented in this study as follows: in Section 2, the numerical tools are presented, while, in Section 3, the test-site along with the measurement devices are described. Section 4 discusses the numerical and experimental results from the measurement campaigns together with a stability analysis and finally a conclusion with outlooks is given in Section 5.

Numerical Model
In the current study, the wind flow is simulated using the software OpenFOAM v6. An unsteady Reynolds Averaged Navier-Stokes (URANS) approach, necessary for the investigation of non-stationary flow features, was considered [34,35]. The governing equations for mass, momentum, and potential temperature, under the Boussinesq approximation, are written relative to a hydrostatic state (index h) in a Cartesian coordinate system as: where u i (i = 1, 2, 3) are the velocity components in a Cartesian framework (x, y, z, where z is aligned with gravity), p and θ are the pressure perturbation from the hydrostatic reference state and potential temperature, µ and µ t are the dynamic molecular and turbulent eddy viscosities, g i is the gravitational acceleration term, and σ h is the turbulent Prandtl number. F Ci and F Di are the source term representing the Coriolis and the drag force exerted by the forest on the wind flow, respectively. The hydrostatic fluid density ρ h is defined in a hydrostatic reference state (subscript 0) as a function of the hydrostatic pressure p h and the temperature T h as: where p 0 is the constant reference pressure set to 1000 hPa, T 0 is the reference temperature equal to 288.5 K, A = 50 K, and R d = 287.05 J kg −1 K −1 [36,37]. The Sogachev k − model is used for the turbulence closure [38]. Additionally, source terms in the turbulent kinetic energy and turbulence dissipation rate equations, using the constants proposed by [39] were added in order to take into account for the drag effect of vegetation on the flow.

Computational Set Up and Mesoscale Forcing
The simulations are conducted on a domain of 10 × 10 × 2.5 km 3 , centered in the escarpment. The spatial resolution consists of 10 m increments in the horizontal directions and the vertical discretization has values between 1.8 m at the ground and 80 m at the top of the domain. The terrain data are based on the digital height model from the Baden-Württemberg Authorities for Spatial Information and Rural Development (LGL) with a spatial resolution of 5 m. In addition, the digitized landscape model from LGL was utilized to describe the land cover on the Earth's surface. Three classes of land uses were retained for this study: urban, ground, and forest ( Figure 1b). In this paper, a deciduous forest with a height of 20 m is considered for the forest patch. The urban and the ground classes were assigned values of 0.50 m and 0.02 m for the aerodynamic roughness length.
The mesoscale flow is simulated using WRF with five nested domains. WRF uses the European Centre for Medium-Range Weather Forecasts analysis data as initial and boundary conditions. The innermost nest has a horizontal grid size of 150 m, a vertical grid size of 15 m close to the ground and is run in a large eddy simulation mode. A one-way nesting method was used for the coupling of WRF and OpenFOAM: the innermost nest provides boundary conditions to the microscale model. Border data include the velocity components, pressure and temperature at 1-minute intervals.
The WRF model provides also the heat flux at the ground. No special temporal and spatial interpolation has been applied as the mesoscale model is considered fine enough near the ground.

Test Site
The field measurements were made at the WINSENT test-site, built on the Swabian Alb in Southern Germany. The main feature of the site is a sloping escarpment up to 200 m in height with respect to the valley (see Figure 1a). The escarpment is characterized by a very steep slope of around 30 • from the Northwest and a more gradual slope at around 15 • further west. A dense vegetation with a tree height between 15 and 25 m is covering the escarpment. To the east of the escarpment, the ground is essentially flat and under agricultural use. The windrose (Figure 1c) measured on a height of 86 m at the test-site shows winds blowing predominantly from the west, and therefore perpendicular to the escarpment, most of the year.

Measurements
Measurements took place on 21 and 22 September 2018. Cup anemometers sampling at 20 Hz are mounted on the tower (48.6652 • latitude, 9.8347 • longitude) along with hygrothermographs, wind vanes, and barometers at various heights between 3 and 100 m ( Table 1). At that time, the sonic anemometers were not yet installed and therefore only the cup anemometer results are presented in this article. To minimise flow distortion due to the tower for the main wind directions, most instruments are mounted on booms pointing in the south-north direction (Figure 2a). An EC station is installed 50 m east of the tower. The raw data were recorded at 20 Hz. However, the 10-min statistics (as prescribed by the international standard IEC 61400-12-1) are used in this study unless otherwise noted. Additionally, in-situ measurements are made using the multipurpose airborne sensor carrier (MASC [24,25,29]). The MASC is an instrumented unmanned aircraft system which can resolve turbulent wind and temperature fluctuations with a resolution of 30 Hz. The wind measurements are performed by measuring the flow speed and flow angles at the aircraft nose with a five-hole flow probe. In order to measure the three-dimensional flow field at the test-site, a vertical grid of racetracks with legs going back and forth over the test-site was performed between 10:55 a.m. and 12:30 p.m. UTC and between 12:55 p.m. and 2:17 p.m. UTC for the first and second day, respectively. This experimental configuration allows for the changes in the flow field to be mapped in the region of the escarpment. An example of straight portions of the flight path selected for the analysis is shown in Figure 2b. Each racetrack consisted of several horizontal straight flight legs, in and against the flow direction. At least three legs (straight and level flight sections of about 1 km length) for each altitude were performed. Table 1 provides a detailed overview of the instruments used.

Results and Discussion
In this section, the simulation results are presented and compared with the tower and flight measurements. The focus in this study is put on the period of the UAS measurement campaigns.

Comparison of Wind and Turbulences Quantities with the Tower and EC Measurements
The simulation results are compared with measurements from cup anemometers and wind vanes installed on the mast. Time series of the 10-min averaged horizontal velocity and wind direction are shown in Figures 3 and 4   It can be seen in Figure 3 and 4 that, during the period of the flight campaign (marked with a grey box), the wind direction is remaining nearly constant between 240 • and 275 • for the first day and 265 • and 285 • for the second one. Only the wind speed shows a drop from 11 m s −1 to 7 m s −1 during the last 30 min for the first date. Except for the wind speed at 11:20 a.m. and 12:10 p.m. UTC on 21 September, the plots indicate a wind speed well simulated during the flight campaigns with values for mean absolute error at 100 m a.g.l. of 2.00 m s −1 and 0.74 m s −1 for the first and second day, respectively.
Turbulence is evaluated in terms of horizontal turbulence intensity (TI) given by the standard deviation of the horizontal wind speed normalized with the mean wind speed. Figure 5 compares the 10-min average simulated horizontal turbulence intensity against the measured data at 10 and 86 m a.g.l. The higher horizontal turbulence intensities are observed at the lower heights, i.e., at 10 m, with values above 40%. For 22 September, the TI values are even higher and can be explained by the definition of TI, which becomes large for small wind speeds. Several experimental [40,41] and numerical [42] studies have revealed that the maximum turbulent energy is produced in the upper half of the forest and carried into the region behind the forest [43,44]. For a westerly flow as in the present work, the 10 m sensor is located about 60 m downstream of the forested escarpment. This explains the high horizontal turbulence intensity values in comparison to the higher sensor at 86 m. The model simulates reasonably well the horizontal TI at the upper level and fails at catching the high values in the lower heights. The vegetative sink introduced by the canopy model [39] in the turbulence equations may be responsible and leads to an underestimate of the turbulent fluctuations. It should be reminded that TI is calculated with the data from a cup anemometer. As mentioned earlier, the sonic anemometers were not yet installed at that time. The influence of the vertical turbulence intensity is therefore not considered in this paper. at 45 m a.g.l. except for some specific times. The wind speed increases progressively with increasing height, with a visible wind speed reduction in the lower part of the profiles, i.e., at 10 m. On the first day, the wind speed at the lower cup anemometer is between 36 and 56% of the one at 100 m. The wind speed reduction is higher on the second day with ratio values ranging between 28 and 35%. Overall, the simulated and observed wind speed agree well in profile shape, except for some specific times.  Figure 7 illustrates the steepness of the terrain with a slope of around 20 • . It also shows the position of the modelled forest in green and in blue the legs flown on both days. Except for the three upper altitudes, the patterns are almost similar between the two flight dates as shown in Figure 7a The scale for the velocity was deliberately changed for both days in order to give a better visualisation of the flow along the UAS paths. The measured velocities show a similar structure for both days: directly after the slope, the velocities are higher due to a local updraft generated by the orography. Behind the escarpment, at 20 m a.g.l., the velocity decreases rapidly to 3 m s −1 . Generally, the simulated wind speed matches the UAS measurements very well, with a slight overprediction at upper levels for the first day, as already observed in the time series in Figure 3.
Not only wind speed but also wind direction is a relevant parameter of wind power output. Several studies have shown that the occurrence of large directional change results in lower turbine performance [45,46]. The measured and simulated wind direction at the WINSENT test-site are presented in Figure 9. Inclination angles are important for each turbine site: wind turbines have to be certified by the standards such as IEC 61400-1, where inflow angles are usually required to be within ±8 • [47]. The future turbine locations should comply with this requirement to ensure that the turbines will withstand the wear and tear during their expected lifetime. The modelled and measured inflow angles compared to each other are presented in Figure 10. Sudden changes in the inclination angle are numerically and experimentally noticed over the escarpment and are due to a local phenomenon induced by local features of the terrain. Positive flow angles of up to 15 • , numerically, and 20 • , experimentally, are reached. The position of this maximum reveals a flow dominated by upward movements and is located over or straight after the escarpment. It is noticeable that, for the second day, a negative angle can experimentally and numerically be seen in the last thirds of the flight path, indicating a flow not yet recovered. The flow remains disrupted even at upper levels (190 and 200 m a.g.l.), where inclination angle becomes smaller but still differs from zero. This reveals that the region of influence from the vegetated escarpment is not limited to low levels but extends to a region higher than the height of the escarpment itself. Previous wind tunnel experiments of escarpments with varying slope angles showed that the region of influence from an escarpment can go up to a height of 2 times the escarpment height [48,49].  As a surprise, the inclination angle is negative (i.e., a vertical wind component from above) at the beginning of the escarpment, even in the model, revealing an incoming flow that is already disturbed. For a westerly wind, the air passes over a small hill, located around 1.5 km west to the test-site ( Figure 1a) and a recirculation can be expected in the valley, leading to negative inclination angles at the end of the valley and still at the beginning of the escarpment. The inclination angles then return to positive values at the upper edge of the escarpment as could be expected by an orographically influenced flow that follows the terrain. To check this hypothesis, the simulated vertical wind speeds in a cross-section along the legs for both days are shown in Figure 11. The simulated values are averaged over the complete flight campaign, i.e., 95 and 82 minutes. On 21 September, a large recirculation zone enclosed between the hill and the escarpment can be seen and corroborates our hypothesis.
The variability of the horizontal wind speed and inclination angle, indicated by the standard deviations, represent a measure of the turbulence intensity. Their values are larger at low altitudes in comparison to higher levels due to increased turbulence near the ground. Looking at the two lowest levels, it is even noticeable that the variability is greater right after the escarpment. This is in accordance with the observation from the meteorological tower in Figure 5. Despite the large differences in the wind condition (high and low wind speeds) for both days, the flow structure seems to be quite similar.

Stability Considerations
The stability can be characterized, e.g., using the potential temperature or the Richardson number. The assessment of atmospheric stability using ground stations and towers is difficult in complex terrain as the flow can be characterized by phenomena such as updraft, downdraft, or canopy flow. A UAS has the ability to sample the atmosphere at higher altitudes in comparison to ground-based instrumentation, and can be a tool for atmospheric stability measurement. In our study, the potential temperature and the bulk Richardson number are used to assess the stability behaviour at the test-site. The bulk Richardson number Ri B [50], which can be computed from the wind and temperature measurement at two heights, is defined as: where T a is the average air temperature of the layer, θ u and θ l are the potential temperature at the upper z u and lower level z l , respectively. U is the horizontal wind speed and g = 9.81 m s −2 is the gravitational acceleration. There is still debate about the critical Richardson Ri c values, but typically Ri c = 0.25 for small ∆z. For Ri B < 0.25, the flow is probably unstable and turbulence is likely to occur. The diurnal pattern of potential temperature over the entire two days by the tower is presented in Figure 12. For 21 September, the tower recorded very small vertical potential temperature gradients (also known as potential temperature lapse rate ∆θ/∆z) during the flight campaign, corresponding to near-neutral atmospheric conditions. During the flight campaign on 22 September, a negative lapse rate (calculated from 2 and 96 m) of −1.8 K/100 m is measured and indicates the occurrence of an unstable atmosphere. The sudden drop in the temperature at 4:00 p.m. UTC on 21 September due to a cold front has led to lower temperatures (around 10 K lower) on 22 September. In Figure 13, the bulk Richardson number at the mast is calculated using the mean wind speed and temperature data (10- The potential temperature profiles from the tower, the MASC and the simulation are compared in Figure 14. Since potential temperature values are different for both days, the differences, that is, the potential temperature minus its minimum value (θ − θ min ) are shown. The MASC measurements are separated in three regions: one over the escarpment, the second part is centred around the tower and the last region is further downstream. The simulation results are averaged over the whole flight duration. The potential temperature at all positions decreases with height. For the first day (Figure 14a), the region of decreasing temperature is visible in the lower 50 m, while, on the second day (Figure 14b), the layer up to 100 m is concerned. Above the escarpment, the simulation shows a potential temperature approaching neutral stratification at heights above 40 m a.g.l. The flight pattern is not following the terrain and therefore the lowest available measurement at this position is at 50 m a.g.l. In the middle region, the tower values (green dots) show a superadiabatic layer next to the ground, with a stronger temperature decrease with altitude on the second day. On that day, the potential temperature decreases by 1.5 K at the lowest 100 m. The MASC exhibits a decrease of 0.5 K between 20 and 50 m a.g.l. The agreement between the different systems (MASC and tower) is satisfactory, given the systematic differences in the measurement techniques, as well as in the vertical and temporal resolutions. The simulation shows a strong potential temperature decrease at the ground corresponding to the position of the first computational cell adjacent to the wall. Above the first cell, the decrease is slower and underestimated compared to measured data from the tower and MASC. This may be the result of a low vertical resolution near the ground but also the surface heat fluxes from the mesoscale model, which are underestimated around mid-day, i.e., during the measurement campaigns.
The bulk Richardson number based on a 20 m averaging windows is computed over the flight legs to assess the impact of the escarpment, with its forest, in enhancing critical conditions for shear instability. For the computations, three levels were chosen: a low level near the ground, a medium one and the highest flown altitudes. Figure 15  To summarize the results, it was found: • Both days have a similar flow structure with a wind direction almost perpendicular to the escarpment. The main differences are the wind speed levels: strong and calm wind condition for the first and second day, respectively. The simulated wind speed and wind direction are in accordance with the tower measurements, except for some specific times. During the flight campaigns, values for mean absolute error of 2.00 m s −1 and 0.74 m s −1 are found for the first and second day, respectively. The tower, positioned 60 m behind the forest, clearly shows the impact of the forest with a wind speed in the lower levels (10 m a.g.l.) strongly reduced. The reduction is higher for 22 September, where the wind speed values at the lower cup are about one third of the one at 100 m, on average.

•
The turbulence is evaluated in terms of horizontal turbulence intensity. The model simulates reasonably well at the upper levels but large discrepancies are observed in lower altitudes compared to the tower measurements. This may be directly linked to the canopy model. Possible future improvements require a deeper investigation on the plant canopy and the turbulence model. The influence of the vertical turbulence intensity is not considered in this paper, as no sonic anemometers were available at that time. Further work with the newly installed sonics is required to provide an accurate representation of the turbulence, including the vertical turbulence intensity.

•
The UAS measurements are used for the model validation. In order to avoid any temporal averaging, a real-time strategy is applied, where the model follows spatially and temporally the aircraft. An accelerated flow is numerically and experimentally found over the escarpment. The model slightly over-predicts the wind speeds at higher levels for the first day but still matches very well the UAS measurements. The flow structure remains the same, despite a range of different velocities.

•
Changes in wind direction at turbine heights can be easily measured by the MASC. At the WINSENT test-site, the wind turned counter-clockwise with height for both days, corresponding to a backing wind.

•
Sonic anemometers mounted on a tower can provide the inclination angles. However these measurements are limited to the tower locations, and there can be significant discrepancies between the inclination angle at the tower location and the one at the turbine location. The MASC can overcome this problem and be used to probe the horizontal extent of an area with negative vertical wind, for example. Upward movements over and straight after the escarpment are observed. The inclination angles at an altitude of about 200 m are smaller but still not equal to 0 • , indicating a flow still influenced by the orography and topography. Inclination angles of 5 • and 2 • , at the future turbine location, are found for the first and second day, respectively. The test-site has been intentionally placed at a location that offers high inclination angles.

•
The stability of the atmosphere can be described based on tower measurements, but the UAS offers the opportunity to sample at higher elevations. The calculated and simulated potential temperature profiles are in accordance with a near neutral ABL on the first day. The second day, dominated by a more convective surface layer, at the test-site is not well simulated. Indeed, the superadiabatic layer next to the ground is numerically underestimated. The temperature values at the ground in the OpenFOAM model are obtained by reading the surface heat fluxes from the WRF model and the low resolution (150 m) may be the reason of this underestimation. In future works, it will be interesting to nudge the CFD model with the EC measurements for the ground temperatures. The bulk Richardson values show a systematic drop behind the escarpment for layers near the ground and reveals thermodynamic instabilities.

Conclusions
In the present study, the flow over the WINSENT test-site was numerically and experimentally investigated for two days. Standard measuring systems such as a meteorological tower or an EC station, but also an unmanned aircraft system are used for the numerical validation. The CFD model, driven by the WRF data, relies on an unsteady RANS model with a canopy model for the forested areas. A simple one-way coupling approach shows good agreement with tower and EC measurements for the wind speed and wind direction. However, no acceptable accuracy for the turbulent intensity could be reached at the lower levels, where high wind shears are occurring, and, therefore, need to be more investigated.
The dataset provided by the airborne measurements demonstrates how an UAS system can be helpful for the ABL investigation by providing a high spatial resolution of the airflow near the escarpment. Numerical validation is performed by following the position of the MASC during the flight campaign, which takes more than an hour. In this way, the variation of the flow is caught and no time averaging is needed. A comparison of the horizontal wind speed, wind direction and inclination angle was able to determine the accuracy of our method. Generally, the simulated wind speed, direction and inclination angle match well with the MASC measurements. The stability characterization shows that the model is predicting well for the first day but need improvement for the second day, dominated by more unstable conditions. In the future, a flight strategy following the terrain should be implemented in order to get more information above the escarpment at lower altitudes. The UAS is still an exotic tool for wind resource assessment. This study provides confidence in using UAS as a tool for meteorological diagnostics or for data assimilation in models. A recently published study shows preliminary results from a WRF model, using assimilation of UAS data collected during the LAPSE-RATE Experiment [51]. The applied method shows improved prediction of wind speed and direction and is encouraging for future UAS data assimilation research.
A good level of accuracy is achievable by considering a CFD model, driven by a WRF model, at relatively low computational costs compared to large eddy simulations (LES) or detached eddy simulations (DES). The method presented in this paper will be extended by coupling the OpenFOAM model to a Delayed Detached Eddy Simulations (DDES) solver [33], where the future wind turbines will explicitly be included. In future, intensive observation periods with additional equipment, such as Lidars, and under stronger thermal stratifications will be performed.