Skip to main content Accessibility help
×
Home

A PRECONDITIONED METHOD FOR THE SOLUTION OF THE ROBBINS PROBLEM FOR THE HELMHOLTZ EQUATION

  • JIANG LE (a1) (a2), HUANG JIN (a1), XIAO-GUANG LV (a1) (a2) and QING-SONG CHENG (a2)

Abstract

A preconditioned iterative method for the two-dimensional Helmholtz equation with Robbins boundary conditions is discussed. Using a finite-difference method to discretize the Helmholtz equation leads to a sparse system of equations which is too large to solve directly. The approach taken in this paper is to precondition this linear system with a sine transform based preconditioner and then solve it using the generalized minimum residual method (GMRES). An analytical formula for the eigenvalues of the preconditioned matrix is derived and it is shown that the eigenvalues are clustered around 1 except for some outliers. Numerical results are reported to demonstrate the effectiveness of the proposed method.

    • Send article to Kindle

      To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

      Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

      Find out more about the Kindle Personal Document Service.

      A PRECONDITIONED METHOD FOR THE SOLUTION OF THE ROBBINS PROBLEM FOR THE HELMHOLTZ EQUATION
      Available formats
      ×

      Send article to Dropbox

      To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

      A PRECONDITIONED METHOD FOR THE SOLUTION OF THE ROBBINS PROBLEM FOR THE HELMHOLTZ EQUATION
      Available formats
      ×

      Send article to Google Drive

      To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

      A PRECONDITIONED METHOD FOR THE SOLUTION OF THE ROBBINS PROBLEM FOR THE HELMHOLTZ EQUATION
      Available formats
      ×

Copyright

Corresponding author

For correspondence; e-mail: jiangle@hhit.edu.cn

Footnotes

Hide All

This research is supported by NSFC (10871034) and Research Grant Z2008038 from HuaiHai Institute of Technology.

Footnotes

References

Hide All
[1]Arnold, A. and Ehrhardt, M., “Discrete transparent boundary conditions for wide angle parabolic equations in underwater acoustics”, J. Comput. Phys. 145 (1998) 611638.
[2]Axelsson, O. and Barker, V., Finite element solution of boundary value problems: theory and computation (Academic Press, Orlando, FL, 1984).
[3]Bayliss, A., Goldstein, C. I. and Turkel, E., “An iterative method for the Helmholtz equation”, J. Comput. Phys. 49 (1983) 443457.
[4]Chan, R. and Ng, M., “Conjugate gradient methods for Toeplitz systems”, SIAM Rev. 38 (1996) 427482.
[5]Erlangga, Y. A., “A robust and efficient iterative method for the numerical solution of the Helmholtz equation”, Ph. D. Thesis, Delft University of Technology, 2005.
[6]Erlangga, Y. A., Vuik, C. and Oosterlee, C. W., “On a class of preconditioners for solving the Helmholtz equation”, Appl. Numer. Math. 50 (2004) 409425.
[7]Gander, M. J. and Nataf, F., “AILU for Helmholtz problems: a new preconditioner based on the analytical parabolic factorization”, J. Comput. Acoust. 9 (2001) 14991506.
[8]Greenbaum, A., Iterative methods for solving linear systems, Volume 17 of Frontiers in Applied Mathematics (SIAM, Philadelphia, PA, 1997).
[9]Guo, C. H., “Incomplete block factorization preconditioning for linear systems arising in the numerical solution of the Helmholtz equation”, Appl. Numer. Math. 19 (1996) 495508.
[10]Haldenwang, P., Labrosse, G., Abboudi, S. and Deville, M., “Chebyshev 3-D spectral and 2-D pseudospectral solvers for the Helmholtz equation”, J. Comput. Phys. 55 (1984) 115128.
[11]Hemmingsson, L. and Otto, K., “Analysis of semi-Toeplitz preconditioners for first-order PDE”, SIAM J. Sci. Comput. 17 (1996) 4764.
[12]Ho, A. C. and Ng, M., “Iterative methods for Robbins problems”, Appl. Math. Comput. 165 (2005) 103125.
[13]Hu, X. and Ling, X., “Preconditioners for elliptic problems via non-uniform meshes”, Appl. Math. Comput. 181 (2006) 11821198.
[14]Joubert, W. D. and Manteuffel, T. A., “Iterative methods for non-symmetric linear systems”, in: Iterative methods for large linear systems (eds Kincaid, D. R. and Hayes, L. J.), (Academic Press, New York, 1990) 348358.
[15]Kechroud, R., Soulaimani, A. and Saad, Y., “Preconditioning techniques for the solution of the Helmholtz equation by the finite element method”, 2003 Workshop in Wave Phenomena in Physics and Engineering: New Models, Algorithms, and Applications, May 18–21, 2003 (eds V. Kumar et al.), (Springer, Berlin, 2003).
[16]Laird, A. L. and Giles, M. B., “Preconditioned iterative solution of the 2D Helmholtz equation”, Technical Report NA 02-12, Computing Laboratory, Oxford University, 2002.
[17]Laub, A. J., Matrix analysis for scientists and engineers (Society for Industrial and Applied Mathematics, Philadelphia, PA, 2005).
[18]Manteuffel, T. A. and Parter, S. V., “Preconditioning and boundary conditions”, SIAM J. Numer. Anal. 27 (1990) 656694.
[19]Monga Made, M. M., “Incomplete factorization-based preconditionings for solving the Helmholtz equation”, Internat. J. Numer. Methods Engrg. 50 (2001) 10771101.
[20]Otto, K., “Construction and analysis of preconditioners for first-order PDE”, Ph. D. Thesis, Department of Scientific Computing, Uppsala University, Uppsala, Sweden, 1993.
[21]Pickering, W. and Harley, P., “FFT solution of the Robbins problem”, IMA J. Numer. Anal. 13 (1993) 215233.
[22]Pickering, W. and Harley, P., “On Robbins boundary conditions, elliptic equations, and FFT methods”, J. Comput. Phys. 122 (1995) 380383.
[23]Plessix, R. E. and Mulder, W. A., “Separation-of-variables as a preconditioner for an iterative Helmholtz solver”, Appl. Numer. Math. 44 (2003) 385400.
[24]Saad, Y., Iterative methods for sparse linear system (PWS Press, New York, NY, 1996).
[25]Urbach, H. P. and Merkx, R. T. M., “Finite element simulation of electromagnetic plane wave diffraction at gratings for arbitrary angles of incidence”, in: Mathematical and numerical aspects of wave propagation phenomena, Proceedings of SIAM First Conference on Mathematical and Numerical Aspects of Wave Propagation (eds Cohen, G., Halpern, L. and Joly, P.), (Society for Industrial and Applied Mathematics, Philadelphia, PA, 1992) 8999.
[26]Van Loan, C. and Pitsianis, N. P., “Approximation with Kronecker products”, in: Linear algebra for large scale and real time applications (eds Moonen, M. S. and Golub, G. H.), (Kluwer Publications, Dordrecht, 1992) 293314.
[27]Vuik, C., Erlangga, Y. A. and Oosterlee, C. W., “Shifted Laplace preconditioners for the Helmholtz equations”, Report 03-18, Delft University of Technology, ISSN 1389-6520, 2003.
MathJax
MathJax is a JavaScript display engine for mathematics. For more information see http://www.mathjax.org.

Keywords

MSC classification

Related content

Powered by UNSILO

A PRECONDITIONED METHOD FOR THE SOLUTION OF THE ROBBINS PROBLEM FOR THE HELMHOLTZ EQUATION

  • JIANG LE (a1) (a2), HUANG JIN (a1), XIAO-GUANG LV (a1) (a2) and QING-SONG CHENG (a2)

Metrics

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.