CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Problems with VAWT simulation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 22, 2014, 12:35
Default Problems with VAWT simulation
  #1
New Member
 
Fernando Ruviaro
Join Date: Mar 2014
Posts: 9
Rep Power: 12
Fernando R. is on a distinguished road
Hi. I am trying to simulate the performance of a vertical axis wind turbine using the CFX. The geometry of the rotor, its rotational speed and the wind conditions are based on the parameters utilized in experimental tests on field.
Since the simulation was based on those experimental tests, I was supposed to obtain similar results. The problem is that, while the real turbine provided more than 10kW of power, the one that I've simulated only consumed power to rotate (generated negative torque). I would like to know what I'm doing wrong.

Information about the domain geometry, mesh and boundary conditions are summarized bellow. Images are attached too.

The rotor has straight blades. To perform the simulation, I've used only half of the rotor (cutting with a plane normal to its axis), being imposed the symmetry boundary condition on the plane of the cut. The simulation was transient.

DOMAIN GEOMETRY
As seen on the Figures 1 to 3, I've created two domains: one enclosing the region of the rotor and a bigger one for the surroundings. The inlet would be on the semicylindrical surface of the surroundings domain. I've chosen it to be semicylindrical to reduce a bit the number of elements.

MESH
I've used a tetragonal, path independent mesh in both domains. In the rotor's I've used a "body of influence" sizing in the region of the blades. Besides, I've applied inflation with the blades surface as boundaries. Figures 4 and 5 show a bit of the mesh.

BOUNDARY CONDITIONS
Surroundings domain
As said above, the "Inlet" was the semicylindrical surface of the surroundings domain. The speed was given by cartesian coordinates. The opposed surface was set as "Outlet", with 0 Pa as average static pressure. On the sides and on the top the surfaces were set as "Free Slip Wall". The bottom surface was set as "Symmetry".
Rotor Domain
The bottom of the rotor domain was set as "Symmetry", and the surface of the rotor (blades, axis, supporting arms) were set as "Smooth Wall". The domain was set as "Rotating" with constant speed.
Interface
The interface between the rotor and surroundings domains (a circle and a cylinder wall) was set as "Fluid Fluid" type. The model was "General Connection" with "Transient Rotor Stator" Frame Change/Mixing Model. The Pitch Change was set as "None".
Attached Images
File Type: jpg Figure_1.jpg (22.6 KB, 20 views)
File Type: jpg Figure_2.jpg (14.2 KB, 16 views)
File Type: jpg Figure_3.jpg (26.6 KB, 21 views)
File Type: gif Figure_4.gif (82.6 KB, 27 views)
File Type: gif Figure_5.gif (73.4 KB, 18 views)
Fernando R. is offline   Reply With Quote

Old   September 22, 2014, 15:53
Default
  #2
Senior Member
 
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 20
singer1812 is on a distinguished road
These are general statements based on your pics alone:

Domain seems small front to back. Mesh seems like it needs work.

Are you using ICEM to get the mesh?
Size transitions seem abrupt, which is what you get using Octree.
Inflation layer doesnt seem adequate.

Try generating volume mesh using Octree using an appropriate expansion, blow away the volume mesh, and generate a delauney with an appropriate expansion. You will get a much smoother transition. Generate the prisms on the delauney.

This geometery also looks pretty simple. You should be able to get a decent hex mesh on it. You can try that also...
singer1812 is offline   Reply With Quote

Old   September 23, 2014, 13:10
Default
  #3
New Member
 
Fernando Ruviaro
Join Date: Mar 2014
Posts: 9
Rep Power: 12
Fernando R. is on a distinguished road
I've attached an image of the half rotor geometry, if it's of interest. The blades has the NACA 0018 profile, with the tips made with a rotation of the profile. The arms are elliptical and the axis is a simple cylinder.
For the mesh I've used the default program that comes when I select the CFX on the Workbench. I've used the path independent method because I thought doing so would reduce the number of elements and reduce the time required to create the mesh.
About the domain, how much would the height, length and width be, relatively to the turbine diameter and half it's height (since I'm simulating half rotor)?
About the suggestions you gave me about creating the mesh, could you tell briefly how can I do it in the default meshing program associated with the CFX?
Attached Images
File Type: jpg Rotor.jpg (30.1 KB, 15 views)
Fernando R. is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Having problems running transient fluent simulation in batch mode obylong FLUENT 1 August 14, 2014 23:58
2D vawt simulation meshing and fluent error daysley FLUENT 5 April 15, 2014 19:51
fluent VAWT simulation ahmed_yousry Main CFD Forum 0 January 26, 2013 06:25
Problems with visualization of simulation in Linux cicagol STAR-CCM+ 3 September 23, 2009 05:51
Fire Simulation Problems alkidos FLUENT 1 August 25, 2009 19:26


All times are GMT -4. The time now is 10:54.