A Fast Propagation Method for the Helmholtz equation

A fast method is proposed for solving the high frequency Helmholtz equation. The building block of the new fast method is an overlapping source transfer domain decomposition method for layered medium, which is an extension of the source transfer domain decomposition method proposed by Chen and Xiang \cite{Chen2013a,Chen2013b}. The new fast method contains a setup phase and a solving phase. In the setup phase, the computation domain is decomposed hierarchically into many subdomains of different levels, and the mapping from incident traces to field traces on all the subdomains are set up bottom-up. In the solving phase, first on the bottom level, the local problem on the subdomains with restricted source is solved, then the wave propagates on the boundaries of all the subdomains bottom-up, at last the local solutions on all the subdomains are summed up top-down. The total computation cost of the new fast method is $O(n^{\frac{3}{2}} \log n)$ for 2D problem. Numerical experiments shows that with the new fast method, Helmholtz equations with half billion unknowns could be solved efficiently on massively parallel machines.

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

[2]  Zhiming Chen,et al.  A Source Transfer Domain Decomposition Method For Helmholtz Equations in Unbounded Domain Part II: Extensions , 2013 .

[3]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[4]  Haijun Wu,et al.  An improved pure source transfer domain decomposition method for Helmholtz equations in unbounded domain , 2015, 1505.06052.

[5]  Laurent Demanet,et al.  The method of polarized traces for the 2D Helmholtz equation , 2014, J. Comput. Phys..

[6]  B. Engquist,et al.  Sweeping preconditioner for the Helmholtz equation: Hierarchical matrix representation , 2010, 1007.4290.

[7]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[8]  M. Gander,et al.  AILU for Helmholtz problems: A new Preconditioner Based on the Analytic Parabolic Factorization.∗ , 2016 .

[9]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

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

[11]  Christophe Geuzaine,et al.  Double sweep preconditioner for optimized Schwarz methods applied to the Helmholtz problem , 2014, J. Comput. Phys..

[12]  Lexing Ying,et al.  Additive Sweeping Preconditioner for the Helmholtz Equation , 2015, Multiscale Model. Simul..

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

[14]  Weiying Zheng,et al.  Convergence of the Uniaxial Perfectly Matched Layer Method for Time-Harmonic Scattering Problems in Two-Layered Media , 2010, SIAM J. Numer. Anal..