The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed...The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed first, which is the solution of displacement fields for elastic half space with circle canyon under output plane harmonic line loading at horizontal surface. Then the integral equation for determining the unknown forces in the problem can be changed into the algebraic one and solved numerically so that crack DSIF can be determined. Last when the medium parameters are altered, the influence on the crack DSIF is discussed partially with the displacement between circle canyon and crack.展开更多
When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly M...When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly Matched Layers （PML） are deduced and their Crank-Nicolson scheme are presented in this paper. We use the second-, sixth-, and tenth-order finite difference and pseudo-spectral algorithms to compute the spatial derivatives. Two numerical models, a homogeneous isotropic medium and a multi-layer model with a cave, are designed to investigate how the absorbing boundary width and the algorithms determine PML effects. Numerical results show that, for PML, the low-order finite difference algorithms have fairly good absorbing effects when the absorbing boundary is thin, whereas, high-order algorithms always have good absorption when the boundary is thick. Finally, we discuss the reflection coefficient and point out its shortcomings, which is why we use the SNR to quantitatively scale the PML effects,展开更多
For linear partial differential equation 〔 2t 2-a 2P( x)〕 m u=f(x,t), where m1,X∈R n,t∈R 1, the author gives the analytic solution of the initial value problem using the operators sh(tP( x) 1/2 )...For linear partial differential equation 〔 2t 2-a 2P( x)〕 m u=f(x,t), where m1,X∈R n,t∈R 1, the author gives the analytic solution of the initial value problem using the operators sh(tP( x) 1/2 )P( x) 1/2 . By representing the operators with integrals, explicit solutions are obtained with an integral form of a given function.展开更多
The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD...The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD method can effectively suppress the numerical dispersion when coarse grids are used.In this paper,the ONAD method is extended to solve the 2-dimensional SH-wave equation in the spherical coordinates.To investigate the accuracy and the efficiency of the ONAD method,we compare the numerical results calculated by the ONAD method and other methods for both the homogeneous model and the inhomogeneous IASP91 model.The comparisons indicate that the ONAD method for solving the SH-wave equation in the spherical coordinates has the advantages of less numerical dispersion,small memory requirement for computer codes,and fast calculation.As an application,we use the ONAD method to simulate the SH-wave propagating between the Earth's surface and the core-mantle boundary(CMB).Meanwhile,we investigate the SH-wave propagating in the mantle through analyzing the wave-field snapshots in different times and synthetic seismograms.展开更多
文摘The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed first, which is the solution of displacement fields for elastic half space with circle canyon under output plane harmonic line loading at horizontal surface. Then the integral equation for determining the unknown forces in the problem can be changed into the algebraic one and solved numerically so that crack DSIF can be determined. Last when the medium parameters are altered, the influence on the crack DSIF is discussed partially with the displacement between circle canyon and crack.
基金supported jointly by the 973 Program (Grant No.2007CB209505)the National Natural Science Fund (Grant No.40704019,40674061)+1 种基金the School Basic Research Fund of Tsinghua University (JC2007030)PetroChina Innovation Fund (Grant No.060511-1-1)
文摘When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly Matched Layers （PML） are deduced and their Crank-Nicolson scheme are presented in this paper. We use the second-, sixth-, and tenth-order finite difference and pseudo-spectral algorithms to compute the spatial derivatives. Two numerical models, a homogeneous isotropic medium and a multi-layer model with a cave, are designed to investigate how the absorbing boundary width and the algorithms determine PML effects. Numerical results show that, for PML, the low-order finite difference algorithms have fairly good absorbing effects when the absorbing boundary is thin, whereas, high-order algorithms always have good absorption when the boundary is thick. Finally, we discuss the reflection coefficient and point out its shortcomings, which is why we use the SNR to quantitatively scale the PML effects,
文摘For linear partial differential equation 〔 2t 2-a 2P( x)〕 m u=f(x,t), where m1,X∈R n,t∈R 1, the author gives the analytic solution of the initial value problem using the operators sh(tP( x) 1/2 )P( x) 1/2 . By representing the operators with integrals, explicit solutions are obtained with an integral form of a given function.
基金supported by National Science Fund of Distinguished Young Scholars of China(Grant No. 40725012)40821002)National Natural Science Foundation of China (Grant No. 41074073)
文摘The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD method can effectively suppress the numerical dispersion when coarse grids are used.In this paper,the ONAD method is extended to solve the 2-dimensional SH-wave equation in the spherical coordinates.To investigate the accuracy and the efficiency of the ONAD method,we compare the numerical results calculated by the ONAD method and other methods for both the homogeneous model and the inhomogeneous IASP91 model.The comparisons indicate that the ONAD method for solving the SH-wave equation in the spherical coordinates has the advantages of less numerical dispersion,small memory requirement for computer codes,and fast calculation.As an application,we use the ONAD method to simulate the SH-wave propagating between the Earth's surface and the core-mantle boundary(CMB).Meanwhile,we investigate the SH-wave propagating in the mantle through analyzing the wave-field snapshots in different times and synthetic seismograms.
