Abstract:In this paper, we propose and analyse a novel third-order low-regularity trigonometric integrator for the semilinear Klein-Gordon equation with non-smooth solution in the $d$-dimensional space, where $d=1,2,3$. The integrator is constructed based on the full use of Duhamel's formula and the employment of a twisted function tailored for trigonometric integrals. Robust error analysis is conducted, demonstrating that the proposed scheme achieves third-order accuracy in the energy space under a weak regularity requirement in $H^{1+\max(\mu,1)}(\mathbb{T}^d)\times H^{\max(\mu,1)}(\mathbb{T}^d)$ with $\mu> \frac{d}{2}$. A numerical experiment shows that the proposed third-order low-regularity integrator is much more accurate than some well-known exponential integrators of order three for approximating the Klein-Gordon equation with non-smooth solutions.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new numerical integration method to effectively handle the solution of the semilinear Klein - Gordon equation (SKGE) under low - regularity conditions. Specifically, the paper proposes a third - order low - regularity trigonometric integrator, aiming to solve the following problems:
1. **High - precision solution under low - regularity conditions**: Traditional time - discretization methods usually require that the initial data have high regularity (i.e., smoothness). However, in many practical applications, the initial data may be non - smooth. For example, when a high - intensity laser pulse excites a Langmuir wave in a plasma, non - smooth initial conditions will be generated. Therefore, it is necessary to develop numerical methods that can still maintain high precision under low - regularity conditions.
2. **Improve computational efficiency and simplify the algorithm**: Although the existing first - order and second - order low - regularity integrators can handle some non - smooth problems, their computational complexity is relatively high or their algorithms are more complex. The method proposed in this paper not only achieves third - order accuracy under low - regularity conditions but also has a very simple scheme, which is convenient for researchers in different disciplines to use.
3. **Verification by numerical experiments**: Through numerical experiments, it is proved that the proposed third - order low - regularity integrator is more accurate than some known third - order exponential integrators when approximating the Klein - Gordon equation with non - smooth solutions.
### Formula summary
The key formulas in this paper include:
- The form of the Klein - Gordon equation:
\[
\begin{cases}
\partial_{tt}u(t, x)-\Delta u(t, x)+\rho u(t, x)=f(u(t, x)), & 0 < t\leq T, \; x\in T^d\subset\mathbb{R}^d, \\
u(0, x)=u_0(x), \quad \partial_t u(0, x)=v_0(x), & x\in T^d,
\end{cases}
\]
where \(d = 1, 2, 3\), \(\rho\geq0\) is a given parameter, and \(f(u)\) is a given nonlinear function.
- The update formula of the proposed third - order low - regularity integrator:
\[
\begin{aligned}
u_{n + 1}&=\cos(h\sqrt{A})u_n+h\operatorname{sinc}(h\sqrt{A})v_n+h^2\Phi_1(h\sqrt{A})f(u_n) \\
&\quad+h^3\Psi_1(h\sqrt{A})f'(u_n)v_n+h^4\Psi_2(h\sqrt{A})F_1(u_n, v_n), \\
v_{n + 1}&=-hA\operatorname{sinc}(h\sqrt{A})u_n+\cos(h\sqrt{A})v_n+h\Phi_2(h\sqrt{A})f(u_n) \\
&\quad+h^2\Phi_1(h\sqrt{A})f'(u_n)v_n+h^3\Psi_1(h\sqrt{A})(F_1(u_n, v_n)+f'(u_n)f(u_n)),
\end{aligned}
\]
where \(F_1(u_n, v_n)=f''(u_n)(v_n^2 - (\nabla u_n)^2)+\rho f(u_n)-\rho f'(u_n)u_n\), and the coefficient functions are defined as:
\[
\begin{aligned}
\Phi_1(h\sqrt{A})&=