Non-standard FDTD for elastic wave simulation: Two-dimensional P-SV case

Geophysical Journal International

Volumn 178, Issue 1, 2009, Pages 282-302

  JafarGandomi, Arash ^a   Takenaka, Hiroshi ^a  

^a KYUSHU UNIVERSITY   (Japan)

Author keywords

Computational seismology; Numerical approximations and analysis; Wave propagation

Indexed keywords

ELASTIC WAVE; FINITE DIFFERENCE METHOD; NUMERICAL MODEL; POISSON RATIO; SEISMOLOGY; WAVE PROPAGATION;

EID: 67650102826     PISSN: 0956540X     EISSN: 1365246X     Source Type: Journal    
DOI: 10.1111/j.1365-246X.2009.04101.x     Document Type: Article

References (62)

1. Bogey, C. & Bailly, C., 2004. A family of low dispersive and low dissipative explicit schemes for flow and noise computations, J. Comput. Phys., 194, 194-214.
2. Bohlen, T., 2002. Parallel 3-D viscoelastic finite difference seismic modeling, Comput. Geosci., 28, 887-899.
3. Bouchon, M. & Aki, K., 1977. Discrete wave-number representation of seismic-source wave fields, Bull. seism. Soc. Am., 67, 259-277.
4. Cole, J.B., 1997. A high-accuracy realization of the Yee algorithm using non-standard finite-differences, IEEE Trans. Microw. Theory Tech., 45, 991-996.
5. Cole, J.B., 2002. High accuracy Yee algorithm based on nonstandard finite-difference: New developments and verifications, IEEE Trans. Antennas Propag., 50(9), 1185-1191.
6. Cole, J.B., 2005. High accuracy nonstandard finite-difference time-domain algorithms for computational electromagnetics: Applications to optics and photonics, in Advances in the Applications of Nonstandard Finite Difference Schemes, pp. 89-189, ed. Mickens, R.E., World Scientific, Singapore.
7. Cruz-Atienza, V. & Virieux, J., 2004. Dynamic rupture simulation of nonplanar faults with a finite-difference approach, Geophys. J. Int., 158, 939-954.
8. De Basabe, J.D. & Sen, M.K., 2007. Grid dispersion and stability criteria of some finite-element methods for acoustic and elastic wave equations, Geophysics, 72, T81-T95.
9. Furumura, T., 2005. Large-scale parallel simulation of seismic wave propagation and strong ground motions for the past and future earthquakes in Japan, J. Earth Sim., 3, 29-38.
10. Furumura, T. & Takenaka, H., 1996. 2.5-D modeling of elastic waves using the pseudospectral method, Geophys. J. Int., 124, 820-832.
11. Furumura, T., Kennett, B.L.N. & Takenaka, H., 1998. Parallel 3-D pseudospectral simulation of seismic wave propagation, Geophysics, 63, 279-288.
12. Geller, R.J. & Takeuchi, N., 1995. A new method for computing highly accurate DSM synthetic seismograms, Geophys. J. Int., 123, 449-470.
13. Geller, R.J. & Takeuchi, N, 1998. Optimally accurate second-order time-domain finite difference scheme for the elastic equation of motion: one-dimensional case, Geophys. J. Int., 135, 48-62.
14. Graves, R., 1996. Simulating seismic wave propagation in 3-D elastic media using staggered-grid finite differences, Bull. seism. Soc. Am., 86, 1091-1106.
15. Hansen, T.M. & Jacobsen, B.H., 2002. Efficient finite difference waveform modeling of selected phases using a moving zone, Comput. Geosci., 29, 819-826.
16. Hartzell, S. et al., 2006. Modeling and validation of a 3D velocity structure for the Santa Clara Valley, California, for seismic-wave simulations, Bull. seism. Soc. Am., 96, 1851-1881.
17. Hayashida, T., Takenaka, H. & Okamoto, T., 1999. Development of 2D and 3D codes of the velocity-stress staggered-grid finite-difference method for modeling seismic wave propagation (in Japanese with English abstract), Sci. Rep., Dept. Earth Planet. Sci., Kyushu Univ., 20, 99-110.
18. Hestholm, S., Moran, M., Ketcham, S., Anderson, T., Dillen, M. & McMechan, G., 2006. Effects of free-surface topography on moving-seismic-source modeling, Geophysics, 71, T159-T166.
19. Hixon, R., 2000. Prefactord small-stencil compact schemes, J. Comput. Phys., 165, 522-541.
20. Hixon, R., & Turkle, E., 2000. Compact implicit Mac-Cormack-type schemes with high accuracy, J. Comput. Phys., 158, 51-70.
21. Holberg, O., 1987. Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena, Geophys. Prospect., 35, 629-655.
22. Igel, H., Mora, P. & Riollet, B., 1995. Anisotropic wave propagation through finite-difference grids, Geophysics, 60, 1203-1216.
23. JafarGandomi, A. & Takenaka, H., 2007. Efficient FDTD algorithm for plane-wave simulation for vertically heterogeneous attenuative media, Geophysics, 72, H43-H53.
24. Jasträm, C. & Behle, A., 1993. Accurate finite-difference operators for modeling the elastic wave equation, Geophys. Prospect., 41, 453-458.
25. Kantartzis, N.V. & Tsiboukis, T.D., 2004. A higher order nonstandard FDTD-PML method for the advanced modeling of complex EMC problems in generalized 3-D curvilinear coordinates, IEEE Trans. Electromagn. Compat., 46(1), 2-11.
26. Kantartzis, N.V. & Tsiboukis, T.D., 2006. Higher Order FDTD Schemes for Waveguide and Antenna Structure, Morgan & Claypool, San Rafel, CA.
27. Kantartzis, N.V., Zygiridis, T.T. & Tsiboukis, T.D., 2002. A nonstandard higher-order FDTD algorithm for 3-D arbitrary and fractal-shaped antenna structure on general curvilinear lattices, IEEE Trans. Magn., 38, 737-740.
28. Kashiwa, T., Kudo, H., Sendo, Y., Ohtani, T. & Kanai, Y., 2002. The phase velocity error and stability condition of the three-dimensional nonstandard FDTD method, IEEE Trans. Magn., 38, 661-664.
29. Klimes, L., 1996. Accuracy of elastic finite differences in smooth media, PAGEOPH, 148, 39-76.
30. Kudo, H., Kashiwa, T. & Ohtani, T., 2002. The nonstandard FDTD method for three-dimensional acoustic analysis and its numerical dispersion and stability condition, Electron. Commun. Japan (Part 3), 85(9), 15-24.
31. Kvasnicka, M. & Zahradník, J., 1996. A combined kinematic-dynamic method for fast computations of the first-arriving waveforms, Report 4, Seismic Waves in Complex 3-D Structures, Dept. Geophys., Charles Univ., Prague, pp. 251-286.
32. Lele, S.K., 1992. Compact finite difference schemes with spectral-like resolution, J. Comput. Phys., 103, 16-42.
33. Levander, A., 1988. Fourth-order finite difference P-SV seismograms, Geophysics, 53, 1425-1436.
34. Madariaga, R., 1976. Dynamics of an expanding circular fault, Bull. seism. Soc. Am., 66, 163-182.
35. Mittet, R., 1994. Implementation of the Kirchhoff integral for elastic waves in staggered-grid modeling schemes, Geophysics, 59, 1894-1901.
36. Mizutani, H., Geller, R.J., & Takeuchi, N., 2000. Comparison of accuracy and efficiency of time-domain schemes for calculating synthetic seismograms, Phys. Earth planet. Inter., 119, 75-97.
37. Moczo, P., Kristek, J. & Halada, L., 2000. 3D fourth-order staggered-grid finite-difference schemes: Stability and grid dispersion, Bull. seism. Soc. Am., 90, 587-603.
38. Moczo, P., Kristek, J., Vavrycuk, V., Archuleta, R. & Halda, L., 2002. 3D heterogeneous staggered-grid finite-difference modeling of seismic motion with volume harmonic and arithmetic averaging of elastic moduli and densities, Bull. seism. Soc. Am., 92, 3042-3066.
39. Moczo, P., Robertsson, J.O. & Eisner, L., 2007. The finite-difference time-domain method for modeling of seismic wave propagation, in Advances in Wave Propagation in Heterogeneous Earth, Advances in Geophysics, Vol. 48, pp. 421-516, eds Wu, R.-S. & Maupin, V., Elsevier-Academic Press, London, UK.
40. Ohtani, T., Taguchi, K. & Kashiwa, T., 2005. New optimization parameters for the nonstandard FDTD method, Microw. Opt. Tech. Lett., 47(2), 161-163.
41. Olsen, K.B., 1994. Simulation of three-dimensional wave propagation in the Salt Lake Basin, PhD thesis. University of Utah, Salt Lake City.
42. Olsen, K.B., Day, S.M. & Bradley, C.R., 2003. Estimation of Q for long-period (>2 s) waves in the Los Angeles basin, Bull. seism. Soc. Am., 93, 627-638.
43. Olsen, K.B. et al., 2006. Strong shaking in Los Angeles expected from Southern San Andreas earthquake, Geophys. Res. Lett., 33, L07305.
44. Panaretos, A.H., Aberle, J.T. & Diaz, R.E., 2007. The effect of the 2-D Laplacian operator approximation on the performance of finite-difference time-domain schemes for Maxwell's equations, J. Comput. Phys., 227, 513-536.
45. Pitarka, A., 1999. 3D finite-difference modeling of seismic motion using staggered grids with nonuniform spacing, Bull. seis. Soc. Am., 89(1), 54-68.
46. Saenger, E. & Bohlen, T., 2004. Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid, Geophysics, 69, 583-591.
47. Saenger, E.H., Gold, N. & Shapiro, S.A., 2000. Modeling the propagation of elastic waves using a modified finite-difference grid, Wave Motion, 31, 77-92.
48. Takenaka, H & Kennett, B.L.N., 1996. A 2.5-D time-domain elastodynamic equation for plane-wave incidence, Geophys. J. Int., 125, 5-9.
49. Takenaka, H., Wang, Y. & Furumura, T., 1999. An efficient approach of the pseudospectral method for modeling of geometrical symmetric seismic wavefield, Earth Planet. Space, 51, 73-99.
50. Takenaka, H., Tanaka, H., Okamoto, T., Kennett, B.L.N., 2003, Quasi-cylindrical 2.5 D wave modeling for large-scale seismic surveys, Geophys. Res. Lett., 30, 2086.
51. Takeuchi, N. & Geller, R.J., 2000. Optimally accurate second order time-domain finite difference scheme for computing synthetic seismograms in 2-D and 3-D media, Phys. Earth planet. Inter., 119, 99-131.
52. Tam, C.K.W. & Webb, J.C., 1993. Dispersion-relation-preserving finite difference schemes for computational acoustics, J. Comput. Phys., 107, 262-281.
53. Tam, C.K.W. & Webb, J.C., 1994. Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation, J. Comput. Phys., 113, 122-133.
54. Toyokuni, G., Takenaka, H., Wang, Y. & Kennett, B.L.N., 2005. Quasi-spherical approach for seismic wave modeling in a 2-D slice of a global earth model with lateral heterogeneity, Geophys. Res. Lett., 32, L09305.
55. Virieux, J., 1984. SH-wave propagation in heterogeneous media: velocity-stress finite-difference method, Geophysics, 49, 1933-1957.
56. Virieux, J., 1986. P-SV wave propagation in heterogeneous media: velocity-stress finite-difference method, Geophysics, 51, 889-901.
57. Wang, S. & Teixeria, F.L., 2003a. A finite-difference time-domain algorithm optimized for arbitrary propagation angles, IEEE Trans. Antennas Propag., 51, 2456-2463.
58. Wang, S. & Teixeria, F.L., 2003b. A three-dimensional angle-optimized finite-difference time-domain algorithm, IEEE Trans. Microw. Theory Tech., 51, 811-817.
59. Wolfram, S., 2003. The Mathematica Book, 5th edn, Wolfram Media, available at http://documents.wolfram.com/mathematica.
60. Wu, W.J., Lines, L.R. & Lu, H.X., 1996. Analysis of higher-order, finite-difference schemes in 3-D reverse-time migration, Geophysics, 61, 845-856.
61. Yang, B. & Balanis, C.A., 2006. An isotropy-improved nonstandard finite-difference time-domain method, IEEE Trans. Antennas Propag., 54, 1935-1942.
62. Zahradník, J., 1995. Simple elastic finite-difference scheme, Bull. seism. Soc. Am., 85, 1879-1887.