CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User 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 Search this Thread Display Modes
Old   August 1, 2013, 06:56
Default PisoFoam and PimpleFoam for my problem can not converge
  #1
Senior Member
 
Join Date: Jun 2011
Posts: 163
Rep Power: 14
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, 33 views)
File Type: jpg 2.jpg (89.8 KB, 36 views)
File Type: jpg 3.jpg (89.3 KB, 27 views)
File Type: jpg 4.jpg (97.3 KB, 23 views)
mechy 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
pimpleFoam vs simpleFoam vs pisoFoam vs icoFoam? phsieh2005 OpenFOAM Running, Solving & CFD 45 March 22, 2021 10:14
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 18:21.