Hostname: page-component-74656f9884-g96k6 Total loading time: 0 Render date: 2026-11-11T02:53:46.687Z Has data issue: false hasContentIssue false
  • English
  • Français

Simulation of Incompressible Free Surface Flow Using the Volume Preserving Level Set Method

Part of: Basic methods in fluid mechanics Equations of mathematical physics and other areas of application Incompressible viscous fluids

Published online by Cambridge University Press:  15 October 2015

Ching-Hao Yu  and
Tony Wen-Hann Sheu
Ching-Hao Yu
Affiliation:
Department of Ocean Science and Engineering, Zhejiang University, Yuhangtang Road, Hangzhou, Zhejiang, P.R.China
Tony Wen-Hann Sheu*
Affiliation:
Department of Engineering Science and Ocean Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei, Taiwan Institute of Applied Mathematical Sciences, National Taiwan University, Taiwan Center of Advanced Study in Theoretical Sciences (CASTS), National Taiwan University, Taiwan
*
*Corresponding author. Email address: twhsheu@ntu.edu.tw (T. W.-H. Sheu)

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the 'Save PDF' action button.

This study aims to develop a numerical scheme in collocated Cartesian grids to solve the level set equation together with the incompressible two-phase flow equations. A seventh-order accurate upwinding combined compact difference (UCCD7) scheme has been developed for the approximation of the first-order spatial derivative terms shown in the level set equation. Developed scheme has a higher accuracy with a three-point grid stencil to minimize phase error. To preserve the mass of each phase all the time, the temporal derivative term in the level set equation is approximated by the sixth-order accurate symplectic Runge-Kutta (SRK6) scheme. All the simulated results for the dam-break, Rayleigh-Taylor instability, bubble rising, two-bubble merging, and milkcrown problems in two and three dimensions agree well with the available numerical or experimental results.

Keywords

Level set equationupwinding combined compact difference schemethree-point grid stencilminimize phase errorsymplectic Runge-Kutta

MSC classification

Secondary: 35Q30: Navier-Stokes equations 76D05: Navier-Stokes equations 76D27: Other free-boundary flows; Hele-Shaw flows 76M20: Finite difference methods

Information

Type
Research Article
Information
Communications in Computational Physics , Volume 18 , Issue 4 , October 2015 , pp. 931 - 956
Copyright
Copyright © Global-Science Press 2015 

References

[1]Hirt, C. W., Amsden, A. A. and Cook, J. L., An arbitrary Lagrangian-Eulerian computing method for all flow speeds, J. Comput. Phys., 135(2) (1997), 203–216.Google Scholar
[2]Boultone-Stone, J. M. and Blake, J. R., Gas bubbles bursting at a free surface, J. Fluid Mech., 254 (1993), 437–466.Google Scholar
[3]Best, J. P., The formation of toroidal bubbles upon the collapse of transient cavities, J. Fluid Mech., 251 (1993), 79–107.Google Scholar
[4]Zhang, Y. L., Yeo, K. S., Khoo, B. C. and Wang, C., 3D jet impact and toroidal bubbles, J. Comput. Phys., 166 (2001), 336–360.Google Scholar
[5]Badalassi, V. E., Ceniceros, H. D. and Banerjee, S., Computation of multiphase systems with phase field models, J. Comput. Phys., 190 (2003), 371–397.Google Scholar
[6]Kim, J., A continuous surface tension force formulation for diffuse-interface models, J. Comput. Phys., 204 (2005), 784–804.Google Scholar
[7]Kim, J., Phase-field models for multi-component fluid flows, Commun. Comput. Phys., 12 (2012), 613–661.CrossRefGoogle Scholar
[8]Ding, H., Spelt, P. D. M. and Shu, C., Diffuse interface model for incompressible two-phase flows with large density ratios, J. Comput. Phys., 226 (2007), 2078–2095.CrossRefGoogle Scholar
[9]Hirt, C. W. and Nichols, B. D., Volume of fluid (VOF) method for the dynamics of free boundaries, J. Comput. Phys., 39 (1981), 201–225.Google Scholar
[10]Chen, S., Johnson, D. B. and Raad, P. E., Velocity boundary conditions for the simulation of free surface fluid flow, J. Comput. Phys., 116 (1995), 262–276.Google Scholar
[11]Caiden, R., Fedkiw, R. and Anderson, C., A numerical method for two phase flow consisting of separate compressible and incompressible regions, J. Comput. Phys., 166 (2001), 1–27.Google Scholar
[12]Li, B. and Shopple, J., An interface-fitted finite element level set method with application to solidification and solvation, Commun. Comput. Phys., 10(1) (2011), 32–56.Google Scholar
[13]Osher, S. and Sethian, J. A., Fronts propagating with curvature dependent speed: algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys., 79 (1988), 12–49.Google Scholar
[14]Sussman, M. and Fatemi, E., An efficient interface preserving level set redistancing algorithm and its application to interfacial incompressible fluid flow, SIAM J. Sci. Comput., 20 (1999), 1165–1191.CrossRefGoogle Scholar
[15]Sussman, M., Smereka, P. and Osher, S., A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114 (1994), 146–159.Google Scholar
[16]Marchandise, E., Remacle, J. F. and Chevaugeon, N., A quadrature-free discontinuous Galerkin method for the level set equation, J. Comput. Phys., 212 (2006), 338–357.Google Scholar
[17]Enright, D., Fedkiw, R., Ferziger, J. and Mitchell, I., A hybrid particle level set method for improved interface capturing, J. Comput. Phys., 183 (2002), 83–116.Google Scholar
[18]Sussman, M. and Puckett, E. G., A coupled level set and volume-of-fluid method for computing 3D and axisymmetric incompressible two-phase flow, J. Comput. Phys., 162 (2000), 301–337.Google Scholar
[19]Losasso, F., Fedkiw, R. and Osher, S., Spatially adaptive techniques for level set methods and incompressible flow, Comput. Fluids, 35 (2006), 995–1010.Google Scholar
[20]Olsson, E. and Kreiss, G., A conservative level set method for two phase flow, J. Comput. Phys., 210 (2005), 225–246.Google Scholar
[21]Osher, S. and Fedkiw, R., Level set methods: Anoverview and some recent results, J. Comput. Phys., 169 (2001), 463–502.Google Scholar
[22]Osher, S. and Fedkiw, R., Level Set Methods and Dynamic Implicit Surfaces, Springer-Verlag: Berlin, 2003.Google Scholar
[23]Sussman, M., Fatemi, E., Smereka, P. and Osher, S., An improved level set method for incompressible two-phase flows, Comput. Fluids, 27 (1998), 663–680.Google Scholar
[24]Sethian, J., Level Set Methods and Fast Marching Methods, Cambridge University Press, Cambridge, U. K, 1999.Google Scholar
[25]Sethian, J. A., Evolution, implementation, and application of level set and fast marching methods for advancing fronts, J. Comput. Phys., 169 (2001), 503–555.Google Scholar
[26]Morrison, P. J., Hamiltonian description of the ideal fluid, Rev. Mod. Phys., 70(2) (1998), 467–521.Google Scholar
[27]Shepherd, T. G., A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids, J. Fluid Mech., 213 (1990), 573–587.CrossRefGoogle Scholar
[28]McLachlan, R. I., Area preservation in computational fluid dynamics, Phys. Letter. A, 264 (1999), 36–44.Google Scholar
[29]Tam, C. K. W. and Webb, J. C., Dispersion-relation-preserving finite difference schemes for computational acoustics, J. Comput. Phys., 107 (1993), 262–281.Google Scholar
[30]Oevel, W. and Sofroniou, M., Symplectic Runge-Kutta schemes II: classification of symmetric method, Univ. of Paderborn, Germany, Preprint, 1997.Google Scholar
[31]Jiang, G.-S. and Peng, D., Weighted ENO schemes for Hamilton-Jacobi equations, SIAM J. Sci. Comput., 21 (2000), 2126–2143.Google Scholar
[32]Shu, C. W. and Osher, S., Efficient implementation of essentially non-oscillatory shock-capturing schemes, J. Comput. Phys., 77 (1988), 439–471.Google Scholar
[33]Bell, J. B., Colella, P., Glaz, H. M., A second-order projection method for the incompressible Navier-Stokes equations, J. Comput. Phys., 85 (1989), 257–283.Google Scholar
[34]Sheu, T. W. H. and Chiu, P. H., A divergence-free-condition compensated method for incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Engrg., 196 (2007), 4479–4494.Google Scholar
[35]Nihei, T. and Ishii, K., A fast solver of the shallow water equations on a sphere using a combined compact difference scheme, J. Comput. Phys., 187 (2003), 639–659.Google Scholar
[36]Kim, J. W. and Lee, D. J., Optimized compact finite difference schemes with maximum resolution, AIAA J., 34(5) (1996), 887–893.Google Scholar
[37]Chu, P. C. and Fan, C., A three-point combined compact difference scheme. J. Comput. Phys., 140 (1998), 370–399.CrossRefGoogle Scholar
[38]Vichnevetsky, R. and Bowles, J. B., Fourier Analysis of Numerical Approximations of Hyperbolic Equations, SIAM, Philadelphia, 1982.Google Scholar
[39]De, A. K. and Eswaran, V., Analysis of a new high resolution upwind compact scheme, J. Comput. Phys., 218 (2006), 398–416.Google Scholar
[40]Zalesak, S. T., Fully multidimensional flux-corrected transport algorithms for fluids, J. Comput. Phys., 31 (1979), 335–362.CrossRefGoogle Scholar
[41]Wang, Z., Yang, J., Koo, B. and Stern, F., A coupled level set and volume-of-fluid method for sharp interface simulation of plunging breaking waves, Int. J. Multiphase Flows, 35 (2009), 227–246.Google Scholar
[42]Koshizuka, S., Tamako, H. and Oka, Y., A particle method for incompressible viscous flow with fluid fragmentation, Comput. Fluid Mech. J., 113 (1995), 134–147.Google Scholar
[43]Kees, C. E., Akkerman, I., Farthing, M. W. and Bazilevs, Y., A conservative level set method suitable for variable-order approximations and unstructured meshes, J. Comput. Phys., 230 (2011), 4536–4558.Google Scholar
[44]Elias, R. N. and Coutinho, A. L. G. A., Stabilized edge-based finite element simulation of 956 free-surface flows, Int. J. Numer. Meth. Fluids, 54 (2007), 965–993.Google Scholar
[45]Tryggvason, G., Numerical simulations of the Rayleigh-Taylor instability, J. Comput. Phys., 75 (1988), 253–382.Google Scholar
[46]Guermond, J. -L. and Quartapelle, L., A projection FEM for variable density incompressible flows, J. Comput. Phys., 165 (2000), 167–188.Google Scholar
[47]Brereton, G. and Korotney, D., Coaxial and oblique coalescence of two rising bubbles, in: Sahin, I., Tryggvason, G. (Eds.), Dynamics of Bubbles and Vortices Near a Free Surface, MD-Vol. 119, ASME, New York, 1991.Google Scholar
[48]Chiu, P. H. and Lin, Y. T., A conservative phase field method for solving incompressible two-phase flows, J. Comput. Phys., 230 (2011), 185–204.Google Scholar
[49]Rieber, M. and Frohn, A., A numerical study on the mechanism of splashing, Int. J. Heat Fluid Flow, 20 (1999), 455–461.Google Scholar
[50]Xiao, F., Ikebata, A. and Hasegawa, T., Numerical simulations of free-interface fluids by a multi integrated moment method, Comput. Struct., 83 (2005), 409–423.Google Scholar
[51]Yokoi, K., A numerical method for free-surface flows and its application to droplet impact on a thin liquid layer, J. Sci. Comput., 35 (2008),372–396.Google Scholar