CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Vortex Street with pisoFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 10, 2011, 10:00
Default Vortex Street with pisoFoam
  #1
New Member
 
Join Date: Dec 2010
Posts: 2
Rep Power: 0
mikahs is on a distinguished road
Hello people...

I am fairly new to CFD and OpenFOAM, so I am not very sure about some results I am getting while simulating a vortex street with a flow around a cylinder.

I was asked to use pisoFoam with laminar flow, and simulate a given mesh (with the cylinder) at Re=20, Re=200 (the vortex street should appear here) and Re=2000.

At Re=20, I had no problems with convergence, with 8000 Iterations my residuals were below 10E-5.

However at Re=200 and at Re=2000 the residuals oscillate asymptotically.

Example: Re=200



Uploaded with ImageShack.us

I would like to know whether this behaviour in my residuals is normal for vortex street calculations, or if I am really having some problems of some other kind.

Ah, and I forgot to say... Courant-Nr<1 ... :/

Code:
ddtSchemes
{
    default         Euler;
}

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

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,epsilon) Gauss limitedLinear 1;
    div(phi,R)      Gauss limitedLinear 1;
    div(R)          Gauss linear;
    div(phi,nuTilda) Gauss limitedLinear 1;
    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(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

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

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
}
Code:
solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0.1;
    }

    pFinal
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          1e-4;
    }

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

}

PISO
{
    nCorrectors     3;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;
}
mikahs is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
karman vortex street help please SSeth STAR-CCM+ 7 January 10, 2011 12:31
3D cylinder and Vortex Street shedding issues josik_1982 FLUENT 2 July 17, 2010 11:19
Karman vortex street ??? anisa FLUENT 1 May 5, 2005 14:04
Karman vortex street Zhipeng FLUENT 5 June 24, 2002 12:00
Kármán vortex street in cavitating flow behind bodies in the cavitation tunnel L. Könözsy Main CFD Forum 0 April 17, 2000 14:16


All times are GMT -4. The time now is 21:22.