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

Sectional Drag and lift coefficient?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 22, 2013, 09:03
Default Sectional Drag and lift coefficient?
  #1
Senior Member
 
Join Date: Dec 2010
Posts: 135
Rep Power: 6
eRzBeNgEl is on a distinguished road
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?
eRzBeNgEl is offline   Reply With Quote

Old   November 22, 2013, 09:48
Default
  #2
Member
 
Kai Zhang
Join Date: Jan 2013
Location: yokohama, Japan
Posts: 83
Rep Power: 4
kkpal is on a distinguished road
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.
kkpal is offline   Reply With Quote

Old   November 22, 2013, 09:53
Default
  #3
Member
 
Kai Zhang
Join Date: Jan 2013
Location: yokohama, Japan
Posts: 83
Rep Power: 4
kkpal is on a distinguished road
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?
kkpal is offline   Reply With Quote

Old   November 25, 2013, 03:23
Default
  #4
Senior Member
 
Hkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 10
haakon will become famous soon enough
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).
haakon is offline   Reply With Quote

Old   November 25, 2013, 05:23
Default
  #5
Senior Member
 
Join Date: Dec 2010
Posts: 135
Rep Power: 6
eRzBeNgEl is on a distinguished road
@kkpal:
Quote:
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.

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.
eRzBeNgEl is offline   Reply With Quote

Old   November 26, 2013, 08:50
Default
  #6
Member
 
Kai Zhang
Join Date: Jan 2013
Location: yokohama, Japan
Posts: 83
Rep Power: 4
kkpal is on a distinguished road
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.
kkpal is offline   Reply With Quote

Old   November 27, 2013, 03:48
Default
  #7
Senior Member
 
Join Date: Dec 2010
Posts: 135
Rep Power: 6
eRzBeNgEl is on a distinguished road
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
eRzBeNgEl is offline   Reply With Quote

Old   November 27, 2013, 22:26
Default
  #8
Member
 
Kai Zhang
Join Date: Jan 2013
Location: yokohama, Japan
Posts: 83
Rep Power: 4
kkpal is on a distinguished road
thanks very much! I will try your method!
kkpal 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
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 STAR-CCM+ 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 CD-adapco 5 July 13, 2004 10:21


All times are GMT -4. The time now is 02:27.