CFD Online Discussion Forums

Hi,

i'm trying to simulate a 2D flow over an airfoil.
I'm using the simpleFoam with a k-epsilon turbulence model, these are my system files:

fvSolution:

Code:

/*--------------------------------*- C++ -*----------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |

|  \\    /   O peration     | Version:  2.0.1                                 |

|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |

|    \\/     M anipulation  |                                                 |

\*---------------------------------------------------------------------------*/

FoamFile

{

    version     2.0;

    format      ascii;

    class       dictionary;

    location    "system";

    object      fvSolution;

}

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



solvers

{

    p

    {

        solver          GAMG;

        tolerance       1e-06;

        relTol          0.1;

        smoother        GaussSeidel;

        nPreSweeps      0;

        nPostSweeps     2;

        cacheAgglomeration true;

        nCellsInCoarsestLevel 10;

        agglomerator    faceAreaPair;

        mergeLevels     1;

    }





    U

    {

        solver          PBiCG;

        preconditioner  DILU;

        tolerance       1e-05;

        relTol          0.1;

    }



    k

    {

        solver          PBiCG;

        preconditioner  DILU;

        tolerance       1e-05;

        relTol          0.1;

    }



    epsilon

    {

        solver          PBiCG;

        preconditioner  DILU;

        tolerance       1e-05;

        relTol          0.1;

    }



    

}



SIMPLE

{

    nNonOrthogonalCorrectors 0;

    pRefCell        10;

    pRefValue       0;



    residualControl

    {

        p               1e-6;

        U               1e-6;

        "(k|epsilon|omega)" 1e-3;

    }

}



relaxationFactors

{

    p               0.3;

    U               0.7;

    k               0.7;

    epsilon         0.7;      



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

fvScheme

Code:

/*--------------------------------*- C++ -*----------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |

|  \\    /   O peration     | Version:  2.0.1                                 |

|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |

|    \\/     M anipulation  |                                                 |

\*---------------------------------------------------------------------------*/

FoamFile

{

    version     2.0;

    format      ascii;

    class       dictionary;

    location    "system";

    object      fvSchemes;

}

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



ddtSchemes

{

    default         steadyState;

}



gradSchemes

{

    default         Gauss linear;

    grad(p)         Gauss linear;

    grad(U)         Gauss linear;

}



divSchemes

{

    default         none;

    div(phi,U)      Gauss upwind;

    div(phi,k)      Gauss upwind;

    div(phi,epsilon) Gauss upwind;

    div((nuEff*dev(T(grad(U))))) 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               ;

}





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

The main problem is, there is no convergence especially for the pressure. I get this message:

Code:

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"

#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"

#2  Uninterpreted: 

#3  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"

#4  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"

#5  Foam::fvMatrix<double>::solve() in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/simpleFoam"

#6  Foam::incompressible::RASModels::kEpsilon::correct() in "/opt/openfoam201/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"

#7  

 in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/simpleFoam"

#8  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"

#9  

 in "/opt/openfoam201/platforms/linuxGccDPOpt/bin/simpleFoam"

Gleitkomma-Ausnahme

Has someone an idea how to solve this problem?

Denis