# Oscillating airfoil for wind energy

 Register Blogs Members List Search Today's Posts Mark Forums Read

March 6, 2014, 17:51
Oscillating airfoil for wind energy
#1
New Member

Join Date: Feb 2014
Posts: 1
Rep Power: 0
Hi,

I am trying to simulate an airfoil oscillating vertically over a 6m displacement. i want to write UDFs to describe different relationships between the pitch of the airfoil and the vertical displacement. first i am trying a sinusoidal relationship so the foil is at max angle of attack at the centre of the displacement and the there is also a spring force pulling it back to the centre of the range. the foil is a symmetrical NACA0012 of chord length 1m in a constant air speed of 10m/s
i aim to find the best pitch motion that will completes a complete oscillation the quickest.
to do this i am using an interpreted six degrees of freedom UDF shown below:

#include "udf.h"

DEFINE_SDOF_PROPERTIES(fmotion5, prop, dt, time, dtime)
{

real x[ND_ND];

real THETA[2];
real vel[ND_ND];

real a ;

#define PI 3.14159265;
real K=10;
real THETAMAX =0.349;
real SPRING;

SPRING=-K*(x[1]);

prop[SDOF_MASS] = 5;
prop[SDOF_IZZ] = 0.00001;

-3<=x[1]<=3;

prop[SDOF_ZERO_TRANS_X]=TRUE;
prop[SDOF_ZERO_TRANS_Z]=TRUE;
prop[SDOF_ZERO_ROT_X]=TRUE;
prop[SDOF_ZERO_ROT_Y]=TRUE;

a = ((vel[1])/(sqrt(vel[1]*vel[1])))*THETAMAX*cos(asin((x[1])/3));
THETA[2]=a;

}

I can only interpret the UDF due to the restrictions on the university computer i am using.

When i hook the UDF in, the airfoil shows no motion and i sometimes get negative cell volumes. Can anyone help me on this?

I am quite new to CFD by the way.

Henry
Attached Images
 mesh.jpg (96.0 KB, 13 views)

 Tags airfoil2d, oscillating airfoil, sixdof, udf and programming

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post robyTKD SU2 Shape Design 6 June 13, 2013 19:13 robyTKD SU2 Shape Design 21 May 29, 2013 09:26 dancfd OpenFOAM Meshing & Mesh Conversion 3 August 26, 2010 18:52 DarrenC Main CFD Forum 5 July 19, 2010 22:33 Joseph CFX 14 April 20, 2010 15:45

All times are GMT -4. The time now is 02:51.