CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Problem obtaining right Cd

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

LinkBack Thread Tools Display Modes
Old   February 19, 2013, 17:19
Exclamation Problem obtaining right Cd
New Member
Join Date: Oct 2012
Posts: 3
Rep Power: 5
matdgonz is on a distinguished road
Hi guys, quite new to openFoam.

I'v been running these 2 cases recently. A cylinder as well as a square cylinder.

I decided on using a value of Re of 100,000 so as to compare with a book.
Im using pisoFoam, k-e model and used initial values of k and e from the equations given here:

Currently running on a single core (dont know who to run in parallel just yet)

Meshing for both are quite dense but i get values of Cd of 0.7 for cylinder (should 1.2) and Cd of 1.7 for sqaure cylinder (should be 2.1)

Could these inaccurate values be a meshing problem? I did increase the mesh and after almost a 6hr simulation, i still got roughly the same Cd for the square cylinder. so im not sure exactly what my next step is.

Control dic:

application pisoFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 100;

deltaT 1;

writeControl timeStep;

writeInterval 5;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

type forces;
functionObjectLibs (""); //Lib to load
patches (cylinder); // change to your patch name
rhoName rhoInf;
rhoInf 1.1614; //Reference density for fluid
CofR (0 0 0); //Origin for moment calculations
outputControl timeStep;
outputInterval 1;
// rhoInf - reference density
// CofR - Centre of rotation
// dragDir - Direction of drag coefficient
// liftDir - Direction of lift coefficient
// pitchAxis - Pitching moment axis
// magUinf - free stream velocity magnitude
// lRef - reference length
// Aref - reference area
type forceCoeffs;
functionObjectLibs ("");
patches (cylinder);
rhoName rhoInf;
rhoInf 1.1614;
CofR (0 0 0);
liftDir (0 1 0);
dragDir (1 0 0);
pitchAxis (0 0 1);
magUInf 1;
lRef 1; //
Aref 0.1; //

outputControl timeStep;
outputInterval 1;


Any help would be much appreciated!
matdgonz is offline   Reply With Quote


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
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43
Gambit - meshing over airfoil wrapping (?) problem JFDC FLUENT 1 July 11, 2011 05:59
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 14:52

All times are GMT -4. The time now is 19:08.