Skip to main content Accessibility help
×
Home
Hostname: page-component-544b6db54f-fg2fv Total loading time: 0.182 Render date: 2021-10-20T20:27:10.104Z Has data issue: true Feature Flags: { "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true, "metricsAbstractViews": false, "figures": true, "newCiteModal": false, "newCitedByModal": true, "newEcommerce": true, "newUsageEvents": true }

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*
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:twhsheu@ntu.edu.tw
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.

Type
Research Article
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), 167184.CrossRefGoogle Scholar
[2]Anderson, C. R., A vortex method for flows with slight density variations, J. Comput. Phys., 61(3) (1985), 417444.CrossRefGoogle Scholar
[3]Boultone-Stone, J. M. and Blake, J. R., Gas bubbles bursting at a free surface, J. Fluid Mech., 254 (1993), 437466.CrossRefGoogle 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), 201225.CrossRefGoogle Scholar
[5]Unverdi, S. O. and Tryggvason, G., A front-tracking method for viscous, incompressible, multi-fluid flows, J. Comput. Phys., 100 (1992), 2537.CrossRefGoogle Scholar
[6]Badalassi, V. E., Ceniceros, H. D. and Banerjee, S., Computation of multiphase systems with phase field models, J. Comput. Phys., 190 (2003), 371397.CrossRefGoogle Scholar
[7]Kim, J., A continuous surface tension force formulation for diffuse-interface models, J. Com-put. Phys., 204 (2005), 784804.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), 20782095.CrossRefGoogle Scholar
[9]Sethian, J. A. and Smereka, P., Level set methods for fluid interfaces, Annu. Rev. Fluid Mech., 35 (2003), 341372.CrossRefGoogle 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), 203216.CrossRefGoogle Scholar
[11]Harlow, F. and Welch, J., Volume tracking methods for interfacial flow calculations, Phys. Fluids, 8 (1965), 2182.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), 1249.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), 338357.CrossRefGoogle 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), 83116.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), 663680.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), 301337.CrossRefGoogle Scholar
[17]Strain, J., Tree methods for moving interfaces, J. Comput. Phys., 151 (1999), 616648.CrossRefGoogle Scholar
[18]Osher, S. and Fedkiw, R., Level Set Methods and Dynamic Implicit Surfaces, Springer-Verlag, Berlin, 2003.CrossRefGoogle 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), 225246.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), 785807.CrossRefGoogle Scholar
[22]Chu, P. C. and Fan, C., A three-point combined compact difference scheme, J. Comput. Phys., 140 (1998), 370399.CrossRefGoogle Scholar
[23]Zingg, D. W., Comparison of high-accuracy finite-difference methods for linear wave propagation, SIAM J. Sci. Comput., 22(2) (2000), 476502.CrossRefGoogle Scholar
[24]Tam, C. K. W. and Webb, J. C., Dispersion-relation-preserving finite difference schemes for computational acoustics, J. Comput. Phys., 107 (1993), 262281.CrossRefGoogle Scholar
[25]Ashcroft, G. and Zhang, X., Optimized prefactored compact schemes, J. Comput. Phys., 190 (2003), 459477.CrossRefGoogle 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), 44794494.CrossRefGoogle Scholar
[27]McLachlan, R. I., Area preservation in computational fluid dynamics, Phys. Lett. A, 264 (1999), 3644.CrossRefGoogle 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), 146159.CrossRefGoogle 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), 662709.CrossRefGoogle Scholar
[31]Li, Y., Wavenumber-extended high-order upwind-biased finite-difference schemes for convective scalar transport, J. Comput. Phys., 133 (1997), 235255.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), 853884.Google Scholar
[33]Zalesak, S. T., Fully multidimensional flux-corrected transport algorithms for fluids, J. Com-put. Phys., 31 (1979), 335362.CrossRefGoogle 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), 312324.CrossRefGoogle Scholar
[35]Tryggvason, G., Numerical simulations of the Rayleigh-Taylor instability, J. Comput. Phys., 75 (1988), 253382.CrossRefGoogle Scholar
[36]Guermond, J. L. and Quartapelle, L., A projection FEM for variable density incompressible flows, J. Comput. Phys., 165 (2000), 167188.CrossRefGoogle Scholar
3
Cited by

Send article to Kindle

To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

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

Send article to Dropbox

To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

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

Send article to Google Drive

To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

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

Reply to: Submit a response

Please enter your response.

Your details

Please enter a valid email address.

Conflicting interests

Do you have any conflicting interests? *