Calculating Lift and Drag in Paraview (paraFoam)

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

 May 24, 2011, 22:11 Calculating Lift and Drag in Paraview (paraFoam) #1 Member   Scott Join Date: Sep 2009 Posts: 41 Rep Power: 9 Hi there, I have a body that I wish to extract the lift and drag forces for in Paraview. I assume that this can be achieved using the calculator but I cannot work it out. Could someone please give me some pointers on this. Basically I need to calculate the normals for each surface, then find the pressure on each surface in the components of X (or Y or Z, depending on if it is Drag, Lift or side force), then intergrate these values for the whole body to get the Force value. If someone has a procedure for something like this then please let me know! Thanks for your help, Scott

 May 25, 2011, 02:49 #2 Member   Scott Join Date: Sep 2009 Posts: 41 Rep Power: 9 I found a solution: Extract Surfaces Generate Surface Normals Use the calculator to multiply pressure by the surface normals (Normals*p) This then prints the values in the data tab at the right in three columns, one for each normal direction. I then integrate these to get total forces on the area using Integrate Variables. Let me know if you see any problems with this approach! Cheers, Scott

 July 18, 2011, 15:17 #3 New Member   Join Date: Aug 2010 Posts: 15 Rep Power: 8 Hey Scott, I am trying to do the same thing and am attempting your approach. I was just wondering if you've made any more progress since you posted. I'd like to collaborate with you if possible.

 July 18, 2011, 17:33 #4 New Member   Join Date: Aug 2010 Posts: 15 Rep Power: 8 Scott, Where (and how) do you calculate the area of each panel? I'm struggling finding a particularly good way of doing this. Right now I'm performing a similar calculation in Fortran, but I would prefer to do this in paraview.

 September 12, 2011, 10:28 #5 Member   Daniel Join Date: Apr 2010 Location: Manchester Posts: 31 Rep Power: 8 Any advances on this? Also, how would you go about including the effect of skin friction? I guess the wallShearStress post-pro function would be a start.

 January 3, 2014, 04:19 lift and drag coeffs in paraview #6 Member   arash mahboubi Join Date: Jun 2013 Location: Iran Posts: 40 Rep Power: 5 Hi How can I calculate lift and drag coefficient in paraview? I've run an oscillating airfoil with pitching motion.I've calculated lift and drag coefficient through libforces in controlDict, but the values ​​are wrong. I want to calculate them in paraview and plot them over the time or angle. Please tell me how to do this?

 March 27, 2014, 12:38 #7 Senior Member     Louis Gagnon Join Date: Mar 2009 Location: Québec, QC, Canada Posts: 228 Rep Power: 10 In 2D, you might want to try something like this in paraview on the surfaces you want to evaluate, 1. filter: generate surface normals (disable splitting, activate compute cell normals) 2. filter: calculator using something like this for moments : Normals_Y*p*yourDensity*(yourCoRx-16) + Normals_X*p*yourDensity*(yourCoRy-coordsY) replacing yourDensity and yourCoR* by appropriate values or simplify it to get forces 3. filter: integrate variables and look for the results of the integration in the table that appears in the end this procedure gives you sum_cells(P*rho*A*norm_y) I think this can be easily adapted to 3D. cheers Ruli, brucecumt, arashfluid and 1 others like this.

 April 16, 2014, 07:13 #8 Member   Daniel Join Date: Apr 2010 Location: Manchester Posts: 31 Rep Power: 8 I've had some success working with cell data (as opposed to point data) on 3D surface meshes, if it's of use: Cell area can be obtained using the Mesh Quality filter and the Area option. Use the Point Data to Cell Data on your pressure data. Then use Generate Surface Normals (may need to Extract Surface first) with the Compute Cell Nomals option enabled. Pop these into a calculator using the Cell Data mode and that should have you on the right track for the pressure side of things. arashfluid likes this.

 February 12, 2016, 08:58 #9 New Member   Join Date: Jan 2016 Posts: 3 Rep Power: 2 Hello louisgag, could you please explain what do you mean by "yourCoR"? Thank you very much.

 February 12, 2016, 08:59 #10 Senior Member     Louis Gagnon Join Date: Mar 2009 Location: Québec, QC, Canada Posts: 228 Rep Power: 10 Center of Rotation Regards, -Louis

February 12, 2016, 09:00
#11
New Member

Join Date: Jan 2016
Posts: 3
Rep Power: 2
Quote:
 Originally Posted by louisgag replacing yourDensity and yourCoR* by appropriate values or simplify it to get forces cheers
Could you please explain what is "yourCoR"?
Thank you very much.

 February 12, 2016, 09:07 #12 New Member   Join Date: Jan 2016 Posts: 3 Rep Power: 2 Thank you for your quick reply. Another question if you don't mind, So, If my center of rotation is (0,0,0), Does this mean that CoRx=0, and CoRy=0? Sorry if the question is too simple, I am new to CFD.

 February 12, 2016, 09:10 #13 Senior Member     Louis Gagnon Join Date: Mar 2009 Location: Québec, QC, Canada Posts: 228 Rep Power: 10 Yes, that is correct. -Louis

October 24, 2016, 06:50
#14
New Member

Join Date: Nov 2015
Location: Turkey
Posts: 17
Rep Power: 3
Quote:
 Originally Posted by louisgag In 2D, you might want to try something like this in paraview on the surfaces you want to evaluate, 1. filter: generate surface normals (disable splitting, activate compute cell normals) 2. filter: calculator using something like this for moments : Normals_Y*p*yourDensity*(yourCoRx-16) + Normals_X*p*yourDensity*(yourCoRy-coordsY) replacing yourDensity and yourCoR* by appropriate values or simplify it to get forces 3. filter: integrate variables and look for the results of the integration in the table that appears in the end this procedure gives you sum_cells(P*rho*A*norm_y) I think this can be easily adapted to 3D. cheers
Hello and thank you very much. But after doing this I still can't see the drag force (Normal_X*p*1.16)

November 18, 2016, 10:20
#15
New Member

bruce
Join Date: May 2012
Posts: 12
Rep Power: 6
Quote:
 Originally Posted by louisgag In 2D, you might want to try something like this in paraview on the surfaces you want to evaluate, 1. filter: generate surface normals (disable splitting, activate compute cell normals) 2. filter: calculator using something like this for moments : Normals_Y*p*yourDensity*(yourCoRx-16) + Normals_X*p*yourDensity*(yourCoRy-coordsY) replacing yourDensity and yourCoR* by appropriate values or simplify it to get forces 3. filter: integrate variables and look for the results of the integration in the table that appears in the end this procedure gives you sum_cells(P*rho*A*norm_y) I think this can be easily adapted to 3D. cheers
Hi louisgag,
May I ask why it's 16 in your equation: Normals_Y*p*yourDensity*(yourCoRx-16) + Normals_X*p*yourDensity*(yourCoRy-coordsY) ?
Another quenstion is about the center of rotation in paraView. My understanding is that the center of the geometry will be the center of rotation. please correct me if I am wrong. It will be very appreciated if you can indicate how to get or display the coordinates of the center of rotation.
Thank you very much in advance!

 November 18, 2016, 10:42 #16 Senior Member     Louis Gagnon Join Date: Mar 2009 Location: Québec, QC, Canada Posts: 228 Rep Power: 10 Hi Bruce, I'm pretty sure that's a typo, and should read "coordsX" thanks for noting it I'll correct it. Regards, -Louis

 November 18, 2016, 22:34 #17 New Member   bruce Join Date: May 2012 Posts: 12 Rep Power: 6 Hi Louis, Thank you very much for your prompt reply. It's clear now. By the way, I found a way to get the center of rotation: simply click the Probe Location button in the tool bar and the 'Center' followed by three values in the Properties tab would be the Center of Rotation used in the equation. If you cannot find the Probe Location button in your tool bar, go to the main menu, then Fliters>Data Analysis>Probe Location, you will get the same thing. This may be very redundant for experienced ParaView user, but I think it should be helpful for newbie of ParaView, like me. Cheers, Bruce louisgag likes this.

 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 morteza08 FLUENT 4 December 13, 2012 16:34 mahbub03 Main CFD Forum 3 May 19, 2011 00:21 Josh CFX 24 May 9, 2011 10:38 zx Main CFD Forum 4 July 27, 2007 23:38 Jamcam Main CFD Forum 0 October 31, 2002 05:39

All times are GMT -4. The time now is 23:58.