CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Post-Processing

Force on a patch

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 14, 2012, 01:33
Default Force on a patch
  #1
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
I need to calculate the force caused by the flow on a patch.

In the first picture you can see the result of the simpleFoam solver. In the second there's the patch I want to calculate Fx on.

I tried in Parafoam to integrate pressure on the surface but maybe I miss something (BTW, is there a tutorial regarding this managements ?)

Another simple question: p for incompressible solver is actually divided by density.
That means that the real pressure at the end has to be multiplied by density, right ?

Daniele
Attached Images
File Type: jpg sfoam-p.jpg (9.4 KB, 26 views)
File Type: jpg sfoam-pisto.jpg (7.2 KB, 21 views)
danvica is offline   Reply With Quote

Old   March 14, 2012, 04:04
Default
  #2
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
Probably I should have posted this under the parafoam forum. Is there any way to move it by me ?

Anyway, I tried the following in Parafoam:

1. Select the above patch.
2. Generate surface normals.
3. Use calculator to define a new var called Px defined as
p*N_x*1000. (where 1000 is density of water).
4. Integrate Px on the surface.

Is it correct ? The obtained value is in Newton, right ?
danvica is offline   Reply With Quote

Old   March 15, 2012, 02:42
Default
  #3
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
Another way, more flexible is to use the libforces library as follow:

Code:
functions
{
forcespisto
{
type forces;
functionObjectLibs ("libforces.so"); //Lib to load
patches (pisto); // change to your patch name
pName p;
Uname U;
rhoName rhoInf;
rhoInf 1000;
CofR (0 0 0); //Origin for moment calculations
outputControl timeStep;
outputInterval 5; 
}
}
I'm interested in the Fx value and it seems quite the same, compared with the one obtained by parafoam (see prev post).
Just a confirmation: is it right ?

Daniele
danvica is offline   Reply With Quote

Old   March 15, 2012, 05:54
Default
  #4
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 154
Rep Power: 7
Aurelien Thinat is on a distinguished road
Hi Daniele,

Yes it is. Just modify the CoR (center of reference) if you want the values of Mx My and Mz.
There is also a "forceCoeffs" function for Cx Cy and Cz calculations.

Aurélien
Aurelien Thinat 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
Using starToFoam clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04
StarToFoam error Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38
CheckMeshbs errors ivanyao OpenFOAM Running, Solving & CFD 2 March 11, 2009 03:34
Import gmsh msh to Foam adorean Open Source Meshers: Gmsh, Netgen, CGNS, ... 24 April 27, 2005 08:19
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


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