[lbordier]

Hi,

I am trying to use dynamicRefineFvMesh within a interDyMFoam computation to refine based on free surface (alpha.water field).

When my computation is in PISO mode, the refinement is done but then, pressure number of iteration increase and forces on body become instable such that finally computation crashes.

When my computation runs in PIMPLE mode, solver crashes just after the first solve of p_rgh :

PIMPLE: iteration 1
Selected 377308 cells for refinement out of 11056721.
Refined from 11056721 to 13697877 cells.
Selected 0 split points out of a possible 377308.
Execution time for mesh.update() = 18.03 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 9.78769e-05, No Iterations 74
time step continuity errors : sum local = 5.99048e-12, global = -1.6006e-14, cumulative = 3.108e-10
smoothSolver:  Solving for alpha.water, Initial residual = 1.8908e-08, Final residual = 8.08454e-11, No Iterations 1
Phase-1 volume fraction = 0.75008  Min(alpha.water) = -1.56267e-05  Max(alpha.water) = 1.16759
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.75008  Min(alpha.water) = -0.00340939  Max(alpha.water) = 1.16624
GAMG:  Solving for p_rgh, Initial residual = 0.0813079, Final residual = 0.000809642, No Iterations 7
[13] #0  Foam::error::printStack(Foam::Ostream&)[7] #0  Foam::error::printStack(Foam::Ostream&) addr2line failed
[13] #1  Foam::sigFpe::sigHandler(int) addr2line failed
[13] #2  ? in "/home/applis/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
[7] #1  Foam::sigFpe::sigHandler(int) addr2line failed
[13] #3   in "/home/applis/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
[7] #2  ??
[13] #4   in "/lib64/libc.so.6"
[7] #3  ?
[13] #5  ?
[13] #6  __libc_start_main? addr2line failed
[13] #7  Foam::UOPstream::write(float)

My fvSolution is as following :

// ************************************************************************* //
[bordier@visu: dtmb5415]$ cat system/fvSolution
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

        MULESCorr       yes;
        nLimiterIter    10;
        alphaApplyPrevCorr  yes;

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

  /*  p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-7;
        relTol          0.01;
    };*/

   p_rgh
    {
        solver          GAMG;

        smoother        DIC;

        tolerance       1e-7;
        relTol          0.01;
        cacheAgglomeration false;
    };

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "pcorr.*"
    {
        $p_rghFinal;
        relTol          0.0001;
    }


    "(U|k|omega).*"
    {
        solver          smoothSolver;

        smoother        symGaussSeidel;
        nSweeps         1;

        tolerance       1e-7;
        relTol          0.1;
        minIter         1;
    };
}

PIMPLE
{
    momentumPredictor   no;

    nOuterCorrectors    25;
    nCorrectors         2; //2
    nNonOrthogonalCorrectors 0;

    maxCo               15;
    maxAlphaCo          5;

    rDeltaTSmoothingCoeff 0.05;
    rDeltaTDampingCoeff 0.5;
    nAlphaSpreadIter    0;
    nAlphaSweepIter     0;
    maxDeltaT           0.005;

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

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

relaxationFactors
{
    fields
    {
        p_rgh 0.4;
        p_rghFinal 0.4;
    }
    equations
    {
        U 0.6;
        UFinal 0.6;
    }
}

cache
{
    grad(U);
}


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

My dynamicMeshDict is :

dynamicRefineFvMeshCoeffs
{
    // How often to refine
    refineInterval  5;

    // Field to be refinement on
    field           alpha.water;

    // Refine field inbetween lower..upper
    lowerRefineLevel 0.001;
    upperRefineLevel 0.999;

    // If value < unrefineLevel unrefine
    unrefineLevel   10;

    // Have slower than 2:1 refinement
    nBufferLayers   2;

    // Refine cells only up to maxRefinement levels
    maxRefinement   1;

    // Stop refinement if maxCells reached
    maxCells        30000000;

    // Flux field and corresponding velocity field. Fluxes on changed
    // faces get recalculated by interpolating the velocity. Use 'none'
    // on surfaceScalarFields that do not need to be reinterpolated.
    correctFluxes
    (
        (phi none)
        (nHatf none)
        (rhoPhi none)
        (alphaPhi none)
        (ghf none)
        ("(alphaPhi-flux(alpha.water))" none)
    );

    // Write the refinement level as a volScalarField
    dumpLevel       true;
}

Does anybody have an idea how to solve this issue ?
Any help is appreciated.