Development of a finite element/discontinuous Galerkin/level set approach for the simulation of incompressible two phase flow. (April 2018)

• Record Type: Journal Article
• Authors: Gao, Puyang; Ouyang, Jie; Zhou, Wen
• Is Part Of: Advances in engineering software. Volume 118(2018)
• Page Start: 45
• Page End: 59
• Publication Date: 2018-04
• Subjects: Incompressible two phase flow -- Level set approach -- Re-initialization -- Finite element method -- Discontinuous Galerkin -- Metal casting; Computer-aided engineering -- Periodicals; Engineering -- Computer programs -- Periodicals; Engineering -- Software -- Periodicals; Periodicals; 620.0028553
• DOI: 10.1016/j.advengsoft.2018.01.006
• Languages: English
• ISSNs: 0965-9978

Abstract

Highlights: A combined finite element/discontinuous Galerkin/level set method is presented to simulate the incompressible two phase flow. here is no need of resorting to any stabilization terms in the whole calculation process. We can avoid reinitializing the level set function in each time step, and carry out the re-initialization procedure after about each 40 time steps. The good mass conservation feature is able to be preserved. The good results of several numerical examples and the related engineering problem validate the versatility, accuracy, robustness and practicability of the resultant approach. Abstract: In this paper, we present a combined finite element/ discontinuous Galerkin/ level set method to simulate the incompressible two phase flow. The level set method is employed to capture the moving interface because of its simplicity and efficiency when dealing with the significant interface deformations. Due to the hyperbolic nature of the level set equation and the level set re-initialization equation, we apply the second order Runge Kutta Discontinuous Galerkin (RKDG) method to get the stable and precise results. Moreover, in order to obtain the accurate velocity field, the hybrid continuous and discontinuous Galerkin approach is utilized to solve the incompressible Navier–Stokes equations. In our combined method, there is no need to re-initialize the level set function in every time step and the re-initialization process is implemented after suitable time steps. The desirable mass conservation property is able to be guaranteed with a mass correction technique in the combined algorithm. In addition, the stabilization terms can be avoided in the whole computational process. Several challenging problems, i.e., a rising bubble, the dam-break flow, the Rayleigh–Taylor instability, and the complex metal casting process in the engineering applications are investigated to evaluate the feasibility, validity and practicability of our approach for solving the problem of the incompressible two phase flow.