Based on existing direct and imaging methods of a staggered finite-difference scheme, an improved algorithm for staggered finite-difference is proposed to implement rugged topographic free boundary conditions. This me...Based on existing direct and imaging methods of a staggered finite-difference scheme, an improved algorithm for staggered finite-difference is proposed to implement rugged topographic free boundary conditions. This method assumes that the free surface can be implemented with horizontal and vertical free surface segments and their corners; the free surface passes through the grid points of shear stress components, instead of the normal stress components. Imaging is carried out for stress components in both horizontal and vertical directions, thus increasing the accuracy. To update particle- velocities, imaging and updating are first performed in the horizontal direction, and then in the vertical direction. The numerical results for elastic flat horizontal free surface with the imaging method and those for flat free surfaces of various slope angles with the proposed method are compared, and are shown to be in good agreement. The advantage of the proposed method is that only the stresses are dealt with in implementing the free surface into the staggered algorithm, which improves computation efficiency.展开更多
基金The first author wishes to thank Prof.John McDonald in Curtin University of Technology of Australia for reviewing this paper.This work was partly supported by the National Natural Science Foundation of China(Grant No.10134020).
文摘Based on existing direct and imaging methods of a staggered finite-difference scheme, an improved algorithm for staggered finite-difference is proposed to implement rugged topographic free boundary conditions. This method assumes that the free surface can be implemented with horizontal and vertical free surface segments and their corners; the free surface passes through the grid points of shear stress components, instead of the normal stress components. Imaging is carried out for stress components in both horizontal and vertical directions, thus increasing the accuracy. To update particle- velocities, imaging and updating are first performed in the horizontal direction, and then in the vertical direction. The numerical results for elastic flat horizontal free surface with the imaging method and those for flat free surfaces of various slope angles with the proposed method are compared, and are shown to be in good agreement. The advantage of the proposed method is that only the stresses are dealt with in implementing the free surface into the staggered algorithm, which improves computation efficiency.
