INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS, vol.14, no.5, 2017 (SCI-Expanded)
We present a fully second order IMplicit/EXplicit (IMEX) time integration technique for solving incompressible multi-phase flow problems. A typical incompressible multiphase flow model consists of the Navier-Stokes equations plus an interface dynamics equation (e.g., the level set equation). Our IMEX strategy is applied to such a model in the following manner. The hyperbolic terms of the Navier-Stokes equations together with the interface dynamics equation are solved explicitly (Explicit Block) making use of the well-understood explicit numerical schemes [Leveque, R. J. [1998] Finite Volume Methods for Hyperbolic Problems, "Texts in Applied Mathematics", (Cambridge University Press); Thomas, J. W. [1999] Numerical Partial Differential Equations II (Conservation Laws and Elliptic Equations), " Texts in Applied Mathematics" (Springer-Verlag, New York)]. On the other hand, the nonhyperbolic (stiff) parts of the flow equations are solved implicitly (Implicit Block) within the framework of the Jacobian-Free Newton Krylov (JFNK) method [Knoll, D. A. and Keyes, D. E. [2004] Jacobian-free Newton Krylov methods: A survey of approaches and applications. J. Comput. Phys. 193, 357397; Saad, Y. [2003] Iterative Methods for Sparse Linear Systems (Siam); Kelley, C. T. [2003] Solving Nonlinear Equations with Newton's Method (Siam)]. In our algorithm implementation, the explicit block is embedded in the implicit block in a way that it is always part of the nonlinear function evaluation. In this way, there exists a continuous interaction between the implicit and explicit algorithm blocks meaning that the improved solutions (in terms of time accuracy) at each nonlinear iteration are immediately felt by the explicit block and the improved explicit solutions are readily available to form the next set of nonlinear residuals. This continuous interaction between the two algorithm blocks results in an implicitly balanced algorithm in that all nonlinearities due to coupling of different time terms are converged with the desired numerical time accuracy. In other words, we obtain a self-consistent IMEX method that eliminates the possible order reductions in time convergence that is quite common in certain types of nonlinearly coupled systems. We remark that an incompressible multi-phase flow model can be a highly nonlinearly coupled system with the involvement of very stiff surface tension source terms. These kinds of flow problems are difficult to tackle numerically.