Error analysis of the Strang splitting for the 3D semilinear wave equation with finite-energy data
Maximilian Ruff
TL;DR
The article analyzes a Strang-splitting time integrator for the 3D semilinear wave equation on $\mathbb{T}^3$ with finite-energy data, establishing sharp convergence rates for cubic and quartic nonlinearities under minimal regularity. The authors develop continuous and discrete Strichartz estimates on the torus, introduce a Fourier-filter to manage nonlinear interactions, and derive detailed semi-discrete error recursions, followed by fully discrete results using a Fourier pseudo-spectral spatial discretization. For $\alpha=3$, they obtain almost second-order convergence in $L^2$ and almost first-order in the energy norm, while for $\alpha=4$ the temporal order is reduced by a half. They also provide fully discrete error bounds combining temporal and spatial discretizations and corroborate the theory with numerical experiments that demonstrate sharpness and the role of the filter. Overall, the work advances low-regularity, structure-preserving time integration for nonlinear wave equations on bounded domains with finite-energy data, and offers practical guidance for implementing high-accuracy simulations on $\mathbb{T}^3$.
Abstract
We study a variant of the Strang splitting for the time integration of the semilinear wave equation under the finite-energy condition on the torus $\mathbb{T}^3$. In the case of a cubic nonlinearity, we show almost second-order convergence in $L^2$ and almost first-order convergence in $H^1$. If the nonlinearity has a quartic form instead, we show an analogous convergence result with an order reduced by 1/2. To our knowledge these are the best convergence results available for the 3D cubic and quartic wave equations under the finite-energy condition. Our approach relies on continuous- and discrete-time Strichartz estimates. We also make use of the integration and summation by parts formulas to exploit cancellations in the error terms. Moreover, error bounds for a full discretization using the Fourier pseudo-spectral method in space are given. Finally, we discuss a numerical example indicating the sharpness of our theoretical results.
