[Michael@UW]

I use overPimpleDyMFoam in OpenFOAM v2006 to simulate the oscillating rotation of two valves.

The cellTypes of some cells are not switched back to 'calculated' but keep 'hole' after the valves leave these cells (see figure attached). It seems the 'cellTypes' are updated correctly during the open process but the issue occurs during the close process.

I use adjustableRunTime and adjustTimeStep to gurantee maxCo to be smaller than 2.0.

I also attached the dictionary files. Hope someone can help me figure out how to solve the problem.

Thanks in advance!

Michael


controlDict

application     overPimpleDyMFoam;


startFrom		latestTime;

deltaT          0.0001028;

startTime		0;

endTime		0.073502;

writeControl    adjustableRunTime;
writeInterval   0.0001028;

purgeWrite      0;

writeFormat     binary;

writePrecision  10;

writeCompression off;

timeFormat      general;

timePrecision   10;

runTimeModifiable yes;
adjustTimeStep	yes;
maxCo 2.0;
maxDeltaT 0.0001028;

fvSchemes

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss upwind;

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

laplacianSchemes
{
    default         Gauss linear corrected;
    laplacian(diffusivity,cellDisplacement)  Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

oversetInterpolation
{
    method          cellVolumeWeight;
}

fluxRequired
{
    default         no;
    pcorr           ;
    p               ;
}

oversetInterpolationSuppressed
{
    grad(p);
    surfaceIntegrate(phiHbyA);
	//cellMotionU;
	//pointMotionU;
}

fvSolution

solvers
{
    "(cellMotionU|cellDisplacement.*)"
    {
        solver          PCG;
        preconditioner  DIC;

        tolerance       1e-9;
        relTol          0;
    }

    p
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-11;
        relTol          0.01;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    pcorr
    {
        $p;
        solver          PCG;
        preconditioner  DIC;
    }

    pcorrFinal
    {
        $pcorr;
        relTol          0;
    }

    "(U|k|epsilon)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-6;
        relTol          0;
    }

    "(U|k|epsilon)Final"
    {
        $U;
        tolerance       1e-6;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor   false;
    correctPhi          false;
    nOuterCorrectors    1;
    nCorrectors         3;
    nNonOrthogonalCorrectors 1;

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

    ddtCorr             true;

    pRefPoint           (-0.000372638 -0.00338073 0);
    pRefValue           0.0;

	checkMeshCourantNo 		yes;
	moveMeshOuterCorrectors true;
}

relaxationFactors
{
    fields
    {
        p      0.3;
        pFinal   1;
    }
    equations
    {
        "U|k|epsilon"     0.3;
        "(U|k|epsilon)Final"   1;
    }
}