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

simpleFoam convergence problems for pipe flow problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2012, 15:40
Default simpleFoam convergence problems for pipe flow problem
  #1
New Member
 
Join Date: Feb 2011
Posts: 1
Rep Power: 0
Mike Graham is on a distinguished road
I've been solving some very simple pipe flow problems using 1.6-ext's simpleFoam with a k-epsilon turbulence model and have had trouble getting good results. As I let my simulation go, in many time steps the p equation solver maxes out the number of iterations required before moving on (I've tried 1000 maximum for GAMG, 2500 for PCG), and eventually the pressures and velocities get huge and the continuity error gets large and the simulation dies due to a SIGFPE.

I've adjusted my p solver and parameters, my grad(p) schemes and various divSchemes, my relaxation parameters, and my initial conditions, but the solution inevitably blows up. If I change my boundary conditions to inlet-pressure/outlet-velocity or pressure/pressure, I don't blow up, but I do not achieve convergence.

Can anyone explain why I'm having trouble or advise me how to proceed?

One case that I ran is posted at http://dl.dropbox.com/u/19693178/examplecase.tar.gz

Here are some of my input files

0/fvSolution
Code:
solvers
{
    p
    {
        solver           GAMG;
        tolerance        1e-6;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;
    };

    U
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    };

    k
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    };

    epsilon
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    };
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
    p               0.3;
    U               0.6;
    k               0.2;
    epsilon         0.2;
}
fvSchemes
Code:
ddtSchemes
{
    default steadyState;
}

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

divSchemes
{
    default         none;
    div(phi,U)      Gauss linearUpwindV Gauss linear;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon)  Gauss upwind;
    div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

fluxRequired
{
    default         no;
    p;
}
0/U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    Inlet
    {
        type            fixedValue;
        value           uniform (-1.68669369 -7.75789138 0);
    }
    Outlet
    {
        type            zeroGradient;
    }
    Sym
    {
        type            symmetryPlane;
    }
    Wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
}
0/p
Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    Outlet
    {
        type            fixedValue;
        value           uniform -29.89;
    }

    Inlet
    {
        type            zeroGradient;
    }

    Wall
    {
        type            zeroGradient;
    }

    Sym
    {
        type            symmetryPlane;
    }
}
0/k
Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 90.583e-6;

boundaryField
{
    Inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    Outlet
    {
        type            zeroGradient;
    }

    Wall
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    Sym
    {
        type            symmetryPlane;
    }
}
0/epsilon
Code:
dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 12.89e-6;

boundaryField
{
    Inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    Outlet
    {
        type            zeroGradient;
    }

    Wall
    {
        type            epsilonWallFunction;
        value           $internalField;
    }

    Sym
    {
        type            symmetryPlane;
    }
}
0/nut
Code:
dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    Inlet
    {
        type            calculated;
        value           uniform 0;
    }

    Outlet
    {
        type            calculated;
        value           uniform 0;
    }

    Wall
    {
        type            nutWallFunction;
        value           uniform 0;
    }

    Sym
    {
        type            symmetryPlane;
    }
}
Mike Graham 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
Multiphase Flow - Convergence Problems bastian26 FLUENT 0 September 19, 2011 07:26
Convergence Problem icoFoam steady flow over an airfoil Lucas OpenFOAM Running, Solving & CFD 5 February 18, 2011 13:46
Convergence problem with simpleFoam solver on tilted channel daev OpenFOAM Running, Solving & CFD 1 December 23, 2010 17:20
Problem in convergence for simulate wind flow around building lalula2 FLUENT 2 June 8, 2009 03:13
Transient flow convergence problems Syed Mohammed Qamar Hasnain CD-adapco 0 March 20, 2007 07:17


All times are GMT -4. The time now is 06:25.