Iterative convergence of boundary-volume integral equation method  

作　　者：Gengxin Yu Liyun Fu 

机构地区：[1]Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China

出　　处：《Earthquake Science》2011年第5期391-400,共10页地震学报（英文版）

基　　金：supported by the National Natural Science Foundation of China(Nos. 41130418 and 40925013);the National Basic Research Program(973 Program)(No.2009CB219403)

摘　　要：The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.

关 键 词：boundary-volume integral equation generalized Lipmann-Schwinger integral equation GMRES method diagonal preconditioner 

分 类 号：P315.01[天文地球—地震学]