Hostname: page-component-586b7cd67f-dsjbd Total loading time: 0 Render date: 2024-12-01T01:28:55.099Z Has data issue: false hasContentIssue false
  • English
  • Français

Numerical Simulation of Free Surface by an Area-Preserving Level Set Method

Published online by Cambridge University Press:  20 August 2015

Tony W. H. Sheu  and
C. H. Yu
Tony W. H. Sheu*
Affiliation:
Department of Engineering Science and Ocean Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei, Taiwan Taida Institute of Mathematical Sciences (TIMS), National Taiwan University Center for Quantum Science and Engineering (CQSE), National Taiwan University
C. H. Yu
Affiliation:
Department of Engineering Science and Ocean Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei, Taiwan
*
Corresponding author.Email:[email protected]
Get access

Abstract

We apply in this study an area preserving level set method to simulate gas/water interface flow. For the sake of accuracy, the spatial derivative terms in the equations of motion for an incompressible fluid flow are approximated by the fifth-order accurate upwinding combined compact difference (UCCD) scheme. This scheme development employs two coupled equations to calculate the first- and second-order derivative terms in the momentum equations. For accurately predicting the level set value, the interface tracking scheme is also developed to minimize phase error of the first-order derivative term shown in the pure advection equation. For the purpose of retaining the long-term accurate Hamiltonian in the advection equation for the level set function, the time derivative term is discretized by the sixth-order accurate symplectic Runge-Kutta scheme. Also, to keep as a distance function for ensuring the front having a finite thickness for all time, the re-initialization equation is used. For the verification of the optimized UCCD scheme for the pure advection equation, two benchmark problems have been chosen to investigate in this study. The level set method with excellent area conservation property proposed for capturing the interface in incompressible fluid flows is also verified by solving the dam-break, Rayleigh-Taylor instability, two-bubble rising in water, and droplet falling problems.

Keywords

65E0576T10Level set methodphase errorupwinding combined compact schemeHamiltoniansymplectic Runge-Kuttaarea conservation property
Type
Research Article
Information
Communications in Computational Physics , Volume 11 , Issue 4 , April 2012 , pp. 1347 - 1371
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Sethian, J. A. and Adalsteinsson, D., An overview of level set methods for etching, deposition and lithography development, IEEE T. Semiconduct. M., 10(1) (1997), 167–184.CrossRefGoogle Scholar
[2]Anderson, C. R., A vortex method for flows with slight density variations, J. Comput. Phys., 61(3) (1985), 417–444.CrossRefGoogle Scholar
[3]Boultone-Stone, J. M. and Blake, J. R., Gas bubbles bursting at a free surface, J. Fluid Mech., 254 (1993), 437–466.Google Scholar
[4]Hirt, C. W. and Nichols, B. D., Volume of fluid method (VOF) for the dynamics of free boundaries, J. Comput. Phys., 39 (1981), 201–225.Google Scholar
[5]Unverdi, S. O. and Tryggvason, G., A front-tracking method for viscous, incompressible, multi-fluid flows, J. Comput. Phys., 100 (1992), 25–37.Google Scholar
[6]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
[7]Kim, J., A continuous surface tension force formulation for diffuse-interface models, J. Com-put. Phys., 204 (2005), 784–804.Google 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]Sethian, J. A. and Smereka, P., Level set methods for fluid interfaces, Annu. Rev. Fluid Mech., 35 (2003), 341–372.Google Scholar
[10]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.CrossRefGoogle Scholar
[11]Harlow, F. and Welch, J., Volume tracking methods for interfacial flow calculations, Phys. Fluids, 8 (1965), 21–82.Google Scholar
[12]Osher, S. and Sethian, J. A., Fronts propagating with curvature dependent speed: algorithms based on Hamilton-Jacobi formulation, J. Comput. Phys., 79 (1988), 12–49.CrossRefGoogle Scholar
[13]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
[14]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.CrossRefGoogle Scholar
[15]Sussman, M., Fatermi, E., Smereka, P. and Osher, S., An improved level set method for incompressible two-fluid flows, Comput. Fluids, 127 (1988), 663–680.Google Scholar
[16]Sussman, M. and Puckett, E., 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.CrossRefGoogle Scholar
[17]Strain, J., Tree methods for moving interfaces, J. Comput. Phys., 151 (1999), 616–648.Google Scholar
[18]Osher, S. and Fedkiw, R., Level Set Methods and Dynamic Implicit Surfaces, Springer-Verlag, Berlin, 2003.Google Scholar
[19]Sethian, J. A., Level Set Methods and Fast Marching Methods, Cambridge University Press, Cambridge, 2003.Google Scholar
[20]Olsson, E. and Kreiss, G., A conservative level set method for two phase flow, J. Comput. Phys., 210 (2005), 225–246.CrossRefGoogle Scholar
[21]Olsson, E., Kreiss, G. and Zahedi, S., A conservative level set method for two phase flow II, J. Comput. Phys., 225 (2007), 785–807.Google Scholar
[22]Chu, P. C. and Fan, C., A three-point combined compact difference scheme, J. Comput. Phys., 140 (1998), 370–399.Google Scholar
[23]Zingg, D. W., Comparison of high-accuracy finite-difference methods for linear wave propagation, SIAM J. Sci. Comput., 22(2) (2000), 476–502.Google Scholar
[24]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
[25]Ashcroft, G. and Zhang, X., Optimized prefactored compact schemes, J. Comput. Phys., 190 (2003), 459–477.Google Scholar
[26]Sheu, T. W. H. and Chiu, P. H., A divergence-free-condition compensated method for incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Eng., 196 (2007), 4479–4494.CrossRefGoogle Scholar
[27]McLachlan, R. I., Area preservation in computational fluid dynamics, Phys. Lett. A, 264 (1999), 36–44.Google Scholar
[28]Sussman, M., Smereka, P. and Osher, S., A level set approach for computing solutions to in-compressible two-phase flow, J. Comput. Phys., 114 (1994), 146–159.Google Scholar
[29]Oevel, W. and Sofroniou, W., Symplectic Runge-Kutta schemes II: classification of symplectic methods, Univ. of Paderborn, Germany, Preprint, 1997.Google Scholar
[30]Zhong, X., High-order finite-difference schemes for numerical simulation of hypersonic boundary-layer transition, J. Comput. Phys., 144(2) (1998), 662–709.Google Scholar
[31]Li, Y., Wavenumber-extended high-order upwind-biased finite-difference schemes for convective scalar transport, J. Comput. Phys., 133 (1997), 235–255.CrossRefGoogle Scholar
[32]Wusi, Y., Lin, C. L. and Patel, V. C., Numerical simulation of unsteady multidimensional free surface motions by level set method, Int. J. Numer. Meth. Fluids, 42 (2003), 853–884.Google Scholar
[33]Zalesak, S. T., Fully multidimensional flux-corrected transport algorithms for fluids, J. Com-put. Phys., 31 (1979), 335–362.Google Scholar
[34]Martin, J. C. and Moyce, W. J., An experimental study of the collapse of liquid columns on a rigid horizontal plane, Philos. T. R. Soc. Lond., 224 (1952), 312–324.Google Scholar
[35]Tryggvason, G., Numerical simulations of the Rayleigh-Taylor instability, J. Comput. Phys., 75 (1988), 253–382.CrossRefGoogle Scholar
[36]Guermond, J. L. and Quartapelle, L., A projection FEM for variable density incompressible flows, J. Comput. Phys., 165 (2000), 167–188.CrossRefGoogle Scholar