[leeminjae]

Hello ladies & gentlemen
I'm working to simulate "flow around cylinder at Re=40"
I generated the mesh using Gambit.
It's a very simple simulation.
the objective of simulation is to get a good matching results of Cd & Cl
But I think my result is not that matching with an another paper.
I don't know what's wrong.
I attached my code below.



And another problem, Co(Courant Number) is around 1.6
Is it O.K ?
The property I used
Ux 0.071576 (m/s)
Diameter 0.01 (m)
Nu 1.7894e-05
Re 40



I'm sorry to attach my file not comfortably.
I don't know how to attach my code with files.

I'm looking forward to taking advice from you.









/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/

forceCoeffs1
{
type forceCoeffs;

functionObjectLibs ( "libforces.so" );

outputControl outputTime;

log yes;

patches ( cy );
pName p;
UName U;
rhoName rhoInf; // Indicates incompressible
log true;
rhoInf 1.225; // Redundant for incompressible
liftDir (0 1 0);
dragDir (1 0 0);
CofR (0 0 0); // Axle midpoint on ground
pitchAxis (0 0 1);
magUInf 0.4;
lRef 0.01; // Wheelbase length
Aref 0.0001; // Estimated


}


// ************************************************** *********************** //
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
down
{
type symmetryPlane;
}
up
{
type symmetryPlane;

}
wall
{
type zeroGradient;

}
back
{
type symmetryPlane;

}
front
{
type symmetryPlane;

}
outlet
{
type fixedValue;
value uniform 0;

}
inlet
{
type zeroGradient;

}
}

// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
down
{
type symmetryPlane;

}
up
{
type symmetryPlane;

}
wall
{
type fixedValue;
value uniform (0 0 0);

}
back
{
type symmetryPlane;

}
front
{
type symmetryPlane;

}
outlet
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform (0.071576 0 0);

}
}
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

libs
(
"libOpenFOAM.so"
"libincompressibleTurbulenceModel.so"
"libincompressibleRASModels.so"
);

application icoFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 150;

deltaT 0.005;

writeControl timeStep;

writeInterval 200;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

functions
{

#include "forceCoeffs"
}


// ************************************************** *********************** //

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default backward;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear orthogonal;
laplacian((1|A(U)),p) Gauss linear orthogonal;
}

interpolationSchemes
{
default linear;
interpolate(HbyA) linear;
}

snGradSchemes
{
default orthogonal;
}

fluxRequired
{
default no;
p ;
}


// ************************************************** *********************** //

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}

U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
}

PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}


// ************************************************** *********************** //

[joshm]

This could be any number of issues. If you post your case as a zip file, people might take a look at it for you.