An improved interface preserving level set method for simulating three dimensional rising bubble. (December 2016)

• Record Type: Journal Article
• Authors: Yu, C.H.; Ye, Z.T.; Sheu, Tony W.H.; Lin, Y.T.; Zhao, X.Z.
• Is Part Of: International journal of heat and mass transfer. Volume 103(2016:Dec.)
• Page Start: 753
• Page End: 772
• Publication Date: 2016-12
• Subjects: Interface preserving level set method -- Bubble rising -- Improved smoothed Heaviside function -- Mass conservation term -- Re-initialization equation -- Surface tension force; Heat -- Transmission -- Periodicals; Mass transfer -- Periodicals; Chaleur -- Transmission -- Périodiques; Transfert de masse -- Périodiques; Electronic journals; 621.4022
• DOI: 10.1016/j.ijheatmasstransfer.2016.07.096
• Languages: English
• ISSNs: 0017-9310

Abstract

Highlights: An improved interface preserving level set method which conserves the mass excellently have been developed to predict the bubble rising flow in the two-phase flow model. The Navier–Stokes equations with surface tension terms are calculated effectively using the level set base density-scaled CSF model. The tangent of hyperbola for interface capturing with the slope weighting (THINC/SW) scheme is also used to validate the abilities of accurately predicting interface for a bubble bursting at a free surface. Abstract: A new interface preserving level set method is developed in three steps to simulate bubble rising problems. In the first step of the solution algorithm, the level set function ϕ is advected by a pure advection equation. An intermediate step is performed to obtain new level set function through an improved smoothed Heaviside function. To keep the new level set function as a distance function and to conserve mass bounded by the interface, in the final solution step a mass correction term is added to the re-initialization equation. This two-phase numerical model is developed underlying the projection method to compute the incompressible Navier–Stokes solutions in collocated grids. In the discretizations of the level set advection equation and the re-initialization equation, the fifth-order weighted essentially non-oscillatory scheme is applied to prevent numerical oscillations occurring around discontinuous interface. The performance of the proposed level set method in conserving mass is compared with conventional level set method applied to solve the single bubble rising problem and the bubble bursting problem at a free surface. Merger of two bubbles is also investigated. Numerical results show that not only the surface tension force can be accurately calculated but also the mass can be conserved excellently using the present level set method.