|
[Sponsors] |
January 15, 2014, 13:01 |
OpenFoam flow around a cylinder
|
#1 |
Member
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 12 |
hi guys!
I am a new comer and i need yours help! 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 |
|
January 15, 2014, 15:25 |
|
#2 |
Senior Member
|
Hi,
you should change one line in fvSchemes: Code:
div((nuEff*dev(grad(U).T()))) Gauss linear; Code:
div((nuEff*dev(T(grad(U))))) Gauss linear; Code:
kFinal { solver PBiCG; preconditioner DILU; tolerance 1e-07; relTol 0; } |
|
January 15, 2014, 17:21 |
|
#3 | |
Member
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 12 |
Quote:
Thank you for your help. That is very helpful for me. Best wishes, Yibo |
||
August 26, 2022, 17:48 |
|
#4 | |
New Member
Gavin.
Join Date: Aug 2022
Posts: 4
Rep Power: 3 |
Quote:
|
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
About Some Concepts:Laminar flow, turbulent flow, steady flow and time-dependent flow | Jing | Main CFD Forum | 8 | October 5, 2018 18:02 |
[Other] OpenFoam Flow over a Circular Cylinder | WolfgangS. | OpenFOAM Meshing & Mesh Conversion | 12 | March 3, 2014 11:53 |
pipe flow OpenFOAM | stix | OpenFOAM Running, Solving & CFD | 1 | March 4, 2013 17:44 |
benchmark: flow over a circular cylinder | goodegg | Main CFD Forum | 12 | January 22, 2013 12:47 |
Flow over a cylinder | Anna | Main CFD Forum | 9 | March 24, 2006 15:32 |