×

High-accuracy finite-difference schemes for solving elastodynamic problems in curvilinear coordinates within multiblock approach. (English) Zbl 1325.74157

Summary: We propose highly accurate finite-difference schemes for simulating wave propagation problems described by linear second-order hyperbolic equations. The schemes are based on the summation by parts (SBP) approach modified for applications with violation of input data smoothness. In particular, we derive and implement stable schemes for solving elastodynamic anisotropic problems described by the Navier wave equation in complex geometry. To enhance potential of the method, we use a general type of coordinate transformation and multiblock grids. We also show that the conventional spectral element method (SEM) can be treated as the multiblock finite-difference method whose blocks are the SEM cells with SBP operators on GLL grid.

MSC:

74S10 Finite volume methods applied to problems in solid mechanics
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
74J05 Linear waves in solid mechanics

References:

[1] Appelö, D.; Petersson, N. A., A stable finite difference method for the elastic wave equation on complex geometries with free surfaces, Commun. Comput. Phys., 5, 84-107 (2009) · Zbl 1364.74016
[2] Cagniard, L., Réflexion et réfraction des ondes séismiques progressives (1939), Gauthier-Villars: Gauthier-Villars Paris · JFM 65.1490.02
[3] Carpenter, M. H.; Gottlieb, D.; Abarbanel, S., The stability of numerical boundary treatments for compact high-order finite-difference schemes, J. Comput. Phys., 108, 2, 272-295 (1994) · Zbl 0791.76052
[4] Despres, B., Uniform asymptotic stability of Strang’s explicit compact schemes for linear advection, SIAM J. Numer. Anal., 47, 3956-3976 (2009) · Zbl 1211.65111
[5] Dovgilovich, L., High-order FD method on curvilinear grids for anisotropic elastodynamic simulations, (75th EAGE Conference & Exhibition (2013)), Extended abstracts
[6] Dovgilovich, L. E.; Sofronov, I. L., High-order finite-difference method of calculating wavefields in anisotropic media, Seismic Technol., 2, 24-30 (2013), (in Russian)
[7] Gassner, G. J., A skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods, SIAM J. Sci. Comput., 35, 3, A1233-A1253 (2013) · Zbl 1275.65065
[8] Gorshkov, A. G.; Medvedskii, A. L.; Rabinskii, L. N.; Tarlakovskii, D. V., Waves in Solids, 1-472 (2004), FizMatLit: FizMatLit Moscow, (in Russian)
[9] Holberg, O., 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 (1987)
[10] de Hoop, A. T., A modification of Cagniard’s method for solving seismic pulse problems, Appl. Sci. Res. Sect. B, 8, 4, 349-356 (1960) · Zbl 0100.44208
[11] Komatitsch, D.; Vilotte, J. P., The spectral element method: an efficient tool to simulate the seismic response of 2D and 3D geological structures, Bull. Seismol. Soc. Am., 88, 368-392 (1998) · Zbl 0974.74583
[12] Kozdon, J. E.; Dunham, E. M.; Nordström, J., Simulation of dynamic earthquake ruptures in complex geometries using high-order finite difference methods, J. Sci. Comput., 55, 1, 92-124 (2013) · Zbl 1278.86009
[13] Kreiss, H.-O.; Scherer, G., Finite element and finite difference methods for hyperbolic partial differential equations, (Mathematical Aspects of Finite Elements in Partial Differential Equations (1974), Academic Press), 195-212 · Zbl 0355.65085
[14] Lamb, H., On the propagation of tremors over the surface of an elastic solid, Philos. Trans. R. Soc. Lond. A, 203, 1-42 (1904) · JFM 34.0859.02
[15] Lekhnitskii, S. G., Theory of Elasticity of an Anisotropic Body, 1-431 (1981), Mir · Zbl 0467.73011
[16] Lisitsa, V.; Vishnevskiy, D., Lebedev scheme for the numerical simulation of wave propagation in 3D anisotropic elasticity, Geophys. Prospect., 58, 4, 619-635 (2010)
[17] Mattsson, K., Summation by parts operators for finite difference approximations of second-derivatives with variable coefficients, J. Sci. Comput., 51, 3, 650-682 (2012) · Zbl 1252.65055
[18] Mattsson, K.; Ham, F.; Iaccarino, G., Stable boundary treatment for the wave equation on second-order form, J. Sci. Comput., 41, 3, 366-383 (2009) · Zbl 1203.65145
[19] Mattsson, K.; Nordström, J., Summation by parts operators for finite difference approximations of second derivatives, J. Comput. Phys., 199, 503-540 (2004) · Zbl 1071.65025
[20] Nordström, J.; Mattsson, K.; Svärd, M., Stable and accurate artificial dissipation, J. Sci. Comput., 21, 1, 57-79 (2004) · Zbl 1085.76050
[21] Petersson, N. A.; Sjögreen, B., An energy absorbing far-field boundary condition for the elastic wave equation, Commun. Comput. Phys., 6, 483-508 (2009) · Zbl 1364.74019
[22] Ryaben’kii, V. S., The method of intrinsic boundary conditions in the theory of difference boundary value problems, Usp. Mat. Nauk, 26, 3, 105-160 (1971) · Zbl 0257.35076
[23] Ryaben’kii, V. S., Method of Difference Potentials and Its Applications, 1-556 (2001), Springer
[24] Saenger, E. H.; Gold, N.; Shapiro, S. A., Modeling the propagation of elastic waves using a modified finite-difference grid, Wave Motion, 31, 77-92 (2000) · Zbl 1074.74648
[25] Samarski, A. A., The Theory of Difference Schemes, 1-786 (2001), CRC Press · Zbl 0971.65076
[26] Sofronov, I. L., Differential part of transparent boundary conditions for certain hyperbolic systems of second-order equations, Dokl. Math., 79, 3, 412-414 (2009) · Zbl 1183.35195
[27] Sofronov, I. L.; Zaitsev, N. A.; Daryin, A.; Voskoboinikova, O.; Dovgilovich, L., Multi-block FD method for 3D geophysical simulation with explicit representation of sub-horizontal interfaces, (74th EAGE Conference & Exhibition (2012)), Extended abstracts
[28] Strand, B., Summation by parts for finite difference approximations for d/dx, J. Comput. Phys., 110, 47-67 (1994) · Zbl 0792.65011
[29] Virieux, J., P-SV wave propagation in heterogeneous media: velocity-stress finite-difference method, Geophysics, 51, 889-901 (1986)
[30] (2013)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.