Volume 24 Number 3
  • E-ISSN: 1365-2478



The manner in which boundary conditions are approximated and introduced into finite difference schemes may have an important influence on the stability and accuracy of the results. The standard von Neumann condition for stability applies only for points which are not in the vicinity of the boundaries. This stability condition does not take into consideration the effects caused by introducing the boundary conditions to the scheme.

Working on elastic media with free stress boundary conditions we found that the boundary approximation gives rise to serious stability problems especially for regions with high Poisson's ratio. In order to detect these effects apriori and to analyse them, we have used a more elaborate procedure for checking the stability of the scheme which takes into consideration the boundary conditions. It is based on studying a locally spaced time propagating matrix which governs the time‐space behavior of a small region of the grid which includes free surface points.

By using this procedure a better insight into the nature of instability caused by the approximations to the boundary conditions was gained which led us to a new stable approximation for the free surface boundary conditions.


Article metrics loading...

Loading full text...

Full text loading...


  1. Alterman, Z., and Kornfeld, P., 1963, Propagation of a pulse within a sphere, J. Acoust. Soc. Am. 35, 1649–1662.
    [Google Scholar]
  2. Alterman, Z., and Karal, F., 1968, Propagation of elastic waves in layered media by finite differences methods, Bull. Seism. Soc. Am. 58, 367–398.
    [Google Scholar]
  3. Alterman, Z., and Loewenthal, D., 1969, Algebraic expressions for the impulsive motion of an elastic half space, Israel J. of Technology7, 495–504.
    [Google Scholar]
  4. Alterman, Z., and Rotenberg, A., 1969, Seismic waves in a quarter plane, Bull. Seism. Soc. Am. 59, 347–368.
    [Google Scholar]
  5. Boore, D. M., 1970, Finite‐difference solution to the equations of elastic wave propagation, with application to Love waves over dipping interface, Ph. D. Thesis, Mass. Inst. Tech.
  6. Cagniard, L., 1962, Reflection and refraction of progressive seismic waves. Translated and revised by E. A. Flinn and C. H. Dix, McGraw‐Hill, New York .
    [Google Scholar]
  7. Garvin, W. W., 1956, Exact transient solution of the buried line source problem, Proc. Roy. Soc. A. 234, 528–541.
    [Google Scholar]
  8. Godunov, S. K., and Ryabenkii, V. S., 1964, Introduction to the theory of difference schemes, Interscience Publishers, New York .
    [Google Scholar]
  9. Ilan, A., Ungar, A., and Alterman, Z., 1975, An improved representation of boundary conditions in finite difference schemes for seismological problems, Geophys. J. Roy. Astr. Soc. 43, 727–747.
    [Google Scholar]
  10. Kreiss, H. O., 1968, Stability theory for difference approximations of mixed initial boundary value problems, Math. Comp. 22, 703–714.
    [Google Scholar]
  11. O'Brien, G. G., Hyman, M. A., and Kaplan, S., 1950, A study of the numerical solution of partial differential equation, J. Math. Physics29, 223.
    [Google Scholar]
  12. Osher, S., 1969, Systems of difference equations with general homogeneous boundary conditions, Trans. Amer. Math. Soc. 137, 177–201.
    [Google Scholar]
  13. Richtmyer, R. D., and Morton, K. W., 1967, Difference methods for initial‐value problems, Interscience Publishers, New York , London , Sydney .
    [Google Scholar]
  • Article Type: Research Article
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