1887
Volume 35 Number 6
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Conventional finite‐difference operators for numerical differentiation become progressively inaccurate at higher frequencies and therefore require very fine computational grids. This problem is avoided when the derivatives are computed by multiplication in the Fourier domain. However, because matrix transpositions are involved, efficient application of this method is restricted to computational environments where the complete data volume required by each computational step can be kept in random access memory.

To circumvent these problems a generalized numerical dispersion analysis for wave equation computations is developed. Operators for spatial differentiation can then be designed by minimizing the corresponding peak relative error in group velocity within a spatial frequency band. For specified levels of maximum relative error in group velocity ranging from 0.03% to 3%, differentiators have been designed that have the largest possible bandwidth for a given operator length.

The relation between operator length and the required number of grid points per shortest wavelength, for a required accuracy, provides a useful starting point for the design of cost‐effective numerical schemes. To illustrate this, different alternatives for numerical simulation of the time evolution of acoustic waves in three‐dimensional inhomogeneous media are investigated. It is demonstrated that algorithms can be implemented that require fewer arithmetic and I/O operations by orders of magnitude compared to conventional second‐order finite‐difference schemes to yield results with a specified minimum accuracy.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.1987.tb00841.x
2006-04-27
2024-04-25
Loading full text...

Full text loading...

References

  1. Alford, R.M., Kelly, K.R. and Boore, D.M.1974. Accuracy of finite‐difference modeling of the acoustic wave equation, Geophysics39, 834–842.
    [Google Scholar]
  2. Brigham, E. O.1974. The Fast Fourier Transform, Prentice‐Hall, Inc.
    [Google Scholar]
  3. Cerjan, C, Kosloff, D., Kosloff, R. and Reshef, M.1985. A non‐reflecting boundary condition for discrete acoustic and elastic wave equations, Geophysics50, 705–708.
    [Google Scholar]
  4. Dablain, M.A.1986. The application of high‐order differencing to the scalar wave equation, Geophysics51, 54–66.
    [Google Scholar]
  5. Gazdag, J.1973. Numerical convective schemes based on accurate computation of space derivatives, Journal of Computational Physics13, 100–113.
    [Google Scholar]
  6. Gazdag, J.1981. Modeling of the acoustic wave equation with transform methods, Geophysics46, 854–859.
    [Google Scholar]
  7. Gibbs, A.J.1969. On the frequency domain responses of causal digital filters, Ph.D. thesis, Department of Electrical Engineering, University of Wisconsin.
  8. Gupta, R.N.1966. Reflection of plane elastic waves from transition layers with arbitrary variation of velocity and density, Bulletin of the Seismological Society of America56, 633–642.
    [Google Scholar]
  9. Hockney, R.W.1984. Performance of parallel computers, in Proceedings of the NATO Workshop on High Speed Computations, J.Kowalik (ed.), NATO ASI Series, F‐7, Springer‐Verlag, Inc.
    [Google Scholar]
  10. Johnson, O.G.1984. Three‐dimensional wave equation computations on vector computers, Proceedings of the IEEE 18–9219, 90–95.
  11. Jordan, T.L.1982. A guide to parallel computation and some Cray‐1 experiences, in Parallel Computations, G.Rodrigue (ed.), Academic Press, Inc.
    [Google Scholar]
  12. Kelly, K.R., Ward, R.W., Treitel, S. and Alford, R.M.1976. Synthetic seismograms: a finite‐difference approach, Geophysics41, 2–27.
    [Google Scholar]
  13. Kosloff, D. and Baysal, E.1982. Forward modeling by a Fourier method, Geophysics47, 1402–1412.
    [Google Scholar]
  14. Kosloff, D., Reshef, M. and Loewenthal, D.1984. Elastic wave calculations by the Fourier method, Bulletin of the Seismological Society of America74, 875–891.
    [Google Scholar]
  15. Kreiss, H.O. and Oliger, J.1972. Comparison of accurate methods for the integration of hyperbolic equations, Tellus24, 199–215.
    [Google Scholar]
  16. Marfurt, K.J.1984. Accuracy of finite‐difference and finite‐element modeling of the scalar and elastic wave equations, Geophysics49, 533–549.
    [Google Scholar]
  17. Mc Clellan, J.H., Parks, T.W. and Rabiner, L.R.1973. A computer program for designing optimum FIR linear phase digital filters, IEEE Transactions on Audio and ElectroacousticsAU‐21, 506–526.
    [Google Scholar]
  18. Moore, J.J., Garbow, B.S. and Hillstrøm, K.E.1980. User guide for MINPACK‐l, Argonne National Laboratory, ANL‐80–74.
    [Google Scholar]
  19. Neves, K.W.1984. Vectorization of scientific software, in Proceedings of the NATO Workshop on High Speed Computations, J.Kowalik (ed.), NATO ASI Series, F‐7, Springer‐Verlag, Inc.
    [Google Scholar]
  20. Orzag, S.A.1980. Spectral methods for problems in complex geometries, Journal of Computational Physics37, 70–92.
    [Google Scholar]
  21. Rabiner, L.R. and Schafer, R.W.1974. On the behavior of minimax relative error FIR digital differentiators, Bell System Technical Journal53, 333–346.
    [Google Scholar]
  22. Rabiner, L.R. and Steiglitz, K.1970. The design of wide‐band recursive and nonrecursive digital differentiators, IEEE Transactions on Audio and Electroacoustics AU‐18, 204–209.
  23. Swarztrauber, P.N.1982. Vectorizing the FFTs, in Parallel Computations, G.Rodrigue (ed.), Academic Press, Inc.
    [Google Scholar]
  24. Trefethen, L.N.1982. Group velocity in finite difference schemes, Society of Industrial and Applied Mathematics Review24, 113–136.
    [Google Scholar]
  25. Wolfe, M.A.1978. Numerical Methods For Unconstrained Optimization, Van Nostrand Reinhold.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.1987.tb00841.x
Loading
  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

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