[ginop]

I am a bit puzzled with this case, and need your opinion.


Here, I am running a transient, incompressible multi-phase simulation with interFoam (v2012). I've already run the same geometry/mesh under different conditions, with success.
I suspect that the issue may be due to the mesh, though the same grid worked under a different oscillation frequency (same amplitude)



This particular case crashes after many timesteps, about 3000. The case is dynamic, using

angularOscillatingDisplacement

Attached is an extract of my log file, including a few timesteps before crashing.
My finite volume files are included below.

Any help is much appreciated.


fvSchemes

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(U)         cellLimited Gauss linear 1;
}

divSchemes
{
    div(rhoPhi,U)   Gauss linearUpwindV grad(U); //Gauss vanLeerV;
    div(phi,U)      Gauss linearUpwindV grad(U);
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss vanLeer; //interfaceCompression;
    div(phi,k)      bounded Gauss linearUpwind  grad(k);
    div(phi,omega)  bounded Gauss linearUpwind  grad(omega);
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
        method meshWave;
 }

fvSolution

solvers
{
    "alpha.water.*"
    {
        nAlphaCorr      3;
        nAlphaSubCycles 1;
        cAlpha          1;
        icAlpha         0;

        MULESCorr       yes;
        nLimiterIter    15;
        alphaApplyPrevCorr  yes;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-10;
        relTol          0;
        minIter         1;
    }

    "pcorr.*"
    {
        solver          GAMG;
        smoother        DIC;
        tolerance       1e-05;
        relTol          0;
    }

    p_rgh
    {
        solver          GAMG;
        smoother        DICGaussSeidel;
        tolerance       1e-09;
        relTol          0.01;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "(U|k|omega)"
    {
        solver          smoothSolver;
        smoother        DILUGaussSeidel;
        tolerance       1e-06;
        relTol          0.01;
        minIter         1;
        maxIter         20;
    }

    "(U|k|omega)Final"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-08;
        relTol          0;
        maxIter         50;
    }

    "cellDisplacement.*"
    {
        solver          GAMG;
        tolerance       1e-8;
        relTol          0;
        smoother        GaussSeidel;
    }

    Phi
    {
        solver          smoothSolver;
        tolerance       1e-08;
        relTol          0;        
        smoother        GaussSeidel;
    }
}

potentialFlow 
{ 
    nNonOrthogonalCorrectors 10; 
    PhiRefCell 0; 
    PhiRefValue 0; 
}

PIMPLE
{
    nNonOrthogonalCorrectors 2;
    nCorrectors         1;
    nOuterCorrectors    25;
    correctPhi          true;
    moveMeshOuterCorrectors true;
    turbOnFinalIterOnly yes;
    pRefCell            0;
    pRefValue           0;

    residualControl
    {
        U
        {
            tolerance  1e-4;
            relTol      0;
        }
        p_rgh
        {
            tolerance  1e-5;
            relTol      0;
        }
    }
}

relaxationFactors
{
    fields
    {
        p_rgh       0.5;
        p_rghFinal  1;
    }
    equations
    {
        ".*"        1;
        U           0.7;
        UFinal      1;
    }
}

cache
{
    grad(U);
}