Skip to main content Accessibility help

A Numerical Method on Eulerian Grids for Two-Phase Compressible Flow

  • Yonghui Guo (a1), Ruo Li (a2) and Chengbao Yao (a3)


We develop a numerical method to simulate a two-phase compressible flow with sharp phase interface on Eulerian grids. The scheme makes use of a levelset to depict the phase interface numerically. The overall scheme is basically a finite volume scheme. By approximately solving a two-phase Riemann problem on the phase interface, the normal phase interface velocity and the pressure are obtained, which is used to update the phase interface and calculate the numerical flux between the flows of two different phases. We adopt an aggregation algorithm to build cell patches around the phase interface to remove the numerical instability due to the breakdown of the CFL constraint by the cell fragments given by the phase interface depicted using the levelset function. The proposed scheme can handle problems with tangential sliping on the phase interface, topological change of the phase interface and extreme contrast in material parameters in a natural way. Though the perfect conservation of the mass, momentum and energy in global is not achieved, it can be quantitatively identified in what extent the global conservation is spoiled. Some numerical examples are presented to validate the numerical method developed.


Corresponding author

*Corresponding author. Email: (Y. H. Guo), (R. Li), (C. B. Yao)


Hide All
[1]Brackbill, J. U., Kothe, D. B., and Ruppel, H. M., FLIP: a low-dissipation, particle-in-cell method for fluid flow, Comput. Phys. Commun., 48(1) 91988), pp. 2538.
[2]Di, Yana, Li, Ruo, Tang, Tao, and Zhang, Pingwen, Level set calculations for incompressible two-phase flows on a dynamically adaptive grid, J. Sci. Comput., 31(1-2) (2007), pp. 7598.
[3]Enright, Douglas, Fedkiw, Ronald, Ferziger, Joel, and Mitchell, Ian, A hybrid particle level set method for improved interface capturing, J. Comput. Phys., 183(1) (2002), pp. 83116.
[4]Hieber, Simone E and Koumoutsakos, Petros, A Lagrangian particle level set method, J. Comput. Phys., 210(1) (2005), pp. 342367.
[5]Koshizuka, Seiichi, Nobe, Atsushi, and Oka, Yoshiaki, Numerical analysis of breaking waves using the moving particle semi-implicit method, Int. J. Numer. Methods Fluids, 26(7) (1998), pp. 751769.
[6]Koumoutsakos, Petros, Multiscale flow simulations using particles, Ann. Rev. Fluid Mech., 37 (2005), pp. 457487.
[7]Monaghan, Joe J., Smoothed particle hydrodynamics, Reports Progress Phys., 68(8) (2005), pp. 1703.
[8]Noh, William F.and Woodward, Paul(Simple Line Interface Calculation). In Proceedings of the Fifth International Conference on Numerical Methods in Fluid Dynamics June 28–July 2, 1976, Twente University, Enschede, pp. 330340, Springer, 1976
[9]Osher, Stanley and Fedkiw, Ronald P., Level set methods: an overview and some recent results, J. Comput. Phys., 169(2) (2001), pp. 463502.
[10]Qiao, Dengjiang, An Introduction to Nuclear Explosion Physics, National Defence Industry Press (in Chinese), 2003.
[11]Rider, William J and Kothe, Douglas B, Stretching and tearing interface tracking methods, in AIAA Computational Fluid Dynamics Conference, 12th, and Open Forum, San Diego, CA, pp. 806816, 1995.
[12]Sbalzarini, I. F., Walther, Jens Honore, Bergdorf, M., Hieber, S. E., Kotsalis, E. M., and Koumoutsakos, P., PPM-a highly efficient parallel particle-mesh libibrary for the simulation of continuum systems, J. Comput. Phys., 215(2) (2006), pp. 566588.
[13]Scardovelli, Ruben and Zaleski, Stéphane, Direct numerical simulation of free-surface and interfacialflow, Ann. Rev. Fluid Mech., 31(1) (1999), pp. 567603.
[14]Sethian, James A., Evolution, implementation, and application of level set and fast marching methods for advancing fronts, J. Comput. Phys., 169(2) (2001), pp. 503555.
[15]Sussman, Mark, A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles, J. Comput. Phys., 187(1) (2003), pp. 110136.
[16]Sussman, Mark, Smereka, Peter, and Osher, Stanley, A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114(1) (1994), pp. 146159.
[17]Tryggvason, Grétar, Bunner, Bernard, Esmaeeli, Asgsghar, Juric, Damir, Al.-Rawahi, N., Tauber, W., Han, J., Nas, S., and Jan, Y. J., a front-tracking method for the computations of multiphase flow, J. Comput. Phys., 169(2) (2001), pp. 708759.
[18]Unverdi, S. O. and Tryggvason, G., A front-tracking method for viscous, incompressible, multi-fluid flows, J. Comput. Phys., 100(1) (1992), pp. 2537.
[19]Walther, Jens Honore, Werder, T., Jaffe, R. L., and Koumoutsakos, P., Hydrodynamic properties of carbon canotubes, Phys. Rev. E, 69(6) (2004), pp. 062201.
[20]Han Young, Yoon, Koshizuka, Seiichi, and Oka, Yoshiaki, Direct calculation of bubble growth, departure, and rise in nucleate pool boiling, Int. J. Multiphase Flow, 27(2) (2001), pp. 277298.


MSC classification

Related content

Powered by UNSILO

A Numerical Method on Eulerian Grids for Two-Phase Compressible Flow

  • Yonghui Guo (a1), Ruo Li (a2) and Chengbao Yao (a3)


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.