A New High-Order Immersed Interface Method for Multi-Phase Flow Simulation

Abstract

Numerical simulation of two-phase flow with imbedded interface of discontinuity is challenging in two major aspects. First, the interface can undergo change, merge and breakup during the course of the simulation. Examples of current successful methods in modeling flow with interface are, among others, the volume-of-fluid method, fronttracking method, and level-set method. Second, flow variables and their derivatives can be discontinuous across the interface. This discontinuity poses severe limitation on the accuracy of commonly used numerical methods. Current available methods in treating the interface discontinuity, such as the immersed boundary method and the ghost fluid method, are mostly first order accurate. Though the immersed interface method of LeVeque and Li (1994) can be globally second order, it is often difficult to apply to complex multi-dimensional flow. This paper presents a new high-order immersed interface method for elliptic equations with imbedded interface of discontinuity. The new method can be arbitrarily high-order accurate, and it can be easily applied to practical two-phase flow because only the physical jump conditions for variables and their first derivatives are needed in the finite difference formulas. In addition, the new interface difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The new interface algorithms of up to accuracy have been tested for one and two-dimensional elliptic equations with imbedded interface. The extension to practical two-phase flow applications will be presented in a future paper.

A fixed-mesh method for incompressible flow–structure systems with finite solid deformations

“…For now, our primary interest is the general investigation of fixed-mesh algorithms, so the first order, easy-to-implement immerse-interface method by Lee and LeVeque [29] is adopted in this study. Extension to higher order schemes is possible [53,56].…”

On high-order shock-fitting and front-tracking schemes for numerical simulation of shock–disturbance interactions

“…Moreover, the ad hoc parameter r* would not be needed. Recently, Zhong [73,78] developed a new immersed interface algorithm that requires jumps only in values of flow variables and their first derivatives while showing arbitrarily high-order accuracy. Since up to first-order jump conditions can be derived from the physics of the problem, this method is more suitable for high-order computations of real problems than other IIM schemes which need higher order derivatives of flow variables to obtain higher order of accuracy.…”

A high-order cut-cell method for numerical simulation of hypersonic boundary-layer instability with surface roughness

“…For multi-phase flows with moving interfaces, an immersed interface method was developed and has been widely used [25][26][27][28][29][30][31][32][33][34][35]. For an immersed interface method, a Cartesian grid is often used where a sharp interface moves freely across the fixed grid lines.…”