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

Finding Lift/Drag/Moment Coefficients on 3D model aircraft

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 25, 2012, 14:29
Default Finding Lift/Drag/Moment Coefficients on 3D model aircraft
  #1
New Member
 
Eric
Join Date: Jan 2012
Posts: 18
Rep Power: 5
EVBUCF is on a distinguished road
This may be a simple question but it is tripping me up. I am solving for the coefficient of lift/drag/moment on a 3D model of a Boeing x-48B and was curious which reference area and reference length to use? I understand the lift direction to be the direction perpendicular to the aircraft's top surface and the drag direction straight into the nose of the aircraft.

For the reference area, am I to use only the area of one wing or both? Since the X-48B is a blended wing body, if I select the area of both wings then that will essentially be the entire aircraft body area.

Also, do I only find the area of the bottom surface of the wing, or the entire surface area top and bottom?

Is this reference area the same for drag as lift?

As far as the reference length to be used for the moment coefficient where do I measure from? The X-48B has swept wings which are blended into the body of the aircraft so this confuses me as well. I attached an image where I labelled 3 possible lengths. Which one should I use? (A,B, or C)

I am using the following function:

functions
{
forces
{
type
forceCoeffs;
functionObjectLibs ( "libforces.so" );
outputControl timeStep;
outputInterval 1;
patches
(
plane
);
pName p;
UName U;
rhoName rhoInf;
log true;
rhoInf 1;
CofR ( 0 0 0 );
liftDir ( 0 1 0 );
dragDir ( 0 0 -1 );
pitchAxis ( 1 0 0 );
magUInf 10;
lRef 1;
Aref 1;
}
}
Attached Images
File Type: png Screenshot.png (55.0 KB, 67 views)
EVBUCF is offline   Reply With Quote

Old   April 25, 2012, 18:02
Default
  #2
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 463
Rep Power: 9
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
You can use whatever you want. For instance you could you your medium aerodynamic chord times your span (c in you figure) but it won't make any different… You can use half of it or you can use 1, the drag and the lift will still be the same, only your coefficients will change according to what you choose as a reference area…
lovecraft22 is offline   Reply With Quote

Old   April 25, 2012, 18:58
Default
  #3
New Member
 
Eric
Join Date: Jan 2012
Posts: 18
Rep Power: 5
EVBUCF is on a distinguished road
So I could just use a reference area and length of 1 to make it simple and then when I solve for Lift it will really be Lift per unit area? So to get the total lift on the aircraft I would have to multiply by the surface area of the aircraft?

Also I want to be sure I am correct with my thinking here. Is the lift direction perpendicular to the wing surface or is it always the vertical axis regardless of the angle of attack? I am doing tests at multiple angles of attack so I want to be certain I have these directions sorted out.
EVBUCF is offline   Reply With Quote

Old   April 26, 2012, 02:54
Default
  #4
Member
 
Join Date: Apr 2012
Location: France
Posts: 72
Rep Power: 5
Rider is on a distinguished road
Quote:
Originally Posted by EVBUCF View Post
So I could just use a reference area and length of 1 to make it simple and then when I solve for Lift it will really be Lift per unit area? So to get the total lift on the aircraft I would have to multiply by the surface area of the aircraft?
Like lovecraft22 has said, "You can use whatever you want".

Drag force => D=0.5*rho*Sref*V^2*Cd.

So,
If you have Cd=0.2 with Sref=2 (m^2) or if you have Cd=0.4 with Sref=1 (m^2), then you will have the same results.

For exemple :
rho=1.225 (kg.m^-3), V=10 (m.s^-1), Cd=0.2, Sref=2 (m^2)
====> D=24.5 N

rho=1.225 (kg.m^-3), V=10 (m.s^-1), Cd=0.4, Sref=1 (m^2)
====> D=24.5 N


Quote:
Originally Posted by EVBUCF View Post
Also I want to be sure I am correct with my thinking here. Is the lift direction perpendicular to the wing surface or is it always the vertical axis regardless of the angle of attack? I am doing tests at multiple angles of attack so I want to be certain I have these directions sorted out.
The lift direction depends on the angle of attack. To be precise, they are in the aerodynamic marker.
Rider is offline   Reply With Quote

Old   April 26, 2012, 04:56
Default
  #5
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 463
Rep Power: 9
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Actually the lift is defined as the aerodynamic force acting perpendicular to the flow direction and the drag is the force that has the same direction as the flow. So, if you rotate your plane around the pitch axis the lift and drag directions won't change.
lovecraft22 is offline   Reply With Quote

Old   April 26, 2012, 05:21
Default
  #6
Member
 
Join Date: Apr 2012
Location: France
Posts: 72
Rep Power: 5
Rider is on a distinguished road
I think we say the same thing.

We have two cases :
If you change the components of the velocity (and don't rotate the geometry), the lift and drag direction must change.
If you rotate the geometry, the lift and drag direction don't change.

With his question, I think that he wants to change the angle of attack like the first case.
Rider is offline   Reply With Quote

Old   April 26, 2012, 06:03
Default
  #7
Senior Member
 
lore
Join Date: Mar 2010
Location: Italy
Posts: 463
Rep Power: 9
lovecraft22 is on a distinguished road
Send a message via Skype™ to lovecraft22
Yes, I thought of that actually… You are right, we are saying the same thing, that's why I said:

Quote:
So, if you rotate your plane around the pitch axis the lift and drag directions won't change
and to summarize they will change if you change your flow direction.
lovecraft22 is offline   Reply With Quote

Old   April 27, 2012, 10:40
Default
  #8
New Member
 
Eric
Join Date: Jan 2012
Posts: 18
Rep Power: 5
EVBUCF is on a distinguished road
Thanks so much for your help guys! You all cleared up any doubts I had. Thanks again!
EVBUCF is offline   Reply With Quote

Old   August 11, 2012, 05:33
Default Re
  #9
Member
 
nauman hashmi
Join Date: Jun 2012
Posts: 54
Rep Power: 5
nauman55 is on a distinguished road
sir can i have the 3d model of x48B i want to study the aerodynamic loads on it. if you can send me the 3d model at naumaneshanhashmi@gmail.com. i would extremely thankful to you.
regards
nauman55 is offline   Reply With Quote

Old   July 20, 2013, 05:09
Default
  #10
Member
 
Join Date: May 2010
Posts: 55
Rep Power: 7
hewei is on a distinguished road
Quote:
Originally Posted by lovecraft22 View Post
Actually the lift is defined as the aerodynamic force acting perpendicular to the flow direction and the drag is the force that has the same direction as the flow. So, if you rotate your plane around the pitch axis the lift and drag directions won't change.
HI is there some way to calculate the multi-element wing force coeffecient?
I write the function like this, but the openfoam just output the flap's coefficient not the total wing. so what should I do? thank you!
functions

{ forces

{

type forceCoeffs;

functionObjectLibs ( "libforces.so" );

outputControl timeStep;

outputInterval 1;

patches

(
wing_main
wing_flap
);


rhoName rhoInf;

log true;

rhoInf 1;

CofR ( 0 0 0 );

liftDir (0 1 0 );

dragDir ( 1 0 0 );

pitchAxis ( 0 0 1 );

magUInf 10;

lRef 1.0;

Aref 1.0;
}
hewei is offline   Reply With Quote

Old   July 20, 2013, 17:01
Default
  #11
New Member
 
reza sadeghi
Join Date: May 2013
Posts: 16
Rep Power: 4
reza1111 is on a distinguished road
Hi,
I'm not sure. Instead of:

Code:
wing_main
wing_flap
try this:

Code:
"wing.*"
Hope this work.

Reza
reza1111 is offline   Reply With Quote

Old   July 27, 2013, 14:03
Default
  #12
Member
 
Join Date: May 2010
Posts: 55
Rep Power: 7
hewei is on a distinguished road
Thank you Reza, I tried ,but failed, there must be sth wrong what I did, anyway thanks for your help.

Wei
Quote:
Originally Posted by reza1111 View Post
Hi,
I'm not sure. Instead of:

Code:
wing_main
wing_flap
try this:

Code:
"wing.*"
Hope this work.

Reza
hewei is offline   Reply With Quote

Old   August 7, 2013, 16:22
Default
  #13
New Member
 
reza sadeghi
Join Date: May 2013
Posts: 16
Rep Power: 4
reza1111 is on a distinguished road
Hi wei,
Try this:
Code:
 

"wing_*"
or
"wing*"
hope this works

the point is that where there are some walls with different names and the part of their names are same, we can refere to the part which is same and * symbole will include every words of the rest.
reza1111 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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
New implemented algebraic Reynolds stress model volker OpenFOAM 13 October 21, 2013 00:01
CFD of Model Aircraft Sashankh FLUENT 9 December 8, 2008 06:53
CFD Analysis of Model Aircraft Sashankh FLUENT 0 December 1, 2008 23:49
MAKING A CFD MODEL OF AN AIRCRAFT, HELP WANTED kjvsingh Main CFD Forum 1 June 30, 2002 09:58


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