An Energy Potential Estimation Methodology and Novel Prototype Design for Building-Integrated Wind Turbines

: ROSEO-BIWT is a new Building-Integrated Wind Turbine (BIWT) intended for installation on the edge of buildings. It consists of a Savonius wind turbine and guiding vanes to accelerate the usual horizontal wind, together with the vertical upward air stream on the wall. This edge effect improves the performance of the wind turbine, and its architectural integration is also beneﬁcial. The hypothetical performance and design conﬁguration were studied for a university building in Eibar city using wind data from the ERA5 reanalysis (European Centre for Medium-Range Weather Forecasts’ reanalysis), an anemometer to calibrate the data, and the actual small-scale behavior in a wind tunnel. The data acquired by the anemometer show high correlations with the ERA5 data in the direction parallel to the valley, and the calibration is therefore valid. According to the results, a wind speed augmentation factor of three due to the edge effect and concentration vanes would lead to a increase in working hours at the rated power, resulting annually in more than 2000 h.


Introduction
In general terms, the market for small wind turbines is currently growing, although the sector of small wind turbines intended for installation in buildings is increasing at a lower rate. According to the World Wind Energy Association (WWEA) [1], the installation of small wind turbines will increase by around 12% annually in the 2015-2020 period. The good economic profitability of small wind turbines and the consistency of technological advancement are determinant factors that explain the growth of the small wind turbine market. On the other hand, in the last several years, research is increasingly being focused on the development of different technologies that help minimize the energy consumption of buildings. This philosophy is known as nZEB (nearly Zero-Energy Building) [2], and it is included in the EU 2010/31/CE directive related to the energy efficiency of buildings. After 2018, every new public building should be constructed in accord with this regulation and, after 2020, every new building should be compliant.
The goal is to maximize energy efficiency and reduce the use of primary energy derived from fossil resources so that the required energy demand can be met by renewable sources. In this sense, mini wind technology, which involves generating energy with wind turbines of 100 kW or less to cover an 1.
They can work as standalone devices, so they can provide energy in isolated locations without a connection to the electric grid. 2.
They work in distributed micro-generation mode, thus minimizing energy losses due to transport and distribution. These devices generate energy at a site that is close to the final user, thus dramatically reducing the need for electric infrastructures.

3.
Furthermore, it can be combined with photovoltaic energy in hybrid installations to enable the optimal use and management of shared electric accumulators.
The recent developments in wind energy for urban environments have inspired different types of Building-Integrated Wind Turbine (BIWT) projects. For example, in London, Strata SE1 is a tall building with 43 floors that will include three wind turbines with diameters of 9 m on the roof of the structure. These wind turbines will be used to meet the building's lighting demand [3].
On a smaller scale, there are a lot of projects that include Horizontal-Axis Wind Turbines (HAWTs) integrated with buildings, as well as Vertical-Axis Wind Turbines (VAWTs). These projects are focused on integrating wind turbines with existing buildings. Thus, these buildings were not previously designed to accelerate air streams, unlike the World Trade Center of Baharein [4] or the mentioned Strata SE1 building. According to this post-integration trend, building-integrated wind turbines are being implemented in strategic locations to capture the acceleration of air streams that are produced because of different geometries. In this sense, the most interesting locations are the upper and lateral edges of a building, especially the former because it is at a reasonable distance from homes.
Nowadays, there are several ongoing projects working to develop an optimal system that harnesses wind energy in urban environments. Most of them have concluded that wind turbines located in obstacle-free environments are not adequate for urban environments because of the urban turbulent flow, which can present a relevant turbulence intensity on the superior edges of the buildings [5,6]. For that reason, HAWT devices, which usually exhibit good performances with laminar flows, perform poorly in urban environments, in addition to their generation of noise as high as 200 dB within a radius of 500 m [7]. Conversely, VAWTs play an essential role in generating wind energy in urban areas since their performance is not much affected by turbulent flows, and they tend to be noiseless [8]. Additionally, the VAWT has a lower cut-in speed than HAWT and a larger or even unlimited cut-off speed, ensuring longer operating times [9][10][11]. Although the power coefficient is lower, the design is simpler and the manufacturing process is easier to carry out.
Along these lines, the existing urban wind energy potential has encouraged researchers to develop a proper methodology for wind energy estimation in urban environments [12]. The use of anemometers at specific locations can be combined with advanced computational simulations of buildings situated in complex urban terrains using CFD (Computational Fluid Dynamics). In this way, wind energy potential estimation using reanalysis and meteorological mesoscale models, which is a well-known offshore and onshore method and also developed by the authors [13,14], can be complemented with different back-end tools.
In this work, the authors present the design of a Savonius drag-driven turbine that is intended for integration into buildings. The proposed turbine is called ROSEO-BIWT, which has been specially designed to work in urban environments. The wind in urban areas is characterized by its turbulence, thus it is important to take advantage of low-speed air streams. The germinal project of ROSEO won the first award in the EDP-RENEWABLE UNIVERSITY CHALLENGE 2017 [15], and the members of the project have now created a university start-up called ROSEO. Although it is typically used as a vertical-axis turbine, ROSEO-BIWT is formed by a Savonius turbine in a horizontal position and concentration vanes that accelerate the air streams by the Venturi effect (see Section 2.2). These types of vanes are usually called PAGVs (Power Augmentation Guiding Vanes) [16][17][18]. The proposed turbine was also designed to be easily architectonically integrated. This was the case for the design proposed by Park et al. [19], in which several Savonius turbines were incorporated into the facade of a building at different heights to take advantage of the vertical currents created by the wind on the walls of the building.
The Savonius wind turbine is a drag-based device, unlike the majority of turbines, which are lift-based. This particular aspect allows for low noise levels and few vibrations, and these factors are very important in building installations [20,21]. The PAGV increases the wind speed as the catching area grows, resulting in a system that is able to start at wind speeds of about 1 m/s, thus ensuring a great number of energy-producing hours. Furthermore, energy generation continues no matter how high the wind speed is.
This paper proceeds as follows: a possible location for the installation, which was established using ERA5 data, is presented. ERA5 is a powerful tool for global atmospheric analysis that is updated in real time (see Section 2.1). The authors also installed an anemometer on the roof of their university to calibrate the wind data for a period of eight months against ERA5 (Section 2.1.2). In this way, an empirical method for the estimation of wind energy potential on buildings with a low computational cost will be developed in subsequent work, as discussed in Section 3.5. Finally, a preliminary small-scale experiment was developed for a wind tunnel with a small Savonius and different configurations of the PAGV (Sections 2.2 and 2.3). The authors finish this work with some relevant conclusions and a future outlook of possible research directions. The qualitative methodology used here can be considered within the scope of analogical reasoning and model construction [22].

Anemometers and ERA5
The university building of Eibar (Engineering School of Gipuzkoa) was selected (longitude: 2.946 • W; latitude: 43.258 • N) to demonstrate a method for formulating a preliminary estimation of energy production using quantile-matching calibration versus a cup anemometer installed on the roof. Figure 1 shows a satellite view of the engineering school and the position of two anemometers installed on two buildings. After eight months of data acquisition, Anemometer 1 showed the best correlations with the ERA5 reanalysis, and its dataset was used for the calibration and energy estimation procedure.
The ERA5 reanalysis, ECMWF's most recent atmospheric reanalysis, covers the second half of the 20th century and this century [23]. For this study and the calibration, 40 years of data were used (from 1979 to 2018), because the complete reanalysis is not yet available. The correlation between ERA5 and the anemometers was computed within their period of intersection (from June 2018 to February 2019). These data include atmospheric and oceanic variables, and are an appropriate tool for estimating wind energy potential [24] offshore and onshore. In this study, ERA5 hourly data with a resolution of 0.3 • × 0.3 • was used.

Quantile-Mapping Calibration
The cup anemometers installed on the buildings of the University of Basque Country in Eibar enabled the development of a preliminary calibration methodology based on quantile-matching techniques that were used previously by the authors for wind energy and wave energy [25][26][27]. In the scientific literature, different calibration or bias correction techniques have been developed and compared for the analysis of several parameters, such as temperature and precipitation (see [28][29][30]). Data from models and reanalysis are compared with observations. In the present study, a simple but effective statistical procedure based on quantile mapping was used.
For this approach, several other terms can be found in the literature: "probability mapping" [31], "quantile-quantile mapping" [32,33], "statistical downscaling" [34], and "histogram equalization" [35]. With this general approach, empirical quantile-mapping bias correction was applied to calibrate ERA5 versus an anemometer in the building. In [36], the same procedure was used for estimating wind energy trends. To summarize, this method of calibration or bias correction is fundamentally statistical, and the idea is to match values with the same quantile in two empirical probability distributions: the one to be calibrated (ERA5), and the one that is the basis for the calibration (anemometer). Figure 2 illustrates the main aspects of this calibration procedure, including the intersection periods and the concept of applying the transference function.
For this paper, the authors obtained an eight-month, 10-min data series, which was filtered every 6 h to match the ERA5 reanalysis for a 10-year period (1-h time resolution, in this case). Thus, there were around 34,500 cases in the anemometer time series and around 87,600 cases in the ERA5 series. Taking 1-h data for both series in the intersection period resulted in 5390 cases, from which the correlation was measured and the subsequent calibration transference function was generated.
Having determined the average wind speed U after calibration on the corresponding facade and considering the typical shape parameter of the Weibull distribution (Rayleigh distribution, k = 2), the corresponding scale parameter can be obtained: (1) Then, the cumulative distribution function and the fraction of time between two wind speeds are determined: and the augmentation factor AF of the PAGV (the ratio between the outlet and inlet free wind speed) can be incorporated into the c parameter [9] because it is proportional to the average wind speed U. The value of AF can only be based on a virtual definition of the outlet velocity of the flux, since the complex interaction between the diffusive flux in the exterior part of the vanes and the motion of the rotor do not permit a simplistic application of the Venturi effect according to the relation between the capture width of the free wind U and the outlet width. However, the optimum tip speed ratio (TSR opt ) at which the power coefficient C p is maximized is directly related to the outlet effective velocity (U out ), because it is well known that TSR opt0 ≈ 1/3 for a drag turbine without augmentation techniques [9]. Due to the Magnus and lift effects in the Savonius rotor, this value can reach 0.4-0.5. Therefore, the augment of TSR opt should be similar to AF considering an effective U out at the position of the rotor in relation to the blade tip speed V tip . Being V A tip the tip speed in the augmented rotor and TSR A the tip speed ratio in the augmented rotor, the hypothesis is that the TSR should be the same for the conventional Savonius and for the augmentation technique if the outlet velocity AF · U is the reference wind speed: However, the TSR with augmentation TSR A should be defined with respect to the free wind speed: Consequently, AF can be also computed using the ratio of the TSR with augmentation versus the TSR without augmentation.
On the other hand, experiment using nozzles by Shika et al. [37] have shown that AF can be 4 or even 5 measuring directly U out at the position of the rotor for U between 0.6 and 0.9 m/s. This relevant increment for low free wind speed is very interesting for our purpose, since reducing significantly the cut-in speed of the rotor. Furthermore, these AFs ensure a great quantity of working hours at rated power, as shown below.

The Location on the Upper Edge of the Building
The effect of wind against buildings has been largely studied by the architectural sector for the purpose of studying the dynamic loads generated by air streams. Because of their work, there is considerable knowledge about the behavior of wind in urban environments. Much of the information has been obtained through experiments with scale models and CFD simulations, similar to the depiction generated by the authors in Figure 3, which was re-created based on the CFD simulation in [38]. Most of the studies that have analyzed the behavior of air streams around buildings agree that the upper edge of the windward face of a building has great wind energy potential. This is because the wind has to surround an object. The effect is even more intense when the building is taller and when the wind direction is perpendicular to the building facade. For example, in a five-story building, the wind velocity increases by 1.2 times at the windward edge [38].
According to the CFD simulations of Balduzzi et al. [5,6], the wind speed increment at the edge can be between 10% and 30%, but the turbulence intensity increases considerable. This is not the worst inconvenience for Savonius turbines, since it is demonstrated that, when turbulence increases, the separation of boundary layer takes place on the lower side of returning blade of the rotor reducing the negative torque [39,40].
Areas of high turbulent intensity create more frequent and stronger gusts [9], but the inertia of a relatively long Savonius rotor (high aspect ratio between the length of the axis and the diameter) can keep the rotation of the turbine without relevant variations. Additionally, it is demonstrated that, as a consequence of blade tips, aerodynamic losses are reduced in Savonius turbines with high aspect ratios [41].

Savonius Turbine
Because of the above-discussed wind behavior, our ROSEO-BIWT's Savonius axis is positioned horizontally along the superior edge of the building. There have been several recent studies on the performance of the Savonius turbine. Mohamed et al. [42] improved the performance using plates to eliminate the negative torque in the returning blade. They carried out tests for a two-bladed and a three-bladed wind turbine, and, in both cases, they improved the power coefficient (C p ) of the wind turbine by up to 27%, with 15% being the typical value.
Apart from these intrinsic improvements, some engineers have developed the mentioned PAGV systems to accelerate air streams. Shikha et al. [43] increased the wind speed by 3.7 times in an experiment using a specific well-studied nozzle. Additionally, Altan et al. [44] studied the influence of the inclination angle of the plates as well as their length. In these experiments, they found that, when the longitude of the PAGV increased, the power also increased. Thus, the important consideration in their study was the relationship between the diameter of the rotor and the length of the PAGV. They even obtained a C p of 38.5%. Other types of PAGVs, referred to as omnidirectional, reached a C p of 48%, implying an increase of 240% relative to a Savonius rotor without a PAGV system.
In terms of longitude and diameter, the Savonius rotor studied in [45] performed best with an aspect ratio of 6:1. Similarly, Park et al. [19] tested different kinds of Savonius rotors, and they discovered that the best design was a six-bladed rotor. Therefore, for our purpose, a similar rotor with these proportions was chosen for the initial test period.

The Final Design
Park et al. [19] developed the idea of using a larger facade surface to generate energy by installing a lot of Savonius rotors at different heights while also using PAGVs to improve the performance of the wind turbines. The system that they proposed is similar to a ventilated facade. It is important to emphasize that they wanted to capture the vertical air streams that are generated on the windward side of the building, as in our case. However, they used parallel vanes in the facade with a small concentration angle; in our case, the upper edge is used to augment the concentration angle and capture not only the vertical stream on the facade but also the horizontal component of the wind. Furthermore, the background of the Savonius rotor is free on the edge of the building: this is an important aspect that is not encountered in turbines installed in the facade. Although the influence of this aspect is out of the scope of this study, it is an obvious aerodynamic advantage.
Another innovation is that the proposed turbine can be installed in existing buildings: it is not a design intended only for new buildings. To summarize, ROSEO-BIWT shows good architectural integration in existing buildings and high economic viability due to the simplicity of the design. Figure 4 shows the ROSEO-BIWT design. This is a schematic perspective that does not take into account the influence of the angle between the two PAGVs; the angle can be adapted for other positions of the vanes. The PAGV areal ratio between the entrance and exit of the air is 4:1, and a similar AF is expected to result from a first simplistic calculus due to the Venturi effect. In any case, as mentioned, the complex interaction between the outlet wind speed and the rotor motion deviates this a priori value of AF = 4.
In their seminal work about a curtain design to increase the performance of a Savonius turbine, Altan et al. [46] determined that the best angles for the capture of wind in their curtain design are 15 • for the superior vane and 45 • for the inferior one that obstructs the negative torque. The authors corroborated the same influence of the inferior vane in the laboratory (see Section 3.1). Figure 5 shows the dimensions for this optimum design with curtains. It is considered a unit of capture width at the inlet, and a geometrical relation of 4:1 for the inlet width (0.25, therefore diameter of 0.50) versus the outlet width. The aspect ratio is six considering the results of Roy and Saha [45]: 3 = 0.5 × 6. These proportions can be established between a capture width of one and two meter; within this size, ROSEO device is manipulable for a worker on the roof in the implementation process and for O&M issues, without the need of a crane .

Experiments in the Wind Tunnel
Although there are results provided by previous studies, in the following sections, the authors describe the general experimental methodology that is being developed. The experimental model construction is proposed by referencing previous findings and design procedures [22]. These are the main steps: 1.
First, according to the literature, the augmentation factor of the wind speed on the edge of the buildings is around 1.2. Wind speed augmentation is the result of the union between the usual horizontal component and the vertical component.

2.
Then, the previous augmentation factor should be multiplied by the new increment AF provided by the vanes. These factors will be measured for different wind speeds in the wind tunnel of the university using a small-scale model of a building with curtain-type vanes (see Figure 5) and a rotor or 2 cm diameter.

3.
A similar experiment will be performed for a real Savonius with one inferior vane and will be critically compared with other studies.

4.
Finally, the Weibull distribution at the location obtained by the previously described calibration methodology will be applied to the measured power curve that includes AF. Thus, the amount of hours at rated power due to this augmentation will be an interesting parameter about energy production. Table 1 describes the main characteristics of the above-mentioned wind tunnel. Figure 6 shows the wind tunnel and the installation of a PAGV and a real Savonius rotor. It should be mentioned that the disposition of the vane below the limit of the rotor's horizontal axis obstructs the negative torque and, simultaneously, accelerates the stream in the opened drag side above the axis. As mentioned, this aerodynamic effect has been properly documented in previous reviews about the performance of the Savonius rotor [39,41,47]. The augmentation may be even higher because of the corner effect of our design. However, until now, these preliminary measurements have only been performed with low values of steady wind speed without considering some important effects, such as the blockage ratio of the tunnel [49,50]. However, the the influence of AF is important at these low wind speeds below the rated power, because it can ensure a sufficient wind speed above the rated wind speed at the outlet of the concentration vanes.
On the other hand, Figure 7 shows the small-building, the guiding vanes with 3:1 inlet/outlet relation, and the six-bladed rotor of 2 cm diameter. Here, the objective is to create an anemometer that is able to capture the wind on the entire outlet area of the vanes. According to our previous experiments, measurements with Pitot tubes result in great fluctuations due to small displacements or inclination deviation in such a narrow area. In this case, the electrical motor is a maxon DCX06M EB KL 6V of 0.529 W. The speed constant is of 3060 min −1 V −1 having a direct way to compute the angular velocity in function of the voltage. Additionally, there is also a constant speed-torque relation of 36,600 min −1 mNm −1 . This characteristic is important because it allows computing the increment of the torque due to the increment of the speed.

Effect of the PAGV in the Real Savonius
AF of around 2 has been corroborated in the experiments of the wind tunnel. This obviously depends on the angle and the length of the vanes and the exact position of the Pitot tube, but this value of AF can be obtained with a suitable disposition of the vanes. However, this measurement is strongly influenced by the exact position and size of the Pitot tube, and future works should study the behavior of AF for higher and more turbulent wind speeds.
Although a maximum outlet/inlet width ratio of 1.4 can be obtained due to the lack of space in the tunnel with the Savonius rotor in the center, the instantaneous power measured for different wind speeds and different angles of the vane in Figure 6 gives a coherent result for this augmentation and subsequent power. Figure 8 shows this behavior with the vane at 30 • , 45 • and 70 • with respect to the horizontal. The pilot test was developed without the vane and the negative-torque wall (NT wall) was applied with the vane in vertical position, obstructing the negative torque's drag. It should be noted that, removing the negative torque, the power is doubled and the other cases (30 • , 45 • and 70 • ) also remove the negative drag.
The results in the curves of Figure 8 show the best working condition for the vane at 45 • , in which the captured power almost triples the pilot test power at each wind speed. Being the geometrical augment relation of 1.4, and 1.4 3 ≈ 3, the power also keeps the typical proportionality relation with U 3 for AF, and, again, a constant AF equal to the geometrical relation is deduced. This fact establishes an important particular case for an hypothetical law that should be demonstrated: for the adequate vane angles, the estimation of power production can be performed using AF × U as the input wind speed for any free wind speed U. The behavior of C p vs. TSR was also studied and is presented in Table 2 for the optima with which AF can be estimated. The presence of the vane increments TSR opt from 0.5 to around 1 with AF ≈ 2. The best case in power augmentation (almost three times) for AF = 2.2 is measured for the vane at 45 • . These results are totally coherent with previous works using different augmentation techniques that obtain TSR opt ≈ 1 compared to a value of 0.5 for a conventional Savonius. In these cases, C p,max is also two or even three times higher thanks to the vanes, deflectors, curtains or other kinds of concentration configurations [44,51]. Additionally, as mentioned above, Altan et al. also already showed for their curtain type augmentation technique that the inferior vane should be at 45 • to optimize the energy capture [52].

Augmentation Factor in the Small-Scale Building Model
In this case, TSR for each wind speed is measured instead of the TSR opt . Consequently, the ratio of the TSRs with (TSR v ) and without (TSR 0 ) the vane should be corrected according to the rotor speed, since the torque is incremented with the speed. Figure 9 shows these results: both TSRs, their ratio, and the corrected ratio that equals AF. This correction factor is established by the increment relation between U and the cut-in wind speed for the pilot experiment (5 m/s). A logical step in the procedure considering the constant speed-torque relation of the DC generator, since both TSRs (therefore, the rotor speeds) are practically linear with respect to U and also fulfill the same increment relation. Thus, AF is between 2.5 and 3, a very relevant result given the fact that the inlet-outlet width relation is 3:1. The 4:1 relation of the initial configuration could not be installed yet, due to the sensible construction details of the small model. Because of this delicate structure, the fabrication process of which has been very laborious, the free wind speed range in the tunnel has been kept below 10 m/s. Figure 9. TSRs, the ratio, and the corrected ratio for the rotor speed with and without the vane.

Wind Rose around the Building
The ERA5 grid around Eibar city is shown in Figure 10 with the ERA5 points in blue. The building on which the anemometer is located is marked in red. The nearest grid point, at a distance of 2.48 km, was chosen to perform the calibration. With the wind rose representing the nearest ERA5 grid point and the anemometer (see Figure 11), it is easy to realize that the ERA5 data have to be calibrated to make an appropriate estimation. ERA5's wind rose shows a strong predominant direction toward the northwest, as it is well-known that the climate of the Basque Country is highly related to the behavior of geostrophic winds [14]. This predominant direction is perpendicular to the valley in Eibar, and it is clearly diminished by the roughness of the terrain and the obstacle of the mountains in the anemometer's data. In fact, Eibar is an industrial city with a population of 20,000 in a deep valley surrounded by mountains that are around 600 m from where the River Deba opens toward the northeast direction.
Thus, this big difference could be explained by the shape and direction of the valley in which Eibar is located. It shows the need for field measurements and indicates that a good calibration methodology must use atmospheric reanalysis to study wind potential in places such as cities and deep valleys, which have high surface roughness. Furthermore, the valley direction determines not only the calibration direction for energy estimation purposes, but it also defines which of the facades of the building should be selected for the implementation of the BIWT. Figure 11. Details of the location of the anemometer and the nearest ERA5 grid point (left). Representation of the corresponding wind roses of ERA5 and anemometer data (right).

Comparison between ERA5 and the Anemometer
These data should be established at a referential height using the log law and the roughness of urban environments [9]. The ERA5 grid-point height is 411 m, thus both datasets should be established at the same height, which is the anemometer's height in this case, since it is the observation.
According to usual considerations in the wind energy sector, the roughness (z 0 ) of the urban terrain is around 1-10 m. Roughness is used to apply the logarithmic law of vertical wind shear, which results in a correction factor between 0.86 and 0.77 for a wind speed at a height of 178 m in ERA5. In terms of the speed of reference at its original height, the authors calculated the following: A correction factor of 0.86 was used prior to the calibration method, which is based on quantile mapping. However, first, the correlation between ERA5 and the anemometer had to be directionally studied, mainly in the direction of the valley line. Furthermore, anemometer data had been previously filtered using advanced filters in meteorology, such as temporal checks, persistence tests, and climate-based range tests [53], which were implemented in the R programming language by the authors [54]. Figure 12 shows a time series of a week in June 2018 when the Pearson's correlation between ERA5 and the anemometer was very high (around 0.95); the wind direction vectors are illustrated above each time point that shows a strong westerly component. The parallel patterns shown by the wind speed series are obvious in the graph. These examples verify the quality of the anemometer's data since their results are comparable to the reliable ERA5 data in the predominant direction line established by the valley (southwest-northeast). In the first days of this week, a cut-off low occurred in the Bay of Biscay, and it caused strong wind and a large amount of precipitation in that area. During the following days, without the influence of the cut-off low, wind moving in the north direction was observed. This is a global-scale synoptic situation that is easier to detect by the ERA5 model than local setups. Therefore, a high correlation between the observed data and ERA5 data was confirmed.
Additionally, the approximation to observation of the corrected ERA5 signal resulting from the application of the log law is clear in the time series using the 0.86 factor, but it is not enough to totally correct the general overestimation presented by ERA5. Although an extreme correction for a high roughness z 0 = 10 m with a factor of 0.77 would strongly reduce this overestimation (Equation (6)), the usual roughness values for urban environments are kept in this graph.
This example is an extraordinary case, but, if all the cases of wind between the south and east were selected during the study period, a good correlation of 0.70 would finally be obtained. This validation was therefore enough to justify the calibration in this directional range, from which the corresponding building facade will be selected for capturing wind energy.

Estimation of the Energy Potential
On the basis of the resource assessment results of the wind potential in buildings and the described wind tunnel experiments, a general methodology is presented that also references previous results from the scientific literature to estimate the annual energy production (AEP) of ROSEO-BIWT: 1.
According to Mertens [38] and the initial experiments with our small-scale building in the wind tunnel, the wind increases its velocity by 20% at the upper edge of a typical building.

2.
The simplest PAGVs have increased the wind speed by four times, with a corresponding increase in C p to a value as high as 0.37 [47]. Although higher values can be obtained with wider entrances, the authors will use an AF of 4 for the estimation, although there is also a 20% augmentation due to the additional architectonic acceleration at the upper edge.

3.
Taking into account the wind rose in Figure 11, the authors only considered the wind data of ERA5 for the valley direction and for our turbine on the corresponding facade.

4.
AF ≈ 3 has been corroborated by our small-scale building with PAGVs for different wind speed values in the wind tunnel. Although the optimum vane angle experiment has not yet been developed, the first test results are consistent with values reported in the literature.

5.
AF ≈ 2 has been corroborated by the real Savonius with the inferior vane. 6.
The analyses of the wind resource in the open direction of the valley and the corresponding facade yield a wind speed histogram or Weibull distribution that can be applied to the power curve of the turbine with AF.

7.
For the first estimation presented here, the working time at rated power due to the increment of wind speed using PAGVs has been computed.
Although the results of the comparison and the calibration of data are very relevant, the objective of this paper is mainly methodological, and a preliminary estimation of AEP should be made using a well-known device. Thus, for the estimation of generated power, a commercial Savonius model (SeaHawk-PACWIND) was used: the rated power is of 1.1 kW, the rated wind speed is 17.9 m/s, the cut-in wind speed is 3.1 m/s, the cut-off without a given limit is a drag device, and the swept area is 0.92 m 2 [55].
When the pure ERA5 wind speed distribution was considered for the best facade, after the quantile-matching calibration using the anemometer data, these are the preliminary results: • The average wind speed is 4.2 m/s, and the shape factor k is around 2, depending on the angle range in the predominant direction, i.e., southwest (see Equations (2) and (1)).

•
The turbine's working hours per year in the interval of rated wind speed (above 17.9 m/s) can be computed if the cumulative density function F(U) (Equation (2)) is applied to c, which results in the following working hours: [1 − F(17.9)] × 365.25 × 24 = 160 (7) • Therefore, AEP is 160 × 1.1 = 170 kWh at the rated power; it is a very small value since the working hours of a profitable turbine should be around 2000 h per year. • However, multiplying the scale parameter c by values between 2 and 4 (AF = 2 is the value obtained in the laboratory using only one inferior vane and 4 the maximum expected value according to the mentioned literature) and keeping the typical value of k = 2, the total augmentation factor AF from the PAGVs and the edge effect increases the AEP and working hours. Figure 13 shows the annual working hours (W H hours ) at rated power in function of the average wind speed U of the site for different factors: AF = 2; 2.5; 3; 3.5; 4.
At low annual average wind speed of 3 m/s, the maximum AF = 4 can produce 2000 h at rated power. At U = 4 m/s, AF between 2.5 and 3 is necessary to ensure the 2000 h. At U = 5 m/s, the minimum AF = 2 obtained with only the inferior vane ( Figure 5) is almost sufficient. At high Us, an AF = 3.5 or 4 implies 75% of the time at rated power. Given the strong directionality presented by the anemometer's wind rose (Figure 11), two ROSEO-BIWTs installed in opposite facades of the building that are perpendicular to the valley would capture almost all of the mentioned hours because winds that are perpendicular to the valley are infrequent.

Conclusions and Future Outlook
An integral methodology with preliminary results is presented for a new type of BIWT. The preliminary results include the energy potential estimation, measurements of small-scale building aerodynamic effects, and the influence of PAGVs. In the future, an AEP increase of 20% via PAGVs at the edge of the buildings must be demonstrated using a real prototype of ROSEO-BIWT at the edge of the building in Eibar. For that, the building in the city of Eibar will be used in the Bizia Lab project of the University of Basque Country, together with the previous wind tunnel experiments for the mentioned small-scale building with PAGVs and short Savonius prototype.
The anemometer was installed on the roof; with the new data provided by ERA5 for the nearest grid point, the identification of the best facade and the corresponding wind distribution were obtained following the methodology described in this paper applied to a longer period. This methodology will be relevant when the one-year period has elapsed. These preliminary results and the methodological discussion developed to date encourage us to implement future refinements of ROSEO-BIWT and the related wind energy estimation methodology.
If the building edge effect and the PAGVs produce a wind speed augmentation of AF, our general mathematical proof for working hours at the rated power shows that the hours without the augmentation can be considerably incremented. This is a very important general rule for turbines with vanes, as shown in Figure 13. Furthermore, the values of AF between two and four are coherent with the literature and the experimental results, even under-valued, since the augment of the free wind in the edge of the building is not considered. The edge effect augment factor of 1.2 documented by the literature and the higher inlet-outlet width relationship other type of PAGVs could increase the overall AF.
Additionally, a novel validation method for anemometers developed by the authors in a recent study for wind farms [56] will be very beneficial since it enables the comparison and combination of data from more than one anemometer installed on the roof of the building. This allows us to consider both the zonal and meridional components in a single comparison score.
Future experiments in the wind tunnel with a small-scale building and the PAGVs will be carried out to obtain the optimum value of the angle between the vanes and the augmentation factor for different wind speeds within the operating range of the turbine. The augmentation factors measured with an interval of 0.5 m/s within this range, together with the measurement of the power curve of the longitudinal profile of the Savonius with the same step, will allow us to apply the corresponding histogram distribution of the corrected and calibrated wind to the augmented power curve. However, it is expected that future energy production results will be similar to the values presented here.
Finally, it should be emphasized that lacking the ability to change the viscosity of the air in the tunnel is an important inconvenience for small-scale building experiments. In the future, a more advanced tunnel with the ability to change the pressure and temperature is necessary, together with a parallel validation of the results using CFD simulations of the edge effect of the building with PAGVs.

Funding:
The research leading to these results was carried out in the framework of the Programme Campus Bizia Lab EHU (Campus Living Lab) with a financial grant from the Office of Sustainability of the Vice-Chancellorship for Innovation, Social Outreach and Cultural Activities of the University of the Basque Country (UPV/EHU). This programme is supported by the Basque Government. We acknowledge also the availability given by the School of Engineering of Gipuzkoa-Eibar in the University of Basque Country, the EDP-Renewable awards in which we obtained the main award in September 2017, the Youth Enterprise Grant of UPV/EHU, and the project GIU17/02 of EHU/UPV. All computations and representations of this work were developed using the programming language R [54].

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: