# What are the points in simulating 2D airfoil??

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

January 28, 2013, 05:08
What are the points in simulating 2D airfoil??
#1
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Hi, guys,

I would like to calculate the aerodynamic dimensionless factors, Cl, Cd, Cm of the airfoils under any situation since I have to develop an unmanned helicopter.

So, for the first step, I have been trying to simulate the 2D airfoil of NACA 0012 with Fluent so far.

But, unfortunately I couldn't derive the reasonable values of the factors, yet. I compared them with the chart of NACA0012 characteristics on 'Helicopter Performance, Stability and Control by Raymond W. Prouty' which I attached here.

The conditions of the first case was 4 degree of the angle of attack with Mach 0.8 using the solvers, SA, RSM and k and epsilon.
But, all of them were failed.
(I could not upload the the case & data file I tried due to the file size.)

Could you give me some advises what's the point to set up the boundary conditions and/or the modeling for simulating 2D airfoil?

Especially, I'm confused with the relation of the near-wall function with the mesh near the airfoil surface.
According to ANSYS help, the mesh should be fine without the wall functions, while it may be coarse if I choose one of the wall functions.
So, do I have to make the mesh coarse if I use RSM or k and epsilon model? Also what is the suitable functions for this case among the wall functions?

Thanks,
Have a good day!

ps) If you have additional tips for the rotor blade which is generally considered its own rotation, could you please let me know them?
Attached Images
 Characteristics of NACA0012 from Prouty book.jpg (76.5 KB, 53 views)

 January 29, 2013, 11:48 #2 Member   Charles Join Date: May 2011 Location: France Posts: 77 Rep Power: 7 Hi, please post a picture of your mesh and the value of Y+ that you get. For a RANS simulation please remember that y+ shall not exceed 300. If you are willing to get precise numerical results I recommande a y+ around 10. Mohd123 likes this.

January 29, 2013, 13:44
#3
New Member

Raj
Join Date: Feb 2011
Posts: 1
Rep Power: 0
Quote:
 Originally Posted by cocobi Hi, guys, I would like to calculate the aerodynamic dimensionless factors, Cl, Cd, Cm of the airfoils under any situation since I have to develop an unmanned helicopter. So, for the first step, I have been trying to simulate the 2D airfoil of NACA 0012 with Fluent so far. But, unfortunately I couldn't derive the reasonable values of the factors, yet. I compared them with the chart of NACA0012 characteristics on 'Helicopter Performance, Stability and Control by Raymond W. Prouty' which I attached here. The conditions of the first case was 4 degree of the angle of attack with Mach 0.8 using the solvers, SA, RSM and k and epsilon. But, all of them were failed. (I could not upload the the case & data file I tried due to the file size.) Could you give me some advises what's the point to set up the boundary conditions and/or the modeling for simulating 2D airfoil? Especially, I'm confused with the relation of the near-wall function with the mesh near the airfoil surface. According to ANSYS help, the mesh should be fine without the wall functions, while it may be coarse if I choose one of the wall functions. So, do I have to make the mesh coarse if I use RSM or k and epsilon model? Also what is the suitable functions for this case among the wall functions? Thanks, Have a good day! ps) If you have additional tips for the rotor blade which is generally considered its own rotation, could you please let me know them?

Generally you have to get closer results for 2d airfoil cases... can you check which boundaries you are using to integrate pressure/skin friction to get total cl, cd, cm...you have to use only airfoil boundary for intgration, not the farfield boundary... are you doing proper axes transformations when converting from cx, cy, cmy to cl, cd, cm

January 29, 2013, 21:23
#4
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Quote:
 Originally Posted by Santos-Dumont Hi, please post a picture of your mesh and the value of Y+ that you get. For a RANS simulation please remember that y+ shall not exceed 300. If you are willing to get precise numerical results I recommande a y+ around 10.
Hi, Santos

The following list is the files I enclosed here.

1. Mesh-wide view
2. Mesh-close view
3. General Check Domain: Check and Report Quality message for the mesh in General.
4. Yplus-SST-4D AoA-Ma08-PFF: Yplus graph, Transition SST solver, Angle of Attack 4 Degree, Mach number 0.8, Boundary Condition-Pressure Far Field
5. Yplus-SST-4D AoA-Ma08-Vel Inlet-Pressure Outlet: Same as the above except the Boundary Condition-Velocity Inlet and Pressure Outlet

Pressure Far Field (PFF) model sets the edge of the domain as the PFF, and Velocity Inlet and Pressure Outlet one sets the front edge of the domain as Velocity Inlet and the rear one as Pressure Outlet.
And I tried several solvers (SA, k-epsilon, k-omega, SST, RSM), but the trend of their result showed similar.

Could you please check if the mesh consition is OK or not?
According to the ANSYS help, y+ should be about 1 or greater than 30. As you can see, both of the Yplus graphs exceed 30 and below 300 that you mentioned. But, there is an obvious difference between the two graphs. The pressure far field model seems to show something like shock or separation after 0.4m on position axis while the vel. inlet model has clear line.
Do you know what makes this difference?

In addition, please let me know how to adjust the Yplus to near 10. As you can see, the value has the order of 2. Do I have to make the mesh more fine or coarse?

Thank you,
For more information, I will attach more images on the next message.

Have a good day,
Attached Images
 Mesh-close view.jpg (53.9 KB, 53 views) Mesh-wide view.jpg (90.0 KB, 48 views) Yplus-SST-4D AoA-Ma08-PFF.jpg (49.6 KB, 45 views) Yplus-SST-4D AoA-Ma08-Vel Inlet-Pressure Outlet.jpg (44.5 KB, 42 views)
Attached Files
 General-Check domain.txt (1.3 KB, 12 views)

Last edited by cocobi; January 29, 2013 at 21:51.

January 29, 2013, 21:42
#5
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
I attached the below list additionally.

1. Ystar-SST-4D AoA-Ma08-PFF (pressure far field)
2. Ystar-SST-4D AoA-Ma08-Vel Inlet-Pressure Outlet
3. Result-Cd-SST-4D AoA-Ma08-PFF: Drag coefficient for PFF model.
4. Result-Cd-SST-4D AoA-Ma08-Vel Inlet-Pressure Outlet: Drag coefficient for Vel Inlet and Pressure Outlet.
5. Mesh nose and tail.

Could you let me know what the difference is between the Yplus and Ystar?

With respect to the Drag coefficients, the result of Pressure Far Field model shows a periodic trend, while the other one shows the converged value(around 0.008) but it is away from the experimental result (0.05).
Like Yplus, the two models have completely different result, even though only different thing is the setting of the edge line of the domain each other.
Could you explain why?

Thank you.
Attached Images
 Ystar-SST-4D AoA-Ma08-PFF.jpg (50.1 KB, 16 views) Ystar-SST-4D AoA-Ma08-Vel Inlet-Pressure Outlet.jpg (44.9 KB, 13 views) Result-Cd-SST-4D AoA-Ma08-PFF.jpg (36.6 KB, 17 views) Result-Cd-SST-4D AoA-Ma08-Vel Inlet-Pressure Outlet.jpg (33.5 KB, 11 views) Mesh-nose and tail.jpg (67.4 KB, 23 views)

January 29, 2013, 22:10
#6
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Quote:
 Originally Posted by cfdraj Generally you have to get closer results for 2d airfoil cases... can you check which boundaries you are using to integrate pressure/skin friction to get total cl, cd, cm...you have to use only airfoil boundary for intgration, not the farfield boundary... are you doing proper axes transformations when converting from cx, cy, cmy to cl, cd, cm
Hi, cfdraj

I tried two different cases in terms of the boundary condition setting.

Case 1: Pressure Far Field set up. In the picture of the Boundary Condition attached, Edge A and B are set as Pressure Far Field and the airfoil surface as wall.

Case 2: Velocity inlet and Pressure Outlet set up. In the picture, Edge A is established as Velocity Inlet and Edge B as Pressure Outlet. The airfoil surface is same as Case 1.

But, the results of them are clearly different as I explained to Santos above.

Also, I enclosed the BC set up dialog windows and monitor set up window of Drag, Lift, Moment.

Thanks,
Have a good day.
Attached Images
 Boundary Condition.jpg (80.3 KB, 19 views) BC-Pressure Far Field set up.jpg (65.5 KB, 31 views) BC-Velocity Inlet set up.jpg (84.0 KB, 26 views) Monitor-drag lift moment.jpg (84.3 KB, 22 views)

January 30, 2013, 05:05
#8
Senior Member

Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 19
Quote:
 Originally Posted by Santos-Dumont If you are willing to get precise numerical results I recommande a y+ around 10.
Why would you do that? As I understand it, y+=10 is one of the worst choices at all, since it is much too high to resolve the viscous layer but also too low to use a wall function. It is right there, where you don't want to have your first grid point.
__________________
The skeleton ran out of shampoo in the shower.

January 30, 2013, 09:31
#9
Member

Charles
Join Date: May 2011
Location: France
Posts: 77
Rep Power: 7
Quote:
 Originally Posted by RodriguezFatz Why would you do that? As I understand it, y+=10 is one of the worst choices at all, since it is much too high to resolve the viscous layer but also too low to use a wall function. It is right there, where you don't want to have your first grid point.
Yes my mistake I agree with you. That is comprehensivly explained in Fluent help, although I never quite like the fact of voluntarely degrading the boundary layer density from 10 to 30 for example. But it seems that it might increase the results accuracy so...

With a regular computer it's hard to reach y+ around 1 so 30 is probably the best choice.

 January 30, 2013, 09:37 #10 Senior Member     Philipp Join Date: Jun 2011 Location: Germany Posts: 1,297 Rep Power: 19 The turbulence (wall) model has two choices: a) calculate the near wall behavior itself, which needs an y+ of about 1 b) use the log-law function, which is valid for y+ > 30 Both are violated for y+=10, so the accuracy generally gets better when you increase (or decrease) a y+ of 10. __________________ The skeleton ran out of shampoo in the shower.

January 30, 2013, 20:01
#11
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Quote:
Hi, Santos,

I am really appreciating for your help.
I am changing the mesh right now, and be careful to set the boundary conditions following your explanation.

I am really clear thanks to your kind explanation.

Just one more question of the sonic speed.
You explained the speed under Mach 0.3 can neglect its compressiblity.
Does that mean I have to consider compressible flow for over Mach 0.3, unconditionally, or because the flow can reach to mach 0.7 which is transonic during travelling the airfoil surface?

Generally speaking, the mach of transonic speed ranges from 0.7 to 1.2 as far as I know, so I thought the flow under 0.7 could neglect its compressibilty.

Of course, for the previous case, it is completely the lack of my knowledge.. ^^;;

Regarding the website, both explorer and chrome cannot connect to 'supercoolprops.com' website. I hope it would be just a IP traffic problem at this moment and then be fine later.

Thank you again.
Have a nice day!

Last edited by cocobi; January 30, 2013 at 20:31.

January 30, 2013, 20:13
#12
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Quote:
 Originally Posted by RodriguezFatz The turbulence (wall) model has two choices: a) calculate the near wall behavior itself, which needs an y+ of about 1 b) use the log-law function, which is valid for y+ > 30 Both are violated for y+=10, so the accuracy generally gets better when you increase (or decrease) a y+ of 10.
Hi, RodriguezFatz,

It is very helpful for me as well as Santos' one.

Just one question I have here.
Do you mean the log-law function stands for a wall function?
If yes, I will be completely clear, but if not.... I would be afraid to be confused..

In any case, you showed me a clear standard when I choose turbulence model. Thank you for that.

Have a good day!

 January 31, 2013, 02:39 #13 Senior Member     Philipp Join Date: Jun 2011 Location: Germany Posts: 1,297 Rep Power: 19 Yes. Maybe there are other wall functions... but log-law is definitly the "default". cocobi likes this. __________________ The skeleton ran out of shampoo in the shower.

January 31, 2013, 04:46
#14
Member

Charles
Join Date: May 2011
Location: France
Posts: 77
Rep Power: 7
Quote:
 Originally Posted by cocobi Just one more question of the sonic speed. You explained the speed under Mach 0.3 can neglect its compressiblity. Does that mean I have to consider compressible flow for over Mach 0.3, unconditionally, or because the flow can reach to mach 0.7 which is transonic during travelling the airfoil surface? Generally speaking, the mach of transonic speed ranges from 0.7 to 1.2 as far as I know, so I thought the flow under 0.7 could neglect its compressibilty.
About the picture I send, you can look here

Page 5 and 6 you have nice wind tunnel pictures showing you the phenomenom.

About your theorical question. Yes generally speaking a flow is transsonic between 0.7 and 1.2. In fact the definition of a transsonic flow is a flow where you have a clear supersonic and subsonic zone. So it indeed depends on the geometry. For example modern aircrafts have a very high transsonic limit (the goal is indeed to reduce the supersonic points to save drag and fuel). The basis of that is the area rule for the fuselage, discovered by the German manufacturer Junkers during the WWII. Using low thickness airfoil also helps because it reduces the speed increase on the upper surface ofthe airfoil.

Anyway, compressibility and transsonic flow are different. The air is always compressible (otherwise aircrafts wouldn't have speed indication under M0.3 and I would be in trouble with my ultralight !) and you get dynamic pressure even at a very low speed.
But in general the physics is simplified under M0.3 and the air is supposed as incompressible (Bernouilli law).

You can use a incompressible setup to slightly higher speed if you want, but if you get close to transsonic you must choose compressible (St-Venant law) otherwise you'll miss the important shock phenomenom.

About the turbulence model, my teachers gave me a simple answer : k-Epsilon for internal flow (turbine, tubes etc...) and Spalart-Allmaras for external flow (aircraft, airfoils...). I like when it's simple !!

January 31, 2013, 20:53
#15
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Quote:
 Originally Posted by RodriguezFatz Yes. Maybe there are other wall functions... but log-law is definitly the "default".
Thank you, RodriquezFatz

Have a good day!

January 31, 2013, 20:55
#16
New Member

Derek Kim
Join Date: Jan 2013
Posts: 15
Rep Power: 6
Quote:
 Originally Posted by Santos-Dumont About the picture I send, you can look here https://www.google.fr/search?q=trans...A8qi0QXkjYCgDA Page 5 and 6 you have nice wind tunnel pictures showing you the phenomenom. About your theorical question. Yes generally speaking a flow is transsonic between 0.7 and 1.2. In fact the definition of a transsonic flow is a flow where you have a clear supersonic and subsonic zone. So it indeed depends on the geometry. For example modern aircrafts have a very high transsonic limit (the goal is indeed to reduce the supersonic points to save drag and fuel). The basis of that is the area rule for the fuselage, discovered by the German manufacturer Junkers during the WWII. Using low thickness airfoil also helps because it reduces the speed increase on the upper surface ofthe airfoil. Anyway, compressibility and transsonic flow are different. The air is always compressible (otherwise aircrafts wouldn't have speed indication under M0.3 and I would be in trouble with my ultralight !) and you get dynamic pressure even at a very low speed. But in general the physics is simplified under M0.3 and the air is supposed as incompressible (Bernouilli law). You can use a incompressible setup to slightly higher speed if you want, but if you get close to transsonic you must choose compressible (St-Venant law) otherwise you'll miss the important shock phenomenom. About the turbulence model, my teachers gave me a simple answer : k-Epsilon for internal flow (turbine, tubes etc...) and Spalart-Allmaras for external flow (aircraft, airfoils...). I like when it's simple !!
Thank you, Santos.

Have a good day!

 February 5, 2013, 14:28 #17 New Member   Join Date: Jan 2013 Posts: 1 Rep Power: 0 Guys, which reference values for area, lenght and depth do you use in the 2D airfoil analysis ?

January 8, 2015, 16:13
#19
Member

Charles
Join Date: May 2011
Location: France
Posts: 77
Rep Power: 7
Quote:
 Originally Posted by sajad 1- in the pressure far field boundary condition you should enter Mach number, is this related to blade tip Mach number and i should enter M=0.8 or this is far field Mach number and i should enter M=0 (since my case is in hover)?
M= 0.8 for the FarField condition, of course

Quote:
 Originally Posted by sajad 2- I don't need flow field far from the blade so i just want to model near blade region but as it's said in fluent user guide when you use pressure far field BC your region radius must be 20 cord length! what can i do about this? if i want to reduce my solution domain(to save computation time) is there any other BC i can use instead of Pressure far field ?(note that flow is compressible in my case)
There is a reason for the Fluent guide to propose the 20 chords empiric rule. The foil will create a pressure variation that will spread several chords away from your airfoil so for your boundary condition you must have the infinite flow without any perturbation caused by the airfoil. You can simulate in a domain that is too small but this will create either divergence or possible false results.
To save computational time, you may use a coarse mesh close from the boundary and a refined mesh near your foil. If you separate and count you'll see that the number of nodes that you need to reach the 20 chord distance between the airfoil and the BC is nothing in comparison with the number of nodes that is needed to simulate only the first millimeter of the boundary layer...

Charles

January 10, 2015, 03:10
#20
New Member

Join Date: Aug 2010
Posts: 17
Rep Power: 8
Quote:
 Originally Posted by Santos-Dumont M= 0.8 for the FarField condition, of course There is a reason for the Fluent guide to propose the 20 chords empiric rule. The foil will create a pressure variation that will spread several chords away from your airfoil so for your boundary condition you must have the infinite flow without any perturbation caused by the airfoil. You can simulate in a domain that is too small but this will create either divergence or possible false results. To save computational time, you may use a coarse mesh close from the boundary and a refined mesh near your foil. If you separate and count you'll see that the number of nodes that you need to reach the 20 chord distance between the airfoil and the BC is nothing in comparison with the number of nodes that is needed to simulate only the first millimeter of the boundary layer... Charles

Hi, Santos

i just don't understand why we should use M=0.8 for far field Mach number. far from the rotor, air is in stagnation so it's velocity and Mach number is zero. where is my mistake?

 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 bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 6 August 26, 2015 09:38 andrew_diablo CFX 2 November 12, 2012 06:28 ANT ANSYS Meshing & Geometry 2 July 23, 2012 15:25 quarkz Main CFD Forum 0 November 12, 2010 19:44 hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36

All times are GMT -4. The time now is 21:00.