CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   OpenFoam validation (

ranas July 7, 2009 05:56

OpenFoam validation

I am validating OpenFOAM on a static 2-D cylinder at Re=200. The radius of the cylinder is 0.5m. The velocity of incoming flow is 1m/s. Mu equals to 0.005.

The full grid is shown here:


I used the following schemes:
default CrankNicholson 0.5;
default Gauss linear;
grad(p) Gauss linear;
default none;
div(phi,U) Gauss linear;
default none;
laplacian(nu,U) Gauss linear corrected;
laplacian(1|A(U),p) Gauss linear corrected;
default linear;
interpolate(HbyA) linear;
default corrected;
default no;

The standard ICCG and BICCG solvers are used untill convergence to resectively 1e-6 and 1e-5 is obtained. The time step was chosen to be 0.01.

Also the standard boundary conditions are applied:
inlet: U fixed and zeroGradient for p
outlet: p fixed and zeroGradient for U
upper/lower boundary: zeroGradient for p and for U

type forceCoeffs;
functionObjectLibs ("");
patches (circle); //change to your patch name
rhoInf 1;
CofR (0 0 0);
liftDir (0 1 0);
dragDir (1 0 0);
pitchAxis (0 0 0);
magUInf 1;
lRef 1;
Aref 1.0;

BrazaM (1985 N um )' result: time average of cd=1. 31; amplitude of cl=0. 65; st=0. 194.
my result: time average of Cd is about 1.9, larger than 1.31; the amplitude of Cl is about 1.0 , more than 0.65; but the st number equals to 0.196 close to the result.

Is someone having any ideas on what's wrong with it? Do I use the proper schemes? Are the bounday conditions correct?


All times are GMT -4. The time now is 12:09.