TY - JOUR

T1 - Approximate momentum conservation for spatial semidiscretizations of semilinear wave equations

AU - Oliver, Marcel

AU - West, Matthew

AU - Wulff, Claudia

PY - 2004/5

Y1 - 2004/5

N2 - We prove that a standard second order finite difference uniform space discretization of the semilinear wave equation with periodic boundary conditions, analytic nonlinearity, and analytic initial data conserves momentum up to an error which is exponentially small in the stepsize. Our estimates are valid for as long as the trajectories of the full semilinear wave equation remain real analytic. The method of proof is that of backward error analysis, whereby we construct a modified equation which is itself Lagrangian and translation invariant, and therefore also conserves momentum. This modified equation interpolates the semidiscrete system for all time, and we prove that it remains exponentially close to the trigonometric interpolation of the semidiscrete system. These properties directly imply approximate momentum conservation for the semidiscrete system. We also consider discretizations that are not variational as well as discretizations on non-uniform grids. Through numerical example as well as arguments from geometric mechanics and perturbation theory we show that such methods generically do not approximately preserve momentum.

AB - We prove that a standard second order finite difference uniform space discretization of the semilinear wave equation with periodic boundary conditions, analytic nonlinearity, and analytic initial data conserves momentum up to an error which is exponentially small in the stepsize. Our estimates are valid for as long as the trajectories of the full semilinear wave equation remain real analytic. The method of proof is that of backward error analysis, whereby we construct a modified equation which is itself Lagrangian and translation invariant, and therefore also conserves momentum. This modified equation interpolates the semidiscrete system for all time, and we prove that it remains exponentially close to the trigonometric interpolation of the semidiscrete system. These properties directly imply approximate momentum conservation for the semidiscrete system. We also consider discretizations that are not variational as well as discretizations on non-uniform grids. Through numerical example as well as arguments from geometric mechanics and perturbation theory we show that such methods generically do not approximately preserve momentum.

UR - http://www.scopus.com/inward/record.url?scp=2942627242&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=2942627242&partnerID=8YFLogxK

U2 - 10.1007/s00211-003-0488-3

DO - 10.1007/s00211-003-0488-3

M3 - Article

AN - SCOPUS:2942627242

SN - 0029-599X

VL - 97

SP - 493

EP - 535

JO - Numerische Mathematik

JF - Numerische Mathematik

IS - 3

ER -