Abstract
We present a simple method for simulating 2-D elastic waves in a model with free‐surface topography of polygonal shape, i.e., a continuous but irregular surface composed of line segments. Our method requires special treatment for each of the six specific cases involving line segments of various slopes as well as transition points between the sloping segments. For brevity, only nonnegatively sloping segments are specifically included. On an inclined free surface, vanishing stress conditions are implemented using a rotated coordinate system parallel to the inclined boundary. At transition points on the topography between line segments, we use a first‐order approximate boundary condition in a locally rotated coordinate system aligned with the bisector of the corner. As in the classical one‐sided explicit approximation scheme widely used for the flat free‐surface case, these extrapolation formulas are accurate to first order in spatial increment. Numerical tests indicate that the present scheme is stable over a range of Poisson’s ratios of practical interest (v > 0.3) for fairly complicated geometric shapes consisting of ridges and valleys with both steep and gentle slopes. Stability for complicated shapes enables us to study realistic problems for which the topography plays a significant role in shaping the wave field and for which analytical solutions are not generally available.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: 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.
