Skip to main content Accessibility help

A Second Order Ghost Fluid Method for an Interface Problem of the Poisson Equation

  • Cheng Liu (a1) and Changhong Hu (a1)


A second order Ghost Fluid method is proposed for the treatment of interface problems of elliptic equations with discontinuous coefficients. By appropriate use of auxiliary virtual points, physical jump conditions are enforced at the interface. The signed distance function is used for the implicit description of irregular domain. With the additional unknowns, high order approximation considering the discontinuity can be built. To avoid the ill-conditioned matrix, the interpolation stencils are selected adaptively to balance the accuracy and the numerical stability. Additional equations containing the jump restrictions are assembled with the original discretized algebraic equations to form a new sparse linear system. Several Krylov iterative solvers are tested for the newly derived linear system. The results of a series of 1-D, 2-D tests show that the proposed method possesses second order accuracy in L norm. Besides, the method can be extended to the 3-D problems straightforwardly. Numerical results reveal the present method is highly efficient and robust in dealing with the interface problems of elliptic equations.


Corresponding author

*Corresponding author. Email addresses: (C. Liu), (C. Hu)


Hide All
[1] Ashton, S. R. and Alexei, F. C., A matlab-based finite difference solver for the Poisson problem with mixed Dirichlet-Neumann boundary conditions, Comput. Phys. Commun., 184 (2012), 783798.
[2] Chen, Z. and Zou, J., Finite element methods and their convergence for elliptic and parabolic interface problems, Numer. Math., 79 (1998), 175202.
[3] Peskin, C. S., Flow patterns around heart valves: A numerical method, J. Comput. Phys., 10 (1972), 252271.
[4] Mohd-Yusof, J., Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries, Center for Turbulence Research Annual Research Briefs, 161 (1997), 317327.
[5] LeVeque, R. J. and Li, Z., The immersed interface method for elliptic equations with discontinuous coefficients and singular sources, SIAM J. Numer. Anal., 31 (1994), 10191044.
[6] Li, Z. and Ito, K., Maximum principle preserving schemes for interface problems with discontinuous coefficients, SIAM J. Sci. Comput., 23 (2001), 339361.
[7] Zhong, X., A new high-order immersed interface method for solving elliptic equations with imbedded interface of discontinuity, J. Comput. Phys., 225 (2007), 10661099.
[8] Fedkiw, R., Aslam, T., Merriman, B. and Osher, S., A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method), J. Comput. Phys., 152 (1999), 457.
[9] Liu, X. D., Fedkiw, R. P. and Kang, M., A boundary condition capturing method for Poissons equation on irregular domains, J. Comput. Phys., 160 (2000), 151178.
[10] Zhou, Y., Zhao, S., Feig, M. and Wei, G., High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources, J. Comput. Phys., 213 (2006), 130.
[11] Zhou, Y. and Wei, G., On the fictitious-domain and interpolation formulations of the matched interface and boundary (MIB) method, J. Comput. Phys., 219 (2006), 228246.
[12] Yu, S. and Wei, G., Three-dimensional matched interface and boundary (MIB) method for treating geometric singularities, J. Comput. Phys., 227 (2007), 602632.
[13] Johansen, H. and Colella, P., A Cartesian grid embedded boundary method for Poisson's equation on irregular domains, J. Comput. Phys., 147 (1998), 6085.
[14] Cisternino, M. and Weynans, L., A parallel second order Cartesian method for elliptic interface problems, Tech. Rep., 7573 INRIA; 2011.
[15] Gibou, F., Fedkiw, R. P., Cheng, L. T. and Kang, M., A second-order-accurate symmetric discretization of the Poisson equation on irregular domains, J. Comput. Phys., 176 (2002), 205227.
[16] Bedrossian, J., Von Brecht, J. H. and Zhu, S. et al., A second order virtual node method for elliptic problems with interfaces and irregular domains, J. Comput. Phys., 229 (2010), 64056426.
[17] Hellrung, J. L., Wang, L. and Sifakis, E. et al., A second order virtual node method for elliptic problems with interfaces and irregular domains in three dimensions, J. Comput. Phys., 231 (2012), 20152048.
[18] Choi, J., Oberoi, R., Edwards, J. and Rosati, J., An immersed boundary method for complex incompressible flows, J. Comput. Phys., 224 (2007), 757784.
[19] Liu, C. and Hu, C. H., An efficient immersed boundary treatment for complex moving object, J. Comput. Phys., 274 (2014), 654680.
[20] Saad, Y., Iterative Methods for Sparse Linear Systems, PWS Publishing Company, Boston, 1996.


MSC classification


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed