Optimized Schwarz Method with Complete Radiation Transmission Conditions for the Helmholtz Equation in Waveguides

In this paper, we present a nonoverlapping optimized Schwarz algorithm with a high-order transmission condition for the Helmholtz equation posed in waveguides. We introduce the new high-order transmission conditions based on the complete radiation boundary conditions (CRBCs) that have been developed for high-order absorbing boundary conditions. We obtain the convergence rate of the algorithm in terms of reflection coefficients of CRBCs, which decrease exponentially with the order of CRBCs. It will be shown that damping parameters involved in the transmission conditions can be selected in an optimal way for enhancing the convergence of the Schwarz algorithm. Also, this algorithm can be employed efficiently for a preconditioner in GMRES implementations based on the substructured form as in [M. J. Gander, F. Magoules, and F. Nataf, SIAM J. Sci. Comput., 24 (2002), pp. 38--60]. Finally, numerical examples confirming the theory will be presented.

[1]  F. Magoulès,et al.  An optimized Schwarz method with two‐sided Robin transmission conditions for the Helmholtz equation , 2007 .

[2]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[3]  Christiaan C. Stolk,et al.  A rapidly converging domain decomposition method for the Helmholtz equation , 2012, J. Comput. Phys..

[4]  Martin J. Gander,et al.  Optimized Schwarz Methods without Overlap for the Helmholtz Equation , 2002, SIAM J. Sci. Comput..

[5]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[6]  Hui Zhang,et al.  Optimized double sweep Schwarz method by complete radiation boundary conditions , 2016, Comput. Math. Appl..

[7]  Lexing Ying,et al.  Sweeping Preconditioner for the Helmholtz Equation: Moving Perfectly Matched Layers , 2010, Multiscale Model. Simul..

[8]  F. Magoulès,et al.  Non-overlapping Schwarz methods with optimized transmission conditions for the H , 2004 .

[9]  Thomas Hagstrom,et al.  Numerical Experiments on a Domain Decomposition Algorithm for Nonlinear Elliptic Boundary Value Problems , 1988 .

[10]  L. Zschiedrich,et al.  Additive Schwarz method for scattering problems using the PML method at interfaces , 2007 .

[11]  Seungil Kim.  Analysis of the convected Helmholtz equation with a uniform mean flow in a waveguide with complete radiation boundary conditions , 2014 .

[12]  Frédéric Nataf,et al.  Optimal Interface Conditions for Domain Decomposition Methods , 1994 .

[13]  Tim Warburton,et al.  Complete Radiation Boundary Conditions: Minimizing the Long Time Error Growth of Local Methods , 2009, SIAM J. Numer. Anal..

[14]  W. Bangerth,et al.  deal.II—A general-purpose object-oriented finite element library , 2007, TOMS.

[15]  Patrick Joly,et al.  Domain Decomposition Method for Harmonic Wave Propagation : A General Presentation , 2000 .

[16]  VLADIMIR DRUSKIN,et al.  Near-Optimal Perfectly Matched Layers for Indefinite Helmholtz Problems , 2015, SIAM Rev..

[17]  Zhiming Chen,et al.  A Source Transfer Domain Decomposition Method for Helmholtz Equations in Unbounded Domain , 2013, SIAM J. Numer. Anal..

[18]  P. Petrushev,et al.  Rational Approximation of Real Functions , 1988 .

[19]  Anne-Sophie Bonnet-Ben Dhia,et al.  Perfectly Matched Layers for the Convected Helmholtz Equation , 2004, SIAM J. Numer. Anal..

[20]  Dan Givoli,et al.  Radiation boundary conditions for time-dependent waves based on complete plane wave expansions , 2010, J. Comput. Appl. Math..

[21]  Robert L. Higdon,et al.  Numerical absorbing boundary conditions for the wave equation , 1987 .

[22]  Martin J. Gander,et al.  Optimized Schwarz Methods , 2006, SIAM J. Numer. Anal..