Three-Dimensional Wind Field Construction and Wind Turbine Siting in an Urban Environment

Three-dimensional urban wind field construction plays an important role not only in the analysis of pedestrian levels of comfort but also in the effectiveness of harnessing wind energy in an urban environment. However, it is challenging to accurately simulate urban wind flow due to the complex land use in urban environments. In this study, a three-dimensional numerical model was developed for urban wind flow construction. To obtain an accurate urban wind field, various turbulence models, including the Reynolds stress model (RSM), k-ω shear stress transport (SST), realizable k-ε, and (Re-Normalisation Group (RNG) k-ε models were tested. Simulation results were compared with experimental data in the literature. The RSM model showed promising potential in simulating urban wind flow. The model was then adopted to simulate urban wind flow for Purdue University Northwest, which is located in the Northwest Indiana urban region. Based on the simulation results, the optimal location was identified for urban wind turbine siting.


Introduction
Three-dimensional urban wind field construction is important for many activities: weather forecasting, natural ventilation system design, pedestrian comfort level analysis, atmospheric dispersion assessments, and urban wind energy harnessing. However, it is challenging to accurately simulate urban wind flow due to complex land usage and terrain conditions, as well as the variety of building shapes, heights, densities, and layout configurations in urban environments. Stathopoulos et al. provided some views on the potential benefits and challenges in urban wind energy [1]; they reviewed the urban wind energy knowledge and discussed some current issues. Recent progress that was made in urban wind energy assessment was also addressed.
Many researchers focus on urban wind field construction and analysis. The straightforward methods are wind tunnel testing and field measurements [2][3][4][5]. Earlier wind tunnel tests were done by Wiren et al. [2], where they conducted experimental studies of wind velocities in passages between and through block-type building models with various types of layouts. Later, Tsang et al. [3] focused on wind tunnel studies of a pedestrian-level wind environment around tall buildings. More recently, Weerasuriya et al. [4] conducted wind tunnel studies to uncover the effects of twisted wind flows on pedestrian-level wind fields in an urban environment. Finally, Li et al. performed wind tunnel studies and tried to find the correlation between urban space quantification and pedestrian-level ventilation [5].
With the rapid development of computational resource and meshing techniques, increasingly more researchers have used computational fluid dynamics along with measured wind tunnel data to simulate and analyze wind flow in an urban environment. Many of the simulations were conducted by using Reynolds-averaged Navier-Stokes (RANS) models or the large-eddy simulation (LES) approach.

Model Validation
The three-dimensional CFD model was first developed to simulate urban wind flow around a single building and two building test cases. The buildings were configured to follow the wind tunnel experiments conducted by Wiren et al. [2].
For the case concerning wind flow around a single building, the detailed building configuration is shown in Figure 1a. The building had a length (L) of 80 m, height (H) of 18 m, and width (B) of 12 m. At the pedestrian level, there is a pedestrian crossing with a length (B) of 12 m, a width (b) of 6 m, and a height (h) of 4 m in the center of the building.

Materials and Methods
Two parts are presented here: the validation study of the various turbulence models adopted for urban wind field simulations and the urban wind flow analysis and turbine siting for Purdue University Northwest, which is located in Northwestern Indiana, a suburban environment.

Model Validation
The three-dimensional CFD model was first developed to simulate urban wind flow around a single building and two building test cases. The buildings were configured to follow the wind tunnel experiments conducted by Wiren et al. [2].
For the case concerning wind flow around a single building, the detailed building configuration is shown in Figure 1a These wind tunnel tests were measured at the height of 2 m at the middle line of the pedestrian channel. The computational domain size was 480 m × 440 m × 108 m. To avoid boundary influences, the distance from the inlet to the front surface and all of the side surfaces of the building to the boundaries of the computational domain was 10H. To make sure the flow reached a fully developed condition, the distance from the outlet to the back surface of the building was 16H.
For the case concerning wind flow around two buildings, the building layout is shown in Figure  2. Two identical buildings were arranged perpendicular to each other. The buildings were separated by a passage with distance b is 4 m. The wind was blowing perpendicular to the front line formed by the two buildings. These wind tunnel tests were measured at the height of 2 m at the middle line of the pedestrian channel. The computational domain size was 480 m × 440 m × 108 m. To avoid boundary influences, the distance from the inlet to the front surface and all of the side surfaces of the building to the boundaries of the computational domain was 10H. To make sure the flow reached a fully developed condition, the distance from the outlet to the back surface of the building was 16H.
For the case concerning wind flow around two buildings, the building layout is shown in Figure 2. Two identical buildings were arranged perpendicular to each other. The buildings were separated by a passage with distance b is 4 m. The wind was blowing perpendicular to the front line formed by the two buildings. Similarly to the previous case, a sufficiently large computational domain was set up, with a height of 108 m, a length of 480 m, and a width of 380 m. The distance from the inlet to the front surface of the building and the building side surfaces to the boundaries of the computational domain was 10H. For this validation study, the downstream of the building was 16H.
For the two validation cases, the inlet velocity was set up using the power law = ( ) . Based  Similarly to the previous case, a sufficiently large computational domain was set up, with a height of 108 m, a length of 480 m, and a width of 380 m. The distance from the inlet to the front surface of the building and the building side surfaces to the boundaries of the computational domain was 10H. For this validation study, the downstream of the building was 16H.
For the two validation cases, the inlet velocity was set up using the power law U U 0 = Z Z 0 α . Based on Wiren's wind tunnel experiments [2], U 0 (10 m/s) and z 0 (2 m) were the reference velocity and height, respectively; α was 0.125, and U is the velocity at vertical distance z from the surface. The pressure outlet was used as the outflow condition, the ground was set as a non-slip boundary condition, and symmetries were set up for all the side and top boundaries.
Different turbulence models were used in the two validation cases. Detailed information is given in the next session for the governing equations, as well as the comparison for the turbulence models adopted in this study.

Turbulence Models
The commercial software ANSYS ® Fluent 19.1 (Ansys Inc., Canonsburg, PA, USA) was used for the simulations. Reynolds-averaged Navier-Stokes models (RANS models) were used to solve the transport equations of the airflow and simulate all turbulence scales. The RANS models are widely used in urban flow simulations because they yield relatively accurate simulation results and are computationally efficient compared with LES models [18]. In this study, various RANS models, including realizable k-ε, RNG k-ε, k-ω SST, and RSM, were implemented to solve the two validation cases. Brief descriptions of each turbulence model are listed below.
The compressible RANS equation for air can be written in the form of a Cartesian tensor using the summation convention: The k-ε models assume that the turbulence viscosity is linked to the turbulence kinetic energy and turbulent frequency via the relation: In this research, turbulence flow was modeled using the realizable k-ε model, RNG k-ε model, k-ω SST model, and Reynolds stress model (RSM) for airflow to compute the turbulent kinetic energy (k) and the turbulent kinetic energy dissipation (ε).
The exact k-ε equations have many terms that have not yet been calculated. In order to facilitate understanding and calculations, the standard k-ε model was proposed by Launder and Spalding, which can be used as a basic method that is applied to engineering flow field calculations. The kinetic energy and turbulence dissipation rates in this model are described as: where u i is the velocity component in the corresponding direction and µ t is the turbulent viscosity computed using k and ε as follows: where C µ , C 1ε , C 2ε , σ k , and σ ε are constants, which are 0.09, 1.44, 1.92, 1.00, and 1.30, respectively. The realizable k-ε model adds a formula to the turbulent viscosity and adds a new transport equation to the dissipation rate, such that the expressions of k and ε are: where: The RNG k-ε model was derived from a strict statistical mathematical technique to solve the problem for the effect of a smaller Reynolds number flow viscosity. The RNG and realizable k-ε models have shown significant improvements over the standard k-ε model in terms of the flow characteristics, including strong streamline curvature, vortex, and rotation [19]. The equations for k and ε are: where: The k-ω shear stress transport (SST) is a two-equation eddy-viscosity model. The inner parts of the boundary layer are solved using the k-ω model, while the k-ε model was implemented in the free shear flow. This model shows good behavior in adverse pressure gradients and separating flow regions. Menter [20] discussed this method in detail. The appropriate transport behavior can be obtained by formulating the vortex viscosity as follows: Note that the production term of ω is given by: This formulation differs from the standard k-ε model. The kinetic energy and the specific dissipation rate equations are: where: The Reynolds stress model (RSM), known as the Reynolds stress transport model, represents the most complete classical turbulence model. The method of closure used is usually called a second-order closure. This model can explain the complex conditions in turbulent flow and the Reynolds stress's effects. The governing equations for the transport and momentum are written as follows: The transport equations of the RSM can be written as [19]:

Urban Wind Turbine Siting
Purdue University Northwest (PNW) is located in Northwest Indiana, a suburban environment. An urban wind flow simulation for PNW was conducted to help identify the optimal location of a VAWT on campus. Figure 3a shows the top view of the campus environment from Google Maps and Figure 3b shows the layout of the campus's main buildings. The main buildings include the Student Union and Library Building (SULB), Gyte building, Porter Hall, Parking Garage building, Power building, Potter building, Classroom Office building (CLO), and Anderson building. All the main buildings were included in the three-dimensional urban wind field geometry models. To reduce the complexity of meshing and further computational efforts, some detailed building structures were neglected.
The computational domain was 740 m in length, 630 m in width, and 100 m in height. It was selected to avoid the boundary effects around the inflow and sides, and to make sure the downstream reached a fully developed flow region. The distance from the inlet to the SULB front surface was 8H, 16H downstream, 5H on each building side, and about 5H above the building's roof. The tallest building height (H) was 17.27 m. Figure 3a shows the top view of the campus environment from Google Maps and Figure 3b shows the layout of the campus's main buildings. The main buildings include the Student Union and Library Building (SULB), Gyte building, Porter Hall, Parking Garage building, Power building, Potter building, Classroom Office building (CLO), and Anderson building. All the main buildings were included in the three-dimensional urban wind field geometry models. To reduce the complexity of meshing and further computational efforts, some detailed building structures were neglected. The computational domain was 740 m in length, 630 m in width, and 100 m in height. It was selected to avoid the boundary effects around the inflow and sides, and to make sure the downstream reached a fully developed flow region. The distance from the inlet to the SULB front surface was 8H, 16H downstream, 5H on each building side, and about 5H above the building's roof. The tallest building height (H) was 17.27 m.
To monitor the on-campus wind flow conditions, three wind sensors were installed on the roof of the SULB. The height of the SULB was 16.43 m, the location of each data sensor is shown in Figure 4. The boundary conditions applied in the PNW campus wind flow simulation are listed in Table 1. To monitor the on-campus wind flow conditions, three wind sensors were installed on the roof of the SULB. The height of the SULB was 16.43 m, the location of each data sensor is shown in Figure 4. The computational domain was 740 m in length, 630 m in width, and 100 m in height. It was selected to avoid the boundary effects around the inflow and sides, and to make sure the downstream reached a fully developed flow region. The distance from the inlet to the SULB front surface was 8H, 16H downstream, 5H on each building side, and about 5H above the building's roof. The tallest building height (H) was 17.27 m.
To monitor the on-campus wind flow conditions, three wind sensors were installed on the roof of the SULB. The height of the SULB was 16.43 m, the location of each data sensor is shown in Figure 4. The boundary conditions applied in the PNW campus wind flow simulation are listed in Table 1. The boundary conditions applied in the PNW campus wind flow simulation are listed in Table 1. The prevailing wind speed and wind direction were determined based on an analysis of year-long wind data [21] from the three data sensors. Recorded wind data from April 2010 to April 2011 were used. After the prevailing wind information was determined, the inflow wind profile was set up, with the details being covered in Section 3.2.
The outlet was a pressure outlet and the other two sides of the domain were also regarded as pressure outlets. The ground was the bottom of the computational domain and was therefore treated as a no-slip wall. The building walls were considered to be smooth no-slip wall boundary conditions. The top boundary of the computational domain was set as a symmetry condition.

Statistical Analysis
To evaluate the wind energy potential at a given site, statistical analysis can be used. The probability distributions and statistical techniques become especially important if only summary data are available. There are two commonly used probability density functions for wind data analysis: the Weibull distribution and the Rayleigh distribution.
The Weibull probability density function is based on two parameters: k and c, with k being a shape factor, while c is a scale factor. Both k and c are functions of the wind speed standard deviation and the average wind velocity. It is not straightforward to obtain k and c. There are several different ways to do the calculation. In this study, the empirical method was used to calculate k and c: The average wind velocity is calculated as: where the Γ(x) is the gamma function, which is defined as Γ(x) = ∞ 0 e −t t x−1 dt and it can be calculated using the recursion formula: In the Weibull distribution, a higher value of k shows a sharper peak of the probability density function curve, indicating that the wind speed variation is less; on the other hand, a lower value of k indicates a bigger variability about the mean wind speed. Generally, the sites with low shape parameters and high mean wind speeds will make them good candidates for wind turbine siting since these features usually indicate that those sites have the highest total mean power density [22]. The Weibull probability density function is defined as: Compared to the Weibull distribution, the Rayleigh distribution is simpler as it only requires the knowledge of the mean wind speed. The probability density function for a Rayleigh distribution is given as: In the Rayleigh distribution, the larger value of the mean wind speed indicates a higher probability at higher wind speeds for specific sites. Sites with a higher mean wind speed are usually associated with good wind energy potential. The two parameters found in the Weibull distribution can better represent a wider variety of wind regimes [23].

Turbulence Model Selection
To compare the performance of different turbulence models when simulating urban wind flow, different turbulence models were adopted for urban wind flow simulations in both the single-building and two-building cases. The tested turbulence models were the RNG k-ε, realizable k-ε model, k-ω SST model, and second-order RSM. The simulation results were compared with experimental results from wind tunnel tests by Wiren et al. [2].  Figure 5b shows the zoomed-in view of the meshes along the building. The boundary layer meshes that were generated are noticeable, along with the fine meshes around the building for enhanced capturing of the boundary layer flow features. Three body of influence regions were selected: around the pedestrian channel with a mesh size of 0.5 m, around the building with a mesh size of 1 m, and another near the upstream and downstream region with a mesh size of 2 m. The rest of the computational domain mesh size was 10 m. In total, 10 inflation layers were used with a 1.2 growth rate. A mesh-independent study was conducted. The medium-sized mesh was adopted in the simulation for different turbulence models as it brought a good balance between the computational accuracy and cost. The total number of elements used in the final simulation for the medium mesh was 2,189,837.
Compared to the Weibull distribution, the Rayleigh distribution is simpler as it only requires the knowledge of the mean wind speed. The probability density function for a Rayleigh distribution is given as: In the Rayleigh distribution, the larger value of the mean wind speed indicates a higher probability at higher wind speeds for specific sites. Sites with a higher mean wind speed are usually associated with good wind energy potential. The two parameters found in the Weibull distribution can better represent a wider variety of wind regimes [23].

Turbulence Model Selection
To compare the performance of different turbulence models when simulating urban wind flow, different turbulence models were adopted for urban wind flow simulations in both the singlebuilding and two-building cases. The tested turbulence models were the RNG k-ε, realizable k-ε model, k-ω SST model, and second-order RSM. The simulation results were compared with experimental results from wind tunnel tests by Wiren et al. [2].  Figure 5b shows the zoomed-in view of the meshes along the building. The boundary layer meshes that were generated are noticeable, along with the fine meshes around the building for enhanced capturing of the boundary layer flow features. Three body of influence regions were selected: around the pedestrian channel with a mesh size of 0.5 m, around the building with a mesh size of 1 m, and another near the upstream and downstream region with a mesh size of 2 m. The rest of the computational domain mesh size was 10 m. In total, 10 inflation layers were used with a 1.2 growth rate. A mesh-independent study was conducted. The medium-sized mesh was adopted in the simulation for different turbulence models as it brought a good balance between the computational accuracy and cost. The total number of elements used in the final simulation for the medium mesh was 2,189,837.   An incompressible transient simulation was conducted. The time-averaged velocity contour at a height of 2 m from the ground is shown in Figure 6. The pedestrian crossing wind speed was higher than its nearby regions due to the reductions in the cross-sectional area because of a passage at the pedestrian level. An incompressible transient simulation was conducted. The time-averaged velocity contour at a height of 2 m from the ground is shown in Figure 6. The pedestrian crossing wind speed was higher than its nearby regions due to the reductions in the cross-sectional area because of a passage at the pedestrian level.   Table 2 show the comparison between the simulation results from different turbulence models and the experimental data. For the single building validation case, the simulation results from all the tested turbulence models showed overall similar trends to the experimental data. The realizable k-ε and the second-order RSM model showcased better matches. The position of the maximum wind speed was also accurately captured compared with the wind tunnel testing results.   Table 2 show the comparison between the simulation results from different turbulence models and the experimental data. For the single building validation case, the simulation results from all the tested turbulence models showed overall similar trends to the experimental data. The realizable k-ε and the second-order RSM model showcased better matches. The position of the maximum wind speed was also accurately captured compared with the wind tunnel testing results. An incompressible transient simulation was conducted. The time-averaged velocity contour at a height of 2 m from the ground is shown in Figure 6. The pedestrian crossing wind speed was higher than its nearby regions due to the reductions in the cross-sectional area because of a passage at the pedestrian level.   Table 2 show the comparison between the simulation results from different turbulence models and the experimental data. For the single building validation case, the simulation results from all the tested turbulence models showed overall similar trends to the experimental data. The realizable k-ε and the second-order RSM model showcased better matches. The position of the maximum wind speed was also accurately captured compared with the wind tunnel testing results.

Validation for Wind in the Two-Building Case
Similarly, urban wind flow for the two-building case was simulated by using the RNG k-ε model, realizable k-ε model, k-ω SST model, and second-order RSM. The computational meshes for the entire domain are shown in Figure 8a and a close-up view of the refined meshes around the two buildings and the passage are shown in Figure 8b. Similar meshing strategies were used as in the single-building validation case. The total element number for the medium mesh was 2,242,748.
Velocity contours at the height of 2 m from the ground are shown in Figure 9. The stagnation region before the building, acceleration regions at the windward sides of the building, the passage, and the wake region after the buildings can all be observed clearly.

Validation for Wind in the Two-Building Case
Similarly, urban wind flow for the two-building case was simulated by using the RNG k-ε model, realizable k-ε model, k-ω SST model, and second-order RSM. The computational meshes for the entire domain are shown in Figure 8a and a close-up view of the refined meshes around the two buildings and the passage are shown in Figure 8b. Similar meshing strategies were used as in the singlebuilding validation case. The total element number for the medium mesh was 2,242,748.
Velocity contours at the height of 2 m from the ground are shown in Figure 9. The stagnation region before the building, acceleration regions at the windward sides of the building, the passage, and the wake region after the buildings can all be observed clearly.

Validation for Wind in the Two-Building Case
Similarly, urban wind flow for the two-building case was simulated by using the RNG k-ε model, realizable k-ε model, k-ω SST model, and second-order RSM. The computational meshes for the entire domain are shown in Figure 8a and a close-up view of the refined meshes around the two buildings and the passage are shown in Figure 8b. Similar meshing strategies were used as in the singlebuilding validation case. The total element number for the medium mesh was 2,242,748.
Velocity contours at the height of 2 m from the ground are shown in Figure 9. The stagnation region before the building, acceleration regions at the windward sides of the building, the passage, and the wake region after the buildings can all be observed clearly.    The simulation results of the velocity at 2 m high along the central line of the passage between the two buildings were compared with experimental data [2], as shown in Figure 10 and Table 3. Except for the results obtained by using the second-order RSM, results from all other turbulence models did not show good agreement with the experimental data. The simulation results of the velocity at 2 m high along the central line of the passage between the two buildings were compared with experimental data [2], as shown in Figure 10 and Table 3. Except for the results obtained by using the second-order RSM, results from all other turbulence models did not show good agreement with the experimental data.

Turbulence Model Max Relative Error
From the two benchmark cases, the RSM continuously showed good simulation results for wind flow in an urban environment and was therefore adopted in the following urban wind flow simulation and wind turbine siting study.

Wind Turbine Siting in an Urban Environment
A wind rose is a convenient method for displaying the distribution of wind direction and distribution of wind velocity in general statistical analysis. It consists of sixteen equidistant radial lines representing the direction, with the length of the line being proportional to both the frequency of the wind and the velocity of the wind from this direction. The wind rose can also show the relative wind velocity.
The data were recorded by the three sensors on the roof of the SULB penthouse. The locations of the three sensors are present in Figure 4. The overall wind rose for April 2010 to April 2011 on the PNW campus is shown in Figure 11. N, S, W, E stand for north, south, west, and east, respectively. The wind direction of the highest frequency was northeast (NE), with a wind speed magnitude of 4.16 m/s. However, when considering the year-long wind direction distribution, the dominant wind direction was from west-northwest (WNW) to south-southeast (SSE). The second most frequent wind direction, west-southwest (WSW) with a wind speed magnitude at 3.17 m/s, was applied as the input wind conditions for the following wind flow simulation.

Turbulence Model Max Relative Error
Realizable k-ε 25.62% k-ω SST 24.65% RNG k-ε 28.53% RSM 5.6% From the two benchmark cases, the RSM continuously showed good simulation results for wind flow in an urban environment and was therefore adopted in the following urban wind flow simulation and wind turbine siting study.

Wind Turbine Siting in an Urban Environment
A wind rose is a convenient method for displaying the distribution of wind direction and distribution of wind velocity in general statistical analysis. It consists of sixteen equidistant radial lines representing the direction, with the length of the line being proportional to both the frequency of the wind and the velocity of the wind from this direction. The wind rose can also show the relative wind velocity.
The data were recorded by the three sensors on the roof of the SULB penthouse. The locations of the three sensors are present in Figure 4. The overall wind rose for April 2010 to April 2011 on the PNW campus is shown in Figure 11. N, S, W, E stand for north, south, west, and east, respectively. The wind direction of the highest frequency was northeast (NE), with a wind speed magnitude of 4.16 m/s. However, when considering the year-long wind direction distribution, the dominant wind direction was from west-northwest (WNW) to south-southeast (SSE). The second most frequent wind direction, west-southwest (WSW) with a wind speed magnitude at 3.17 m/s, was applied as the input wind conditions for the following wind flow simulation. The average wind speed distribution for each month over the year is shown in Figure 12. April 2011 had the highest monthly average wind speed. The Weibull and Rayleigh distribution results for each month over the year are shown in Figure  13. A higher probability at a higher mean wind speed can be observed for April 2011 in both the Weibull and Rayleigh distributions. The average wind speed distribution for each month over the year is shown in Figure 12. April 2011 had the highest monthly average wind speed. The average wind speed distribution for each month over the year is shown in Figure 12. April 2011 had the highest monthly average wind speed. The Weibull and Rayleigh distribution results for each month over the year are shown in Figure  13. A higher probability at a higher mean wind speed can be observed for April 2011 in both the Weibull and Rayleigh distributions.  These statistical results were used for the wind turbine power extraction evaluation, as shown in Session 3.2.2.

Urban Wind Flow Analysis
The Eddy GT wind turbine from Urban Green Energy (New York, NY, USA) was available to be installed on the top of the PNW campus buildings' roofs. With the advantages of being lightweight and producing a low amount of noise, this small-scale vertical axis wind turbine (VAWT) has been used for residential and commercial applications in various regions in the United States. The rated power for this wind turbine is 1 kW.
The detailed properties and performance information of that wind turbine, as stated by the manufacturer, are given in Table 4. The total height of the turbine was 4.70 m. These statistical results were used for the wind turbine power extraction evaluation, as shown in Section 3.2.2.

Urban Wind Flow Analysis
The Eddy GT wind turbine from Urban Green Energy (New York, NY, USA) was available to be installed on the top of the PNW campus buildings' roofs. With the advantages of being lightweight and producing a low amount of noise, this small-scale vertical axis wind turbine (VAWT) has been used for residential and commercial applications in various regions in the United States. The rated power for this wind turbine is 1 kW.
The detailed properties and performance information of that wind turbine, as stated by the manufacturer, are given in Table 4. The total height of the turbine was 4.70 m. In order to identify an ideal siting location for the wind turbine, three-dimensional wind fields for PNW were constructed. Based on the validation results, RSM was selected for the PNW urban wind flow simulation. As explained before, the inlet velocity profile was defined by using the 1/7 power law as follows: U = U 0 (z/19) 1/7 , U 0 = 3.17 m/s. Three different sets of mesh were used for the mesh-independent study, where the element numbers were 2,673,927, 7,826,295, and 12,095,887 for coarse (surface element size 0.8 m, body element size 5 m), medium (surface element size 0.5 m, body element size 5 m), and fine mesh (surface element size 0.4 m, body element size 3 m), respectively. The simulation results obtained from the medium mesh and fine mesh did not show much difference. The medium mesh was adopted in the simulation as it brought a good balance between computational accuracy and cost. Figure 14 shows the velocity contour at a height level of 21 m. It represents the average height of the selected turbine (half of the rotor height) on different buildings' roofs.  In order to identify an ideal siting location for the wind turbine, three-dimensional wind fields for PNW were constructed. Based on the validation results, RSM was selected for the PNW urban wind flow simulation. As explained before, the inlet velocity profile was defined by using the 1/7 power law as follows: U = U0(z/19) 1/7 , U0 = 3.17 m/s. Three different sets of mesh were used for the mesh-independent study, where the element numbers were 2,673,927, 7,826,295, and 12,095,887 for coarse (surface element size 0.8 m, body element size 5 m), medium (surface element size 0.5 m, body element size 5 m), and fine mesh (surface element size 0.4 m, body element size 3 m), respectively. The simulation results obtained from the medium mesh and fine mesh did not show much difference. The medium mesh was adopted in the simulation as it brought a good balance between computational accuracy and cost. Figure 14 shows the velocity contour at a height level of 21 m. It represents the average height of the selected turbine (half of the rotor height) on different buildings' roofs. At a height level of 21 m, the SULB showed good wind potential. Even though the CLO and Anderson buildings were higher than the SULB, at the height level of 21 m, the wind flow above the CLO and Anderson buildings were at the lower part of the boundary layer flow region, which showed lower wind potential compared to the same height as the SULB. Among all the buildings on At a height level of 21 m, the SULB showed good wind potential. Even though the CLO and Anderson buildings were higher than the SULB, at the height level of 21 m, the wind flow above the CLO and Anderson buildings were at the lower part of the boundary layer flow region, which showed lower wind potential compared to the same height as the SULB. Among all the buildings on PNW's campus, the SULB can be seen as a good candidate for siting the turbine when comparing wind energy potentials.
A side view of the wind flow around the SULB is shown in Figure 15. When the wind blows over the edges of a building, the flow will separate at the leading edge of the roof at a certain angle, referred to as the skew angle. The flow was accelerated because of the separation effects. Furthermore, serving as the foundation of the wind turbine, the special penthouse structure on top of the SULB edge can also increase the wind potential due to the skew angle effects; therefore, an optimal location for wind turbine siting was proposed. PNW's campus, the SULB can be seen as a good candidate for siting the turbine when comparing wind energy potentials. A side view of the wind flow around the SULB is shown in Figure 15. When the wind blows over the edges of a building, the flow will separate at the leading edge of the roof at a certain angle, referred to as the skew angle. The flow was accelerated because of the separation effects. Furthermore, serving as the foundation of the wind turbine, the special penthouse structure on top of the SULB edge can also increase the wind potential due to the skew angle effects; therefore, an optimal location for wind turbine siting was proposed. After considering the simulation results, support capability of the roof structure, existing objects on top of the building roof, available budget, and overall aesthetic look, a small-scale VAWT was installed on the roof of the SULB. Figure 16a shows the final selection of the wind turbine siting location on top of the SULB, while Figure 16b shows a picture of the final installation.

Wind Turbine Performance Analysis
In an ideal situation, the wind turbine should be installed in a site where the rated wind speed has a high probability. Based on the wind data analysis, the wind energy potential on the PNW campus was not very high, as the annual (from April 2010 to April 2011) averaged wind speed was 2.67 m/s, which is lower than the cut-in wind speed. The Eddy GT Wind Turbine will not fully perform under this given wind condition. Over the year (April 2010 to April 2011), there were seven months that had an average wind speed higher than the cut-in wind speed. Wind power was 15ft 11ft Figure 15. Contour of the wind velocity around the SULB.
After considering the simulation results, support capability of the roof structure, existing objects on top of the building roof, available budget, and overall aesthetic look, a small-scale VAWT was installed on the roof of the SULB. Figure 16a shows the final selection of the wind turbine siting location on top of the SULB, while Figure 16b shows a picture of the final installation. PNW's campus, the SULB can be seen as a good candidate for siting the turbine when comparing wind energy potentials. A side view of the wind flow around the SULB is shown in Figure 15. When the wind blows over the edges of a building, the flow will separate at the leading edge of the roof at a certain angle, referred to as the skew angle. The flow was accelerated because of the separation effects. Furthermore, serving as the foundation of the wind turbine, the special penthouse structure on top of the SULB edge can also increase the wind potential due to the skew angle effects; therefore, an optimal location for wind turbine siting was proposed. After considering the simulation results, support capability of the roof structure, existing objects on top of the building roof, available budget, and overall aesthetic look, a small-scale VAWT was installed on the roof of the SULB. Figure 16a shows the final selection of the wind turbine siting location on top of the SULB, while Figure 16b shows a picture of the final installation.

Wind Turbine Performance Analysis
In an ideal situation, the wind turbine should be installed in a site where the rated wind speed has a high probability. Based on the wind data analysis, the wind energy potential on the PNW campus was not very high, as the annual (from April 2010 to April 2011) averaged wind speed was 2.67 m/s, which is lower than the cut-in wind speed. The Eddy GT Wind Turbine will not fully perform under this given wind condition. Over the year (April 2010 to April 2011), there were seven months that had an average wind speed higher than the cut-in wind speed. Wind power was

Wind Turbine Performance Analysis
In an ideal situation, the wind turbine should be installed in a site where the rated wind speed has a high probability. Based on the wind data analysis, the wind energy potential on the PNW campus was not very high, as the annual (from April 2010 to April 2011) averaged wind speed was 2.67 m/s, which is lower than the cut-in wind speed. The Eddy GT Wind Turbine will not fully perform under this given wind condition. Over the year (April 2010 to April 2011), there were seven months that had an average wind speed higher than the cut-in wind speed. Wind power was expected to be generated during these periods. Due to some constraints, the wind turbine was selected even though it did not meet all the desires. In conclusion, the selection of a wind turbine is a complex decision and should be made with consideration of the wind potential, availability of resources, and budget constraints.
After the installation, a more specific analysis that focused on the wind turbine power generation was conducted.
The Eddy GT wind turbine power output curve shown in Figure 17 indicates the amount of power that could be extracted from the wind turbine at different wind speeds. The extracted power can be calculated using PT(V w ) = P r is the wind speed, P r is the wind turbine rated power, V 1 is the cut-in wind speed, and V r is the rated wind speed) when the wind speed is greater than the cut-in speed and smaller than the rated speed. If the wind speed is greater than or equal to the rated wind speed and less than the cut-out speed, the rated power will be extracted from the wind turbine. expected to be generated during these periods. Due to some constraints, the wind turbine was selected even though it did not meet all the desires. In conclusion, the selection of a wind turbine is a complex decision and should be made with consideration of the wind potential, availability of resources, and budget constraints. After the installation, a more specific analysis that focused on the wind turbine power generation was conducted.
The Eddy GT wind turbine power output curve shown in Figure 17 indicates the amount of power that could be extracted from the wind turbine at different wind speeds. The extracted power can be calculated using ( ) = 3 − 1 3 3 − 1 3 ( is the wind speed, is the wind turbine rated power, 1 is the cut-in wind speed, and is the rated wind speed) when the wind speed is greater than the cut-in speed and smaller than the rated speed. If the wind speed is greater than or equal to the rated wind speed and less than the cut-out speed, the rated power will be extracted from the wind turbine. Based on the power output curve of the wind turbine and the wind velocity probability density on the installation site, the estimated power generation of the wind turbine can be calculated. As shown in Figure 18, the total power generated over the year (April 2010 to April 2011) was 202 kWh using the Rayleigh distribution compared to a slightly higher value of 204 kWh obtained using the Weibull distribution. Due to the low wind energy potential in the selected urban site, the extracted power was much lower than the energy extraction at the rated power.  Based on the power output curve of the wind turbine and the wind velocity probability density on the installation site, the estimated power generation of the wind turbine can be calculated. As shown in Figure 18, the total power generated over the year (April 2010 to April 2011) was 202 kWh using the Rayleigh distribution compared to a slightly higher value of 204 kWh obtained using the Weibull distribution. Due to the low wind energy potential in the selected urban site, the extracted power was much lower than the energy extraction at the rated power. expected to be generated during these periods. Due to some constraints, the wind turbine was selected even though it did not meet all the desires. In conclusion, the selection of a wind turbine is a complex decision and should be made with consideration of the wind potential, availability of resources, and budget constraints. After the installation, a more specific analysis that focused on the wind turbine power generation was conducted.
The Eddy GT wind turbine power output curve shown in Figure 17 indicates the amount of power that could be extracted from the wind turbine at different wind speeds. The extracted power can be calculated using ( ) = 3 − 1 3 3 − 1 3 ( is the wind speed, is the wind turbine rated power, 1 is the cut-in wind speed, and is the rated wind speed) when the wind speed is greater than the cut-in speed and smaller than the rated speed. If the wind speed is greater than or equal to the rated wind speed and less than the cut-out speed, the rated power will be extracted from the wind turbine. Based on the power output curve of the wind turbine and the wind velocity probability density on the installation site, the estimated power generation of the wind turbine can be calculated. As shown in Figure 18, the total power generated over the year (April 2010 to April 2011) was 202 kWh using the Rayleigh distribution compared to a slightly higher value of 204 kWh obtained using the Weibull distribution. Due to the low wind energy potential in the selected urban site, the extracted power was much lower than the energy extraction at the rated power.

Discussion
A three-dimensional CFD model was developed to simulate wind flow in urban environments. Different turbulence models were adopted to simulate urban wind flow. The simulation results were compared with wind tunnel testing data. The results showed that among all the models tested (RSM, k-ω SST, realizable k-ε, and RNG k-ε models), RSM gave the best matches with experimental data in both benchmark cases. Therefore, RSM was selected to simulate the urban wind flow and wind turbine siting for Purdue University Northwest, which is located in an urban area in Northwestern Indiana.
A one-year wind data analysis result was used to set up the inlet wind profile. Statistical wind data analysis was performed using Rayleigh and Weibull distributions, and the results were used as the wind turbine energy production estimation for the selected wind turbine.
3-D suburban wind characters were studied and analyzed based on wind fields around the SULB. Flow over a building was separated at the roof windward leading edge at a certain skew angle. Wind flow became accelerated due to the flow separation, which may enhance the wind energy potential in that area. Besides serving as the foundation of rooftop wind turbines, buildings may also play an active role in helping to gain more wind energy output in an urban environment.
Furthermore, when considering the wind turbine siting on top of the building roof, a taller building has the benefit of reaching a higher wind power density. However, in order to take advantage of this feature, the selected turbine needs to reach a certain height. Otherwise, due to the lower part of the boundary layer flow effects, a relatively low building may be a better choice for a low-height wind turbine siting.
Small-scale wind turbines are usually lightweight, have low vibrations, and reduced noise. They are widely used in urban environments. Many factors need to be considered for urban wind turbine siting, such as the selected site wind potential, the probability of the wind speed compared to the cut-in wind speed and rated wind speed, the height of the wind turbine vs. the height of the buildings, the building roof structure, and existing objects. The final decision should fit within the budget limits.