
[Sponsors] 
November 22, 2013, 09:03 
Sectional Drag and lift coefficient?

#1 
Senior Member
Join Date: Dec 2010
Posts: 135
Rep Power: 7 
Hi guys,
I know how to get an overall drag and lift over whole domain, What I need is also a sectional drag& lift output, in example: a 3D Cylinder is 5 heigth units high. OF gives me CD for the surface "cyl". And what I want to achieve is also a drag and lift for Z=1,2,3 Any idea? 

November 22, 2013, 09:48 

#2 
Senior Member
Join Date: Jan 2013
Posts: 111
Rep Power: 5 
I have an idea yet have not put into practice.
The drag and lift forces are calculated from pressure and friction forces around the cylinder, one can easily get the pressure distribution around the cylinder at desire Z by utilizing the sample utility. As for the friction force, with my immature knowledge about fluid dynamics I think one can calculate from the wall shear force since the utility wallGradU is available in OpenFoam. These are just some rough thoughts, I'm not sure if it could work. 

November 22, 2013, 09:53 

#3 
Senior Member
Join Date: Jan 2013
Posts: 111
Rep Power: 5 
by the way I'm also doing LES of flow around cylinders at Re=1000 by pimpleFoam, however, as hard as I try, I could not obtain a satisfactory result in terms of drag coefficient. My LES model is smagorinsky with Cs=0.1 and y+ near the cylinder wall is smaller than 1. The Cd in my simulation is about 1.15 but from references it should be around 1.05.
Any advice? 

November 25, 2013, 03:23 

#4 
Senior Member
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 11 
I have slight troubles understanding what you want, but is this it: http://www.openfoam.org/version2.2.0...processing.php (see img.)?
If yes, then look at the appropriate motorBike tutorial (in the incompressible/simpleFoam folder). 

November 25, 2013, 05:23 

#5  
Senior Member
Join Date: Dec 2010
Posts: 135
Rep Power: 7 
@kkpal:
Quote:
Thanks for that idea, I found a workaround with paraview and got my sectional lift LES Drag Coefficient Problem? I am doing kind of a DNS and my results are matching perfectly. Used a modified icoFoam Solver 2nd order spatial and 4th order in time. Are you using a structured mesh? Is your case compressible or incompressible? Did you set the right projected reference Length and area in your controldict for calculation of cd? Did you also reconstructed CD value with paraFoam and compared it with OpenFoam Output? @hakoon Thanks for your answer, but this is not the answer for my problem. I solved by a workaround. 

November 26, 2013, 08:50 

#6 
Senior Member
Join Date: Jan 2013
Posts: 111
Rep Power: 5 
Later I found out Re=1000 could not be a good validation for LES, I tried Re=3900 and the results worked out well. There are some discussion about Re=1000 in this thread
Simulation Flow Around cylinder 3D Could you please elaborate how did you get your Cd and Cl from paraview? I'm very interested in that. 

November 27, 2013, 03:48 

#7 
Senior Member
Join Date: Dec 2010
Posts: 135
Rep Power: 7 
Hi kkpal,
These are my paraviewsteps in short: 1)Extracted Surface 2)Generate Surface Normals 3)Calculated: cp = normals*p; >Pressure Drag cf=normals*mu*wallGradU_y > Friction Drag 4) Build Integral Drag=S(cp_X+cf_X)dA*2/rho*v*v*A Hope this helps, but there are some differences to OpenFoam Output > could be caused by building derivates or integration failures 

November 27, 2013, 22:26 

#8 
Senior Member
Join Date: Jan 2013
Posts: 111
Rep Power: 5 
thanks very much! I will try your method!


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
lift & drag coefficient on airfoil  n. natik  FLUENT  8  March 31, 2015 19:02 
airfoil lift and drag coefficient  amir_14  FLUENT  5  January 1, 2013 09:30 
Drag and Lift coefficient for diffusor  rego  STARCCM+  3  May 7, 2012 18:05 
Automotive test case  vinz  OpenFOAM Running, Solving & CFD  98  October 27, 2008 09:43 
drag and lift coefficient  Noé  CDadapco  5  July 13, 2004 10:21 