CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Negative drag coefficient...

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

Like Tree2Likes
  • 2 Post By fletc900

Reply
 
LinkBack Thread Tools Display Modes
Old   January 20, 2011, 13:32
Default Negative drag coefficient...
  #1
New Member
 
Bruno F.
Join Date: Feb 2010
Location: Boulder, CO
Posts: 17
Rep Power: 7
fletc900 is on a distinguished road
Greetings.

I'm trying to model a NACA 0006 airfoil to compare it experimental data (from Theory of Wing Sections, I. Abbott & V. Doenhoff). I'm trying to replicate in this case the flow for angle of attack of 4 degrees (so low alpha), and a Reynolds number of 3e6. The flow is incompressible and moves along the +x direction, and I'm running it using simpleFOAM and the Spalart-Allmaras model for turbulence.

The lift coefficient matches the experimantal data (~0.4), but the predicted drag coefficient is negative (i.e, producing thrust), and the solution seems to have converged after 500 iterations. The wierd thing is that I get a positive Cd value (as expected) when the angle of attack is zero, but as soon as I give it an angle of attack, the Cd becomes negative. I must be doing something wrong on the settings. I'm including the entire validation case as *.gz file. Any clues anyone?

Regards,
Bruno
Attached Files
File Type: gz naca0006.tar.gz (5.5 KB, 49 views)
fletc900 is offline   Reply With Quote

Old   January 20, 2011, 23:45
Default
  #2
New Member
 
Bruno F.
Join Date: Feb 2010
Location: Boulder, CO
Posts: 17
Rep Power: 7
fletc900 is on a distinguished road
Solved. Dumb mistake. I had set the angle of attack by specifying <u,v> components of velocity (horizontal and vertical components) while the airfoil was perfectly horizontal. I still had my liftDir and dragDir pointing along the (0,1,0) and (1,0,0) directions which is wrong. This will give you the normal and axial forces on the airfoil, not the lift and drag. Two fixes would be:

1 - Just rotate liftDir and dragDir axial components by whatever angle of attack you are using.

2 - Or use the following equation to calculate it directly:


Cl=Cn cos(alpha) - Ca sin(alpha)
Cd=Ca cos(alpha) + Cn sin(alpha)

Where Ca and Cn are the axial and normal forces that are computed when using dragDir along (1,0,0) and liftDir (0,1,0) with libforces.

Bruno
masood135 and jasonbot like this.
fletc900 is offline   Reply With Quote

Old   April 9, 2011, 06:41
Default
  #3
Member
 
Aamer Shahzad
Join Date: Mar 2010
Posts: 58
Rep Power: 7
aamer is on a distinguished road
Hello Bruno.....

lets say that i am translating an airfoil in x direction and the airfoil has an angle of attack of 40 degree. lets assume that instead of using u and v components of velocity, i have rotated my whole grid by 40 degree (and thus airfoil orients by 40 degree), in such case what should be used for Cl and Cd.... should it be (0,1,0) and (1,0,0) respectively now?
aamer is offline   Reply With Quote

Old   April 9, 2011, 16:38
Default
  #4
New Member
 
Bruno F.
Join Date: Feb 2010
Location: Boulder, CO
Posts: 17
Rep Power: 7
fletc900 is on a distinguished road
Good afternoon Aamer. Yes, if you rotate the mesh by 40 deg, and maintain the flow only along the x-direction, then the lift will be along (0 1 0) and the drag will be along (1 0 0). Just remember, by definition, that the lift is always perpendicular to the flow and the a drag is always in the same direction of the flow. Let me know if this is not clear...

Cheers
Bruno
fletc900 is offline   Reply With Quote

Old   April 10, 2011, 04:51
Default
  #5
Member
 
Aamer Shahzad
Join Date: Mar 2010
Posts: 58
Rep Power: 7
aamer is on a distinguished road
Dear Bruno.....

if the airfoil is moving in still air, then there is no flow around it.... so how can we find out the direction of flow then.....
if airfoil is moving in -x direction (in still air), do we have to assume that the flow direction will be in +x direction....?
aamer is offline   Reply With Quote

Old   April 10, 2011, 11:03
Default
  #6
New Member
 
Bruno F.
Join Date: Feb 2010
Location: Boulder, CO
Posts: 17
Rep Power: 7
fletc900 is on a distinguished road
Good morning Aamer.

It doesn't matter if the airfoil is moving or if it is static with wind blowning at the same velocity upon it. The same relative wind velocity is encountered for the airfoil in both cases. Perhaps this link will help you better understand: http://www.allstar.fiu.edu/aero/flight33.htm


Cheers
Bruno
fletc900 is offline   Reply With Quote

Old   April 10, 2011, 12:10
Default
  #7
Member
 
Aamer Shahzad
Join Date: Mar 2010
Posts: 58
Rep Power: 7
aamer is on a distinguished road
Thanks Bruno.... for the explanation and the link
aamer is offline   Reply With Quote

Old   October 28, 2012, 20:28
Default
  #8
Member
 
R. P.
Join Date: Jul 2010
Location: Brazil
Posts: 63
Rep Power: 6
Rophys is on a distinguished road
Hi everyone,

I am trying to calculate the drag, lift and pitching-moment coefficients for a flow over a probe. The code structure below works fine when the flow is perpendicular to the probe. However, when I changed the angle of attack I got a rude difference for lift and pitching-moment when compared with experimental data (the drag values is OK).

I guess, this difference is related to the vectors (liftDir (0 1 0); dragDir (1 0 0); pitchAxis (0 0 1)). I this this vectors values is just valid for a parallel flow, isn't it ?

Anybody knows how to change this vectors according to a certain angle of attack ? My mesh still the same and I just rotate my geometry in 10 degree.

Thank you very much.

Rophys

forces
{
type forces;
enabled true;
functionObjectLibs ( "libforces.so" );
outputControl outputTime;
patches (probe);
directForceDensity true;
fDName fDMean;
CofR (0.02 0 0);
log on;
}

forceCoeffs
{
type forceCoeffs;
functionObjectLibs ("libforces.so");
patches (probe);
outputControl outputTime;
fDName fDMean;
rhoInf 1.73e-5;
CofR (0.02 0 0);
liftDir (0 1 0);
dragDir (1 0 0);
pitchAxis (0 0 1);
magUInf 1503.1;
lRef 0.05;
Aref 9.817e-4;
}
Rophys is offline   Reply With Quote

Old   January 28, 2015, 13:14
Default
  #9
Senior Member
 
Alhasan's Avatar
 
Hasan KJ
Join Date: Dec 2011
Location: Bristol, United Kingdom
Posts: 172
Rep Power: 5
Alhasan is on a distinguished road
Quote:
Originally Posted by fletc900 View Post
Solved. Dumb mistake. I had set the angle of attack by specifying <u,v> components of velocity (horizontal and vertical components) while the airfoil was perfectly horizontal. I still had my liftDir and dragDir pointing along the (0,1,0) and (1,0,0) directions which is wrong. This will give you the normal and axial forces on the airfoil, not the lift and drag. Two fixes would be:

1 - Just rotate liftDir and dragDir axial components by whatever angle of attack you are using.

2 - Or use the following equation to calculate it directly:


Cl=Cn cos(alpha) - Ca sin(alpha)
Cd=Ca cos(alpha) + Cn sin(alpha)

Where Ca and Cn are the axial and normal forces that are computed when using dragDir along (1,0,0) and liftDir (0,1,0) with libforces.

Bruno
Hello Everyone,
Can someone tell me what are the normal and axial forces he is talking about here the Ca and Cn ??? is it the Cd and Cl computed from when using dragDir along (1,0,0) and liftDir (0,1,0) or something else because when I use the Cd and Cl calculated from computations in the place of Ca and Cn to get the actual Cd and Cl I am getting a list of negative and positive values for a range of angle of attacks from -5 to 20

Edit: Answer: yes the Ca and Cl are the computed lift and drag using dragDir along (1,0,0) and liftDir (0,1,0) and the reason I was getting a mixture of positive and negative values is because I was calculating in Excel and I did not use Radians. you have to use RADIANS when using sin and cos to get the right answer.


Thanks,
Hasan K.J
__________________
"Real knowledge is to know the extent of one's ignorance." - Confucius

Last edited by Alhasan; February 1, 2015 at 09:48.
Alhasan is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
multi wing element - negative coefficient of drag? Zweeper FLUENT 10 March 11, 2010 13:20
Automotive test case vinz OpenFOAM Running, Solving & CFD 98 October 27, 2008 09:43
why my Drag coefficient value is negative value Liu CFX 7 June 28, 2008 07:09
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00


All times are GMT -4. The time now is 11:38.