Stable and Accurate Second-Order Formulation of the Shifted Wave Equation
Ken Mattsson 1*, Florencia Parisi 21 Department of Information Technology, Uppsala University, P.O. Box 337, S-751 05 Uppsala, Sweden.
2 FaMAF, Universidad Nacional de Cordoba, Medina Allende S/N Ciudad Universitaria, 5000 Cordoba, Argentina.
Received 20 August 2008; Accepted (in revised version) 4 May 2009
Available online 16 July 2009
High order finite difference approximations are derived for a one-dimensional model of the shifted wave equation written in second-order form. The domain is discretized using fully compatible summation by parts operators and the boundary conditions are imposed using a penalty method, leading to fully explicit time integration. This discretization yields a strictly stable and efficient scheme. The analysis is verified by numerical simulations in one-dimension. The present study is the first step towards a strictly stable simulation of the second-order formulation of Einstein's equations in three spatial dimensions.AMS subject classifications: 35L05, 35L20, 65N06, 65N12, 83C05
Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: High-order finite difference methods, wave equation, numerical stability, second derivatives, Einstein's equations.
Email: firstname.lastname@example.org (K. Mattsson), email@example.com (F. Parisi)