Upwind schemes for the wave equation in second-order form

We develop new high-order accurate upwind schemes for the wave equation in second-order form. These schemes are developed directly for the equations in second-order form, as opposed to transforming the equations to a first-order hyperbolic system. The schemes are based on the solution to a local Riemann-type problem that uses d'Alembert's exact solution. We construct conservative finite difference approximations, although finite volume approximations are also possible. High-order accuracy is obtained using a space-time procedure which requires only two discrete time levels. The advantages of our approach include efficiency in both memory and speed together with accuracy and robustness. The stability and accuracy of the approximations in one and two space dimensions are studied through normal-mode analysis. The form of the dissipation and dispersion introduced by the schemes is elucidated from the modified equations. Upwind schemes are implemented and verified in one dimension for approximations up to sixth-order accuracy, and in two dimensions for approximations up to fourth-order accuracy. Numerical computations demonstrate the attractive properties of the approach for solutions with varying degrees of smoothness.

[1]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[2]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[3]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[4]  With Invariant Submanifolds,et al.  Systems of Conservation Laws , 2009 .

[5]  Leslie Greengard,et al.  An Integral Evolution Formula for the Wave Equation , 2000 .

[6]  Heinz-Otto Kreiss,et al.  Difference Approximations of the Neumann Problem for the Second Order Wave Equation , 2004, SIAM J. Numer. Anal..

[7]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[8]  Patrick Joly,et al.  Construction and analysis of higher order finite difference schemes for the 1D wave equation , 2000 .

[9]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[10]  William E. Schiesser,et al.  Linear and nonlinear waves , 2009, Scholarpedia.

[11]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[12]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[13]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

[14]  Gunilla Kreiss,et al.  A note on the effect of artificial viscosity on solutions of conservation , 1996 .

[15]  William J. Rider,et al.  On sub-linear convergence for linearly degenerate waves in capturing schemes , 2008, J. Comput. Phys..

[16]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[17]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[18]  Bruno Després,et al.  Finite volume transport schemes , 2008, Numerische Mathematik.

[19]  Leslie Greengard,et al.  Strongly consistent marching schemes for the wave equation , 2003 .

[20]  Leslie Greengard,et al.  High order marching schemes for the wave equation in complex geometry , 2004 .

[21]  W. Henshaw,et al.  Composite overlapping meshes for the solution of partial differential equations , 1990 .

[22]  R. C. Swanson.  On Some Numerical Dissipation Schemes , 1998 .

[23]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

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

[25]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[26]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[27]  H. Kreiss.  Stability theory for difference approximations of mixed initial boundary value problems. I , 1968 .

[28]  William D. Henshaw,et al.  A High-Order Accurate Parallel Solver for Maxwell's Equations on Overlapping Grids , 2005, SIAM J. Sci. Comput..

[29]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[30]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[31]  Jeffrey W. Banks,et al.  Numerical methods for solid mechanics on overlapping grids: Linear elasticity , 2010, J. Comput. Phys..

[32]  Bruno Després,et al.  Uniform Asymptotic Stability of Strang's Explicit Compact Schemes for Linear Advection , 2009, SIAM J. Numer. Anal..

[33]  H. Kreiss,et al.  Stability Theory of Difference Approximations for Mixed Initial Boundary Value Problems. II , 1972 .

[34]  Bruno Desprésaff n.  Finite volume transport schemes , 2008 .