[Haroun]

Dear Foamers,

I have a case that I want to solve using a LRR and a SSG turbulence model. The LRR case is running fine but I keep getting a time step continuity error for the SSG model.

• My mesh should be all right (as I had prior error due to this, I worked it again and now it should be fine).
• I tried running the case from 0s and from 6s (taking the '6' folder solved by the LRR model).
• I tried changing my gradSchemes and my divSchemes.
• I tried to change the solver and its preconditionner.
• I tried to change the tolerances and the relaxation factors.
• I tried to decrease the delta T and switching on and off the adjustTimeStep option.
• I tried to set a first order time scheme (Euler) instead of a second (crankNicolson)

They all gives me the same error message :

time step continuity errors : sum local = 1.293e+65, global = 5.57821e+54, cumulative = 1.54365e+55
[2] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #2  ? in /lib64/libpthread.so.0
[2] #3  Foam::tmp<Foam::GeometricField<Foam::innerProduct<Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>::value_type, Foam::Vector<double> >::type, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::surfaceInterpolationScheme<Foam::Vector<double> >::dotInterpolate<Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) at ??:?
[2] #4  Foam::surfaceInterpolationScheme<Foam::Vector<double> >::dotInterpolate(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
[2] #5  ? at gaussLaplacianSchemes.C:?
[2] #6  Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::gammaSnGradCorr(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[2] #7  Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacian(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[2] #8  Foam::fv::laplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacian(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[2] #9  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, Foam::SymmTensor<double> >(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[2] #10  Foam::RASModels::SSG<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
[2] #11  ? at ??:?
[2] #12  __libc_start_main in /lib64/libc.so.6
[2] #13  ? at ??:?
[n075:76036] *** Process received signal ***
[n075:76036] Signal: Floating point exception (8)
[n075:76036] Signal code:  (-6)
[n075:76036] Failing at address: 0x11fc85d00012904
[n075:76036] [ 0] /lib64/libpthread.so.0(+0xf630)[0x7fa7b8f19630]
[n075:76036] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x7fa7b8f194fb]
[n075:76036] [ 2] /lib64/libpthread.so.0(+0xf630)[0x7fa7b8f19630]
[n075:76036] [ 3] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZN4Foam26surfaceInterpolationSchemeINS_6VectorIdEEE14dotInterpolateINS_14GeometricFieldIS2_NS_13fvsPatchFieldENS_11surfaceMeshEEEEENS_3tmpINS5_INS_12innerProductINT_10value_typeES2_E4typeES6_S7_EEEERKSB_RKNS5_IS2_NS_12fvPatchFieldENS_7volMeshEEERKNS9_INS5_IdS6_S7_EEEE+0x7ee)[0x7fa7bf57097e]
[n075:76036] [ 4] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZNK4Foam26surfaceInterpolationSchemeINS_6VectorIdEEE14dotInterpolateERKNS_14GeometricFieldIS2_NS_13fvsPatchFieldENS_11surfaceMeshEEERKNS4_IS2_NS_12fvPatchFieldENS_7volMeshEEE+0x72)[0x7fa7bf572732]
[n075:76036] [ 5] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(+0x165143d)[0x7fa7bfa5343d]
[n075:76036] [ 6] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSchemeIdNS_10SymmTensorIdEEE15gammaSnGradCorrERKNS_14GeometricFieldINS_6VectorIdEENS_13fvsPatchFieldENS_11surfaceMeshEEERKNS5_IdNS_12fvPatchFieldENS_7volMeshEEE+0x1e8)[0x7fa7bfa6c6d8]
[n075:76036] [ 7] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSchemeIdNS_10SymmTensorIdEEE12fvmLaplacianERKNS_14GeometricFieldIS3_NS_13fvsPatchFieldENS_11surfaceMeshEEERKNS5_IdNS_12fvPatchFieldENS_7volMeshEEE+0x1a8)[0x7fa7bfa6f778]
[n075:76036] [ 8] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so(_ZN4Foam2fv15laplacianSchemeIdNS_10SymmTensorIdEEE12fvmLaplacianERKNS_14GeometricFieldIS3_NS_12fvPatchFieldENS_7volMeshEEERKNS5_IdS6_S7_EE+0x60)[0x7fa7bfa69830]
[n075:76036] [ 9] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libfvOptions.so(_ZN4Foam3fvm9laplacianIdNS_10SymmTensorIdEEEENS_3tmpINS_8fvMatrixIT_EEEERKNS_14GeometricFieldIT0_NS_12fvPatchFieldENS_7volMeshEEERKNS9_IS6_SB_SC_EE+0x134)[0x7fa7be032fb4]
[n075:76036] [10] /home/018860125/OpenFOAM/OpenFOAM-common/platforms/linux64GccDPInt64Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels3SSGINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x36a)[0x7fa7bc4f055a]
[n075:76036] [11] pimpleFoam[0x429542]
[n075:76036] [12] /lib64/libc.so.6(__libc_start_main+0xf5)[0x7fa7b8b5e555]
[n075:76036] [13] pimpleFoam[0x42ce3c]
[n075:76036] *** End of error message ***

Here's my control dict file :

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2106                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     pimpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         30;

deltaT          1e-2;

writeControl    adjustableRunTime;

writeInterval   0.25;

purgeWrite      0;

writeFormat     ascii;

writePrecision   6;

writeCompression on;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

adjustTimeStep  no;

maxCo           0.6;

maxDeltaT       1.0e-3;



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

fvScheme :

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

ddtSchemes
{
    default         Euler; //CrankNicolson 0.7; 
}

gradSchemes
{
    default         leastSquares; //Gauss linear limited 0.5; 
}

divSchemes
{
    default         none;

    div(phi,U)       Gauss linearUpwindV grad(U);
    div(phi,h)       Gauss linearUpwind default;
    div(phi,e)       Gauss linearUpwind default;
    div(phi,k)       Gauss upwind default;
    div(phi,epsilon) Gauss upwind default;
    div(phi,R)       Gauss linearUpwind default;
    div(phiv,p)      Gauss linearUpwind default;
    div(phi,K)       Gauss linearUpwind default;
    div(phi,Ekp)     Gauss linearUpwind default;
    div(R)           Gauss linear;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div((nu*dev2(T(grad(U)))))    Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.5;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.5;
}

wallDist
{
    method meshWave;
}

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

fvSolution :

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2106                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p//p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-4;
        relTol          0.1;
    }

    pFinal//p_rghFinal
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-4;
        relTol          1e-8;
    }

    "(U|h|e|k|epsilon|R)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-3;
        relTol          0.1;
    }

    "(U|h|e|k|epsilon|R)Final"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-2;
        relTol          1e-8;
    }
}

PIMPLE
{
    momentumPredictor yes;
    nOuterCorrectors 50;
    nCorrectors     3;
    nNonOrthogonalCorrectors 2;
    pRefCell    0;
    pRefValue   1.0e5;
    turbOnFinalIterOnly no;

    outerCorrectorResidualControl
    {
        p
        {
            tolerance   1e-4;
            relTol      0.1;
        }

        U
        {
            tolerance   1e-4;
            relTol      0.1;
        }

        h
        {
            tolerance   1e-3;
            relTol      0.1;
        }

        T
        {
            tolerance   1e-3;
            relTol      0.01;
        }

        "(k|epsilon|R.*)"
        {
            tolerance   1e-10;
            relTol      0.01;
        }
    }
}

relaxationFactors
{
    equations
    {
        U               0.7;
        k               0.7;
        epsilon         0.3;
        R               0.7;
        p               0.3;
    }
}

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

I read on an other post that it can be due to the boundary condition. I checked it and everything seems correct to me, if the LRR case can run with those BCs I don't see why the SSG couldn't ?

[sr786]

Hi Haroun,

Have you managed to solve this issue. I am having similar problems with SSG model for transient solver (I'm using interFOAM with PIMPLE) and I initialised my solution from a converged and stable SST K-Omega simulation