CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   LTSInterFoam Crashing (https://www.cfd-online.com/Forums/openfoam-solving/126421-ltsinterfoam-crashing.html)

JFM November 17, 2013 01:15

LTSInterFoam Crashing
 
Good day All

I am attempting to create a LTSInterFoam model as the interFoam models I currently have developed are taking an inordinate amount of time to finish (two months so far on 8/16 nodes with approximately 500,000 cells). Anyways the error output I am getting is:
Code:

Time = 1

Flow time scale min/max = 0.0223363, 1
Smoothed flow time scale min/max = 0.0223363, 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191047  Min(alpha1) = -93.8154  Max(alpha1) = 26.8419
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191234  Min(alpha1) = -2410  Max(alpha1) = 8705.97
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191417  Min(alpha1) = -808793  Max(alpha1) = 223717
[0] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2  in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[0] #8 
[0]  in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[0] #9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #10 
[0]  in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[Ubuntu:29885] *** Process received signal ***
[Ubuntu:29885] Signal: Floating point exception (8)
[Ubuntu:29885] Signal code:  (-6)
[Ubuntu:29885] Failing at address: 0x3e8000074bd
[Ubuntu:29885] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fc8032b04a0]
[Ubuntu:29885] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fc8032b0425]
[Ubuntu:29885] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fc8032b04a0]
[Ubuntu:29885] [ 3] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xd2) [0x7fc8043953c2]
[Ubuntu:29885] [ 4] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_RNS1_IS8_EESD_h+0xf1d) [0x7fc804397cfd]
[Ubuntu:29885] [ 5] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x4d3) [0x7fc804398c43]
[Ubuntu:29885] [ 6] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7fc805e86e17]
[Ubuntu:29885] [ 7] LTSInterFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x152) [0x4562a2]
[Ubuntu:29885] [ 8] LTSInterFoam() [0x42c93a]
[Ubuntu:29885] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fc80329b76d]
[Ubuntu:29885] [10] LTSInterFoam() [0x43145d]
[Ubuntu:29885] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 29885 on node Ubuntu exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

I have been using the following fvSolutions file:
Code:

solvers
{
        pcorr
        {
                solver                                PCG;

                preconditioner
                {
                        preconditioner                  GAMG;

                        smoother                                DICGaussSeidel;
                        agglomerator                        faceAreaPair;
                        mergeLevels                                1;
                        nCellsInCoarsestLevel        10;
                        cacheAgglomeration                true;

                        tolerance                                1e-6;
                        relTol                                        0;
                };

                tolerance                1e-6;
                relTol                        0;
        };

        p_rgh
        {
        solver                                  GAMG;

        smoother                                GaussSeidel;
        agglomerator                                  faceAreaPair;
        mergeLevels                            1;
        nCellsInCoarsestLevel                10;
        cacheAgglomeration                        true;

        tolerance                              1e-7;
        relTol                                  0.01;
    };

    p_rghFinal
    {
        $p_rgh;
        tolerance      1e-7;
        relTol          0;
    }

        U
        {
                solver                        smoothSolver;        //PBiCG;
                smoother                GaussSeidel;        //preconditioner        DILU;
                tolerance                1e-06;
                relTol                        0;
                nSweeps                        1;
        }

        UFinal
        {
                solver                        smoothSolver;        //$U;
                smoother                GaussSeidel;
                tolerance                1e-6;
                relTol                        0;
                nSweeps                        1;
        }

/*        "U.*"
        {
                solver                smoothSolver;
               
                smoother        GaussSeidel;
                nSweeps                1;
               
                tolerance        1e-8;
                relTol                0;
        };
*/

}

PIMPLE
{
    momentumPredictor        no;        // yes;

    nCorrectors                        2;        // 1;
    nNonOrthogonalCorrectors 0;

    nAlphaCorr      1;
    nAlphaSubCycles                3; //1;
    cAlpha          1;
    correctPhi                no;

    maxCo                                0.75;
    maxAlphaCo                        0.75;
    nAlphaSweepIter                1;

    rDeltaTSmoothingCoeff        0.1;
    rDeltaTDampingCoeff                1;
    maxDeltaT                                1;
}

relaxationFactors
{
    fields
    {
    }
    equations
    {
            "U.*"                                1;
    }
}


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

And fvSchemes file:
Code:

ddtSchemes
{
    default        localEuler rDeltaT;
}

gradSchemes
{
    default        Gauss linear;
}

divSchemes
{
    div(rho*phi,U)  Gauss linearUpwind grad(U);
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
//        div(phi,k)      Gauss upwind;
//        div(phi,omega)  Gauss upwind;
    div((muEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default        Gauss linear corrected;
}

interpolationSchemes
{
    default        linear;
}

snGradSchemes
{
    default        corrected;
}

fluxRequired
{
    default        no;
    p_rgh;
    pcorr;
    alpha1;
}


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

I would appreciate if someone could have a look at these and let me know what I am missing. I suspect this is an issue with alpha?

Regards
JFM


All times are GMT -4. The time now is 10:17.