Abstract:ABSTRACT The curvilinear-grid finite-difference method (FDM), which uses curvilinear coordinates to discretize the nonplanar interface geometry, is extended to simulate acoustic and seismic-wave propagation across the fluid–solid interface at the sea bottom. The coupled acoustic velocity-pressure and elastic velocity-stress formulation that governs wave propagation in seawater and solid earth is expressed in curvilinear coordinates. The formulation is solved on a collocated grid by alternative applications of forward and backward MacCormack finite difference within a fourth-order Runge–Kutta temporal integral scheme. The shape of a fluid–solid interface is discretized by a curvilinear grid to enable a good fit with the topographic interface. This good fit can obtain a higher numerical accuracy than the staircase approximation in the conventional FDM. The challenge is to correctly implement the fluid–solid interface condition, which involves the continuity of tractions and the normal component of the particle velocity, and the discontinuity (slipping) of the tangent component of the particle velocity. The fluid–solid interface condition is derived for curvilinear coordinates and explicitly implemented by a domain-decomposition technique, which splits a grid point on the fluid–solid interface into one grid point for the fluid wavefield and another one for the solid wavefield. Although the conventional FDM that uses effective media parameters near the fluid–solid interface to implicitly approach the boundary condition conflicts with the fluid–solid interface condition. We verify the curvilinear-grid FDM by conducting numerical simulations on several different models and compare the proposed numerical solutions with independent solutions that are calculated by the Luco-Apsel-Chen generalized reflection/transmission method and spectral-element method. Besides, the effects of a nonplanar fluid–solid interface and fluid layer on wavefield propagation are also investigated in a realistic seafloor bottom model. The proposed algorithm is a promising tool for wavefield propagation in heterogeneous media with a nonplanar fluid–solid interface.

An Improved Acoustic/elastic Interface Approach for 2D Staggered Grid Finite-Difference Modeling of Rayleigh Waves in the Presence of Surface Topography

Numerical Investigation of Rayleigh-wave Propagation on Topography Surface

Modeling Of Rayleigh-Wave Propagation In Three-Dimensional Media With Topography Using Staggered-Grid Finite Difference Scheme

Numerical simulation of elastic wave based on the staggered grid finite difference method

An Improved Vacuum Formulation for 2D Finite-Difference Modeling of Rayleigh Waves Including Surface Topography and Internal Discontinuities

An Efficient and Accurate Strategy of Numerical Simulation of 2D Rayleigh Waves under Arbitrary Undulating Surface Conditions

Oscillatory Responses of a Spatially Random Seabed in the Vicinity of Structures: Validation and Application of the Integrated CFD-SFEM

Simulate the elastic wavefields in media with an irregular surface topography based on staggered grid finite difference

Numerical Investigation of Implementation of Air-Earth Boundary by Acoustic-Elastic Boundary Approach

Viscoelastic wave finite-difference modeling in the presence of topography with adaptive free-surface boundary condition

The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order, Staggered-Grid, Finite Difference Scheme

Acoustic Reverse Time Migration and Perfectly Matched Layer in Boundary-Conforming Grids by Elliptic Method

A More Accurate Boundary Spectral Method Based on Subdivision Surface

Boundary Conditions and Case Analysis of Numerical Modeling of Rayleigh Wave

Strongly imposing the free surface boundary condition for wave equations with finite difference operators

A Three-dimensional Immersed Boundary Method for Accurate Simulation of Acoustic Wavefields with Complex Surface Topography

An Overset Mesh-Free Finite-Difference Method for Seismic Modeling Including Surface Topography

An Efficient Immersed Free Surface Boundary Method for 3-D Scalar Seismic Waves Finite-Difference Modeling in Presence of Topography

A non-balanced staggered-grid finite-difference scheme for the first-order elastic wave-equation modeling

3D Seismic-Wave Modeling with a Topographic Fluid–Solid Interface at the Sea Bottom by the Curvilinear-Grid Finite-Difference Method

2-D Poroelastic Wave Modelling with a Topographic Free Surface by the Curvilinear Grid Finite-Difference Method