The present article is concerned with the implementation of a recent semi-analytical method referred to as fractional reduced differential transform method （FRDTM） for computation of approximate solution of time-fra...The present article is concerned with the implementation of a recent semi-analytical method referred to as fractional reduced differential transform method （FRDTM） for computation of approximate solution of time-fractional gas dynamics equation （TFGDE） arising in shock fronts. In this approach, the fractional derivative is described in the Caputo sense. Four numeric experiments have been carried out to confirm the validity and the efficiency of the method. It is found that the exact or a closed approximate analytical solution of a fractional nonlinear differential equations arising in allied science and engineering can be obtained easily. Moreover, due to its small size of calculation contrary to the other analytical approaches while dealing with a complex and tedious physical problems arising in various branches of natural sciences and engineering, it is very easy to implement.展开更多
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 present article is concerned with the implementation of a recent semi-analytical method referred to as fractional reduced differential transform method （FRDTM） for computation of approximate solution of time-fractional gas dynamics equation （TFGDE） arising in shock fronts. In this approach, the fractional derivative is described in the Caputo sense. Four numeric experiments have been carried out to confirm the validity and the efficiency of the method. It is found that the exact or a closed approximate analytical solution of a fractional nonlinear differential equations arising in allied science and engineering can be obtained easily. Moreover, due to its small size of calculation contrary to the other analytical approaches while dealing with a complex and tedious physical problems arising in various branches of natural sciences and engineering, it is very easy to implement.
基金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.
