Wave propagation in poroelastic media: A velocity‐stress staggered‐grid finite‐difference method with perfectly matched layers

2003 ◽  
Author(s):  
Dong‐Hoon Sheen ◽  
Kagan Tuncay ◽  
Chang‐Eob Baag ◽  
Peter J. Ortoleva
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Staggered Grid ◽  
Perfectly Matched Layers ◽  
Poroelastic Media ◽  
Matched Layers

Numerical Modeling of Wave Propagation in Poroelastic Media with a Staggered-Grid Finite-Difference Method

2013 ◽  
Vol 275-277 ◽  
pp. 612-617
Author(s):  
Wen Sheng Zhang ◽  
Li Tong
Keyword(s):  
Numerical Modeling ◽  
Wave Propagation ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Staggered Grid ◽  
Poroelastic Medium ◽  
Numerical Computations ◽  
Poroelastic Media ◽  
Biot’S Equations

In this paper, wave propagation in poroelastic medium is simulated with a staggered-grid finite-difference method. The formulation is discretized based on the second-order Biot’s equations rather than the corresponding velocity-stress form. In order to eliminate boundary reflections, the PML method is applied. Numerical computations are implemented and the results show the correctness and effectiveness of the schemes presented in this paper.

A High-Order Finite-Difference Method on Staggered Curvilinear Grids for Seismic Wave Propagation Applications with Topography

2021 ◽  
Author(s):  
Ossian O’Reilly ◽  
Te-Yang Yeh ◽  
Kim B. Olsen ◽  
Zhifeng Hu ◽  
Alex Breuer ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Vertical Axis ◽  
Cartesian Grid ◽  
Staggered Grid ◽  
Graphic Processing Units ◽  
Curvilinear Grids ◽  
Curvilinear Grid

ABSTRACT We developed a 3D elastic wave propagation solver that supports topography using staggered curvilinear grids. Our method achieves comparable accuracy to the classical fourth-order staggered grid velocity–stress finite-difference method on a Cartesian grid. We show that the method is provably stable using summation-by-parts operators and weakly imposed boundary conditions via penalty terms. The maximum stable timestep obeys a relationship that depends on the topography-induced grid stretching along the vertical axis. The solutions from the approach are in excellent agreement with verified results for a Gaussian-shaped hill and for a complex topographic model. Compared with a Cartesian grid, the curvilinear grid adds negligible memory requirements, but requires longer simulation times due to smaller timesteps for complex topography. The code shows 94% weak scaling efficiency up to 1014 graphic processing units.

Sign in / Sign up

Export Citation Format

Share Document