1. 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. Blocken et al. [
6] reviewed wind-tunnel and computational fluid dynamics (CFD) techniques to address their accuracy for wind comfort assessment. A comparison between the RANS model and LES for urban wind simulations was provided. Rodi used various
k-
ε models, Reynolds stress models, and LES models for urban wind simulations [
7], where the results were compared with experimental data. The performance, cost, and potential of solving complex flow problems using those turbulence models were also compared. Rodi pointed out that the use of RANS calculations underestimated the turbulence fluctuations in general, while the use of LES gave better simulation details of the turbulence flow with the cost of a large increase in computing time. In the LES model, the unresolved small-scale turbulent motion was modeled, which was less influenced by the boundary conditions. Johnson and Hunter simulated urban wind flow using the
k-
ε model and revealed the broad flow features [
8]. Compared with the wind tunnel experiments, a general agreement between the numerical results and experimental data was demonstrated. Willis investigated and tested the CFD validation guidelines in an urban wind simulation [
9], while Song and He conducted a numerical simulation using an LES model when evaluating pedestrian wind flow [
10]; in their studies, different geometry and wind conditions were adopted and simulation results were compared with wind tunnel experimental data.
Different variations of RANS models were widely used by many others for simulating urban wind flow: Santiago and Martilli used the standard
k-
ε model when simulating wind flow over a regular array of cubes [
11]; Zheng et al. conducted the pedestrian-level wind environment study via wind tunnel tests and realizable
k-
ε simulations [
12]; Gnatowska et al. used the realizable
k-
ε model to simulate flow fields in an urban environment, where the results were compared with Particle Image Velocimetry (PIV) data in a wind tunnel and a detailed description of the velocity field in the built-up environment was obtained [
13]; Akashi et al. employed the Re-Normalisation Group (RNG)
k-
ε model to simulate the urban wind flow under different building layout effects [
14]; Hu et al. applied the RNG
k-
ε model during a quantitative investigation of the building density effects on residential wind environments [
15]; and Liu et al. used the RNG
k-
ε model to simulate the wind distribution in an urban community with a full-scale geometrical model coupled with a microscale model [
16].
Among these modeling works, LES simulations can reproduce more accurate details of the turbulence flow, but because of its large computational cost, most researchers choose to use a RANS model for urban wind flow simulations instead. In urban wind flow analysis, understanding the flow features within the passages is more important than capturing the detailed turbulence features near wall regions. The RANS models satisfy these requirements and are more computationally affordable than an LES. Reiter compared the standard
k-
ε model, realizable
k-
ε model, and RSM model in urban wind simulations [
17], but some other popular RANS models, such as the RNG
k-
ε model and
k-
ω shear stress transport (SST) model were not included. The comparisons between different RANS models for urban wind flow simulations are limited in the literature.
There are two main parts to this paper. First, various popular RANS turbulence models, including the RNG k-ε model, realizable k-ε model, k-ω SST model, and RSM, were used to simulate wind flow in an urban environment. The simulation results were compared with wind tunnel experimental data. Among all those models, the RSM showed the most promising results for urban flow simulations, with the simulation results closely matching the experimental data. Second, wind turbine siting in an urban environment was investigated for Purdue University Northwest campus. By implementing the selected turbulence model coupled with a year-long wind data analysis, an optimal location was selected for a small vertical axis wind turbine (VAWT).
2. 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.
2.1. 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.
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 on Wiren’s wind tunnel experiments [
2],
U0 (10 m/s) and
z0 (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.
2.2. 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
is the velocity component in the corresponding direction and
is the turbulent viscosity computed using
k and
ε as follows:
where
,
,
,
, 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]:
2.3. 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.
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 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.
2.4. 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
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].
4. 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.