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

PisoFoam and PimpleFoam for my problem can not converge

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 1, 2013, 05:56
Default PisoFoam and PimpleFoam for my problem can not converge
  #1
Senior Member
 
Join Date: Jun 2011
Posts: 141
Rep Power: 6
mechy is on a distinguished road
I want to solve flow around a cylinder and thin plate
the thin plate is placed horizontally after cylinder

I have test many different settings in fvSchemes and fvSolutions
but all of them give uncorrected results


I will be so grateful if anybody can help me


Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default Euler;
}

gradSchemes
{
    default           leastSquares ;
//    grad(p)         leastSquares ;
//    grad(U)         leastSquares ;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;//Gauss linearUpwind grad(U);//
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,omega)  Gauss limitedLinear 1;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
     div((nuEff*dev(grad(U).T())))  Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;//Gauss linear limited 0.5;//
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
//    pcorr           ;
    p;
}

// ************************************************************************* //
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.3                                   |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version         2.0;
    format          ascii;
    class           dictionary;
    object          fvSolution;
}

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

solvers
{
    p
    {
        solver           GAMG;
        tolerance        1e-5;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       1;
        nPostSweeps      2;
        cacheAgglomeration true;
        directSolveCoarsest true;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 400;
        mergeLevels      1;
        minIter          1;
    }
    pFinal
    {
        solver           GAMG;
        tolerance        1e-9;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       1;
        nPostSweeps      2;
        cacheAgglomeration true;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 400;
        mergeLevels      1;
        minIter          1;
    }

    "(U|k|omega)"
    {
//        solver          PBiCG;
//        preconditioner  DILU;

        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-5;
        relTol           0;
        nSweeps          1;     
        minIter          1;
    }
    "(UFinal|kFinal|omegaFinal)"
    {
//        solver          PBiCG;
//        preconditioner  DILU;

        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-6;
        relTol           0.0;
        nSweeps          1;   
        minIter          1;
    }
    cellMotionU
    {
        solver          GAMG;
        tolerance       1e-6;
        relTol          1e-3;
        smoother        GaussSeidel;
        cacheAgglomeration true;
        directSolveCoarsest true;
        nCellsInCoarsestLevel 20;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }
}

PIMPLE
{
    nOuterCorrectors    3;
    nCorrectors         2;
    nNonOrthogonalCorrectors 3;

    pRefCell            0;
    pRefValue           0;
    correctPhi          no;    
//    momentumPredictor yes;

}

PISO
{
    nOuterCorrectors 2;
    nCorrectors     2;   
    nNonOrthogonalCorrectors 3;
    pRefCell        0;
    pRefValue       0;

//    ddtPhiCorr no;
}


relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.7;
        k               0.7;
        omega           0.7;
    }
}
potentialFlow
{
    nNonOrthogonalCorrectors 0;
}
cache
{
    grad(U);
}

// ************************************************************************* //
Attached Images
File Type: jpg 1.jpg (88.0 KB, 16 views)
File Type: jpg 2.jpg (89.8 KB, 19 views)
File Type: jpg 3.jpg (89.3 KB, 14 views)
File Type: jpg 4.jpg (97.3 KB, 13 views)
mechy 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
pimpleFoam vs simpleFoam vs pisoFoam vs icoFoam? phsieh2005 OpenFOAM Running, Solving & CFD 29 June 13, 2014 09:20
I got complete different results for pisofoam and pimplefoam? lakeat OpenFOAM Running, Solving & CFD 6 December 27, 2012 05:25


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