[cjwhodge]

Dear All,


Does anyone have any ideas why the water surface is following the deformation in the mesh rather than recalculating and following the physically correct water surface level?
I've attached the fvSchemes and fvSolutions. Any hints on what may be wrong would be much appreciated!

fvScheme:

ddtSchemes
{
    default         Euler;
// default 	  CrankNicolson 0.9; 
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    div(rhoPhi,U)   Gauss limitedLinearV 1; 
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;

}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p_rgh;
    pcorr;
    alpha.water;
}
// ************************************************************************* //

fvSolution:

solvers
{
    "alpha.water.*"
    {
        nAlphaCorr      2;
        nAlphaSubCycles 3;               
        cAlpha          1.0
  
        icAlpha         0; 
        MULESCorr       yes; 
        nLimiterIter    8;    
        alphaApplyPrevCorr  true;

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

    "pcorr.*"
    {
        solver          PCG; 
        preconditioner
        {
            preconditioner  GAMG; 
            tolerance       1e-7;
            relTol          0;
            smoother        DICGaussSeidel;
      

            nPreSweeps      0;
            nPostSweeps     2;
            nFinestSweeps   2;

            cacheAgglomeration true;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1; 
           
        }

        tolerance       1e-7;
        relTol          0;
        maxIter         100;
    }

    p_rgh
    {
        solver          GAMG;  
        tolerance       1e-7;
        relTol          0.05;
        smoother        DIC;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    p_rghFinal
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance       1e-8;
            relTol          0;

            smoother        DICGaussSeidel;
            nPreSweeps      2;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration true;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }

        tolerance       1e-8;
        relTol          0;
        maxIter         20;
    }

    "(U|k|epsilon)"  
    {
        solver          PBiCG;          
        preconditioner  DILU;
        smoother        GaussSeidel;
        tolerance       1e-06;
        relTol          0.1;
        nSweeps         1;
    }

    "(U|k|epsilon)Final"
    {
        solver          smoothSolver; 
        smoother        GaussSeidel;
        tolerance       1e-07;
        relTol          0;
        nSweeps         1;
    }
}


PIMPLE
{
    momentumPredictor   no;  
    nCorrectors         3;	 
    nOuterCorrectors    150; 

 
    nNonOrthogonalCorrectors 1;

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

        p_rgh
        {
            tolerance  5e-4;
            relTol     0;
        }
     }

    correctPhi          yes;
    moveMeshOuterCorrectors yes;
}


relaxationFactors
{
    fields
    {
        p_rgh      0.7;              
        p_rghFinal   1;
        // 0.3 ~ 22 iterations 
    }

    equations
    {
        "U|k|epsilon"        0.3;   
        "(U|k|epsilon)Final"   1;
    }
}


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