求解Helmholtz方程的快速算法（英文）  

A Fast Propagation Method for the Helmholtz Equation

作　　者：冷伟[1] 

机构地区：[1]中国科学院科学与工程计算国家重点实验室,北京100190

出　　处：《工程数学学报》2015年第5期726-742,共17页Chinese Journal of Engineering Mathematics

基　　金：The National 863 Project of China(2012AA01A309)

摘　　要：本文提出了求解Helmholtz方程的一个新的快速算法.该算法是建立在有重叠区域的区域分解算法之上的.该算法首先对求解区域进行层次的区域分解,然后建立了各层次的子区域上的入射波到出射波的映射,最后通过层次的传播波的信息,得到Helmholtz方程的解.该方法具有计算复杂度小、适合大规模并行计算的优点,数值实验表明,该方法能够有效的并行求解有上亿自由度的二维Helmholtz方程.A fast method is proposed for solving the high frequency Helmholtz equation. The building block of the new fast method is an overlapping domain decomposition method for layered medium. In the new fast method, the computation domain is firstly decomposed hierarchically into many subdomains on different levels. Then the mapping from incident waves to out-going waves on all the subdomains are set up. Finally, the wave propagates on the subdomain boundaries on different levels to reach the solution to the Helmholtz equation. The new fast method is of low complexity, and suitable for parallel computing. Numerical experiments show that with the new fast method, 2D Helmholtz equations with half billion unknowns could be solved efficiently on massively parallel machines.

关 键 词：HELMHOLTZ方程 有限差分 快速算法 区域分解 完全匹配层 

分 类 号：O241.82[理学—计算数学] O241.6[理学—数学]