×

zbMATH — the first resource for mathematics

A high-order discontinuous Galerkin method for wave propagation through coupled elastic-acoustic media. (English) Zbl 1427.74071
Summary: We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports \(h\)-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic-acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.

MSC:
74J10 Bulk waves in solid mechanics
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
74S99 Numerical and other methods in solid mechanics
65Y05 Parallel numerical computation
Software:
SPECFEM3D; p4est
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Ainsworth, M., Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods, Journal of computational physics, 198, 1, 106-130, (2004) · Zbl 1058.65103
[2] Ainsworth, M.; Monk, P.; Muniz, W., Dispersive and dissipative properties of discontinuous Galerkin finite element methods for the second-order wave equation, Journal of scientific computing, 27, 1-3, 5-40, (2006) · Zbl 1102.76032
[3] Bao, H.; Bielak, J.; Ghattas, O.; Kallivokas, L.F.; O’Hallaron, D.R.; Shewchuk, J.R.; Xu, J., Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers, Computer methods in applied mechanics and engineering, 152, 1-2, 85-102, (1998) · Zbl 0961.74056
[4] Appelö, D.; Petersson, N.A., A stable finite difference method for the elastic wave equation on complex geometries with free surfaces, Communications in computational physics, 5, 1, 84-107, (2009) · Zbl 1364.74016
[5] Lombard, B.; Piraux, J.; Gélis, C.; Virieux, J., Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves, Geophysical journal international, 172, 252-261, (2008)
[6] Lombard, B.; Piraux, J., Numerical treatment of two-dimensional interfaces for acoustic and elastic waves, Journal of computational physics, 195, 1, 90-116, (2004) · Zbl 1087.76079
[7] R.W. Hoberecht, A finite volume approach to modeling injury mechanisms of blast-induced traumatic brain injury, Master’s thesis, University of Washington, 2009.
[8] Voinovich, P.; Merlen, A.; Timofeev, E.; Takayama, K., A Godunov-type finite-volume scheme for unified solid-liquid elastodynamics on arbitrary two-dimensional grids, Shock waves, 13, 3, 221-230, (2003) · Zbl 1063.76064
[9] Komatitsch, D.; Barnes, C.; Tromp, J., Wave propagation near a fluid-solid interface: a spectral element approach, Geophysics, 65, 2, 623-631, (2000)
[10] Chaljub, E.; Capdeville, Y.; Vilotte, J., Solving elastodynamics in a fluid-solid heterogeneous sphere: a parallel spectral element approximation on non-conforming grids, Journal of computational physics, 187, 2, 457-491, (2003) · Zbl 1060.86003
[11] Chaljub, E.; Valette, B., Spectral element modelling of three-dimensional wave propagation in a self-gravitating Earth with an arbitrarily stratified outer core, Geophysical journal international, 158, 1, 131-141, (2004)
[12] Komatitsch, D.; Tsuboi, S.; Tromp, J., The spectral-element method in seismology, (), 205-228
[13] Riviere, B.; Wheeler, M.F., A discontinuous Galerkin method applied to nonlinear parabolic equations, (), 231-244 · Zbl 0946.65078
[14] Grote, M.J.; Schneebeli, A.; Schötzau, D., Discontinuous Galerkin finite element method for the wave equation, SIAM journal on numerical analysis, 44, 6, 2408-2431, (2006) · Zbl 1129.65065
[15] Chung, E.T.; Engquist, B., Optimal discontinuous Galerkin methods for wave propagation, SIAM journal on numerical analysis, 44, 5, 2131-2158, (2006) · Zbl 1158.65337
[16] De Basabe, J.D.; Sen, M.K.; Wheeler, M.F., The interior penalty discontinuous Galerkin method for elastic wave propagation: grid dispersion, Geophysical journal international, 175, 1, 83-93, (2008)
[17] S. Delcourte, L. Fezoui, N. Glinsky-Olivier, A high-order discontinuous Galerkin method for the seismic wave propagation, in: ESAIM: Proceedings, vol. 27, 2009, pp. 70-89. doi:10.1051/proc/2009020. · Zbl 1167.86305
[18] Dumbser, M.; Käser, M., An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes-II. the three-dimensional isotropic case, Geophysical journal international, 167, 1, 319, (2006)
[19] Käser, M.; Dumbser, M., A highly accurate discontinuous Galerkin method for complex interfaces between solids and moving fluids, Geophysics, 73, 3, T23-T35, (2008)
[20] Hesthaven, J.S.; Warburton, T., Nodal discontinuous Galerkin methods: algorithms, analysis, and applications, Texts in applied mathematics, vol. 54, (2008), Springer · Zbl 1134.65068
[21] LeVeque, R., Finite volume methods for hyperbolic problems, (2002), Cambridge University Press · Zbl 1010.65040
[22] Kopriva, D.A., Metric identities and the discontinuous spectral element method on curvilinear meshes, Journal of scientific computing, 26, 3, 301-327, (2006) · Zbl 1178.76269
[23] Kopriva, D.A., Implementing spectral methods for partial differential equations, (2009), Springer · Zbl 1172.65001
[24] Dumbser, M.; Käser, M., An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes-II. the three-dimensional isotropic case, Geophysical journal international, 167, 1, 319, (2006)
[25] Gonzalez, O.; Stuart, A.M., A first course in continuum mechanics, Cambridge texts in applied mathematics, (2008), Cambridge University Press
[26] Feng, K.-A.; Teng, C.-H.; Chen, M.-H., A pseudospectral penalty scheme for 2D isotropic elastic wave computations, Journal of scientific computing, 33, 3, 313-348, (2007) · Zbl 1127.74048
[27] Deville, M.; Fischer, P.; Mund, E., High-order methods for incompressible fluid flow, Cambridge monographs on applied and computational mathematics, vol. 9, (2002), Cambridge University Press · Zbl 1007.76001
[28] S.S. Collis, M. Heinkenschloss, Analysis of the streamline upwind/Petrov Galerkin method applied to the solution of optimal control problems, Technical Report TR02-01, Department of Computational and Applied Mathematics, Rice University, Houston, TX 77005-1892, 2002. <http://www.caam.rice.edu/heinken>.
[29] Kopriva, D.A.; Woodruff, S.L.; Hussaini, M.Y., Computation of electromagnetic scattering with a non-conforming discontinuous spectral element method, International journal for numerical methods in engineering, 53, 1, 105-122, (2002) · Zbl 0994.78020
[30] Kopriva, D.A.; Gassner, G., On the quadrature and weak form choices in collocation type discontinuous Galerkin spectral element methods, Journal of scientific computing, 44, 136-155, (2010) · Zbl 1203.65199
[31] Kennedy, C.A.; Carpenter, M.H.; Lewis, R.M., Low-storage, explicit runge – kutta schemes for the compressible navier – stokes equations, Applied numerical mathematics, 35, 3, 177-219, (2000) · Zbl 0986.76060
[32] Kaufman, A.A.; Levshin, A.L., Acoustic and elastic wave fields in geophysics, III, Methods in geochemistry and geophysics, vol. 39, (2005), Elsevier B.V.
[33] Achenbach, J.D., Wave propagation in elastic solids, North-holland series in applied mathematics and mechanics, vol. 16, (1973), North-Holland Publishing Company · Zbl 0268.73005
[34] Vinh, P.C.; Ogden, R.W., On formulas for the Rayleigh wave speed, Wave motion, 39, 3, 191-197, (2004) · Zbl 1163.74454
[35] C. Burstedde, L.C. Wilcox, O. Ghattas, : scalable algorithms for parallel adaptive mesh refinement on forests of octrees, SIAM Journal on Scientific Computing, submitted for publication. · Zbl 1230.65106
[36] Burstedde, C.; Ghattas, O.; Gurnis, M.; Tan, E.; Tu, T.; Stadler, G.; Wilcox, L.C.; Zhong, S., Scalable adaptive mantle convection simulation on petascale supercomputers, ()
[37] Burstedde, C.; Ghattas, O.; Gurnis, M.; Isaac, T.; Stadler, G.; Warburton, T.; Wilcox, L.C., Extreme-scale AMR, ()
[38] Dziewonski, A.M.; Anderson, D.L., Preliminary reference Earth model, Physics of the Earth and planetary interiors, 25, 4, 297-356, (1981), doi:10.1016/0031-9201(81)90046-7
[39] Carrington, L.; Komatitsch, D.; Laurenzano, M.; Tikir, M.M.; Michéa, D.; Goff, N.L.; Snavely, A.; Tromp, J., High-frequency simulations of global seismic wave propagation using SPECFEM3D GLOBE on 62K processors, ()
[40] Aki, K.; Richards, P.G., Quantitative seismology, (2002), University Science Books
[41] Performance applications programming interface (PAPI). <http://icl.cs.utk.edu/papi/>.
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.