Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm
Yanbin Wang1, Yang Luo2, Yanfang Qin3, Xianbing Zhang1
1 Department of Geophysics, School of Earth and Space Sciences, Peking University, Beijing 100871, China 2 Department of Geosciences, Princeton University, 312 Guyot Hall, Princeton, NJ 08544, USA 3 Equipe de Ge´osciences Marines, Institut de Physique du Globe de Paris, 4 Place Jussieu, 75252 Paris Cedex 05, France
Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm
Yanbin Wang1, Yang Luo2, Yanfang Qin3, Xianbing Zhang1
1 Department of Geophysics, School of Earth and Space Sciences, Peking University, Beijing 100871, China 2 Department of Geosciences, Princeton University, 312 Guyot Hall, Princeton, NJ 08544, USA 3 Equipe de Ge´osciences Marines, Institut de Physique du Globe de Paris, 4 Place Jussieu, 75252 Paris Cedex 05, France
摘要We present a parallel hybrid algorithm based on pseudospectral method (PSM) and finite difference method (FDM) for two-dimensional (2-D) global SHwavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corresponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation.
Abstract:We present a parallel hybrid algorithm based on pseudospectral method (PSM) and finite difference method (FDM) for two-dimensional (2-D) global SHwavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corresponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation.
基金资助:the National Natural Science Foundation of China (Granted Nos. 41174034 and 40874020)
通讯作者:
Yanbin Wang
E-mail: ybwang@pku.edu.cn
引用本文:
Yanbin Wang, Yang Luo, Yanfang Qin, Xianbing Zhang. Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm[J]. 《地震学报》英文版, 2014, 27(3): 277-284.
Yanbin Wang, Yang Luo, Yanfang Qin, Xianbing Zhang. Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm. Earthquake Science, 2014, 27(3): 277-284.