# forceCoeffs and rhoInf dependency

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Search this Thread Display Modes
 July 5, 2013, 16:29 forceCoeffs and rhoInf dependency #1 Member   Florian Join Date: Nov 2009 Posts: 59 Rep Power: 11 Hello, I have an incompressible simpleFoam solution of a flow field. I use forceCoeffs to calculate Cd etc. on some patches. Code: ```coeffsCar { type forceCoeffs; functionObjectLibs ( "libforces.so" ); outputControl timeStep; outputInterval 10; patches \$carPatches; pName p; UName U; rhoName rhoInf; log true; rhoInf 3; liftDir (0 0 1); dragDir (1 0 0); pitchAxis (0 1 0); CofR (0.56725 0 -0.1274); magUInf 45; lRef 1.5; Aref 0.4;``` I tried it on exactly the same case, one time set rhoInf to 3, another time to the actual value of 1.4. Both forceCoeffs got me exactly the same results, regardless of the value of rhoInf. Is there some bug in OF 2.1.1 or a bug in my understanding of physics? Force on a patch is F=p*A, drag coefficient is Cd = F / (0.5 rhoInf U^2). Why is there no influence of rho? Thanks!

 July 24, 2013, 13:51 #2 Senior Member   Joachim Join Date: Mar 2012 Location: Providence, RI Posts: 144 Rep Power: 10 Hi Florian, if the solver is incompressible, the pressure computed by OpenFOAM is p/rho and not p. Hence, you don't need to know the density to get your coefficient (just divide your pressure force by 0.5*Ue^2 A). You can write anything you want for the density in forceCoeffs, it is not read when computing the coefficients. Regards, Joachim

 July 25, 2013, 10:04 #3 Member   Florian Join Date: Nov 2009 Posts: 59 Rep Power: 11 Yeah, that became clear after posting this. But I'm still riddled why OF requires setting a value of rhoInf. In the source code, as far as I can tell, it's determined by units if p' = p/rho (incompressible) or p' = p. Regards...

 August 9, 2018, 03:55 #4 New Member   zein elserfy Join Date: May 2018 Posts: 6 Rep Power: 3 hi Joachim, I want to ask you how to calculate the value of Cd and Cl ? by calculating it by dividing pressure force by (0.5*v*v*A),so this mean that results of openfoam for cd and cl are not correct ?

 Thread Tools Search this Thread Search this Thread: Advanced Search 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 junkie71189 OpenFOAM Running, Solving & CFD 31 March 11, 2017 21:56 srini_esi OpenFOAM 0 October 19, 2012 06:34 appa OpenFOAM Running, Solving & CFD 1 June 7, 2012 05:04 bigbang OpenFOAM 0 July 13, 2011 09:38

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

 Contact Us - CFD Online - Privacy Statement - Top