[thintn222]

Hi,all Foamers

I'm running a simulation about bubble rising using MPPICInterFoam solver, coupling LPT to VOF method . Everything is going well when i test with a simple domain . But when i tried it with a real domain about 90000 hexahedra of cells, It freezing at time step 4.503 because of the value of Courant number is so high , about XX e+12 .

I tried to decrease the deltaT in the ControlDict to be 1e-8 even 1e-9 but it still not work.

But when i decrease the number of cell from 90000->11250 . It worked well.

So, what is the problem ? How i can run this simulation at 90000 cell .

I attached my fvSolution, fvSchemes and checkMesh belowed.

Any help would be highly appreciated ,

Thank you so much for reading !


checkMesh

Mesh stats
    points:           97061
    faces:            276900
    internal faces:   263100
    cells:            90000
    faces per cell:   6
    boundary patches: 3
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     90000
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology
    inlet               36       49       ok (non-closed singly connected)
    outlet              900      961      ok (non-closed singly connected)
    walls               12864    12936    ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-0.075 -0.075 0) (0.075 0.075 0.5)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-1.25066e-16 -1.31008e-17 1.48292e-18) OK.
    Max cell openness = 1.35525e-16 OK.
    Max aspect ratio = 1 OK.
    Minimum face area = 2.5e-05. Maximum face area = 2.5e-05.  Face area magnitu                                                                                                                                                             des OK.
    Min volume = 1.25e-07. Max volume = 1.25e-07.  Total volume = 0.01125.  Cell                                                                                                                                                              volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.11023e-13 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

fvSchemes

ddtSchemes
{
      default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(rhoPhi,U)               Gauss upwind;
    div(phi,alpha)              Gauss vanLeer;
    div(phirb,alpha)            Gauss linear;
    div(alphaRhoPhic,k)         Gauss upwind;
    div(alphaRhoPhic,epsilon)   Gauss upwind;
    div((((alphac*rho)*nuEff)*dev2(T(grad(U))))) Gauss linear;
    div(phiGByA,kinematicCloud:alpha) 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 2;
        cAlpha          1;

        MULESCorr       yes;
        nLimiterIter    2;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-7;
        relTol          0;
        maxIter         100;
    }

    pcorr
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-5;
        relTol          0;
    }

    p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-07;
        relTol          0.05;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

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

    kinematicCloud:alpha
    {
        solver          GAMG;
        tolerance       1e-06;
        relTol          0.1;
        smoother        GaussSeidel;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }
}

PIMPLE
{
    momentumPredictor   no;
    nOuterCorrectors    1;
    nCorrectors         2;
    nNonOrthogonalCorrectors 2;
    pRefCell      0;
    pRefValue     0;

    residualControl
    {
          "(p|U)"
          {
           tolerance 1e-4;
           relTol 0;
          }
    }
   // convergence 1e-08;
}

relaxationFactors
{
    fields
    {
       p  1;
    }
    equations
    {
       U 1;
       k 1;
       epsilon 0.5;
       R       0.5;
       nuTilda 0.5;
       omega   0.5;
    }
}

Once again , thank you so much for reading !