- **OpenFOAM Running, Solving & CFD**
(*https://www.cfd-online.com/Forums/openfoam-solving/*)

- - **Problem obtaining right Cd**
(*https://www.cfd-online.com/Forums/openfoam-solving/113441-problem-obtaining-right-cd.html*)

Problem obtaining right CdHi 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: http://en.wikipedia.org/wiki/Turbulence_kinetic_energy 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; functions { forces { type forces; functionObjectLibs ("libforces.so"); //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; } forceCoeffs { // 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 ("libforces.so"); 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! |

All times are GMT -4. The time now is 21:37. |