×

Optimized double sweep Schwarz method by complete radiation boundary conditions. (English) Zbl 1359.65291

Summary: We present an optimized double sweep nonoverlapping Schwarz method for solving the Helmholtz equation in semi-infinite waveguides. The domain is decomposed into nonoverlapped layered subdomains along the axis of the waveguide and local wave propagation problems equipped with complete radiation conditions for high-order absorbing boundary conditions are solved forward and backward sequentially. For communication between subdomains, Neumann data of local solutions in one domain are transferred to the neighboring subdomain in the forward direction and Dirichlet data are exploited in the backward direction. The complete radiation boundary conditions enable us to not only minimize reflection coefficients for most important modes in an optimal way but also find Neumann data without introducing errors that would be produced if finite difference formulas were used for computing Neumann data. The convergence of the double sweep Schwarz method is proved and numerical experiments using it as a preconditioner are presented to confirm the convergence theory.

MSC:

65N55 Multigrid methods; domain decomposition for boundary value problems involving PDEs
65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation

Software:

deal.ii
PDF BibTeX XML Cite
Full Text: DOI

References:

[1] Gander, M. J.; Magoulès, F.; Nataf, F., Optimized Schwarz methods without overlap for the Helmholtz equation, SIAM J. Sci. Comput., 24, 1, 38-60, (2002), (electronic) · Zbl 1021.65061
[2] Gander, M. J., Optimized Schwarz methods, SIAM J. Numer. Anal., 44, 2, 699-731, (2006) · Zbl 1117.65165
[3] Boubendir, Y.; Antoine, X.; Geuzaine, C., A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation, J. Comput. Phys., 231, 2, 262-280, (2012) · Zbl 1243.65144
[4] Gander, M. J.; Halpern, L.; Magoulès, F., An optimized Schwarz method with two-sided Robin transmission conditions for the Helmholtz equation, Internat. J. Numer. Methods Fluids, 55, 2, 163-175, (2007) · Zbl 1125.65114
[5] Kim, S.; Zhang, H., Optimized Schwarz method with complete radiation transmission conditions for the Helmholtz equation in waveguides, SIAM J. Numer. Anal., 53, 3, 1537-1558, (2015) · Zbl 1327.78025
[6] Chen, Z.; Gander, M. J.; Zhang, H., On the relation between optimized Schwarz methods and source transfer, (Domain Decomposition Methods in Science and Engineering XXII, Lect. Notes Comput. Sci. Eng., vol. 104, (2016), Springer Berlin), 217-225 · Zbl 1339.65230
[7] Gander, M. J.; Nataf, F., An incomplete LU preconditioner for problems in acoustics, J. Comput. Acoust., 13, 3, 455-476, (2005) · Zbl 1189.76362
[8] Engquist, B.; Ying, L., Sweeping preconditioner for the Helmholtz equation: moving perfectly matched layers, Multiscale Model. Simul., 9, 2, 686-710, (2011) · Zbl 1228.65234
[9] Chen, Z.; Xiang, X., A source transfer domain decomposition method for Helmholtz equations in unbounded domain, SIAM J. Numer. Anal., 51, 4, 2331-2356, (2013) · Zbl 1285.65082
[10] Stolk, C. C., A rapidly converging domain decomposition method for the Helmholtz equation, J. Comput. Phys., 241, 0, 240-252, (2013) · Zbl 1349.65426
[11] Vion, A.; Geuzaine, C., Double sweep preconditioner for optimized Schwarz methods applied to the Helmholtz problem, J. Comput. Phys., 266, 171-190, (2014) · Zbl 1296.65169
[12] Magoulès, F.; Iványi, P.; Topping, B. H.V., Non-overlapping Schwarz methods with optimized transmission conditions for the Helmholtz equation, Comput. Methods Appl. Mech. Engrg., 193, 45-47, 4797-4818, (2004) · Zbl 1112.74444
[13] S. Kim, Analysis of complete radiation boundary conditions for the Helmholtz equation in perturbed waveguides, manuscript.
[14] S. Kim, T. Hagstrom, Complete radiation boundary conditions for the Helmholtz equation I: waveguides, manuscript. · Zbl 1412.65187
[15] Higdon, R. L., Absorbing boundary conditions for difference approximations to the multidimensional wave equation, Math. Comp., 47, 176, 437-459, (1986) · Zbl 0609.35052
[16] Higdon, R. L., Numerical absorbing boundary conditions for the wave equation, Math. Comp., 49, 179, 65-90, (1987) · Zbl 0654.65083
[17] Courant, R.; Hilbert, D., Methods of mathematical physics, vol. 1, (1953), Wiley-Interscience New York
[18] Goldstein, C. I., A finite element method for solving Helmholtz type equations in waveguides and other unbounded domains, Math. Comp., 39, 160, 309-324, (1982) · Zbl 0493.65046
[19] Bécache, E.; Bonnet-Ben Dhia, A.-S.; Legendre, G., Perfectly matched layers for the convected Helmholtz equation, SIAM J. Numer. Anal., 42, 1, 409-433, (2004) · Zbl 1089.76045
[20] Bendali, A.; Guillaume, P., Non-reflecting boundary conditions for waveguides, Math. Comp., 68, 225, 123-144, (1999) · Zbl 0907.35127
[21] Nataf, F.; Rogier, F.; de Sturler, E., Optimal interface conditions for domain decomposition methods, tech. rep. internal report no. 301, (1994), CMAP, Ecole Polytechnique
[22] Petrushev, P.; Popov, V., (Rational Approximation of Real Functions, Encyclopedia of Mathematics, vol. 28, (1987), Cambridge University Press Cambridge) · Zbl 0644.41010
[23] Bangerth, W.; Hartmann, R.; Kanschat, G., Deal.II—a general-purpose object-oriented finite element library, ACM Trans. Math. Software, 33, 4, 24, (2007) · Zbl 1365.65248
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.