CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   PisoFoam and PimpleFoam for my problem can not converge (http://www.cfd-online.com/Forums/openfoam/121610-pisofoam-pimplefoam-my-problem-can-not-converge.html)

mechy August 1, 2013 05:56

PisoFoam and PimpleFoam for my problem can not converge
 
4 Attachment(s)
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);
}

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



All times are GMT -4. The time now is 04:58.