CFD Online Discussion Forums - Open Channel Flow - Interfoam

Hello,

I am running a simple open channel flow simulation using interFoam. At the outlet, the max phase 1 volume fraction is greater than 1. After some time of running (+1 min flow time) the solution goes unstable and begins to diverge. Any ideas on what causes that and if there is a way to keep the alpha.water bounded between 0 and 1?

Boundary Conditions and fv files are below:

U

Code:

boundaryField

{

    inlet

    {

        type            variableHeightFlowRateInletVelocity;

        flowRate        constant 100;

        alpha           alpha.water;

        value           uniform (0 0 0);

    }

    outlet

    {

        type            pressureInletOutletVelocity;

        value           uniform (0 0 0);

    }

    concrete

    {

        type            fixedValue;

        value           uniform (0 0 0);

    }

    bed

    {

        type            fixedValue;

        value           uniform (0 0 0);

    }

    wall

    {

        type            slip;

    }

    top

    {

        type            pressureInletOutletVelocity;

        value           uniform (0 0 0);

    }

}

p_rgh:

Code:

boundaryField

{

    inlet       {type            zeroGradient;}



    outlet

    {

        type            fixedValue;

        value           uniform 2673221;

    }



    wall        {type            zeroGradient;}

    bed         {type            zeroGradient;}

    concrete    {type            zeroGradient;}



    top

    {

        type            totalPressure;

        gamma           1;

        p0              uniform 0;

        value           uniform 0;

    }



    defaultFaces

    {

        type            empty;

    }

}

alpha.water:

Code:

boundaryField

{

    inlet

    {

        type            variableHeightFlowRate;

        lowerBound      0;

        upperBound      1;

        value           uniform 0;

    }



    outlet

    {

        type            inletOutlet;

        inletValue      uniform 0;

        value           uniform 1;

    }



    wall

    {

        type            zeroGradient;

    }



    bed

    {

        type            zeroGradient;

    }



    concrete

    {

        type            zeroGradient;

    }





    top

    {

        type            inletOutlet;

        inletValue      uniform 0;

        value           uniform 0;

    }



    defaultFaces

    {

        type            empty;

    }

}

fvSolution:

Code:

solvers

{

    "alpha.water.*"

    {

        nAlphaCorr      2;

        nAlphaSubCycles 2;

        alphaOuterCorrectors yes;

        cAlpha          1;



        MULESCorr       yes;

        nLimiterIter    3;



        solver          smoothSolver;

        smoother        symGaussSeidel;

        tolerance       1e-8;

        relTol          0;

    }



    pcorr GAMG

    {

        tolerance        1e-7;

        relTol           0.0;



        smoother         DIC;//GaussSeidel;

        nPreSweeps       0;

        nPostSweeps      2;

        nFinestSweeps    2;



        cacheAgglomeration true;

        nCellsInCoarsestLevel 10;

        agglomerator     faceAreaPair;

        mergeLevels      1;

        maxIter 10;

    };



    p_rgh GAMG

    {

        tolerance        1e-7;

        relTol           0.0;



        smoother         DIC;//GaussSeidel;

        nPreSweeps       0;

        nPostSweeps      2;

        nFinestSweeps    2;



        cacheAgglomeration true;

        nCellsInCoarsestLevel 10;

        agglomerator     faceAreaPair;

        mergeLevels      1;

    };

    p_rghFinal GAMG

    {

        tolerance        1e-8;

        relTol           0.0;



        smoother         DIC;//GaussSeidel;

        nPreSweeps       0;

        nPostSweeps      2;

        nFinestSweeps    2;



        cacheAgglomeration true;

        nCellsInCoarsestLevel 10;

        agglomerator     faceAreaPair;

        mergeLevels      1;

    };







    "(U|k|omega).*"

    {

        solver          smoothSolver;

        smoother        symGaussSeidel;

        nSweeps         1;

        tolerance       1e-6;

        relTol          0.1;

    };

}





PIMPLE

{

    momentumPredictor no;

    nOuterCorrectors 1;

    nCorrectors     1;

    nNonOrthogonalCorrectors 2;

}



relaxationFactors

{

    fields

    {

    }

    equations

    {

        ".*" 0.1;

    }

}

fvSchemes:

Code:

ddtSchemes

{

    default         Euler;

}



gradSchemes

{

  default cellMDLimited Gauss linear 1;

  grad(k) faceLimited leastSquares 1;

  grad(omega) faceLimited leastSquares 1;

}



divSchemes

{

    div(rhoPhi,U)  Gauss upwind;

    div(phi,alpha)  Gauss vanLeer;

    div(phirb,alpha) Gauss interfaceCompression 1;

    div(phi,p_rgh)  Gauss upwind;

    div(phi,k)      Gauss upwind;//vanLeer;

    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;

    div(phi,omega) Gauss upwind;

}



laplacianSchemes

{

    default         Gauss linear uncorrected;

}



interpolationSchemes

{

    default         linear;

}



snGradSchemes

{

    default         uncorrected;

}



fluxRequired

{

    default         no;

    p_rgh;

    pcorr;

    gamma;

}



wallDist

{

    method meshWave;

}