1887

Abstract

Summary

We present a multi-frontal hierarchically semi-separable solver to perform forward modeling of the 3D Helmholtz acoustic problem. Our frequency-domain solver combines two efficient approaches. First, it uses an optimal 27-point finite-difference scheme to decrease numerical dispersion and reduce required discretization of the model in terms of points per wavelength from 15 to about 4. Second, it uses a supernodal multi-frontal method based on low-rank approximation and hierarchically semi-separable (HSS) structure to improve performance, decrease memory usage and make it practical for realistic size 3D models required by full-waveform inversion. We validate our new solver using a 3D synthetic model and benchmark it against other codes.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201413101
2015-06-01
2024-04-26
Loading full text...

Full text loading...

References

  1. Chandrasekaran, S., Dewilde, P., Gu, M. and Somasunderam, N.
    [2010] On the numerical rank of the off-diagonal blocks of Schur complements of discretized elliptic PDEs. SIAM J. Matrix Anal., 2261–2290.
    [Google Scholar]
  2. Collino, F. and Tsogka, C.
    [2001] Application of the PML absorbing layer model to the linear elasto-dynamic problem in anisotropic heterogeneous media. Geophysics, 66(1), 294–307.
    [Google Scholar]
  3. George, J.A.
    [1973] Nested dissection of a regular finite element mesh. SIAM J. Numer. Anal., 10, 345–363.
    [Google Scholar]
  4. Jo, C., Shin, C. and Suh, J.
    [1996] An optimal 9-point finite-difference frequency-space 2-D scalar wave extrapolator. Geophysics, 61, 529–537
    [Google Scholar]
  5. Xia, J.
    [2012] Robust and efficient multifrontal solver for large discretized PDEs, High-Perform. In: Berry, M.W. et al. (Eds)Sci. Comput. Springer, 199–217.
    [Google Scholar]
  6. Rjasanow, S.
    [2002] Adaptive Cross Approximation of Dense Matrices. Proc. Int. Association for Boundary Element Methods, May28–30, 1–12.
    [Google Scholar]
  7. Wang, S., de Hoop, M.V., Xia, J. and Li, X.S.
    [2012] Massively parallel structured multifrontal solver for time-harmonic elastic waves in 3D anisotropic media. Proceedings of the Project Review, Geo-Mathematical Imaging Group (Purdue University, West Lafayette IN), 1, 97–121.
    [Google Scholar]
  8. Hackbusch, W.
    [1999] A Sparse Matrix Arithmetic Based on H-Matrices. Part I: Introduction to H-Matrices. Computing, 62(2), 89–108.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201413101
Loading
/content/papers/10.3997/2214-4609.201413101
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error