CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

CFD Validation of S809 airfoil.

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 5, 2011, 03:39
Default CFD Validation of S809 airfoil.
Senior Member
Raashid Baig
Join Date: Mar 2010
Location: Bangalore, India
Posts: 135
Rep Power: 16
cfd_newbie is on a distinguished road
I want to perform CFD validation of S809 experimental data using commercial CFD software CFD++ (Somers, D. M., 1989, “Design and Experimental Results for the S809 Airfoil,” Airfoils, Inc., State College,PA)
for the following two reasons -
1. CFD results have not be matched correctly till date. Please correct me if I am wrong here. I have read only two papers in this field -
i) Wolfe, W.P. and Ochs, S.S. “CFD Calculations of S809 Aerodynamic Characteristics”, AIAA Paper AIAA-97-0972, 1997
ii) Langtry, R. B., Gola, J., and Menter, F. R., “Predicting 2D Airfoil and 3D Wind Turbine Rotor Performance using a Transition Model for General CFD Codes,” 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, January 2006, AIAA-2006-395

2. Since this is essentially a 2D problem, I can try different turbulence models and higher order methodologies like LNS, DES, DDES (by extruding the 2D geometry) without much computational overhead.
This will be a validation of various options available for turbulence modelling as well.

I need suggestions from fellow users of cfd-online, if I am thinking on the correct lines. Weather I have to take into account some other factors ?

One problem is that while referring to the original document from Somers about the experiment I couldn't find the the values of Pressure, Temperature and Density etc. which I can use as my boundary values at inlet and outlet in my CFD simulation. Can someone help me point out where I can get these quantities ?

Any suggestion from other users who have worked on S809 airfoil earlier will be very helpful.

cfd_newbie is offline   Reply With Quote

Old   July 6, 2011, 02:26
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19
Martin Hegedus is on a distinguished road
In general the input isn't really dependent on the true values of pressure, temperature, etc. as long as you match Mach number and Reynolds number since you'll be comparing to CL and CD. So choose standard sea level pressure and density. Then pick a velocity to get the correct Mach number. Then pick a viscosity to get the correct Reynolds numbers. The true value of temperature becomes more relevant as the Mach number increases because Sutherland's formula is based on (T0+C)/(T+C) where C is Sutherland's constant.

Onto the S809. It's more involved than most seem to be aware of, or at least admit to. And this is based on my own opinion which is backed up by some CFD work I did on it.

The S809 is very thick. Normal airfoil data and intuition do not apply to it.

There are two regions, pre stall and post stall. I believe the first stall starts to occur around 9.5 degrees.

The TU Delft data shows that there is only a small difference between tripped and untripped data. So I'm skeptical about some of the papers I've seen which are playing with turbulence models to try to match the data post stall. Pre stall, fine. Post stall, no. If you want to check out different turbulence models you should only compare to pre stall, if you are modeling it without the WT walls. In my opinion.

When you look at the data, there appears to be TWO lift curve slopes, prestall and post stall. And the airfoil undergoes TWO stalls.

The reports only state the the flow is 2D prestall. One can not assume that the flow is 2D post stall. Again, in my opinion.

So this is what I discovered when I modeled the airfoil with the WT walls. Once the airfoil stalls, a vortex develops at the junction of the wing and wall. Once the vortex develops it knocks down the lift of the airfoil. This is why there is a dip at 12.5 degrees. As the angle of attack increases past this point, the lift increases until the wing stalls a second time.

Unfortunately, the vortex effect is dependent on the thickness of the boundary layer at the wall. You can have many different results by starting the WT boundary layer at different upstream locations. The thinner the boundary layer, the less likely the vortex will be occur. And it is somewhat discreet. It's either there, or not as a function of axial station. What a pain.

You must also include the top and bottom walls for pre stall values. The top and bottom does effect the lift curve slope.

Unfortunately I wasn't able to finish the test case. The 3D analysis with WT walls takes a lot of grid points and I didn't have the resources available.

Some information about the Delft WT: (This came from Nando Timmer, MSc, Wind Energy Research Group, Faculty of Aerospace Engineering, Delft University of Technology)
"1) The low-speed tunnel has a test section of 1.25 (height) x 1.80 (width)x 2.60 m (length) with solid walls.
2) The axial pressure gradient is virtually 0 as a result of a slightly bigger test section exit as entrance, so compensating for the t.s. boundary layer."

Good Luck!
Martin Hegedus is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
CFD Salary CFD Main CFD Forum 17 January 3, 2017 17:09
[ICEM] Having several dificulties with 3d airfoil meshing in ICEM CFD Teja ANSYS Meshing & Geometry 14 February 23, 2015 13:01
CFD Validation John222 CFX 3 May 8, 2011 08:22
Airfoil y+ for Validation Purposes asd FLUENT 3 June 26, 2007 09:27
public CFD Code development Heinz Wilkening Main CFD Forum 38 March 5, 1999 11:44

All times are GMT -4. The time now is 17:52.