Comparison of the Interface-Capturing Behavior between LSM and MCLS

Abstract:

Numerical simulation of the moving interface is widely investigated in the fields of two-phase flow, crystal growth, molten metal casting and so on. To simulate accurately and quickly motion state of the moving interface, two types of interface-capturing method, LSM (level set method) and MCLS (mass-conserving LSM coupled with VOF (volume of fluid)), were comparatively analysed and studied. Two kinds of moving interface, rotation motion of Zaleska slotted disk interface and shearing deformation of the circle interface, were simulated. The level set advection equation was solved with IATVD (the characteristic integral-averaging finite volume method), RK-WENO (the third order Runge-Kutta method with WENO (weighted essentially non-oscillatory)), or the low order method, respectively. The level set function was re-initialized with modified Godunov method. Computation was implemented in different size mesh. When it is independent with mesh size, results show that the high order method, IATVD and RK-WENO, can get higher computing accuracy in LSM without topological geometry change. But computing time cost of RK-WENO is much greater than that of IATVD. For the low order method, LSM captures the moving interface in bad, however, MCLS’s numerical result is very good. But computing time cost of MCLS is much and much greater than that of LSM. For interface-capturing with topological geometry change, these two high order methods all can not capture the moving interface in good in LSM. Nevertheless, the resulting interfaces all are same as original interface in MCLS. These indicate that an appropriate numerical method and interface-capturing method can be chosen to accurately and quickly capture the related moving interface.

• Periodical: Key Engineering Materials (Volume 730)
• Pages: 533-540
• DOI: https://doi.org/10.4028/www.scientific.net/KEM.730.533
• Online since: February 2017
• Authors: Zhi Lin Liu, Yong Huang*
• Keywords: IATVD, LSM, MCLS, Molten Metal Casting, Moving Interface, RK-WENO

References

[1] S. Osher, J. A. Sethian, Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys. 79 (1988) 12-49.

DOI: 10.1016/0021-9991(88)90002-2

[2] C. W. Hirt, B. D. Nichols, Volume of fluid (VOF) method for the dynamics of free boundaries, J. Comput. Phys. 39 (1981) 201-225.

DOI: 10.1016/0021-9991(81)90145-5

[3] J. A. Sethian, A fast marching level set method for monotonically advancing fronts, Proc. Nat. Acad. Sci. 93 (1996) 1591-1595.

DOI: 10.1073/pnas.93.4.1591

[4] D. Adalsteinsson, J. A. Sethian, A fast level set method for propagating interfaces, J. Comput. Phys. 118 (1995) 269-277.

DOI: 10.1006/jcph.1995.1098

[5] D. Peng, et al., A PDE-based fast local level set method, J. Comput. Phys. 155 (1999) 410-438.

[6] M. Sussman, et al., An adaptive level set approach for incompressible two-phase flows, J. Comput. Phys. 148 (1999) 81-124.

[7] S. Osher, C. Shu, High-order essentially nonoscillatory schemes for Hamilton-Jacobi equations, J. Num. Anal. 28 (1991) 907-922.

DOI: 10.1137/0728049

[8] X. Liu, S. Osher, T. Chan, Weighted essentially non-oscillatory schemes, J. Comput. Phys. 115 (1994) 200-212.

DOI: 10.1006/jcph.1994.1187

[9] G. Jiang, D. Peng, Weighted ENO schemes for Hamilton-Jacobi equations, J. Sci. Comput. 21 (2000) 2126-2143.

DOI: 10.1137/s106482759732455x

[10] H. X. Li, et al., Numerical simulation of interface movement in gas-liquid two-phase flows with level set method, Nucl. Pow. Eng. 03 (2005) 242-248.

[11] M. Sussman, E. G. Puckett, A coupled level set and volume-of-fluid method for computing 3D and axisymmetric incompressible two-phase flows, J. Comput. Phys. 162 (2000) 301-337.

DOI: 10.1006/jcph.2000.6537

[12] G. Son, N. Hur, A coupled level set and volume-of-fluid method for the buoyancy-driven motion of fluid particles, Num. Heat Transf.: Part B: Fundament. 42 (2002) 523-542.

DOI: 10.1080/10407790260444804

[13] G. Son, Efficient implementation of a coupled level-set and volume-of-fluid method for three-dimensional incompressible two-phase flows, Num. Heat Transf.: Part B: Fundament. 43 (2003) 549-565.

DOI: 10.1080/713836317

[14] Q. Yang, et al., An investigation on coupled level-set and VOF algorithm in interface capturing, Aeronaut. Comput. Tech. 04 (2012) 14-19.

[15] S. P. Van der pijl, A. Segal, C. Vuik, A mass-conserving Level-Set method for modeling of multi-phase flows, Int. J. Num. Meth. fluids, 47 (2005) 339-361.

DOI: 10.1002/fld.817

[16] M. Herrmann, Refined level set grid method for tracking interfaces, Annual Research Briefs, Center for Turbulence Research, NASA Ames/Stanford University, 2005, pp.3-18.

[17] R. X. Liu, et al., Tracking and reconstruction methods for moving-interfaces, Appl. Math. Mech. 25 (2004) 279-290.