CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

center of pressure

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

Like Tree4Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 3, 2014, 06:56
Question
  #21
Senior Member
 
Anna Tian's Avatar
 
Meimei Wang
Join Date: Jul 2012
Posts: 494
Rep Power: 12
Anna Tian is on a distinguished road
Quote:
Originally Posted by albert_mathews View Post
Although this thread is old, I found it very useful. For those like me looking for more details, I worked on it a bit and found that the following equations work to find the centre of force on an object in a flow.

If I have made any mistakes, please correct me since, in that case, I would like to know how to do it correctly.

In CFX Post, Expressions tab, right click, create a new expression, I chose COFy as name. For a 2D surface called CarBody, the following calculates the height of the drag force applied to the object.
areaInt(y*(Pressure+Wall Shear))@CarBody/areaInt(Pressure+Wall Shear)@CarBody
Likewise for for lift/downforce, I chose COFz.
areaInt(Z*(Pressure+Wall Shear))@CarBody/areaInt(Pressure+Wall Shear)@CarBody

I'm also looking for the way to extract COP location in CFX. I run exactly the same case by Fluent and CFX. The other values obtained (e.g. lift&drag) are almost exactly the same. Fluent can export COP directly. I use your CEL to obtain COP in CFX.

Fluent and CFX give very different value! So I guess your CEL formula shall have problem.

I can't figure out what the problem is. Anyone could help on that?
__________________
Best regards,
Meimei
Anna Tian is offline   Reply With Quote

Old   September 3, 2014, 11:18
Default
  #22
New Member
 
Albert Mathews
Join Date: Jan 2014
Posts: 2
Rep Power: 0
albert_mathews is on a distinguished road
Anna,

thx for looking into this. question, in Fluent, is the COP calculation integrating pressure and viscous forces or just pressure?

regards
albert_mathews is offline   Reply With Quote

Old   September 3, 2014, 19:30
Default
  #23
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,538
Rep Power: 128
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
There is no need to use functions like "areaInt(Pressure+Wall Shear)@CarBody". The function force()@CarBody is the one to use.

Also, the basic derivation of centre of pressure is T=d cross F where T is the torque vector defined around your origin point, F is the force vector and d is the displacement of the force to give the torque. This means your COP is at -d. When you evaluate this it gives you a set of simultaneous equations to solve to get the value of d from T and F. T is evaluated from torque_x/y/z() and F is evaluated from force_x/y/z.

In your case you may be able to simplify this because you know soemthing about the COP. For instance you know it is on the plane in a 2D simulation.
ghorrocks is offline   Reply With Quote

Old   September 3, 2014, 21:41
Question
  #24
Senior Member
 
Anna Tian's Avatar
 
Meimei Wang
Join Date: Jul 2012
Posts: 494
Rep Power: 12
Anna Tian is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
There is no need to use functions like "areaInt(Pressure+Wall Shear)@CarBody". The function force()@CarBody is the one to use.

Also, the basic derivation of centre of pressure is T=d cross F where T is the torque vector defined around your origin point, F is the force vector and d is the displacement of the force to give the torque. This means your COP is at -d. When you evaluate this it gives you a set of simultaneous equations to solve to get the value of d from T and F. T is evaluated from torque_x/y/z() and F is evaluated from force_x/y/z.

In your case you may be able to simplify this because you know soemthing about the COP. For instance you know it is on the plane in a 2D simulation.

Thanks. I'll calculate in this way. Btw, what's wrong with the formula given as below?

areaInt(y*(Pressure+Wall Shear))@CarBody/areaInt(Pressure+Wall Shear)@CarBody
__________________
Best regards,
Meimei
Anna Tian is offline   Reply With Quote

Old   September 3, 2014, 23:44
Default
  #25
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,538
Rep Power: 128
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
areaInt(Pressure+Wall Shear)@CarBody will give you the magnitude of the force vector. I ahev not checked your equation in detail, but for your equation to work you will need to have assumptions about the position on the COP.

But my main problem is why use these functions when a direct function [force()] already exists? It is a little pedantic, but I prefer to keep my stuff simple as possible.
ghorrocks is offline   Reply With Quote

Old   September 4, 2014, 10:20
Default
  #26
Senior Member
 
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 17
singer1812 is on a distinguished road
Glenn:

The method you describe provides a non-unique solution. The force vector acts in a line and you can have many solutions to the set of equations along that line.

That method works well if you are looking for the center of pressure acting on a surface or a plane (or one where you know one of the points). This allows you to fix one of the points and then the other 2 points become unique.

It also can work for an airfoil. You can assume that the center of pressure acts along the cord and then the other 2 points become unique.

Alberts equations are, by definition, the center of pressure.
Batis likes this.
singer1812 is offline   Reply With Quote

Old   September 4, 2014, 10:32
Default
  #27
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,538
Rep Power: 128
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
You are correct, thanks for the correction and clarification.
ghorrocks is offline   Reply With Quote

Reply

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
Pressure Rise Error emueller CFX 0 May 5, 2009 12:08
center of pressure location for car praveenreddy FLUENT 2 December 3, 2007 12:33
variation of center of pressure with angle of atta Virag Mishra CFX 0 August 27, 2007 11:27
Center of Pressure Astrid CFX 0 December 21, 2006 05:07
center of pressure with fluent giorgio FLUENT 0 June 29, 2004 13:36


All times are GMT -4. The time now is 13:57.