MPI- and CUDA- implementations of modal finite difference method for P-SV wave propagation modeling

Authors

  • Reza Khajavi Earthquake Research Center, Ferdowsi University of Mashhad, Iran
Abstract:

Among different discretization approaches, Finite Difference Method (FDM) is widely used for acoustic and elastic full-wave form modeling. An inevitable deficit of the technique, however, is its sever requirement to computational resources. A promising solution is parallelization, where the problem is broken into several segments, and the calculations are distributed over different processors. For the present FD routines, however, such parallelization technique inevitably needs domain-decomposition and inter-core data exchange, due to the coupling of the governing equations. In this study, a new FD-based procedure for seismic wave modeling, named as ‘Modal Finite Difference Method (MFDM)” is introduced, which deals with the simulation in the decoupled modal space; thus, neither domain-decomposition nor inter-core data exchange is anymore required, which greatly simplifies parallelization for both MPI- and CUDA implementations over CPUs and GPUs. With MFDM, it is also possible to simply cut off less-significant modes and run the routine for just the important ones, which will effectively reduce computation and storage costs. The efficiency of the proposed MFDM is shown by some numerical examples.

Upgrade to premium to download articles

Sign up to access the full text

Already have an account?login

similar resources

mpi- and cuda- implementations of modal finite difference method for p-sv wave propagation modeling

among different discretization approaches, finite difference method (fdm) is widely used for acoustic and elastic full-wave form modeling. an inevitable deficit of the technique, however, is its sever requirement to computational resources. a promising solution is parallelization, where the problem is broken into several segments, and the calculations are distributed over different processors. ...

full text

Mimetic Finite Difference Modeling Of 2D Elastic P-SV Wave Propagation

Recently, efforts have been made by Kristek et al.(2002) toward the implementation of a zero-traction boundary condition for an elastic medium by using high-order staggered-grid finite-difference modeling and avoiding symmetry conditions, vacuum formulations, or, other approaches that require grid points located beyond the physical boundary (ghost points). In this work, a new set of numerical d...

full text

Modal Fields Calculation Using the Finite Difference Beam Propagation Method

A method is described to construct modal fields for an arbitrary oneor two-dimensional refractive index structure. An arbitrary starting field is propagated along a complex axis using the slowly varying envelope approximation (SVEA). By choosing suitable values for the step-size, one mode is maximally increased in amplitude on propagating, until convergence has been obtained. For the calculatio...

full text

Compact finite difference modeling of 2-D acoustic wave propagation

We present two fourth-order compact finite difference (CFD) discretizations of the velocity-pressure formulation of the acoustic wave equation in 2-D rectangular grids. The first method uses standard implicit CFD on nodal meshes and requires solving tridiagonal linear systems along each grid line, while the second scheme employs a novel set of mimetic CFD operators for explicit differentiation ...

full text

Wave propagation across acoustic / Biot’s media: a finite-difference method

Numerical methods are developed to simulate the wave propagation in heterogeneous 2D fluid / poroelastic media. Wave propagation is described by the usual acoustics equations (in the fluid medium) and by the low-frequency Biot’s equations (in the porous medium). Interface conditions are introduced to model various hydraulic contacts between the two media: open pores, sealed pores, and imperfect...

full text

P-SC/ wave propagation in heterogeneous media: Velocity-stress finite-difference method

I present a finite-difference method for modeling P-SV wave propagation in heterogeneous media. This is an extension of the method I previously proposed for modeling SH-wave propagation by using velocity and stress in a discrete grid. The two components of the velocity cannot be defined at the same node for a complete staggered grid: the stability condition and the P-wave phase velocity dispers...

full text

My Resources

Save resource for easier access later

Save to my library Already added to my library

{@ msg_add @}


Journal title

volume 2  issue 2

pages  185- 202

publication date 2016-07-01

By following a journal you will be notified via email when a new issue of this journal is published.

Hosted on Doprax cloud platform doprax.com

copyright © 2015-2023