# American Institute of Mathematical Sciences

2012, 17(4): 1155-1174. doi: 10.3934/dcdsb.2012.17.1155

## Some new finite difference methods for Helmholtz equations on irregular domains or with interfaces

 1 Lilly Corporate Center, DC 4108, Eli Lilly and Company, Indiana, IN 46285, United States 2 Center For Research in Scientiﬁc Computation & Department of Mathematics, North Carolina State University, Raleigh, NC 27695-8205

Received  December 2010 Revised  September 2011 Published  February 2012

Solving a Helmholtz equation $\Delta u + \lambda u = f$ efficiently is a challenge for many applications. For example, the core part of many efficient solvers for the incompressible Navier-Stokes equations is to solve one or several Helmholtz equations. In this paper, two new finite difference methods are proposed for solving Helmholtz equations on irregular domains, or with interfaces. For Helmholtz equations on irregular domains, the accuracy of the numerical solution obtained using the existing augmented immersed interface method (AIIM) may deteriorate when the magnitude of $\lambda$ is large. In our new method, we use a level set function to extend the source term and the PDE to a larger domain before we apply the AIIM. For Helmholtz equations with interfaces, a new maximum principle preserving finite difference method is developed. The new method still uses the standard five-point stencil with modifications of the finite difference scheme at irregular grid points. The resulting coefficient matrix of the linear system of finite difference equations satisfies the sign property of the discrete maximum principle and can be solved efficiently using a multigrid solver. The finite difference method is also extended to handle temporal discretized equations where the solution coefficient $\lambda$ is inversely proportional to the mesh size.
Citation: Xiaohai Wan, Zhilin Li. Some new finite difference methods for Helmholtz equations on irregular domains or with interfaces. Discrete & Continuous Dynamical Systems - B, 2012, 17 (4) : 1155-1174. doi: 10.3934/dcdsb.2012.17.1155
##### References:
 [1] J. Adams, P. Swarztrauber and R. Sweet, FISHPACK: Efficient FORTRAN subprograms for the solution of separable elliptic partial differential equations., Available from: \url{http://www.netlib.org/fishpack/}., (). [2] J. B. Bell, P. Colella and H. M. Glaz., A second-order projection method for the incompressible Navier-Stokes equations,, J. Comput. Phys., 85 (1989), 257. doi: 10.1016/0021-9991(89)90151-4. [3] D. L. Brown, R. Cortez and M. L. Minion, Accurate projection methods for the incompressible Navier-Stokes equations,, J. Comput. Phys., 168 (2001), 464. doi: 10.1006/jcph.2001.6715. [4] D. Calhoun, A Cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions,, J. Comput. Phys., 176 (2002), 231. doi: 10.1006/jcph.2001.6970. [5] D. L. Chopp, Some improvements of the fast marching method,, SIAM J. Sci. Comput., 23 (2001), 230. doi: 10.1137/S106482750037617X. [6] S. Deng, "Immersed Interface Method for Three Dimensional Interface Problems and Applications,", Ph.D thesis, (2001). [7] S. Deng, K. Ito and Z. Li, Three-dimensional elliptic solvers for interface problems and applications,, J. Comput. Phys., 184 (2003), 215. doi: 10.1016/S0021-9991(02)00028-1. [8] W. E and J.-G. Liu, Projection method. I. Convergence and numerical boundary-layers,, SIAM J. Numer. Anal., 32 (1995), 1017. doi: 10.1137/0732047. [9] T. Hou, Z. Li, S. Osher and H. Zhao, A hybrid method for moving interface problems with application to the Hele-Shaw flow,, J. Comput. Phys., 134 (1997), 236. doi: 10.1006/jcph.1997.5689. [10] J. Hunter, Z. Li and H. Zhao, Reactive autophobic spreading of drops,, J. Comput. Phys., 183 (2002), 335. doi: 10.1006/jcph.2002.7168. [11] Z. Li, "The Immersed Interface Method: A Numerical Approach for Partial Differential Equations with Interfaces,", Ph.D thesis, (1994). [12] Z. Li, A fast iterative algorithm for elliptic interface problems,, SIAM J. Numer. Anal., 35 (1998), 230. doi: 10.1137/S0036142995291329. [13] Z. Li and K. Ito, Maximum principle preserving schemes for interface problems with discontinuous coefficients,, SIAM J. Sci. Comput., 23 (2001), 339. doi: 10.1137/S1064827500370160. [14] Z. Li and K. Ito, "The Immersed Interface Method. Numerical Solutions of PDEs Involving Interfaces and Irregular Domains,", Frontiers in Applied Mathematics, 33 (2006). [15] Z. Li, S. R. Lubkin and X. Wan, An augmented IIM-level set method for Stokes equations with discontinuous viscosity,, in, 15 (2007), 193. [16] Z. Li and B. Soni, Fast and accurate numerical approaches for Stefan problems and crystal growth,, Numerical Heat Transfer, 35 (1999), 461. doi: 10.1080/104077999275848. [17] Z. Li, H. Zhao and H. Gao, A numerical study of electro-migration voiding by evolving level set functions on a fixed cartesian grid,, J. Comput. Phys., 152 (1999), 281. doi: 10.1006/jcph.1999.6249. [18] M. N. Linnick and H. F. Fasel, A high-order immersed interface method for simulating unsteady incompressible flows on irregular domains,, J. Comput. Phys., 204 (2005), 157. doi: 10.1016/j.jcp.2004.09.017. [19] S. Osher and R. Fedkiw, "Level Set Methods and Dynamic Implicit Surfaces,", Applied Mathematical Sciences, 153 (2003). [20] D. Russell and Z. J. Wang, A Cartesian grid method for modeling multiple moving objects in 2D incompressible viscous flow,, J. Comput. Phys., 191 (2003), 177. doi: 10.1016/S0021-9991(03)00310-3. [21] J. A. Sethian, "Level Set Methods and Fast Marching Methods. Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science,", 2nd edition, 3 (1999). [22] X. Wan, "Numerical Simulation Methods for Biological Tissue Interactions,", Ph.D thesis, (2007). [23] W.-J. Ying and C. S. Henriquez, A kernel-free boundary integral method for elliptic boundary value problems,, J. Comput. Phy., 227 (2007), 1046. doi: 10.1016/j.jcp.2007.08.021. [24] P. M. de Zeeuw, Matrix-dependent prolongations and restrictions in a blackbox multigrid solver,, J. Comput. Appl. Math., 33 (1990), 1. doi: 10.1016/0377-0427(90)90252-U.

show all references

##### References:
 [1] J. Adams, P. Swarztrauber and R. Sweet, FISHPACK: Efficient FORTRAN subprograms for the solution of separable elliptic partial differential equations., Available from: \url{http://www.netlib.org/fishpack/}., (). [2] J. B. Bell, P. Colella and H. M. Glaz., A second-order projection method for the incompressible Navier-Stokes equations,, J. Comput. Phys., 85 (1989), 257. doi: 10.1016/0021-9991(89)90151-4. [3] D. L. Brown, R. Cortez and M. L. Minion, Accurate projection methods for the incompressible Navier-Stokes equations,, J. Comput. Phys., 168 (2001), 464. doi: 10.1006/jcph.2001.6715. [4] D. Calhoun, A Cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions,, J. Comput. Phys., 176 (2002), 231. doi: 10.1006/jcph.2001.6970. [5] D. L. Chopp, Some improvements of the fast marching method,, SIAM J. Sci. Comput., 23 (2001), 230. doi: 10.1137/S106482750037617X. [6] S. Deng, "Immersed Interface Method for Three Dimensional Interface Problems and Applications,", Ph.D thesis, (2001). [7] S. Deng, K. Ito and Z. Li, Three-dimensional elliptic solvers for interface problems and applications,, J. Comput. Phys., 184 (2003), 215. doi: 10.1016/S0021-9991(02)00028-1. [8] W. E and J.-G. Liu, Projection method. I. Convergence and numerical boundary-layers,, SIAM J. Numer. Anal., 32 (1995), 1017. doi: 10.1137/0732047. [9] T. Hou, Z. Li, S. Osher and H. Zhao, A hybrid method for moving interface problems with application to the Hele-Shaw flow,, J. Comput. Phys., 134 (1997), 236. doi: 10.1006/jcph.1997.5689. [10] J. Hunter, Z. Li and H. Zhao, Reactive autophobic spreading of drops,, J. Comput. Phys., 183 (2002), 335. doi: 10.1006/jcph.2002.7168. [11] Z. Li, "The Immersed Interface Method: A Numerical Approach for Partial Differential Equations with Interfaces,", Ph.D thesis, (1994). [12] Z. Li, A fast iterative algorithm for elliptic interface problems,, SIAM J. Numer. Anal., 35 (1998), 230. doi: 10.1137/S0036142995291329. [13] Z. Li and K. Ito, Maximum principle preserving schemes for interface problems with discontinuous coefficients,, SIAM J. Sci. Comput., 23 (2001), 339. doi: 10.1137/S1064827500370160. [14] Z. Li and K. Ito, "The Immersed Interface Method. Numerical Solutions of PDEs Involving Interfaces and Irregular Domains,", Frontiers in Applied Mathematics, 33 (2006). [15] Z. Li, S. R. Lubkin and X. Wan, An augmented IIM-level set method for Stokes equations with discontinuous viscosity,, in, 15 (2007), 193. [16] Z. Li and B. Soni, Fast and accurate numerical approaches for Stefan problems and crystal growth,, Numerical Heat Transfer, 35 (1999), 461. doi: 10.1080/104077999275848. [17] Z. Li, H. Zhao and H. Gao, A numerical study of electro-migration voiding by evolving level set functions on a fixed cartesian grid,, J. Comput. Phys., 152 (1999), 281. doi: 10.1006/jcph.1999.6249. [18] M. N. Linnick and H. F. Fasel, A high-order immersed interface method for simulating unsteady incompressible flows on irregular domains,, J. Comput. Phys., 204 (2005), 157. doi: 10.1016/j.jcp.2004.09.017. [19] S. Osher and R. Fedkiw, "Level Set Methods and Dynamic Implicit Surfaces,", Applied Mathematical Sciences, 153 (2003). [20] D. Russell and Z. J. Wang, A Cartesian grid method for modeling multiple moving objects in 2D incompressible viscous flow,, J. Comput. Phys., 191 (2003), 177. doi: 10.1016/S0021-9991(03)00310-3. [21] J. A. Sethian, "Level Set Methods and Fast Marching Methods. Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science,", 2nd edition, 3 (1999). [22] X. Wan, "Numerical Simulation Methods for Biological Tissue Interactions,", Ph.D thesis, (2007). [23] W.-J. Ying and C. S. Henriquez, A kernel-free boundary integral method for elliptic boundary value problems,, J. Comput. Phy., 227 (2007), 1046. doi: 10.1016/j.jcp.2007.08.021. [24] P. M. de Zeeuw, Matrix-dependent prolongations and restrictions in a blackbox multigrid solver,, J. Comput. Appl. Math., 33 (1990), 1. doi: 10.1016/0377-0427(90)90252-U.
 [1] Jian Hao, Zhilin Li, Sharon R. Lubkin. An augmented immersed interface method for moving structures with mass. Discrete & Continuous Dynamical Systems - B, 2012, 17 (4) : 1175-1184. doi: 10.3934/dcdsb.2012.17.1175 [2] Champike Attanayake, So-Hsiang Chou. An immersed interface method for Pennes bioheat transfer equation. Discrete & Continuous Dynamical Systems - B, 2015, 20 (2) : 323-337. doi: 10.3934/dcdsb.2015.20.323 [3] So-Hsiang Chou. An immersed linear finite element method with interface flux capturing recovery. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2343-2357. doi: 10.3934/dcdsb.2012.17.2343 [4] Zhongyi Huang. Tailored finite point method for the interface problem. Networks & Heterogeneous Media, 2009, 4 (1) : 91-106. doi: 10.3934/nhm.2009.4.91 [5] Qiang Du, Manlin Li. On the stochastic immersed boundary method with an implicit interface formulation. Discrete & Continuous Dynamical Systems - B, 2011, 15 (2) : 373-389. doi: 10.3934/dcdsb.2011.15.373 [6] Shi Jin, Xu Yang, Guangwei Yuan. A domain decomposition method for a two-scale transport equation with energy flux conserved at the interface. Kinetic & Related Models, 2008, 1 (1) : 65-84. doi: 10.3934/krm.2008.1.65 [7] S. L. Ma'u, P. Ramankutty. An averaging method for the Helmholtz equation. Conference Publications, 2003, 2003 (Special) : 604-609. doi: 10.3934/proc.2003.2003.604 [8] Sheng Xu. Derivation of principal jump conditions for the immersed interface method in two-fluid flow simulation. Conference Publications, 2009, 2009 (Special) : 838-845. doi: 10.3934/proc.2009.2009.838 [9] Shi Jin, Min Tang, Houde Han. A uniformly second order numerical method for the one-dimensional discrete-ordinate transport equation and its diffusion limit with interface. Networks & Heterogeneous Media, 2009, 4 (1) : 35-65. doi: 10.3934/nhm.2009.4.35 [10] Daniele Boffi, Lucia Gastaldi. Discrete models for fluid-structure interactions: The finite element Immersed Boundary Method. Discrete & Continuous Dynamical Systems - S, 2016, 9 (1) : 89-107. doi: 10.3934/dcdss.2016.9.89 [11] Tao Lin, Yanping Lin, Weiwei Sun. Error estimation of a class of quadratic immersed finite element methods for elliptic interface problems. Discrete & Continuous Dynamical Systems - B, 2007, 7 (4) : 807-823. doi: 10.3934/dcdsb.2007.7.807 [12] Zhenlin Guo, Ping Lin, Guangrong Ji, Yangfan Wang. Retinal vessel segmentation using a finite element based binary level set method. Inverse Problems & Imaging, 2014, 8 (2) : 459-473. doi: 10.3934/ipi.2014.8.459 [13] Siu-Long Lei. Adaptive method for spike solutions of Gierer-Meinhardt system on irregular domain. Discrete & Continuous Dynamical Systems - B, 2011, 15 (3) : 651-668. doi: 10.3934/dcdsb.2011.15.651 [14] Anita T. Layton, J. Thomas Beale. A partially implicit hybrid method for computing interface motion in Stokes flow. Discrete & Continuous Dynamical Systems - B, 2012, 17 (4) : 1139-1153. doi: 10.3934/dcdsb.2012.17.1139 [15] Jing Xu, Xue-Cheng Tai, Li-Lian Wang. A two-level domain decomposition method for image restoration. Inverse Problems & Imaging, 2010, 4 (3) : 523-545. doi: 10.3934/ipi.2010.4.523 [16] Xueyong Wang, Yiju Wang, Gang Wang. An accelerated augmented Lagrangian method for multi-criteria optimization problem. Journal of Industrial & Management Optimization, 2017, 13 (5) : 1-9. doi: 10.3934/jimo.2018136 [17] Ben A. Vanderlei, Matthew M. Hopkins, Lisa J. Fauci. Error estimation for immersed interface solutions. Discrete & Continuous Dynamical Systems - B, 2012, 17 (4) : 1185-1203. doi: 10.3934/dcdsb.2012.17.1185 [18] Moulay Rchid Sidi Ammi, Ismail Jamiai. Finite difference and Legendre spectral method for a time-fractional diffusion-convection equation for image restoration. Discrete & Continuous Dynamical Systems - S, 2018, 11 (1) : 103-117. doi: 10.3934/dcdss.2018007 [19] Ömer Oruç, Alaattin Esen, Fatih Bulut. A unified finite difference Chebyshev wavelet method for numerically solving time fractional Burgers' equation. Discrete & Continuous Dynamical Systems - S, 2019, 12 (3) : 533-542. doi: 10.3934/dcdss.2019035 [20] Liuyang Yuan, Zhongping Wan, Jingjing Zhang, Bin Sun. A filled function method for solving nonlinear complementarity problem. Journal of Industrial & Management Optimization, 2009, 5 (4) : 911-928. doi: 10.3934/jimo.2009.5.911

2017 Impact Factor: 0.972