The pseudospectral method for solving vorticity equations on spherical surface is discussed. An interpolation procedure, which is different from the usual onesl is proposed. Based on such an interpolation, the pseudos...The pseudospectral method for solving vorticity equations on spherical surface is discussed. An interpolation procedure, which is different from the usual onesl is proposed. Based on such an interpolation, the pseudospectral scheme is constructed. Its generalized stability and convergence are analyzed rigorously. The theoretical analysis and computational skills can also be applied to other nonlinear partial differential equations defined on spherical surface.展开更多
We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical ...We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical coordinates and the model is taken through a great circle of the Earth. Spatial derivatives in the wave equations are calculated in the wavenumber domain by multiplication, and the transformation between spatial and wavenumber domains is performed via fast Fourier transformation. Because of the high accuracy and high speed of the Fourier pseudospectral method, the scheme enables us to calculate a short-wavelength global SH-wavefield with accurate waveforms and arrival times for models with heterogeneities that can be approximated as azimuthally symmetric. Comparing with two-dimensional simulation methods based on an axisymmetric model, implementing the seismic source in the present scheme is more convenient. We calculated the global SH-wavefield for the preliminary reference Earth model to identify the generation, reflection and refraction of various seismic phases propagating in the Earth. Applications to a heterogeneous global model with low-velocity perturbation above the core-mantle boundary were conducted to analyze the effect of lateral heterogeneity on global SH-wave propagation.展开更多
文摘The pseudospectral method for solving vorticity equations on spherical surface is discussed. An interpolation procedure, which is different from the usual onesl is proposed. Based on such an interpolation, the pseudospectral scheme is constructed. Its generalized stability and convergence are analyzed rigorously. The theoretical analysis and computational skills can also be applied to other nonlinear partial differential equations defined on spherical surface.
基金supported by National Natural Science Foundation of China (Grant Nos. 40874020, 40474012 and 40821062)National R&D Special Fund for Public Welfare Industry (Grant No. 20070804)
文摘We present a scheme to simulate SH-wave propagation in a whole-Earth model with arbitrary lateral heterogeneities employing the Fourier pseudospectral method. Wave equations are defined in two-dimensional cylindrical coordinates and the model is taken through a great circle of the Earth. Spatial derivatives in the wave equations are calculated in the wavenumber domain by multiplication, and the transformation between spatial and wavenumber domains is performed via fast Fourier transformation. Because of the high accuracy and high speed of the Fourier pseudospectral method, the scheme enables us to calculate a short-wavelength global SH-wavefield with accurate waveforms and arrival times for models with heterogeneities that can be approximated as azimuthally symmetric. Comparing with two-dimensional simulation methods based on an axisymmetric model, implementing the seismic source in the present scheme is more convenient. We calculated the global SH-wavefield for the preliminary reference Earth model to identify the generation, reflection and refraction of various seismic phases propagating in the Earth. Applications to a heterogeneous global model with low-velocity perturbation above the core-mantle boundary were conducted to analyze the effect of lateral heterogeneity on global SH-wave propagation.
