
[Sponsors] 
October 28, 2012, 20:00 
Drag and lift force using OpenFoam2.1.1

#41 
Member
R. P.
Join Date: Jul 2010
Posts: 72
Rep Power: 9 
Hi everyone,
I am trying to calculate the drag, lift and pitchingmoment coefficients for a flow over a probe. The code structure below works fine when the flow is parallel to the probe. However, when I changed the angle of attack I got a rude difference for lift and pitchingmoment when compared with experimental data (the drag values is OK). I guess, this difference is related to the vectors (liftDir (0 1 0); dragDir (1 0 0); pitchAxis (0 0 1)). I this this vectors is just valid for a parallel flow, isn't it ? Anybody knows how to change this vectors according to a certain angle of attack ? Thank you very much. Rophys forces { type forces; enabled true; functionObjectLibs ( "libforces.so" ); outputControl outputTime; patches (probe); directForceDensity true; fDName fDMean; CofR (0.02 0 0); log on; } forceCoeffs { type forceCoeffs; functionObjectLibs ("libforces.so"); patches (probe); outputControl outputTime; fDName fDMean; rhoInf 1.73e5; CofR (0.02 0 0); liftDir (0 1 0); dragDir (1 0 0); pitchAxis (0 0 1); magUInf 1503.1; lRef 0.05; Aref 9.817e4; } 

October 29, 2012, 04:29 

#42 
Member
Join Date: Apr 2012
Location: France
Posts: 72
Rep Power: 7 
Hi,
We have two cases : If you change the components of the velocity (and don't rotate the geometry), the lift and drag direction must change. If you rotate the geometry, the lift and drag direction don't change. For exemple, if you change your flow direction (angle of attack) like this : (V*cos(alpha) 0 V*sin(alpha)). You must have to change the direction like this : liftDir(sin(alpha) 0 cos(alpha)) dragDir (cos(alpha) 0 sin(alpha)) I hope that I have answered to your question 

October 29, 2012, 05:20 

#43 
Member
R. P.
Join Date: Jul 2010
Posts: 72
Rep Power: 9 
Hi Rider,
Thanks for your response. It was clear enough. In my case, the mesh still the same and I just rotate the geometry (the liftDir and dragDir remained the same, liftDir (0 1 0); dragDir (1 0 0)). Even doing this the lift and pitchingmoment is not matching with experimental results. It is something around 2030% difference. The drag force is matching quite well, as well the Ca and Cn. If anybody have some suggestions in how to solve this, I'll appreciate. Cheers. 

October 29, 2012, 06:10 

#44 
Member
Join Date: Apr 2012
Location: France
Posts: 72
Rep Power: 7 
Your domain is enough biggest to capture all the tail?
What is your y+ value ? 

October 29, 2012, 06:34 

#45 
Member
R. P.
Join Date: Jul 2010
Posts: 72
Rep Power: 9 
I am using the OpenFOAM but I am not using a NavierStokes solver. It is called DSMC and for it, I don't have problems with y+


October 29, 2012, 06:37 

#46 
Member
Join Date: Apr 2012
Location: France
Posts: 72
Rep Power: 7 
Ok, it's your company solver or a free solver ?


October 29, 2012, 06:46 

#47 
Member
R. P.
Join Date: Jul 2010
Posts: 72
Rep Power: 9 
Yes it is free....DSMC solver is part of the OpenFoam. It is used, for example, to simulate reentry problems in very high altitudes where the NS is no long valid (is too rarefied to solve NS equations). It take into account the molecular movement instead to solve differential equation but when used for "low altitude" the results should be quite close to CFD results.


October 29, 2012, 07:06 

#48 
Member
Join Date: Apr 2012
Location: France
Posts: 72
Rep Power: 7 
Ok.
Thanks. 

May 18, 2013, 07:52 
why do you remove Uname and pName ,and add rhoName rhoInf?

#49 
New Member
jiujiumin
Join Date: May 2013
Posts: 9
Rep Power: 6 
i am a new user in openFoam ,I want to know why it do not use the U and p in caculating the forces and forcecoeffs? And i also want to know why it always has the error :From function void forces::read(const dictionary& dict)
in file forces/forces.C at line 277 Could not find U, p or rho in database. Deactivating forces. Because I have the U and p files in the 0 file,so i am so confused about it . 

January 3, 2014, 04:38 
lift and drag coeffs in paraview

#50 
Member

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? 

October 31, 2016, 17:29 
Could not find rho

#51 
New Member
Join Date: Oct 2016
Posts: 22
Rep Power: 3 
Hey,
I want to calculate the drag and lift coefficients of an 2D airfoil. But when I use the force and force coeffs script I always get an error saying "could not find rho". Without the drag and lift it works fine and rho is defined in transport properties. I am using openfoam 4.1. Anybody have an idea or had a similar problem. Thanks! 

November 3, 2016, 08:41 

#52 
New Member
Join Date: Oct 2016
Posts: 22
Rep Power: 3 
works with OpenFoam 3. Hence there should be a bug in OpenFoam4.


December 11, 2017, 03:10 
calculating forces

#53 
Senior Member
A. Min
Join Date: Mar 2015
Posts: 158
Rep Power: 4 
Hi all
I saw in forces code that shear stress is calculated by: shear stress =  mu*gama in code: rho()*nu*dev(twoSymm(fvc::grad(U))) what is the minus in this formula? Now I want to calculate forces for viscoelastic fluid. should I write: total shear stress = tau  mu*gama that tau is polymeric stress. Thanks 

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 
Fluent Good Lift coefficient BAD drag coefficient  Rif  Main CFD Forum  4  March 9, 2010 11:52 
Drag coefficient & Lift coefficient??  mehrdadeng  CFX  15  December 9, 2009 22:49 
Automotive test case  vinz  OpenFOAM Running, Solving & CFD  98  October 27, 2008 09:43 
drag and lift coefficient  Noé  Siemens  5  July 13, 2004 10:21 