CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

OpenFoam flow around a cylinder

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree3Likes
  • 3 Post By alexeym

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 15, 2014, 13:01
Question OpenFoam flow around a cylinder
  #1
Member
 
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 12
Scabbard is on a distinguished road
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
Scabbard is offline   Reply With Quote

Old   January 15, 2014, 15:25
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

you should change one line in fvSchemes:

Code:
div((nuEff*dev(grad(U).T()))) Gauss linear;
to

Code:
div((nuEff*dev(T(grad(U))))) Gauss linear;
Also you need to add following to fvSolutions (just after dictionary for k):

Code:
    kFinal
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-07;
        relTol          0;
    }
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.
jimbean, Scabbard and utkunun like this.
alexeym is offline   Reply With Quote

Old   January 15, 2014, 17:21
Default
  #3
Member
 
Join Date: Dec 2013
Location: Newcastle
Posts: 54
Rep Power: 12
Scabbard is on a distinguished road
Quote:
Originally Posted by alexeym View Post
Hi,

you should change one line in fvSchemes:

Code:
div((nuEff*dev(grad(U).T()))) Gauss linear;
to

Code:
div((nuEff*dev(T(grad(U))))) Gauss linear;
Also you need to add following to fvSolutions (just after dictionary for k):

Code:
    kFinal
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-07;
        relTol          0;
    }
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
Scabbard is offline   Reply With Quote

Old   August 26, 2022, 17:48
Default
  #4
New Member
 
Gavin.
Join Date: Aug 2022
Posts: 4
Rep Power: 3
GavinRe is on a distinguished road
Quote:
Originally Posted by Scabbard View Post
Is there another way for me to see this link?
GavinRe is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


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