Hostname: page-component-cd9895bd7-8ctnn Total loading time: 0 Render date: 2024-12-27T09:42:02.806Z Has data issue: false hasContentIssue false
  • English
  • Français

Numerical Simulation of Time-Harmonic Waves in Inhomogeneous Media using Compact High Order Schemes

Published online by Cambridge University Press:  20 August 2015

Steven Britt ,
Semyon Tsynkov  and
Eli Turkel
Steven Britt*
Affiliation:
Department of Mathematics, North Carolina State University, Box 8205, Raleigh, NC 27695, USA
Semyon Tsynkov*
Affiliation:
Department of Mathematics, North Carolina State University, Box 8205, Raleigh, NC 27695, USA
Eli Turkel*
Affiliation:
School of Mathematical Sciences, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
*
Email:dsbritt@ncsu.edu
Corresponding author.Email:tsynkov@math.ncsu.edu
Email:turkel@post.tau.ac.il
Get access

Abstract

In many problems, one wishes to solve the Helmholtz equation with variable coefficients within the Laplacian-like term and use a high order accurate method (e.g., fourth order accurate) to alleviate the points-per-wavelength constraint by reducing the dispersion errors. The variation of coefficients in the equation may be due to an inhomogeneous medium and/or non-Cartesian coordinates. This renders existing fourth order finite difference methods inapplicable. We develop a new compact scheme that is provably fourth order accurate even for these problems. We present numerical results that corroborate the fourth order convergence rate for several model problems.

Keywords

65N0678A4878M20Helmholtz equationvariable coefficientshigh order accuracycompact finite differences
Type
Research Article
Information
Communications in Computational Physics , Volume 9 , Issue 3 , March 2011 , pp. 520 - 541
Copyright
Copyright © Global Science Press Limited 2011

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Agranovich, M. S., Spectral properties of diffraction problems, In Voitovich, N. N., Kat-senelenbaum, B. Z. and Sivov, A. N., editors, Method of Generalized Resonances in Diffraction Theory, Nauka, Moscow, 1977, [in Russian].Google Scholar
[2]Agranovich, M. S., Katsenelenbaum, B. Z., Sivov, A. N. and Voitovich, N. N., Generalized Method of Eigenoscillations in Diffraction Theory, WILEY-VCH Verlag Berlin GmbH, Berlin, 1999, Translated from the Russian Manuscript by Nazaikinskii, Vladimir.Google Scholar
[3] Ivo Babžka, M. and Sauter, Stefan A., Is the pollution effect of the FEM avoidable for the Helmholtz equation considering high wave numbers?, SIAM Rev., 42(3) (2000), 451–484 (electronic), reprint of SIAM J. Numer. Anal., 34(6) (1997), 2392–2423, [MR1480387 (99b:65135)].Google Scholar
[4]Baruch, G., Fibich, G. and Tsynkov, S., High-order numerical method for the nonlinear Helmholtz equation with material discontinuities in one space dimension, J. Comput. Phys., 227 (2007), 820–850.Google Scholar
[5]Baruch, G., Fibich, G. and Tsynkov, S., A high-order numerical method for the nonlinear Helmholtz equation in multidimensional layered media, J. Comput. Phys., 228 (2009), 3789–3815.CrossRefGoogle Scholar
[6]Bayliss, A., Goldstein, C. I. and Turkel, E., On accuracy conditions for the numerical computation of waves, J. Comput. Phys., 59(3) (1985), 396–404.Google Scholar
[7]Steven Britt, D., Tsynkov, Semyon and Turkel, Eli, A compact fourth order scheme for the Helmholtz equation in polar coordinates, J. Sci. Comput., 2010, accepted for publication.Google Scholar
[8]Caruthers, J. E., Steinhoff, J. S. and Engels, R. C., An optimal finite difference representation for a class of linear PDE’s with application to the Helmholtz equation, J. Comput. Acoust., 7(4) (1999), 245–252.Google Scholar
[9]Davis, Timothy A., Direct methods for sparse linear systems, Volume 2 of Fundamentals of Algorithms, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2006.Google Scholar
[10] Iain Duff, S. and Pralet, Stéphane, Strategies for scaling and pivoting for sparse symmetric indefinite problems, SIAM J. Matrix. Anal. Appl., 27(2) (2005), 313–340 (electronic).Google Scholar
[11]Erlangga, Y. and Turkel, E., Iterative schemes for high order compact discretizations to the exterior Helmholtz equation, M2AN, 2009, submitted for publication.Google Scholar
[12]Fibich, G. and Tsynkov, S. V., High-order two-way artificial boundary conditions for nonlinear wave propagation with backscattering, J. Comput. Phys., 171 (2001), 632–677.Google Scholar
[13]Fibich, G. and Tsynkov, S. V., Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions, J. Comput. Phys., 210(1) (2005), 183–224.Google Scholar
[14]Hagemann, M. and Schenk, O., Weighted matchings for preconditioning symmetric indefinite linear systems, SIAM J. Sci. Comput., 28(2) (2006), 403–420 (electronic).CrossRefGoogle Scholar
[15]Harari, I. and Turkel, E., Accurate finite difference methods for time-harmonic wave propagation, J. Comput. Phys., 119(2) (1995), 252–270.Google Scholar
[16]Kreiss, Heinz-Otto and Petersson, N. Anders, A second order accurate embedded boundary method for the wave equation with Dirichlet data, SIAM J. Sci. Comput., 27(4) (2006), 1141–1167 (electronic).CrossRefGoogle Scholar
[17]Kreiss, Heinz-Otto, Petersson, N. Anders and Yström, Jacob, Difference approximations for the second order wave equation, SIAM J. Numer. Anal., 40(5) (2002), 1940–1967 (electronic).Google Scholar
[18]Nabavi, M., Siddiqui, K. and Dargahi, J., A new 9-point sixth-order accurate compact finite-difference method for the Helmholtz equation, J. Sound. Vibrat., 307 (2007), 972–982.Google Scholar
[19]Nehrbass, J. W., Jevtic, J. O. and Lee, R., Reducing the phase error for finite-difference methods without incresing the order, IEEE Trans. Antenna. Propagat., 46 (1998), 1194–1201.CrossRefGoogle Scholar
[20]Nilsson, Stefan, Petersson, N. Anders, Sjögreen, Björn and Kreiss, Heinz-Otto, Stable difference approximations for the elastic wave equation in second order formulation, SIAM J. Numer. Anal., 45(5) (2007), 1902–1936.CrossRefGoogle Scholar
[21]Ryaben’kii, V. S., Method of Difference Potentials and Its Applications, Volume 30 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 2002.Google Scholar
[22]Singer, I. and Turkel, E., High-order finite difference methods for the Helmholtz equation, Comput. Methods. Appl. Mech. Engrg., 163(1-4) (1998), 343–358.Google Scholar
[23]Singer, I. and Turkel, E., Sixth-order accurate finite difference schemes for the Helmholtz equation, J. Comput. Acoust., 14(3) (2006), 339–351.Google Scholar
[24]Smirnov, V. I., Course of Higher Mathematics, Vol. IV, Part 1, Nauka, Moscow, Sixth revised edition, 1974, [in Russian].Google Scholar
[25]Smirnov, V. I., Course of Higher Mathematics, Vol. IV, Part 2, Nauka, Moscow, Sixth revised edition, 1981, [in Russian].Google Scholar
[26]Sutmann, G., Compact finite difference schemes of sixth order for the Helmholtz equation, J. Comput. Appl. Math., 203(1) (2007), 15–31.Google Scholar
[27]Tsynkov, S. V., On the definition of surface potentials for finite-difference operators, J. Sci. Comput., 18(2) (2003), 155–189.Google Scholar