An enhanced FIVER method for multi-material flow problems with second-order convergence rate

详细信息    查看全文

• 作者：Alex Main^c ; Xianyi Zeng^c ; Philip Avery^a ; Charbel Farhat^a ; ^b ; ^c ; ¹ ; ^{cfarhat@stanford.edu}
• 关键词：Embedded boundary method ; Finite volume method ; Immersed boundary method ; Large density jump ; Multi-material ; Multi-phase ; Total variation bounded ; Two-phase Riemann solver
• 刊名：Journal of Computational Physics
• 出版年：2017
• 出版时间：15 January 2017
• 年：2017
• 卷：329
• 期：Complete
• 页码：141-172
• 全文大小：3042 K
• 卷排序：329

文摘

The finite volume (FV) method with exact two-material Riemann problems (FIVER) is an Eulerian computational method for the solution of multi-material flow problems. It is robust in the presence of large density jumps at the fluid–fluid interfaces, and the presence of large structural motions, deformations, and even topological changes at the fluid–structure interfaces. To achieve simplicity in implementation, it approximates each material interface by a surrogate surface which conforms to the control volume boundaries. Unfortunately, this approximation introduces a first-order error of the geometric type in the solution process. In this paper, it is first shown that this error causes the original version of FIVER to be inconsistent in the neighborhood of material interfaces and degrades its global order of spatial accuracy. Then, an enhanced version of FIVER is presented to rectify this issue, restore consistency, and achieve for smooth problems the desired global convergence rate. To this effect, the original definition of a surrogate material interface is retained because of its attractive simplicity. However, the solution at this interface of a two-material Riemann problem is enhanced with a simple reconstruction procedure based on interpolation and extrapolation. Next, the extrapolation component of this procedure is equipped with a limiter in order to achieve nonlinear stability for non-smooth problems. In the one-dimensional inviscid setting, the resulting FIVER method is also shown to be total variation bounded. Focusing on the context of a second-order FV semi-discretization, the nonlinear stability and second-order global convergence rate of this enhanced FIVER method are illustrated for several model multi-fluid and fluid–structure interaction problems. The potential of this computational method for complex multi-material flow problems is also demonstrated with the simulation of the collapse of an air bubble submerged in water and the comparison of the computed results with corresponding experimental data.