Editor’s note: this article is an edited version of the case study provided by:
Matheus C. Fernandes,
Case Western Reserve University Dept. of Mechanical Engineering
David H. Matthiesen PhD
Case Western Reserve University, Dept. of Material Sciences and Engineering, Cleveland, Ohio,
david.matthiesen@case.edu
The Case Western Reserve University Wind Turbine, a NorthWind 100kW research turbine, partially powers the adjacent Veale Convocation Center. This turbine was erected in November 2010 and is useful as an example of an urban wind turbine due to its location. It is in the center of the University campus in which the surrounding building heights vary from 20 to 40 meters high. To analyze and visualize the wind flow from every possible direction, a Computational Fluid Dynamics (CFD) model using COMSOL Multiphysics was created. By using local wind direction measurements, the analysis of the model focused mainly on the prevailing wind directions. These wind measurements were captured by cup anemometers, which were located at three different heights (28, 31, 40 meters from ground) and locations. The model includes dimensions of the buildings as well as their footprint location relative to the turbine as shown in Figure 1. However, the model does not include buildings that are not immediately surrounding the turbine, nor does it include the buildings complex shapes. Figure 1 Google aerial view of Case Western Reserve University wind turbine and surrounding buildings.
Governing equations
This model uses the Navier-Stokes iterations and the continuity equation as implemented in the COMSOL Multiphysics laminar flow package. The boundaries of the model are set to at least 60m from the top or the sides of any building in order to avoid any tunnel effect created by the fictitious boundaries. Also, in order to calculate the type of fluid flow, the Reynolds number calculation was used.
Methods
Geometry
To generate the model, a footprint sketch consisting of each buildings location and dimension was drawn using the 2D work plane feature of COMSOL Multiphysics. Each building was then extruded to a bit larger than its dimensional height. The reason for broadening the buildings dimensions in the model is to overshoot any effects the buildings may have on the turbine at hub height. If the buildings with these dimensions do not impact the turbine, then the smaller dimensions will not impact the turbine either.
To generate the model universe, an arbitrary square was also drawn on the 2-D work plane and extruded to a proper height. This square was centered in the model in order to be easily rotated using the rotation angle feature of COMSOL Multiphysics. The rotation feature allows for easy change in the direction of the induced wind flow permitting the analysis on the effect of the buildings on the wind turbine as a function of wind direction. The turbine is constructed by creating a union of regular blocks, cylinders and a vertical extruded sketch to create the shape of the blades. 3.2 Boundary Conditions Once the model universe is created, all walls are set to slip conditions except for the ground which is set to slip condition, one face is set to inflow and the opposing face is set as the outflow. After creating these boundaries the cubes material is set to air at room temperature (27°C). Once all settings are defined for the model universe, each building is cut from the cube using the difference feature of COMSOL Multiphysics. This action recreates the bottom boundary of the model by taking into account the shapes and locations of the buildings. The wind turbine geometry is then also cut using a similar method as to the buildings.
Laminar flow
After the geometry was fully created, the wind was uniformly induced through the face at a steady rate with zero initial values. This model uses the laminar flow feature of COMSOL Multiphysics and therefore uses time invariant stationary solvers.
Numerical model
Using the three years of wind direction and velocity data acquired from the three cup anemometers, the wind was found to prevail from the West North West direction. This wind varies slightly also in a few other directions depending on the season.
Model validation
To validate the model, real wind measurements made by the cup anemometers were compared to nodes created at the same location within the model. Given four specific different wind directions, the model was compared between its nodes and the actual wind data collected by the anemometers.
Wind velocity surface of CWRU turbine area; air flow induced from West (270°prevailing) 15 m/s; note the yellow boundary is below the blades.
General discussion
From the CFD results created by COMSOL Multiphysics, there is a very small effect due to the buildings in most cases, especially when the wind is induced from the North West direction, the prevailing direction.
In the above image, readers can see the colors at hub height that there is no effect on the wind turbine.
Wind velocity isosurface of turbine area; airflow induced at East (90°non-prevailing) 15 m/s; note the yellow slightly going over the blades. The turbine is visible about in the middle of the image.
Laminar flow validation
For this project the COMSOL Multiphysics basic package was used which only solves laminar flow problems. For air at room temperature with velocity of 15 m/s the Reynolds number was calculated to be around 1.1×106. However, a test to validate using laminar flow was conducted by looking at the difference in boundary layer as a function of the Reynolds number. In order to conduct this experiment, the Reynolds number was decreased by increasing the dynamic viscosity of air by keeping all else constant. The dynamic viscosity was increased by an order of magnitude each time resulting in the Reynolds number increasing an order of magnitude.
The COMSOL model was run seven times and the boundary layers were analyzed by using the surface view feature. The difference between boundary layers was minimal. There was a 0.4 meters change in the boundary layer when the Reynolds number was changed from 1.1×106 to 11.0.
Conclusions
- The model had reasonable agreement with the wind measurements.
- The model visually demonstrated the effects of the buildings on the wind profile, thus on the turbine.
- For the prevailing wind direction the wind did not have an effect on the turbine at hub height. From an unlikely wind direction (North East), the turbine was slightly affected by the MCCO smoke stack.
Comsol Multiphysics
Filed Under: News
