OpenFoam flow around a cylinder
 

hi guys!

I am a new comer and i need yours help! :o

I like to simulate a turbulent flow over a circular cylinder and I have find an example from the OpenFoam WIKI.

Here is the link
http://openfoamwiki.net/index.php/Be...coftac_ufr2-02

However , when I run this case, it shows an error keyword div((nuEff*dev(T(grad(U))))) is undefined in dictionary "DATA/0115/ufr2-02/system/fvSchemes.divSchemes"

How should I solve this problem? and also I can not run this case in parallel it shows that keyword div((nuEff*dev(T(grad(U))))) is undefined in dictionary "IOstream.divSchemes"

here is The fvSchemes.

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / 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;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,k) Gauss limitedLinear 0.95;
div(phi,B) Gauss limitedLinear 1;
div(phi,nuTilda) Gauss limitedLinear 1;
div(B) Gauss linear;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DBEff,B) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

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

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}


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

Thank you for your help.

Best wishes

Yibo

Hi,

you should change one line in fvSchemes:

to

Also you need to add following to fvSolutions (just after dictionary for k):

At least you'll be able to run the case after these modifications, guess you'll need to modify the case further to get meaningful results.

Hi，alexeym,

Thank you for your help. That is very helpful for me.

Best wishes,
Yibo:)

Is there another way for me to see this link?