CFD Online Discussion Forums - Reporting bugs for pimple loop in VOF solvers

Hi foamers,

I would like to report a bug with regard to the pimple implementation in compressibleInterFoam and similar solvers.

Code:

#include "alphaEqnsSubCycle.H"

            // correct interface on first PIMPLE corrector

            if (pimple.corr() == 1)

            {

                interface.correct();

            }

should be replaced with

Code:

// correct interface on first PIMPLE corrector

            if (pimple.corr() == 1)

            {

                #include "alphaEqnsSubCycle.H"

                interface.correct();

            }

This way the number of iteration for updating interface is consistent with number of iteration for solving alphaEqn. In this case, they are only updated in the first pimple loop. I have tested both implementation that the second one is much much more stable.

Hope this helps.