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

Question to liftDrag

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 29, 2007, 23:17
Default Srinath Madhavan, I figured
  #21
Member
 
Alessandro Spadoni
Join Date: Mar 2009
Location: Atlanta, GA
Posts: 65
Rep Power: 9
gtg627e is on a distinguished road
Srinath Madhavan,

I figured out how to compute the drag for turbulent cases; however, I have doubts about the results obtained from liftDrag.C.

I followed the instructions in the following thread:

http://www.cfd-online.com/OpenFOAM_D...s/1/5067.html?


The reason why I am doubful is because I have a simple naca0012 at 0 angle of attack, and OpenFOAM predicts a total drag of 0.00314, which is well below the experimental value of 0.0058 from "theory of wing sections." Later this evening, I tried the same case at 4-deg angle of attack, and now I get a total drag of 0.008 and a lift coefficient of 0.05. For this case, the drag should be approximately 0.006 and the coefficient of lift shuld be about 0.5.

I ran checkMesh and got no errors.
Everything is incompressible, inlet V = 50 m/s
chord = 1 m
k-epsilon model:

inlet k = 0.375
inlet epsilon 0.03

Everything is reasonable, the analysis converges reasonably well, but the force coefficients are off. As Vincent suggested in

http://www.cfd-online.com/OpenFOAM_D...s/1/5067.html?

I will compare pressure coefficent, velocity distribution, etc. For the pressure coefficient, however, I will have to do a little coding.

It is possible that my analysis is ill-posed, but it converges well, and both velocity and pressure seem very reasonable.

Thank you,

Alessandro
gtg627e is offline   Reply With Quote

Old   December 8, 2007, 08:51
Default Alessandro, I understand your
  #22
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
Alessandro, I understand your problem. However, the minute turbulence models enter the picture, I switch off the need for accurate experimental predictions in my head. A bad habit, but nevertheless quite true most of the time. Anyways, I would be interested in knowing whether you have tried the same case with Fluent or CFX and if so what dimensionless force coefficients (lift, drag etc.) do they predict?
msrinath80 is offline   Reply With Quote

Old   March 17, 2008, 12:03
Default Could someone verify that the
  #23
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
Could someone verify that the turbulent drag calculation using liftDrag is correct......we are finding differences....anyone?

Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   March 20, 2008, 14:36
Default Hi Frank & others, I am som
  #24
Senior Member
 
Srinath Madhavan (a.k.a pUl|)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13
msrinath80 is on a distinguished road
Hi Frank & others,

I am somewhat confused as to why most posts in this forum suggest that modifications are required to the liftDrag utility if one has to obtain lift/drag predictions when using turbulence models. If the function of a turbulence model is to simulate (or more accurately, emulate) turbulent behavior in the fluid and if this is reasonably well done by the turbulence model in question, then turbulent behavior should automatically show up in the result (i.e. in the primitive variables u,v,w,p). Then all that remains is using these results in the integration of shear and pressure forces over the walls in question, which should directly give the correct lift and drag coefficients. Turbulent viscosity should not enter the picture at this stage as it is merely a means to get to the correct (hopefully) u,v,w and p. So when lift/drag forces are calculated, one should only use the actual density/viscosity of the fluid in question and not the turbulent viscosity?

To put it in a nutshell, if the velocity and pressure fields are accurate in space (and in time for transient problems), then the liftDrag utility should give the correct forces without any modifications.

I would appreciate if someone can correct me if I am wrong.
msrinath80 is offline   Reply With Quote

Old   March 20, 2008, 16:47
Default Hi Srinath, In fact, this i
  #25
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
Hi Srinath,

In fact, this is exactly what I meant by my previous question. A turbulence model is used to give a better solution of u,v,w,p when the flow is turbulent. So, in fact the boundary layer becomes fuller (compared to laminar flow) which has its influence on the viscous drag, just like the laminar case. In my view, the turbulent viscosity term is not needed in the drag/lift calculation......

Anyone else??

Frank

btw, I did not wrote this utility, but made some modifications such that integration in icoFoam i.e. becomes easier....
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   March 25, 2008, 05:52
Default Works like this: tauw = nu
  #26
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
Works like this:

tauw = nu * dU/dy_wall right?

Except we do not have dU/dy_wall, we have dU/dy_0 between the first node and the wall. Thus we calculate wall shear using a wall function. To keep things compact in the momentum equation we stick the wall function in the turbulent viscosity so that:

tauw = nuEff_wall * dU/dy_0

So the use of turbulent viscosity to calculate the wall shear is purely an artifact used by OpenFOAM to store the influence of the wall function - nuEff_wall.
eugene is offline   Reply With Quote

Old   March 25, 2008, 06:14
Default Allright, so in the turbulent
  #27
Senior Member
 
Frank Bos
Join Date: Mar 2009
Location: The Netherlands
Posts: 340
Rep Power: 10
lr103476 is on a distinguished road
Allright, so in the turbulent case the viscous force calculation contains some turbulent contribution.....via wall function.

Then, why is there still need for a turb force contribution,present in liftDrag, derived from turbulence->R().

Frank
__________________
Frank Bos
lr103476 is offline   Reply With Quote

Old   March 25, 2008, 06:46
Default Because (for example in kEpsil
  #28
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 13
eugene is on a distinguished road
Because (for example in kEpsilon) turbulence->R() is calculated from

2/3*I*k - nut_*2*symm(grad(U))

The wall boundary of nut_ contains the effect of the wall function. (See "#include "wallViscosity.H" line 88 kEpsilon.C)
eugene is offline   Reply With Quote

Old   March 25, 2008, 07:04
Default Formally and for consistency w
  #29
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,781
Rep Power: 22
hjasak will become famous soon enough
Formally and for consistency we should have a correction for snGrad(U) on the wall to get the correct drag though.

Does anyone feel like changing them all?

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   September 19, 2014, 03:38
Default Drag coefficient amplitude
  #30
New Member
 
Join Date: Sep 2014
Posts: 10
Rep Power: 3
Konstantine is on a distinguished road
Hello everyone,

I need some assistance on understanding why in the drag coefficient time diagram, i am getting the correct frequency according to the Re number but the amplitude is totally different from past studies in papers. more specifically it should oscillate around 1.44 and mine oscillates around 3. Also in other cases with different mesh it drops to 0.3. although frequency is correct. Does anybody know about amplitude dependency on either numerical or physical factors that could help me understand why. I am using openfoam and i have also used another code (non commercial) that shows the same thing.

Regards
Konstantine
Konstantine 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
LiftDrag tool nuovodna OpenFOAM Running, Solving & CFD 45 September 2, 2009 17:56
LiftDrag for 141 ryan_m OpenFOAM Running, Solving & CFD 2 August 24, 2009 21:26
LiftDrag coefficient in LES fabian_korn OpenFOAM Post-Processing 1 September 22, 2008 02:34
LiftDrag utility question msrinath80 OpenFOAM Running, Solving & CFD 8 March 28, 2008 11:55
LiftDrag utility not available guggi OpenFOAM Running, Solving & CFD 1 August 2, 2006 12:36


All times are GMT -4. The time now is 04:44.