Convergence of Summation-by-Parts Finite Difference Methods for the Wave Equation
Tóm tắt
When using a finite difference method to solve a time dependent partial differential equation, the truncation error is often larger at a few grid points near a boundary or grid interface than in the interior. In computations, the observed convergence rate is often higher than the order of the large truncation error. In this paper, we develop techniques for analyzing this phenomenon, and particularly consider the second order wave equation. The equation is discretized by a finite difference operator satisfying a summation by parts property, and the boundary and grid interface conditions are imposed weakly by the simultaneous approximation term method. It is well-known that if the semi-discretized wave equation satisfies the determinant condition, that is the boundary system in Laplace space is nonsingular for all Re
$$(s)\ge 0$$
, two orders are gained from the large truncation error localized at a few grid points. By performing a normal mode analysis, we show that many common discretizations do not satisfy the determinant condition at
$$s=0$$
. We then carefully analyze the error equation to determine the gain in the convergence rate. The result shows that stability does not automatically imply a gain of two orders in the convergence rate. The precise gain can be lower than, equal to or higher than two orders, depending on the boundary condition and numerical boundary treatment. The accuracy analysis is verified by numerical experiments, and very good agreement is obtained.
Tài liệu tham khảo
Abarbanel, S., Ditkowski, A., Gustafsson, B.: On error bounds of finite difference approximations to partial differential equations—temporal behaviour and rate of convergence. J. Sci. Comput. 15, 79–116 (2000)
Appelö, D., Kreiss, G.: Application of a perfectly matched layer to the nonlinear wave equation. Wave Motion 44, 531–548 (2007)
Carpenter, M.H., Gottlieb, D., Abarbanel, S.: Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes. J. Comput. Phys. 111, 220–236 (1994)
Fernández, D.C.D.R., Hicken, J.E., Zingg, D.W.: Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations. Comput. Fluids 95, 171–196 (2014)
Grote, M.J., Schneebeli, A., Schötzau, D.: Discontinuous Galerkin finite element method for the wave equation. SIAM J. Numer. Anal. 44, 2408–2431 (2006)
Gustafsson, B.: The convergence rate for difference approximations to mixed initial boundary value problems. Math. Comp. 29, 396–406 (1975)
Gustafsson, B.: High Order Difference Methods for Time Dependent PDE. Springer, Berlin (2008)
Gustafsson, B., Kreiss, H.O., Oliger, J.: Time-Dependent Problems and Difference Methods. Wiley, New Jersey (2013)
Hagstrom, T., Hagstrom, G.: Grid stabilization of high-order one-sided differencing II: second-order wave equations. J. Comput. Phys. 231, 7907–7931 (2012)
Kreiss, H.O., Oliger, J.: Comparison of accurate methods for the integration of hyperbolic equations. Tellus 24, 199–215 (1972)
Kreiss, H.O., Scherer, G.: Finite element and finite difference methods for hyperbolic partial differential equations, Mathematical aspects of finite elements in partial differential equations. Symp. Proc. 33, 195–212 (1974)
Mattsson, K., Almquist, M.: A solution to the stability issues with block norm summation by parts operators. J. Comput. Phys. 253, 418–442 (2013)
Mattsson, K., Ham, F., Iaccarino, G.: Stable and accurate wave-propagation in discontinuous media. J. Comput. Phys. 227, 8753–8767 (2008)
Mattsson, K., Ham, F., Iaccarino, G.: Stable boundary treatment for the wave equation on second-order form. J. Sci. Comput. 41, 366–383 (2009)
Mattsson, K., Nordström, J.: Summation by parts operators for finite difference approximations of second derivatives. J. Comput. Phys. 199, 503–540 (2004)
Mattsson, K., Nordström, J.: High order finite difference methods for wave propagation in discontinuous media. J. Comput. Phys. 220, 249–269 (2006)
Nissen, A., Kreiss, G., Gerritsen, M.: Stability at nonconforming grid interfaces for a high order discretization of the Schrödinger equation. J. Sci. Comput. 53, 528–551 (2012)
Nissen, A., Kreiss, G., Gerritsen, M.: High order stable finite difference methods for the Schrödinger equation. J. Sci. Comput. 55, 173–199 (2013)
Petersson, N.A., Sjögreen, B.: Stable grid refinement and singular source discretiztion for seismic wave simulations. Comm. Comput. Phys. 8, 1074–1110 (2010)
Svärd, M.: A note on \(L^{\infty }\) bounds and convergence rates of summation-by-parts schemes. BIT. Numer. Math. 54, 823–830 (2014)
Svärd, M., Nordström, J.: On the order of accuracy for difference approximations of initial-boundary value problems. J. Comput. Phys. 218, 333–352 (2006)
Svärd, M., Nordström, J.: Review of summation-by-parts schemes for initial-boundary-value problems. J. Comput. Phys. 268, 17–38 (2014)
Taylor, N.W., Kidder, L.E., Teukolsky, S.A.: Spectral methods for the wave equation in second-order form. Phys. Rev. D 82, 024037 (2010)
Virta, K., Mattsson, K.: Acoustic wave propagation in complicated geometries and heterogeneous media. J. Sci. Comput. 61, 90–118 (2014)