Optimal blended spectral-element operators for acoustic wave modeling

Spectral-element methods, based on high-order polynomials, are among the most commonly used techniques for computing accurate simulations of wave propagation phenomena in complex media. However, to retain computational efficiency, very high order polynomials cannot be used and errors such as numerical dispersion and numerical anisotropy cannot be totally avoided. In the present work, we devise an approach for reducing such errors by considering modified discrete wave operators. We analyze consistent and lumped operators together with blended operators (weighted averages of consistent and lumped operators). Furthermore, using the operator-blending approach and a novel dispersion analysis method, we develop optimal spectral-element operators that have increased numerical accuracy, without resorting to very high order operators. The new operators are faster and computationally more efficient than consistent operators. Our approach is based on the tensor product decomposition of the element matrices into 1D f...

[1]  Alexander Graham,et al.  Kronecker Products and Matrix Calculus: With Applications , 1981 .

[2]  K. Marfurt.  Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[3]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[4]  Isaac Fried,et al.  Finite element mass matrix lumping by numerical integration with no convergence rate loss , 1975 .

[5]  B. Fornberg.  The pseudospectral method: Comparisons with finite differences for the elastic wave equation , 1987 .

[6]  Patrick Joly,et al.  Higher-order finite elements with mass-lumping for the 1D wave equation , 1994 .

[8]  W. A. Mulder,et al.  Higher-order triangular and tetrahedral finite elements with mass lumping for solving the wave equation , 1999 .

[9]  J. Z. Zhu,et al.  The finite element method , 1977 .

[10]  M. Jensen,et al.  HIGH CONVERGENCE ORDER FINITE ELEMENTS WITH LUMPED MASS MATRIX , 1996 .

[11]  W. A. Mulder,et al.  Spurious modes in finite-element discretizations of the wave equation may not be all that bad , 1999 .

[12]  T. Belytschko,et al.  Dispersion analysis of finite element semidiscretizations of the two‐dimensional wave equation , 1982 .

[13]  M. Christon.  The influence of the mass matrix on the dispersive nature of the semi-discrete, second-order wave equation , 1999 .

[14]  Géza Seriani,et al.  Numerical simulation of interface waves by high‐order spectral modeling techniques , 1992 .

[15]  Géza Seriani,et al.  LOW AND HIGH ORDER FINITE ELEMENT METHOD: EXPERIENCE IN SEISMIC MODELING , 1994 .

[16]  Géza Seriani,et al.  3-D large-scale wave propagation modeling by spectral element method on Cray T3E multiprocessor , 1998 .

[17]  Henk A. van der Vorst,et al.  Alternatives to the Rayleigh Quotient for the Quadratic Eigenvalue Problem , 2003, SIAM J. Sci. Comput..

[18]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[19]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[20]  P. Pinsky,et al.  Complex wavenumber Fourier analysis of the p-version finite element method , 1994 .

[21]  Rosemary A. Renaut,et al.  Accuracy, Resolution, and Stability Properties of a Modified Chebyshev Method , 2002, SIAM J. Sci. Comput..

[22]  Slimane Adjerid,et al.  A posteriori finite element error estimation for second-order hyperbolic problems , 2002 .

[23]  Géza Seriani,et al.  Spectral element method for acoustic wave simulation in heterogeneous media , 1994 .

[24]  H. Saunders.  Book Reviews : The Finite Element Method (Revised): O.C. Zienkiewicz McGraw-Hill Book Co., New York, New York , 1980 .

[25]  Shen R. Wu,et al.  Lumped mass matrix in explicit finite element method for transient dynamics of elasticity , 2006 .

[26]  Ashley F. Emery,et al.  Accuracy in modeling the acoustic wave equation with Chebyshev spectral finite elements , 1997 .

[27]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[28]  Géza Seriani,et al.  High-order spectral element method for elastic wave modeling: 62nd Ann , 1991 .

[29]  Géza Seriani,et al.  A Parallel Spectral Element Method for Acoustic Wave Modeling , 1997 .