Analysis of the dissipation and dispersion properties of the multi-domain Chebyshev pseudospectral method

The accuracy of the multi-domain Chebyshev pseudospectral method is investigated for wave propagation problems by examining the properties of the method in the wavenumber space theoretically in terms of dispersion and dissipation errors. For a number of ( N + 1 ) points in the subdomains used in the literature, with N typically between 8 to 32, significant errors can be obtained for waves discretized by more than π points per wavelength. The dispersion and dissipation errors determined from the analysis in the wavenumber space are found to be in good agreement with those obtained in test cases. Accuracy limits based on arbitrary criteria are proposed, yielding minimum resolutions of 7.7, 5.2 and 4.0 points per wavelength for N = 8 , 16 and 32 respectively. The numerical efficiency of the method is estimated, showing that it is preferable to choose N between 16 and 32 in practice. The stability of the method is also assessed using the standard fourth-order Runge-Kutta algorithm. Finally, 1-D and 2-D problems involving long-range wave propagation are solved to illustrate the dissipation and dispersion errors for short waves. The error anisotropy is studied in the 2-D case, in particular for a hybrid Fourier-Chebyshev configuration.

[1]  Gregor Gassner,et al.  A Comparison of the Dispersion and Dissipation Errors of Gauss and Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods , 2011, SIAM J. Sci. Comput..

[2]  D. Gottlieb,et al.  Numerical analysis of spectral methods , 1977 .

[3]  Bengt Fornberg,et al.  A practical guide to pseudospectral methods: Introduction , 1996 .

[4]  D. Gottlieb,et al.  Spectral methods for hyperbolic problems , 2001 .

[5]  Klaus Holliger,et al.  A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition , 2013, J. Comput. Phys..

[6]  John Thuburn,et al.  Dispersion analysis of the spectral element method , 2012 .

[7]  C. Bogey,et al.  A family of low dispersive and low dissipative explicit schemes for flow and noise computations , 2004 .

[8]  Abdollah A. Afjeh,et al.  Multidimensional optimization of finite difference schemes for Computational Aeroacoustics , 2008, J. Comput. Phys..

[9]  Qing Huo Liu,et al.  The 2.5-D multidomain pseudospectral time-domain algorithm , 2003 .

[10]  Chris Lacor,et al.  Dispersion and dissipation properties of the 1D spectral volume method and application to a p-multigrid algorithm , 2007, J. Comput. Phys..

[11]  J. Forssén,et al.  The extended Fourier pseudospectral time-domain method for atmospheric sound propagation. , 2010, The Journal of the Acoustical Society of America.

[12]  David A. Kopriva.  Spectral solution of acoustic wave-propagation problems , 1990 .

[13]  Zhi Qi,et al.  Theory of Bessel Functions of High Rank , 2015 .

[14]  L. Trefethen.  Spectral Methods in MATLAB , 2000 .

[15]  Woojae Seong,et al.  Pseudospectral time-domain method for pressure-release rough surface scattering using a surface transformation and image method , 2009, Numerical Algorithms.

[16]  Ameneh Taleei,et al.  A Chebyshev pseudospectral multidomain method for the soliton solution of coupled nonlinear Schrödinger equations , 2011, Comput. Phys. Commun..

[17]  Jan S. Hesthaven,et al.  Regular Article: Spectral Collocation Time-Domain Modeling of Diffractive Optical Elements , 1999 .

[18]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[19]  Xin Zhang,et al.  A Fourier Pseudospectral Method for Some Computational Aeroacoustics Problems , 2006 .

[20]  David I. Gottlieb,et al.  Optimal Decomposition of the Domain in Spectral Methods for Wave-Like Phenomena , 2000, SIAM J. Sci. Comput..

[21]  C. Bogey,et al.  Low-dissipation and low-dispersion fourth-order Runge–Kutta algorithm , 2006 .

[22]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[23]  Qing Huo Liu,et al.  The PSTD algorithm: A time-domain method requiring only two cells per wavelength , 1997 .

[24]  Ben-yu Guo,et al.  Fourier-Chebyshev pseudospectral method for two-dimensional vorticity equation , 1993 .

[25]  F Levassort,et al.  Modeling of piezoelectric transducers with combined pseudospectral and finite-difference methods. , 2008, The Journal of the Acoustical Society of America.

[26]  W. Pernice,et al.  Pseudo-spectral time-domain simulation of the transmission and the group delay of photonic devices , 2008 .

[27]  Edip Baysal,et al.  Forward modeling by a Fourier method , 1982 .

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

[29]  Steven G. Johnson,et al.  A Modified Split-Radix FFT With Fewer Arithmetic Operations , 2007, IEEE Transactions on Signal Processing.

[30]  M. Hornikx,et al.  THE EXTENDED FOURIER PSEUDOSPECTRAL TIME-DOMAIN (PSTD) METHOD FOR FLUID MEDIA WITH DISCONTINUOUS PROPERTIES , 2010 .

[31]  M. Y. Hussaini,et al.  An Analysis of the Discontinuous Galerkin Method for Wave Propagation Problems , 1999 .

[32]  Christopher K. W. Tam,et al.  Radiation Boundary Condition and Anisotropy Correction for Finite Difference Solutions of the Helmholtz Equation , 1994 .

[33]  Alfred Behle,et al.  Multi-domain Chebyshev-Fourier method for the solution of the equations of motion of dynamic elasticity , 1992 .

[34]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[35]  J. Carcione,et al.  A 2D Chebyshev differential operator for the elastic wave equation , 1996 .

[36]  Chris Lacor,et al.  On the Stability and Accuracy of the Spectral Difference Method , 2008, J. Sci. Comput..