Analysis of Accelerating Devices for Enclosure Wind Turbines

Copyright: © 2017 Sedaghat A, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. *Corresponding author: Ahmad Sedaghat, School of Engineering, Australian College of Kuwait, Mishref, Kuwait, E-mail: a.sedaghat@ack.edu.kw Received: March 10, 2017: Accepted: August 04, 2017: Published: August 07, 2017


Introduction
From the classical fluid mechanics it is well known that nozzles accelerate and diffusers decelerate steady internal flows with constant mass flow rate. In such devices, the outlet to inlet velocities will be proportion to the ratio of inlet and outlet areas; hence, there would be appreciable conversion factor for these fluid devices. For example, in design of wind tunnels, nozzles and diffus-From 1981 forward, attention was also drawn more in mathematical and computational modelling of shrouded wind turbines. Fletcher [4] has used blade element computational method on diffuser augmented wind turbines including wake rotation and blade Reynolds number effects. His computational results matched well with some experiments on power coefficient and axial velocities. Maximum power augmentation of four times was reported using solidity factor of 0.1 ≤ σ ≤ 0.15 depending on the aerofoil section used for the wind turbine blades. Nagai and Irabu [5] have also developed the momentum based theory for diffuser augmented wind turbine in Japan.
Helmy [6] further developed theoretically shrouded wind turbines assisted with rear edge flanged diffuser. His studies revealed that power coefficient was doubled and he suggested optimal dimension of diffuser and the flange. In parallel in UK and Denmark, Phillips, et al. [7] used CFD modelling and the development of the diffuser augmented wind turbine and Hansen, et al. [8] studied effect of placing a diffuser around a wind turbine. Bet and Grassmann [9] studied computationally two diffuser rings around the wind turbine rotor and reported that by creating a field of low pressure behind the wind turbine the power coefficient was doubled. design wind turbines by devices to accelerate wind velocity at rotor plane. Because of fundamental differences of external flows over nozzles and diffusers compared with internal flows, a large number of literatures focuses on diffuser or shrouded wind turbines. The orientation of accelerating devices within shrouded wind turbines were mainly horizontal for improvement of horizontal axis wind turbines; but, no study were conducted to study the flow accelerating devices in vertical orientation. With recent growth interest in using wind catchers in architects of modern buildings and applications of renewable energies, there is opportunity to develop wind generators within wind catcher architecture.

Backgrounds in shrouded wind turbines
Research works exploiting applications of wind turbine inside shrouded aero generator may be dated back as early as 1976 when Igra [1] built and tested two wind turbines in wide range of axial and angular inlet velocities for several numbers of blades. His investigation also included wind turbine design in enclosures which indicated that his simple model can be used for design of shrouded wind turbines. His experiments for a load factor below 0.3 showed good agreement with his model for a range of solidity factor (1.57 ≤ σ ≤ 6.28). Eventually he observed that for low load factors (below 0.3) the number of blades above 6 will not improve power performance. Igra [1] also summarized that for higher load factors the viscous effects should be incorporated within his model and for higher Reynolds number than narrow band of his experiments (10 5 ≤ Re ≤ 2.6 × 10 5 ) from similarity laws the power performance can be better. A relatively constant power capture for wide range of a velocity factor reported. Igra [2] reported research and development of shrouded wind turbines with significant power augmentation using a compact shroud and up to 80% extra power improvement by using ring shape flaps while additional 25% enhancement were also observed using appropriate bleeding from exterior of shroud. He also presented design of the axial wind turbine for the proposed shroud. A pilot wind turbine with 1 hp power at 5 m/s with 3 meters in diameter were built and its field test results were compared with models ( Figure 1). He reported measurements of wind power which shows augmentation of two times compared with shroud less wind turbine. The power augmentation was related to two important factors: the in viscid wake recovery causing subatmospheric pressure at the shroud exit plane and the kinetic energy conversion into pressure in the rear part of diffuser. From economical point of view also the given concept proved to be more efficient than conventional wind turbines.
The new idea of using diffuser assisted wind turbine spread in US, Japan, and Europe. In US, Gilbert Foreman  tex formation and interaction occur inside the diffuser. Their optimal design claimed to produce 4 to 5 times more power compared with their bare wind turbine (500 W). Both wind tunnel tests and field tests performed as expected. Ohya, et al. [14,15] introduced a new diffuser design with a broad-ring brim at the exit in several compact dimensions. Two to five times wind power augmentation was achieved using new compact designs due to low-pressure region mainly produced by strong vortex formation behind the broad brim. This brings more mass flow rate to the wind turbine inside the diffuser shroud.
In another work, Matsushima, et al. [16] studied a similar diffuser-flanged wind turbine using field testing and also computation of the bare diffuser-flanged without wind turbine. From computation, they reported the maximum increase of wind speed by a factor of 1.7 of the free stream velocity using appropriate diffuser length and expansion angle. They also reported fastest air speed near the inlet of diffuser. The field test results also confirmed that the wind power increased by up to 2.4 times compared with the bare wind turbine. They however argued that their numerical results overestimated the maximum wind speed ratio and output power of 1.7 and 5, respectively, compared with measured values. Therefore, they suggested that their wind turbine had negative effects when introduced to the diffuser or the frequent and rapid changes in wind speed and wind direction has caused the discrepancy leading lower performance of the wind turbine generator than computed optimal one.
In UK, Wang, et al. [17,18] reported on aerodynamic design of a small urban wind turbine with scoop. CFD simulations of the wind turbine and detailed design for blades, nosecones and nacelles were given. Estimation of the annual power output of this wind turbine calculated using hourly wind data of Scottish Weather Station was compared with other commercial turbines and reported as rather competitive. The final design of the scoop wind turbine boosted wind speed and power output by a factor of 1.5 and 2.2 times, respectively, compared with the bare wind turbine of same swept area confirmed also by wind tunnel tests. The results also indicated that power capture improved particularly at locations where average wind speed was lower and wind was more turbulent. Grant, et al. [19] also considered urban wind power using a roof-mounted ducted wind turbine to use pressure difference by wind flow around a building in UK. A simple mathematical model presented and compared with outcomes from tests at model and prototype scale. It was concluded that ducted turbines have, for in the most promising applications, power coefficients well in excess of the conventional Betz limit. Another advantage was low visual impact and safety.
In Taiwan, Hua and Cheng [20] designed a vertical Later in Japan, Ohya, et al. [10,11] have developed a diffuser-flanged shroud system around small wind turbines. The flanged-diffuser shroud job was to collect and accelerate the free stream wind. The flange at the exit of the diffuser shroud generated a low-pressure region in the wake region of the diffuser by vortex formation and entraining more air into the wind turbine inlet. The shrouded wind turbine was optimized for which wind power augmentation of about 4-5 times is reported compared with the bare wind turbine.
Abe, et al. [12] have carried out an extensive experimental and numerical study on the same wind turbine with a flanged diffuser of Ohya, et al. [10,11] (Figure 2). They observed four times augmentation of wind turbine power coefficient using both hot-wire anemometry measurements and computations. In addition, they noticed when the power performance of diffuser-flanged was normalized by a mean velocity behind blades, both the bare and diffuser-shrouded wind turbines presented similar performances. Thus, they conclude that in both situations the wind turbine operates similarly and the power augmentation was mainly related to the acceleration of approaching air by the flanged diffuser. Wind tunnel experiments and filed measurements indicated that the output performance was as expected.
In 2008, Ohya, et al. [13] have developed their optimized diffuser shroud wind turbine with a broad-ring flange ( Figure 2) at the exit of diffuser. Their emphasis was on placing the exit flange dimension and position so that a low-pressure region generated behind of the wind turbine in the wake of rotor where a strong tip vor- steadiness effects. Their results show that the interaction with upstream bowed wakes has a beneficial effect on casing end wall flow behaviour resulting in weakening secondary flows. The combined unsteady effects on turbine performance depend on the balance of unsteady end wall mixing and end wall secondary flows and the negative impact of these effects should be alleviated to achieve enhanced performance.

Review of wind catchers
Wind catchers (Badgirs) are important infrastructure part of the Iranian traditional architecture for supplying natural conditioned air to buildings in both dry hot and humid hot weather conditions for thousands of years. The wind catchers not only used for air ventilation of buildings but also used for cooling in water cisterns (Ab-anbar) [25,26]. The natural ventilation of buildings is the tradition of ancient Persia buildings to allow wind entry from tall towers for cooling and conditioning air above a basement pond into rooms and hall of building. The wind catchers are widely used in architecture of desert type buildings in the central and the southern parts of Iran and are known as Bawdgir [27][28][29][30]. The wind catchers configurations are varied among one sided, two sided, four sided and eight sided towers. The design of the wind catcher is varied from one city to another city. For example, the design of wind catchers in the city of Yazd are taller with four to eight sided configuration compared with the nearby city of Maybod with lower height and one sided types [31]. The reason for this difference may be related to the fact that Yazd city is located among mountains with weak winds while Maybod is located in a flat desert with frequent strong dust storms. In general, wind catchers with short and one sided designs or tall and with four to eight sided designs are widely adapted to suit air conditioning of building in different meteorological conditions. The wind catchers are traditionally built for heights from just 2 meters to 20 meters or higher on top of buildings using wood-reinforced masonry constructed by clay and sun dried clay bricks. Some authors suggest that taller towers capture higher wind speeds with lower dust entry [27,30,32,33]. Figure 3 shows a typical four and eight sided wind catcher in Dowlat-Abad house and garden in the city of Yazd. The house was built in the 18 th century and was once the home of the Persian regent, Karim Khan Zand [34].
Recently, many researches were conducted to assess operation and performance of wind catchers in Iran. Montazeri [35] has investigated wind catchers with different number of opening using smoke visualization and computational techniques. His results suggest that the bucket-shape ducted wind turbine. They argued that conventional ducted wind turbines required large inlets to absorb more mass flow rate; however, the task should be increasing the speed of wind. They optimized geometry of the ducted wind turbine using object-oriented optimization and CFD simulations. The optimal shape of their duct appeared to be an unconventional nozzle, which increased wind speed by 60%. The results of their field tests of their prototype ducted wind turbine showed power augmentation of 1.6 to 1.8 times the bare wind turbine. Chen, et al. [21] also developed a shrouded, small, horizontal-axis wind turbine installing on moving vehicles. They studied the effects of flanged diffusers on rotor performance of a small (30 cm rotor diameter) wind turbine with different rotor solidities (20-60%) and wind speeds (10-20 m/s). Wind tunnel experiments confirmed that the flanged diffuser may significantly increase the power output of the wind turbine for higher solidity of 30% to 40% and low wind speeds. Their study also showed that a small wind turbine has the characteristics of low torque and high rotational speed, and high rotor solidity for maximum power output compared to a conventional large wind turbine.
In Australia, Kosasih and Tondelli [22] studied shrouded micro-wind turbine by examining diffuser shape and geometries, blade airfoils, and wind condition of a mounting site. The effect of different configurations namely: straight diffuser, nozzle-diffuser combination, and diffuser-brim combination were investigated. Test results confirmed that straight diffuser improved wind turbine performance by 60%, the nozzle-diffuser enhancement of 63%, and brimmed diffuser also improved performance. Increasing the diffuser length (L/D) did not affect the optimum Cp of the wind turbine but shifted the performance curve and the optimum Cp to higher tip-speed-ratio. But adding brim (flange) at the exit plane of the diffuser increased the performance, Cp as well as reduced the cut in speed and shifted optimum to higher value. Jafari and Kosasih [23] also reported experimental and Computational Fluid Dynamic (CFD) studies for diffuser shrouded an AMPAIR300 wind turbine which demonstrated significant power augmentation. Their studies showed strong dependency on the shape and dimensions of the diffuser. Different diffuser shapes and dimensions were investigated to understand their effects on power augmentation. From the CFD results, sub-atmospheric back pressure is found to be the most influential factor in power augmentation. It was argued that flow separation on the diffuser surface led to significant reduction in the pressure recovery coefficient and so the overall power augmentation.
Gao, et al. [24] have computationally studied internal unsteady flow in a 1.5-stage shrouded turbine with straight and bowed vanes. Their aim was to control of casing end wall secondary flows and losses due to un-  Figure 4). The innovative feature of INVELOX is by combining wind tunnel and wind catcher systems to capture wind flow through a vertical enclosure intake and hence eliminating the need to yaw control. The flow accelerates through a shrouded Venturi and exhaust into atmosphere by a diffuser.
CFD simulations and field testing were conducted to evaluate performance of this wind turbine. The results of their studies showed that it is possible to capture, accelerate and concentrate the wind. The total average wind power augmentation of INVELOX over 8 days was reported about 314% compared with bare wind turbine. The cut-in wind speed of 1 m/s is promising for development of enclosure wind turbines integrated within wind catcher designs for low speed regions in hot deserts.
The Allaei's idea seems to stem on the concept of wind catcher designs which used to generate natural cooling comfort through drift of wind into building. Wind catchers considered as the tradition of Persian architecture which was developed and used for many centuries. However, there is a fundamental difference between crucial parameter in performance of the studied wind catchers is the number of openings. Kalantar [36] research in the city of Yazd showed that the evaporative cooling is the main cooling mechanism in this dry hot city and he suggested equipping the wind catchers with the water vaporization system to reduce considerably air temperature inside buildings. Bahadori, et al. [37] investigated two new designs for wind catchers, with wetted column and wetted surfaces, and compared their performances with a conventional wind-catcher in the city of Yazd. The two new designs have improved temperature reduction and increased relative humidity compared with traditional one.
Kazemi Esfeh, et al. [38] applied smoke visualization to assess wind catcher models with flat, inclined and curved roofs. His study showed superior performance using the wind catcher with a curved roof.

Development of an enclosure wind turbine
Most recently in US, Allaei and Andreopoulos [39] have introduced the new concept of enclosure wind tur-  found that the best concentrator should locate 5 cm in front of the wind turbine and it should be a nozzle shape with inlet to outlet ratio of 6.
The aim of present study is to exploit external flows over vertical accelerating devices such as nozzles and diffusers for application in such wind generators. The paper is outlined as to give a background study on shroud-diffuser wind power augmentation followed by theoretical modelling of the accelerating devices. Then, the results of mathematical modelling are compared with results of CFD simulations and some available experimental data on horizontal axis accelerating devices. Next, computational results are presented for vertical nozzle and diffusers at wide range of wind speeds. Finally, conclusions and guidelines are provided for future development of vertical accelerating devices in enclosure wind generating systems.

Mathematical Modelling
Literature discussed above has merely reflected on simple models for wind turbine integrated with shroud diffusers. Some authors suggested that the main cause of power augmentation was due to shroud diffuser and others considered interaction between wind turbine and shroud was responsible for this. Ohya, et al. [13] examined three bare hollow circular nozzle, pipe, and diffuser geometries, with no wind turbine and area ratio of 1/4, 1, and 4, respectively, and L/D = 7.7 in wind tunnel experiments as shown in Figure 5. Air flows with smoke strip visualization technique were used through inside and outside of the three model devices. Results indicated that air flows avoided the nozzle-type model instead air better flows into the diffuser-type model because air was sucked into diffuser due to the back suction side. It was also shown that wind speed gradually increased with Allaei's model as shown in Figure 4 and the traditional wind catchers. As seen in Figure 3, the wind is drawn into a big space from a uniform cross section canal. This is more likely act as a diffuser than a nozzle. The openings of the wind catcher are narrow and much smaller compared with room space where wind is drawn. Therefore, the above literature of research on diffusers and the proven concept of wind catchers may suggest that a enclosure wind turbine should utilize a diffuser to collect and expand air into the designs such as INVELOX and then the concentrated wind should be accelerated in a contraction nozzle and an expanding diffuser (or Venturi) similar to design of blower wind tunnels to augment wind speed. The advantage of such design compare with Ohya's diffuser assisted horizontal wind turbine is that the wind can enter into the structure from any directions without the need for yaw mechanism control. Design of enclosure wind turbine can also integrate with wind catchers to be incorporated in architecture of modern buildings with natural air conditioning systems [31,40].
Bussel [41] have reviewed 50 years of Diffuser Augmented Wind Turbines (DAWTs). Bussel concluded that using 1D momentum theory power extracted from a DAWT is the same as a bare wind turbine. He reported that power augmentation beyond 3 has never achieved. However, Wong, et al. [42] reviewed various flow augmentation methods for the purpose of utilizing with Vertical Axis Wind Turbines (VAWTs) and reported large magnification of 9 times output power. Chen, et al. [43] has reviewed innovative wind turbine concepts. They believe the main challenges of new innovative wind turbines are not only augmentation of power but also economical and environmental friendly aspects of designs. Shonhiwa and Makaka [44] has reviewed the concept of Concentrator Augmented Wind Turbine (CAWT). They Energy equation for air flows inside pipe between point 1 and point 2 is expressed by [45]: From continuity equation, the volume flow rates at inlet and outlet of the pipe must be equal; hence, for constant cross sectional area of the pipe, the inlet and outlet average velocities are equal; i.e. 1 2 V = V , although the kinetic energy coefficients are not equal because of different velocity profiles at inlet and outlet for this undeveloped fluid flow; i.e. α 1 ≠ α 2 . At inlet however, the inlet velocity is assumed nearly uniform, V 1c ≈ V 0 , for which α 1 = 1, and at the outlet, the velocity profile is as- ; α 2 = 1.2245. From Ohya, et al. [13] experiments, the outlet pressure can be nearly approximated by the free stream pressure; i.e. P 2 ≈ P 0 ; hence, equations (1) and (2) can be simplified as: Equating right hand sides of (3) and (4), the relation between inlet and free stream velocities are obtained as L/D; hence, it was concluded that the diffuser structure is most effective for collecting and accelerating the wind. A remarkable increase in wind speed ratio was obtained at inlet of diffuser, between 1.6 to 2.4, if a long diffuser over L/D = 3 is used; however, it was preferable to have a shorter diffuser body with a L/D of less than 2 for their shroud wind turbine so that they examined several short diffuser-type structures with end flanges to achieve high velocity ratio at inlet.
Ohya, et al. [13] presented experimental wind velocity and static pressure distribution on the central axis of the hollow circular devices. These data will be used for validation of our simple mathematical model and also to verify our CFD computations. Below, the simple mathematical models are presented for horizontal devices. The same idea can be explored when using the devices in vertical orientation; however, it is avoided here due to lack of experimental data; instead, more reliable CFD simulation are carried out to compare if there is any appreciable difference between the accelerating devices with horizontal and vertical orientation.

Horizontal pipe flows
Here, simple mathematical models based on energy equation by some head losses are developed for external flows over pipe, nozzle, and diffuser at locations 0, 1, and 2 as shown in Figure 5.
In general, the external flows over and inside the horizontal pipe is highly complicated due to inlet and outlet interaction with free stream flow, developing inner and outer boundary layers outside and inside pipe, entrance pipe undeveloped flow, turbulent flows inside and outside pipe. To simplify this, it is assumed that the inlet effect is to raise entrance pressure so that pressure gradient along inside pipe allows for flowing air through the pipe. Hence, air velocity at inlet and through pipe will be lower than free stream velocity. By writing energy equation from point 0 to point 1, including inlet minor loss with   [13] suggested that a high pressure zone created at the nozzle inlet avoided air to easily enter the nozzle and this limited mass flow rate passing through the nozzle. It was observed that air velocity at the outlet of nozzle did not increased from the value at upstream value and hence the nozzle was not considered as suitable device to accelerate air speed. The inlet velocity was also considerably reduced compared with upstream velocity. From these observations, air velocity and pressure coefficient at inlet and outlet of the nozzle can be estimated using continuity and energy equations. Writing energy equation from point 0 to point 1 for the horizontal nozzle including the nozzle inlet loss (k 1 ) as follows [45]: Energy equation for air flows inside the nozzle between point 1 and point 2 can be expressed by [41]: The nozzle loss coefficient k 2 for a typical nozzle may be found in reference [45]. From continuity equation, the volume flow rates at inlet and outlet of the pipe must be equal; hence, for the area ratio of 4 between inlet and outlet of the nozzle, the inlet and outlet average velocities are given by, i.e.
Ohya, et al. [13] experiments, the outlet centreline velocity is nearly equals to upstream velocity, i.e. V 2c ≈ V 0 , for which α 2 = 1 Similar to the pipe flows, the kinetic energy coefficients at inlet and outlet are different so that by assuming nearly uniform flow at inlet, we can write For the air velocity of V 0 = 5 m/s, the Reynolds number value is Re = 4.1 × 10 4 , hence by using Moody diagram, the friction coefficient for a smooth pipe is obtained as: f = 0.021. For a typical loss coefficient of k 1 = 0.02 (from reference [45]), although the loss coefficients should be obtained from experiments, the inlet average velocity is calculated from (5) as: 1 0 = 0.84 V V , and consequently the outlet centreline velocity is obtained as: V 2c = 1.03 V 0 . This is in good agreement with experiments as shown in Figure 6. Although for each case, the coefficients of the head losses should be determined from experiments.
In order to determine pressure coefficient, for the inlet and outlet centreline of the pipe, equation (4) is rewritten as follow to obtain pressure coefficient: For the pipe with the given values above, C p1 = 0.237, C p2 = 0. This slightly overestimates experiment as shown in Figure 7.

Horizontal nozzle flows
As shown in Figure 5 found in reference [45]. Taking k 1 = 0.4, k 2 = 0.02, A 2 /A 1 = 4 although the loss coefficients should be obtained from experiments, the average outlet velocity is determined as:  Figure 6).
Inserting the above findings into equation (7) and manipulating to find the inlet pressure coefficient, the following is found: Substituting the above values, the inlet pressure coefficient is determined as C p1 = -2.02 . The pressure coefficients at inlet and outlet are shown in Figure 7 which nearly matches with experiments.

CFD Simulations
In order to develop the vertical enclosure wind turbine, the components such as pipe, nozzle, diffuser, and Venturi should be studied when oriented in vertical direction. From the mathematical models developed in section 2, it may be discovered that the potential heads due to elevation between inlet and outlet play an important role. However, it is difficult to test such devices in vertical arrangement because most of widely used wind tunnels are designed to test objects with horizontal wind. Hence, the objective of this section is to evaluate performance of such devices in vertical wind using Computational Fluid Dynamics (CFD) simulations.
In this work, CFD simulations are performed for solving air flows over the Ohya's nozzle and diffuser outlined in section 2. The purpose is to verify and validate CFD solutions with the available wind tunnel test data [13]. The commercial CFD software, Fluent package [46], is used in this study to solve the Reynolds average Navier-Stokes (RANS) equations with (k-ω)-SST turbulence model for the three dimensional Ohya's nozzle/diffuser model. The computational domain consists of unstructured mesh with 452,000 mesh elements. The mesh inside the nozzle/diffuser arrangement was constructed using hexagonal elements as shown by Figure  6. Outside the nozzle/diffuser arrangement, the hybrid pipe flows, the velocity profile is assumed turbulent with . Taking n = 7 for turbulent flows, the average velocity in outlet is given by: Inserting the above findings into equation (7) and manipulating to find the inlet pressure coefficient, the following is found: Substituting k 1 = 0.02 [45] and 1 , although the loss coefficients should be obtained from experiments, the value of pressure coefficient at inlet is found to be C p1 = 0.94. Similarly by substituting the above given data in equation (8) and manipulating the outlet pressure coefficient is obtained as follows: Substituting above values and k 2 = 0.02 [45] in equation (10), the outlet value of pressure coefficient is obtained as C p1 = 0.268. The pressure coefficients at inlet and outlet are shown in Figure 7 which exactly matches with experiments.

Horizontal diffuser flows
As shown in Figure 5, the external and internal flow over an open diffuser is considered here. A mathematical model is developed based on energy equation considering inlet and interior head losses. Experiments by Ohya, et al. [13] suggested that a low pressure zone developed at the diffuser inlet which accelerated air flow inside the diffuser increasing mass flow rate passing through the diffuser. It was observed that air velocity at the inlet considerably increased by a factor of 1.6 to 2.4 times of the upstream velocity. From these observations, air velocity and pressure coefficient at inlet and outlet of the diffuser can be estimated using continuity and energy equations. The same energy equations as the nozzle given in equations (7) and (8) can be used for the diffuser. Assuming the outlet pressure as atmospheric, P 2 = P 0 , and using continuity equation, , the equations (7) and (8) are combined to find outlet velocity as follows: The loss coefficient k 1 for a typical diffuser may be diffuser model, the CFD results agrees very well half way through the diffuser with the wind tunnel test data; but, it fails considerably the reminder of diffuser and downstream of flow because of strong viscous interaction and possible separation zone with vortex formation in the diffuser and its downstream. It appeared that the CFD results and our simple analytical model agree well for both inlet and outlet of diffuser flow. Therefore, the simple external flow over diffuser may cast an interesting example for CFD model validations to address some physical phenomena underneath flow features of such typical devices. From Figure 8 as expected, the pressure coefficient has exactly matched up with experimental data for both nozzle and diffuser flows. This indicates that pressure field has perfectly obtained using CFD simulation. However, the mismatch in velocity field should be addressed to the flow viscosity and shear stresses obtained from the studied turbulence model. This can be separately investi-mesh element was used. The boundary conditions utilised in the current study are depicted in Figure 6. The solution is deemed converged if the difference between two consecutive iteration was less the 10 -4 for all involved parameters (velocities, k,ω). The Semi Implicit Method for Pressure Linked Equation (SIMPLE) algorithm was employed for the calculation of pressure and velocity fields. Third order MUSCAL discretization scheme was used with momentum, turbulent kinetic energy, and specific dissipation and energy equations. However, the second order upwind discretization scheme was used with the pressure mode. Furthermore, the values of Y+ at the walls of nozzle/diffuser ranged from 20-85.
The CFD results for the centreline velocity and the pressure coefficient distribution are shown in Figure 7 and vertical direction of wind. The interesting finding from the results in Figure 9 is that the nozzle perform as discussed in previous section by retarding flow and hence increasing pressure at inlet and then regain almost the same upstream velocity at the outlet. This clearly shows that there is fundamental difference between internal nozzle flows which accelerate upstream velocity by the nozzle area ratio. However for the open nozzle flows with internal and external flows, the external flows dictate that the velocity retard at inlet and then regain to the external upstream velocity value.
Changing the upstream air velocity from 1 m/s to 5 m/s, the effects on the centreline velocity and the pressure coefficient are shown in Figure 10 and Figure 11; respectively, for the vertical nozzle. The nozzle inlet is lo-gated and is out of scope of present study.
For the purpose of present study, two devices including circular nozzle and diffuser in vertical arrangement are examined using FLUENT software [29]. The nozzle and diffuser used here are with the Dimensions of: Length of L = 0.85 m, smaller diameter of D = 0.25 mso L/D = 3.4, and the area ratio of A 2 /A 1 = 6. The different geometries are adopted here for a practical application for which shorter lengths of the devices are required.

External flows over nozzle
Here, CFD simulations around the vertical nozzle in external flows are presented. As shown in Figure 9, computational results for the nozzle suggest that there is no appreciable differences between lateral (horizontal) or  there. The general structure of pressure contours around the nozzle remains unchanged except the high pressure values magnify by increasing wind speed ( Figure 11).

External flows over diffuse
CFD simulations around the vertical diffuser are investigated here in external flows. The objective is to examine by changing the upstream air velocity from 1 m/s to 5 m/s what will be the effects on the centreline velocity and the pressure coefficient. Similar to the nozzle problem, the diffuser inlet is located at X/L = 0 and its outlet at X/L = 1. The centreline velocity and pressure coefficient are shown in Figure 12 and Figure 13, respectively. From Figure 12, it is observed that inlet velocity of diffuser increases by a factor of 1.2 to 1.3, for this particular design, by changing wind velocity from 1 m/s to 5 m/s. However unlike Ohya's diffuser, the outlet velocity of the cated at X/L = 0 and its outlet at X/L = 1. From Figure 10, it is observed that the nozzle outlet velocity has slightly increased from the upstream wind velocity by increasing wind speeds from 1 m/s to 5 m/s. This is in agreement with what was observed in experimental works in wind tunnels. However, the most interesting observation here is the pressure coefficient as shown in Figure 11. By increasing wind velocity from 1 m/s to 5 m/s, both inlet and outlet pressure coefficient increase considerably which may suggest that air avoids entering the nozzle due to formation of a high pressure zone at inlet. From CFD simulations, it can be seen that by increasing wind velocity the accelerated zone forms at exterior part of inlet causing less air mass flow rate to pass through the interior part of nozzle. Moreover, high pressure zones forms at inlet and outlet of nozzle pushing more air to pass from exterior part of nozzle due to low pressure  diffuser has not recovered to the free stream velocity. The strong suction at outlet also caused a reverse flow region just downstream of diffuser shown in Figure 12 by negative velocity zone from X/L = 1.26 to X = 1.84. Considering the pressure coefficient as shown in Figure 13, there is considerable decrease of pressure at inlet causes air sucked into the interior of the diffuser. This is interestingly in agreement with experiments of Ohya, et al. [13], although there is vacuum at outlet of diffuser too which increases by increasing wind speed.

Conclusions
In this study, the review of accelerating devices in wind power generation was given. From several studies in wind turbines with shroud diffuser, it is observed that it may be more practical to use diffuser assisted wind collectors. The concept was also proved practical in the literature of wind catchers widely used in architecture of buildings. In this study the following are found: 1. Some of literature suggests that the bare accelerator devices are responsible for wind velocity and wind power augmentation. But, other part of literature considers the interaction between a wind turbine and an accelerator device such as diffuser is the cause for flow augmentation.
2. To clear these points and also to have better understanding of air flows around accelerator devices, some simple mathematical models introduced which compares well with the available experimental data in cases of pipe and nozzle flows; however, there is a mismatch between diffuser outlet data with experiments perhaps due to high viscous flow interaction there.
3. There is not much known on performance of accelerating devices when are used in vertical enclosure wind turbines. To examine the function of different accelerating devices such as nozzles and diffusers in vertical orientation, CFD simulation were performed which indicated that there is no appreciable differences of using these devices in horizontal or vertical arrangements.
4. CFD solutions over nozzle flows shows high pressure augmentation at inlet which may prohibit air for entry to the nozzle and so reduces mass flow rate. The device may be utilized after air collection by a diffuser to accelerate wind speed to desirable speed for operation of wind turbine.
5. CFD simulations of diffuser devices also confirmed the generation of a suction zone which improves considerably intake of air into the diffuser. 6. From these observations, it is understand that to collect and accelerate wind in enclosure wind turbines,