In the Eulerian two-fluid model for two-phase flows, the standard formulation for the momentum equations becomes singular as the volume fractions tend to zero. Some papers suggest that resulting velocities can fluctuate widely in the regions where volume fractions tend to vanish.

How do you deal with the problem of singularity in your code?

Or could you give me some effective suggestions for the problem?

Thanks!

As far as I know, the VOF can often lead to such a problem due to the numerical diffusion. See the paper written by R.Abrall.(?).

I am usin level set method to distinguish between fluids. It inherently has a function to keep the level set reasonable over all the time.

There are a number of way for multi-fuids flow simulation. I am not sure which one is best.

To jinwon

I think the level set method is perfer to the VOF method and is more prevalent.

Thanks for commending the paper of R.Abrall.

Do you give me some more information about the paper of R.Abrall, for example the title of the paper?

I really hope that you can give me some suggestion for my problem.

Thanks!

Now, I don't have the list. But you can easily search it through googling. Just type as

"compressible multifluid"