Compact High Order Accurate Schemes for the Three Dimensional Wave Equation

We construct a family of compact fourth order accurate finite difference schemes for the three dimensional scalar wave (d’Alembert) equation with constant or variable propagation speed. High order accuracy is of key importance for the numerical simulation of waves as it reduces the dispersion error (i.e., the pollution effect). The schemes that we propose are built on a stencil that has only three nodes in any coordinate direction or in time, which eliminates the need for auxiliary initial or boundary conditions. These schemes are implicit in time and conditionally stable. A particular scheme with the maximum Courant number can be chosen within the proposed class. The inversion at the upper time level is done by FFT for constant coefficients and multigrid for variable coefficients, which keeps the overall complexity of time marching comparable to that of a typical explicit scheme.

[1]  A. Brandt.  Rigorous quantitative analysis of multigrid, I: constant coefficients two-level cycle with L 2 -norm , 1994 .

[2]  A. Bayliss,et al.  On accuracy conditions for the numerical computation of waves , 1985 .

[3]  Semyon Tsynkov,et al.  A Compact Fourth Order Scheme for the Helmholtz Equation in Polar Coordinates , 2010, J. Sci. Comput..

[4]  Jean-Pierre Bérenger,et al.  Perfectly Matched Layer (PML) for Computational Electromagnetics , 2007, PML for Computational Electromagnetics.

[5]  C. Greif,et al.  Closed-form multigrid smoothing factors for lexicographic Gauss–Seidel , 2012 .

[6]  Hui Liang,et al.  Stability of theta-schemes in the numerical solution of a partial differential equation with piecewise continuous arguments , 2010, Appl. Math. Lett..

[7]  I. Babuska,et al.  Dispersion and pollution of the FEM solution for the Helmholtz equation in one, two and three dimensions , 1999 .

[8]  Semyon Tsynkov,et al.  Numerical Simulation of Time-Harmonic Waves in Inhomogeneous Media using Compact High Order Schemes , 2011 .

[9]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[10]  S. Armfield,et al.  A compact fourth-order spatial discretisation applied to the Navier-Stokes equations , 2016 .

[11]  Semyon Tsynkov,et al.  A method of boundary equations for unsteady hyperbolic problems in 3D , 2018, J. Comput. Phys..

[12]  Sébastien Imperiale,et al.  Introduction and study of fourth order theta schemes for linear wave equations , 2013, J. Comput. Appl. Math..

[13]  Semyon Tsynkov,et al.  Optimization of Acoustic Source Strength in the Problems of Active Noise Control , 2003, SIAM J. Appl. Math..

[14]  Stefan A. Sauter,et al.  Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers? , 1997, SIAM Rev..

[15]  Yuanyuan Wang,et al.  A novel unsplit perfectly matched layer for the second-order acoustic wave equation. , 2014, Ultrasonics.

[16]  Semyon Tsynkov,et al.  Numerical solution of the wave equation with variable wave speed on nonconforming domains by high-order difference potentials , 2018, J. Comput. Phys..

[17]  M. Grote,et al.  Efficient PML for the Wave Equation , 2010, 1001.0319.

[18]  E. Turkel,et al.  Operated by Universities Space Research AssociationAccurate Finite Difference Methods for Time-harmonic Wave Propagation* , 1994 .

[19]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[20]  D. Fishelov,et al.  A new fourth-order compact scheme for the Navier-Stokes equations in irregular domains , 2016, Comput. Math. Appl..

[21]  B. Zeghmati,et al.  An efficient parallel high-order compact scheme for the 3D incompressible Navier–Stokes equations , 2017 .

[22]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[23]  Jinhua Yu,et al.  An efficient complex‐frequency shifted‐perfectly matched layer for second‐order acoustic wave equation , 2014 .

[24]  N. K. Kulman,et al.  Method of difference potentials and its applications , 2001 .

[25]  Yiping,et al.  COMPACT FOURTH-ORDER FINITE DIFFERENCE SCHEMES FOR HELMHOLTZ EQUATION WITH HIGH WAVE NUMBERS , 2008 .

[26]  S. B. Atienza-Samols,et al.  With Contributions by , 1978 .

[27]  Semyon Tsynkov,et al.  Compact 2D and 3D sixth order schemes for the Helmholtz equation with variable wave number , 2013, J. Comput. Phys..

[28]  I. Singer,et al.  High-order finite difference methods for the Helmholtz equation , 1998 .

[29]  Semyon Tsynkov,et al.  A High Order Compact Time/Space Finite Difference Scheme for the Wave Equation with Variable Speed of Sound , 2018, J. Sci. Comput..