Abstract
We extend the Robertson’s image method (Robertsson, 1996) to a 3D finite-difference time-domain (FDTD) wave propagation in the presence of surface topography. Robertson method is accurate for a horizontal free surface and easily implemented on a standard staggered grid (SSG), but has some artefact for a topographic surface if it is staircase represented. We use a surface-conforming mapping to transfer the physical topographic domain (x,y,z) to a regular Cartesian domain (ξ,κ,η), which is equivalent to a curvilinear coordinate system. We take the velocity-stress first-order equation form for the (anisotropic) elastic wave equation, and discretize the transferred Cartesian domain on a fully staggered grid (FSG). Following the Robertsson procedure, we put all the particle velocity zero above the free surface, the velocity derivative with respect to η at the free surface is replaced with velocity derivatives with respect to ξ and κ, and anti-symmetrize the traction component normal to the curvilinear surface. We present the numerical results to show the efficiency of this extension.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.