Abstract
In this paper we introduce a class of second-order exponential schemes for the time integration of semilinear wave equations. They are constructed such that the established error bounds only depend on quantities obtained from a well-posedness result of a classical solution. To compensate missing regularity of the solution the proofs become considerably more involved compared to a standard error analysis. Key tools are appropriate filter functions as well as the integration-by-parts and summation-by-parts formulas. We include numerical examples to illustrate the advantage of the proposed methods.
Similar content being viewed by others
1 Introduction
In this paper we are interested in solving abstract wave equations of the form
in some Hilbert space \(H\) where L is a positive, self-adjoint operator and G is a sufficiently regular nonlinearity (e.g., Fréchet-differentiable). Such equations arise in many physical models. A prominent example is the cubic wave equation
posed on some interval \(I \subseteq {\mathbb {R}}\) and equipped with appropriate initial and boundary conditions.
Our aim is to construct and investigate time integration schemes for (1.1) under physically realistic assumptions, in particular finite energy conditions. Hence, the solution will in general be of low regularity and we thus restrict ourselves to second-order schemes. Clearly, a standard time integrator (e.g., a Runge–Kutta scheme or an exponential integrator) can only be applied to an abstract evolution equation if it is unconditionally stable due to the unbounded operator L.
In the finite dimensional case (\(\dim H< \infty \)), unconditionally stable integrators (in the sense that L does not cause any restriction on the time step) for this equation were already considered in [8, 12, 17, 24]. Such exponential (or trigonometric) integrators were shown to be second-order convergent while only assuming a finite-energy condition. This was somewhat surprising since usually, second-order (exponential) schemes need two bounded time derivatives of the solution in the error analysis. The key ingredient are certain matrix functions that act as filters. The effect of these filters is that they remove resonances in the local error, which, in contrast to a standard error analysis, enforce cancellation effects in the global error. In fact one can prove that local and global error are of the same order if the filters are chosen appropriately.
Recently, in [2, 3] we presented a completely new technique to prove related results for ordinary differential equations by reformulating a trigonometric integrator as a Strang splitting applied to a modified problem. Using ideas from [20, 21], a specific representation of the local error was derived, which allowed us to separate terms of order three (which can be treated in a standard way) and the leading local error term, which is of order two only. A carefully adapted Lady Windermere’s fan argument is employed to treat these terms in the global error accumulation.
In this paper we prove error bounds for different classes of exponential integrators applied to an evolution equation (1.1) in a unified way. More precisely, we characterize the structure of the defects and the properties of filter functions which allow second-order convergence under a finite-energy condition in different abstract frameworks (i.e., in different function spaces), which define the assumptions on L, G, and the initial data. Within this framework we can handle various boundary conditions.
We point out that our results are not restricted to globally Lipschitz continuous functions G, but also apply to locally Lipschitz ones that satisfy certain growth conditions. Our analysis thus covers the class of nonlinearities for which the existence of a classical solution can be guaranteed. In particular, this includes polynomial nonlinearities up to a certain degree which is determined by the spatial dimension and the corresponding Sobolev embeddings. In the one-dimensional case such equations with periodic boundary conditions and arbitrary high polynomial degree have been studied in [9] for nonlinearities with Lipschitz properties on a whole scale of Sobolev spaces. However, this rich structure is not available in our general framework and, in contrast to our work, well-posedness cannot be guaranteed.
Further work on exponential integration schemes for the time-integration of semilinear wave equations was conducted in [1] where a sine-Gordon equation is studied and the difficulties arise in the proper treatment of a single constant \(c\rightarrow \infty \) which induces high oscillations in time. In the paper [10] the approach from [9] was extended and in the one-dimensional case a quasilinear wave equation with periodic boundary conditions was studied. However, they assume smooth coefficients and high regularity for the analysis. Exponential splitting schemes for linear evolution equations have been analyzed in [15]. The error estimates depend on commutator bounds that are not available in our scenario. Finally, we point out that we do not address long-term behavior as in [4, 6].
The paper is organized as follows. In Sect. 2, we give an informal overview over the methods of interest, the main concepts, and the main results and also present numerical examples illustrating the main results of our work. In particular, the necessity of using averaging techniques in the regime of low-regularity is shown.
The informal overview is made rigorous in Sect. 3, where we introduce the analytic framework and a functional calculus which allows us to define the operator-valued filters and ensures well-posedness of the problem as well as of the schemes.
We further state the assumptions on the operator L, the nonlinearity G, and the initial data that on the one hand will guarantee the well-posedness of (1.1) and on the other hand allow to carry out the error analysis.
In Sect. 4 we characterize filter functions which allow to prove that the exact solution of the original problem and the solution of the averaged problem only differ up to terms of order \(\tau ^2\), where \(\tau >0\) denotes the step size. Section 5 provides a characterization of numerical methods in terms of the structure of their defects, which are necessary to derive error bounds.
Finally, Sects. 6 and 7 contain our main results the error bounds for one-step and for multistep methods, respectively.
2 Informal overview of methods, concepts and results
Before we present the analytical framework necessary to formulate our results rigorously, we first give an informal overview of the methods of interest, the main concepts, and the main results. In the finite dimensional case \(\dim H < \infty \) (which is not the situation of interest in this paper), all the approximations presented are well-defined and the statements valid. However, for evolution equations posed in appropriate function spaces, this is no longer true unless additional assumptions are imposed. Since some of them are rather technical, we postpone them to Sect. 3.
2.1 Problem statement: Second-order differential equation
Let L be a linear, self-adjoint, and positive-definite operator on \(H\) and \(G: [0,t_{\mathrm {end}}] \times H\rightarrow H\). We consider the differential equation
and assume that the solution q satisfies the finite-energy condition
where \(\langle \cdot ,\cdot \rangle _H\) denotes the inner product on H. In first-order formulation, the differential equation can be written as
with
and inner product
Obviously, A is skew-adjoint with respect to \(\langle \cdot ,\cdot \rangle \) and hence has a purely imaginary point spectrum.
2.2 Methods
In the following we shortly present four different types of methods to discretize equation (2.2) in time with a constant stepsize \(\tau >0\).
2.2.1 Strang splitting
The exact flows \(\varphi _\tau ^A\) and \(\varphi _\tau ^{f}\) of the two subproblems
are given explicitly by
We consider the Strang splitting in the variants \(\bigl (A,f,A\bigr )\) and \(\bigl (f,A,f\bigr )\) given by
respectively. Note that the \(\bigl (f,A,f\bigr )\)variant (2.3b) is equivalent to a trigonometric integrator without filter functions, see, e.g., [13, XIII.2.2].
2.2.2 Corrected Lie Splitting
Next we consider the second-order corrected Lie splitting given by
with the correction term
It is inspired by a fourth-order method of this type proposed in [22, 4.9.3 (c)]. Note that in the linear case, where \(f(t,u) = F u \), the correction term reduces to the commutator
Hence, one can consider (2.4) as an approximation to the method
which was considered in [26, (3.37)].
2.2.3 Exponential Runge–Kutta methods
General two-stage exponential Runge–Kutta methods are of the form
where \(c_2 \in (0,1]\) is a given quadrature node. Recall that the \(\varphi \)-functions are defined as
If the coefficient functions \(b_1,b_2\) satisfy
the method is second-order convergent for parabolic problems, see [18, Theorem 4.3.]. Popular choices are \(c_2 = \frac{1}{2} \), \(b_1 = 0\) or \(c_2 = 1 \), \(b_2(z) = \varphi _2(z)\). All our results also apply to the symmetric, but implicit exponential Runge–Kutta scheme from [5, Example 2.1] and to ERKN methods, e.g., those considered in [28]. The necessary modifications are straightforward so that we omit the details.
2.2.4 Exponential multistep methods
The two-step exponential multistep method from [19, (2.7)]
is derived from the variation-of-constants formula for the exact solution of (1.1) by approximating the nonlinearity f in the integral term by an interpolation polynomial using the last two approximations \(u_{n-1}, u_n\).
In a similar manner we consider a method that was used in [7, (B 4)], namely
For \(A = 0 \) it reduces to an explicit Nyström method, cf. method (1.13’) in [14].
2.3 Averaged differential equation
Let \(\chi = \phi ,\psi :i {\mathbb {R}}\rightarrow {\mathbb {R}}\) be even (i.e., \(\chi (-z) = \chi (z)\)) and analytic functions satisfying \(\chi (0) = 1\). Then we define
and an averaged nonlinearity
Using the block diagonal operators
we consider the averaged differential equation
The averaging is done such that the solution \({\widetilde{u}}\) of (2.8) also satisfies a finite-energy condition (2.1) (with a modified constant \({\widetilde{K}}\), which is independent of \(\tau \) and n, cf., Lemma 4.2 below). In Theorem 4.1, we provide sufficient conditions on \(\psi , \phi \) such that
where \(\left\Vert \cdot \right\Vert \) denotes the norm induced by \(\langle \cdot ,\cdot \rangle \).
Remark 2.1
We only consider fixed stepsizes in this paper. Introducing variable stepsizes would require the investigation of a whole family of averaged problems (2.8). In addition, our analysis below would involve a much higher technical effort.
2.4 Averaged methods
The main idea is to apply one of the numerical methods to the averaged equation (2.8) instead of the original one (2.2). Equivalently, one could modify the numerical scheme in an appropriate way using filter functions. This is illustrated in Fig. 1.
Since the solutions of (2.2) and (2.8) only differ by terms of order \(\tau ^2\), one might hope for second-order accuracy if the method is of order two at least. In fact we will later see that in the case of evolution equations, this intuition is not always justified, i.e., order reduction might appear. The main goal in this paper is to characterize the numerical methods, the assumptions on L and G, and the choice of the filter functions which lead to second-order error bounds.
2.5 Main results
Our main results, which are detailed in Theorem 6.2 for exponential one-step methods and in Sect. 7 for exponential multistep methods, are the following error bounds.
-
(a)
The Strang splitting, the exponential Runge–Kutta, and the exponential multistep methods applied to the original equation (2.2) satisfy
$$\begin{aligned} \Vert u(t_n) -u_n \Vert \le C_1 \tau \,. \end{aligned}$$ -
(b)
All methods of Sect. 2.2 applied to the averaged equation (2.8) with appropriate filters \(\phi ,\psi \) satisfy
$$\begin{aligned} \Vert u(t_n) -u_n \Vert \le C_2 \tau ^2 \,. \end{aligned}$$
The constants \(C_1, C_2\) only depend on the initial value \(u_0\), the finite energy K, properties of G, and \(t_{\mathrm {end}}\), but not on n and \(\tau \).
The strategy to prove these bounds is to split the error into two terms, namely
The first term is bounded by Theorem 4.1, the second by Theorem 6.1 or Corollaries 7.1, and 7.2, respectively. A crucial step is to show that the averaged solution inherits the regularity of the original solution, which is done in Lemma 4.2.
2.6 Numerical example
In this section we illustrate the effect of averaging within numerical methods by approximating the solution of a variant of the sine-Gordon equation given on the torus \( {\mathbb {T}}= {\mathbb {R}}/ ( 2 \pi {\mathbb {Z}})\) by
with \(t \in [0,1]\) and \(m_i,m_a \in L^{\infty }({\mathbb {T}})\). Note that for \(q \in L^2({\mathbb {T}})\) and
we have G(q) in \(L^2({\mathbb {T}})\), but even for \(q \in H^1({\mathbb {T}})\) we cannot expect \(G(q) \in H^\epsilon ({\mathbb {T}})\) for any \(\epsilon >0\). Hence, the analysis of [9, 10] does not apply to such non-smooth nonlinearities. For the spatial discretization, we used a Fourier spectral method in order to control the regularity of the solution. The initial values \((q_0,v_0) \in H^1({\mathbb {T}}) \times L^2({\mathbb {T}})\) are constructed such that
for \(\epsilon = 10^{-6}\), see [16] for details.
In Fig. 2 we computed the approximate solution with the Strang splitting variant (2.3a), i.e., \(\bigl (A,{\widetilde{f}},A\bigr )\), with filters (blue, dots)
and without filters, i.e., \(\phi =\psi =1\), (red, crosses) with \(N= 2^j\), \(j=9,10,11\), spatial grid points. The codes are available from the authors on request. We observe order reduction of the non-averaged scheme to order one in the stiff regime, while in the non-stiff regime, the two errors of both schemes are quite close. The non-stiff regime is characterized by time steps \(\tau \) for which \(\varphi _1(\tau A)\) is invertible for all \(\tau < \tau _0\). Since \(\left\Vert A \right\Vert \approx N/2\), this is true for \( \tau _0 \approx 4 \pi / N\). For abstract evolution equations, only the stiff regime is relevant, i.e., the limit \(N\rightarrow \infty \).
3 Analytical framework
We fix some notation for the rest of the paper. For Hilbert spaces X, Y, \(\langle \cdot ,\cdot \rangle \,=\,\langle \cdot ,\cdot \rangle _X \) denotes the scalar product on X and \( {\mathcal {B}}(X,Y)\) the set of all bounded operators \(T:X\rightarrow Y\) equipped with the standard operator norm \(\Vert T \Vert _{Y\leftarrow X}\). Further, \(C^k (X,Y)\) is the space of all k-times Fréchet-differentiable functions from X to Y. We write \(W^{k,p}(\varOmega )\), \(k\in {\mathbb {N}}_0\), \(1\le p \le \infty \), for the Sobolev space of order k with all (weak) derivatives in \(L^p(\varOmega )\) and abbreviate \(H^k(\varOmega ) :=W^{k,2}(\varOmega )\). For multi-indices \(\alpha ,\beta \in {\mathbb {N}}^\ell \) we write \(\alpha \le \beta \) if \(\alpha _i \le \beta _i\) for all \(i=1,\dots ,\ell \).
3.1 Second-order equation
Let H be a real, separable Hilbert space and \(L :{\mathcal {D}}(L) \subseteq H \rightarrow H\) be a positive, self-adjoint operator with compact resolvent. We consider the abstract second-order evolution equation (1.1) in H. To reformulate it as a first-order system we use the intermediate space \(V = {\mathcal {D}}(L^{1/2})\) with
with dense and compact embeddings, in particular, there is a constant \(C_{\rm emb}\) such that
We exemplify the abstract framework considered in the rest of the paper by a class of semilinear wave equations.
Example 3.1
We consider the semilinear evolution equation (1.1) in the following setting:
-
(a)
\(\emptyset \ne \varOmega \subseteq {\mathbb {R}}^d\) is a convex, bounded Lipschitz domain with \(d\in \{1,2,3\}\).
-
(b)
\(L = -\mathop {\mathrm {div}} ( {\mathbf {A}} \nabla )\) with uniformly positive definite \({\mathbf {A}}\in L^{\infty }(\varOmega )^{d\times d}\).
-
(c)
For \(g :[0,t_{\mathrm {end}}] \times \varOmega \times {\mathbb {R}}\rightarrow {\mathbb {R}}\) there is some \(\alpha = ( \alpha _t , \alpha _x, \alpha _y ) \in {\mathbb {N}}^3\) such that all partial derivatives \( \partial ^\beta g \), \(\beta \le \alpha \), exist, are continuous in t and y and bounded in x.
-
(d)
There is \(\gamma >1\) and a constant \(C_\text{g}>0\) such that for all \((t,x,y) \in [0,t_{\mathrm {end}}] \times \varOmega \times {\mathbb {R}}\) we have
$$\begin{aligned} \begin{aligned} | g(t,x,y) | , |\partial _t g(t,x,y) |&\le C_\text{g} \bigl (1 + |y|^\gamma \bigr ) \,, \\ |\partial _y g(t,x,y) |&\le C_\text{g} \bigl (1 + |y|^{\gamma -1} \bigr ) \,. \end{aligned} \end{aligned}$$(3.2a)For the corrected Lie Splitting (2.4) we assume in addition
$$\begin{aligned} |\partial _{yy} g(t,x,y) | & \le C_\text{g} \bigl (1 + |y|^{\gamma -1} \bigr ). \end{aligned} $$(3.2b)For \((t,x) \in [0t_{\mathrm {end}}] \times \varOmega \) and \(q \in V\) we define
$$\begin{aligned} G(t,q) (x) :=g(t,x,q(x)) . \end{aligned}$$
In Table 1 at the end of this section, details for different choices of H are presented for these problems.
In the following we recall sufficient conditions on the nonlinearity G to guarantee well-posedness of the equation and to establish the error analysis presented in Sects. 4, 5, 6, and 7.
Assumption 3.2
(Well-posedness) For G we have \(G \in C^1([0,t_{\mathrm {end}}] \times V, H) \), i.e., G is Fréchet-differentiable with Fréchet-derivative \(J_{G}(t,q) \in {\mathcal {B}}\bigl ([0,t_{\mathrm {end}}] \times V, H\bigr ) \) for all \(q \in V, t\in [0,t_{\mathrm {end}}]\).
The most subtle assumption is given now. It states the necessary regularity for G evaluated at a sufficiently smooth function.
Assumption 3.3
(Regularity of G evaluated at a smooth function) For \(q \in C^1([0,t_{\mathrm {end}}] , V) \cap C([0,t_{\mathrm {end}}] , {\mathcal {D}}(L))\) we have
The next assumption states bounds of G and \(J_{G}\). We point out that the dependency of the constants arising from different radii is crucial for the error analysis.
Assumption 3.4
(Regularity of G) There are constants \(C = C(r)\) such that for given \(r_V,r_L>0\) and q with \(\left\Vert q \right\Vert _{V} \le r_V\), \(\left\Vert q \right\Vert _{{\mathcal {D}}(L)} \le r_L\), \(p \in V\), and \(t\in [0,t_{\mathrm {end}}]\) the following inequalities are satisfied:
For the corrected Lie Splitting (2.4) we assume in addition for \(\left\Vert p_i \right\Vert _{V} \le r_{V}, \; i= 1, \; 2,\)
Remark 3.5
Note that shifting G to \(G + c I\) for some \(c\in {\mathbb {R}}\) does not affect the validity of Assumptions 3.2 to 3.4. Hence, we can also treat positive semidefinite operators L by applying a shift.
In Table 1 the main example 3.1 is specified more precisely. We collected three examples where we stated for a given Hilbert space \(H\) the dimension d of the domain \(\varOmega \) and additional assumptions on the data such that Assumptions 3.2 to 3.4 are satisfied. All examples are posed with homogeneous Dirichlet boundary conditions. By possibly shifting L, we can also treat Neumann, Robin, or periodic boundary conditions.
Higher order Sobolev spaces \(H= H^k(\varOmega )\), \(k\ge 2\), can be handled as well but the spaces and conditions for the operators and parameters become more complicated.
Remark 3.6
Note that from Assumption 3.2 and the chain rule one can only conclude that
instead of (A1).
-
(a)
In Example 3.1 the additional regularity \(q \in C([0,t_{\mathrm {end}}] , {\mathcal {D}}(L))\) is sufficient to verify the Assumption (A1).
-
(b)
For \(G \in C^1([0,t_{\mathrm {end}}] \times V, V)\) the chain rule immediately yields Assumption (A1). However, in Example 3.1 with \(H = H^{-1}(\varOmega )\) and \(V = L^2(\varOmega )\), this would imply that G is already affine-linear, see [11, Sect. 3]. Hence, not even the function \(q \,\mapsto\, \sin (q) \) would be covered by the analysis.
3.2 First-order equation
We consider the first-order formulation (2.2) of equation (1.1) on the separable Hilbert space \(X= V\times H\). The skew-adjoint operator A is given on its domain \({\mathcal {D}}(A) = {\mathcal {D}}(L) \times V\). Hence, A is the generator of a unitary group \(\left( e^{t A} \right) _{t\in {\mathbb {R}}}\). We call u a classical solution of (2.2) on \([0,t^*)\) if u solves (2.2), \(u(0) = u_0\), and
for any \(t_{\mathrm {end}}< t^*\). The Assumptions 3.2 to 3.4 are translated into this setting by means of the following three lemmas. The first one provides a classical solution of (2.2) by standard semigroup theory. All statements in the lemmas directly follow from the special structure of f and the assumptions in Sect. 3.1.
Lemma 3.7
(Well-posedness) Let G satisfy Assumption 3.2. Then \(f:[0,t_{\mathrm {end}}] \times X\rightarrow X\) defined in (2.2) satisfies \(f\in C^1([0,t_{\mathrm {end}}] \times X, X)\) with Fréchet derivative \(J_{f}\bigl ( t, u \bigr ) \in {\mathcal {B}}\left( [0,t_{\mathrm {end}}] \times X, X\right) \) for all \(u \in X\) and \(t\in [0,t_{\mathrm {end}}]\).
The following lemma shows differentiability of f in the stronger \({\mathcal {D}}(A)\) norm.
Lemma 3.8
(Regularity of \(f\) evaluated at a smooth function) Let G satisfy Assumption 3.3and u satisfy (3.3). Then we have
The next Lemma contains two Lipschitz properties of f which easily follow from the corresponding bound on the derivative. They are crucial for the forthcoming error analysis.
Lemma 3.9
(Regularity of \(f\)) Let G satisfy Assumption 3.4. Then there are constants \(C = C(r)\) such that for given \(r_X,r_A>0\) and \(u_i\) with \( \left\Vert u_i \right\Vert \le r_X \), \( \left\Vert u_i \right\Vert _{{\mathcal {D}}(A)} \le r_A \), \(i=1,2\), \(v \in X\), and \(t\in [0,t_{\mathrm {end}}]\) the following inequalities are satisfied:
For the corrected Lie Splitting (2.4) we further have for \( \left\Vert v_i \right\Vert \le r_X \), \(i=1,2,\)
Lemma 3.7 guarantees local well-posedness of (2.2), see [23, Thm. 6.1.5]. Our error analysis only requires assumptions on the data, which then implies the following regularity of the solution.
Proposition 3.10
Let Assumption 3.2be satisfied and take an initial value \(u_0\,\in \,{\mathcal {D}}(A)\). Then there exists a time \(t^* > 0\) and a classical solution of (2.2) on \([0,t^* )\) satisfying (3.3). Hence, for every \(0<t_{\mathrm {end}}<t^*\) there exists a constant \(K>0\) with
In the following we refer to (3.4) as the generalized finite-energy condition.
Remark 3.11
-
(a)
Note that for \(u=\bigl (q, q' \bigr )\) in the situation of Example 3.1 with \(H = H^{-1}(\varOmega )\), the generalized finite energy condition implies
$$\begin{aligned} \left\Vert Au (t) \right\Vert ^2 = \left\Vert q \right\Vert _{{\mathcal {D}}(L)} ^2 + \left\Vert q' \right\Vert _{V} ^2 = \big \Vert {\mathbf {A}}^{1/2} \, \nabla q \big \Vert _{L^2} ^2 + \left\Vert q' \right\Vert _{L^2} ^2 \le K^2 \,. \end{aligned}$$This corresponds to the finite energy condition used in [8, 12, 17, 24].
-
(b)
The bound (3.4) also implies
$$\begin{aligned} \left\Vert u' (t) \right\Vert ^2 = \left\Vert q' \right\Vert _{V} ^2 + \left\Vert q'' \right\Vert _{H} ^2 \le K^2 \,. \end{aligned}$$
3.3 Filter
From the compact resolvent property of L and the compact embeddings we can infer that also A has a compact resolvent. Hence, A admits an orthonormal basis of eigenvectors
where \( M\subseteq {\mathbb {N}}\) and \(\lambda _k \in {\mathbb {R}}\). Any \(x\in X\) can thus be represented as
with the equivalence
This enables us to define the following functional calculus on the set
see [25, Theorem 5.9]. It leads to the following properties of operator functions.
Theorem 3.12
Let \(A:{\mathcal {D}}(A) \rightarrow H\) be a skew-adjoint operator on a separable Hilbert space \(X\) with compact resolvent. Then the map \(\varPsi _A :{\mathcal {C}}_b \left( i {\mathbb {R}}\right) \rightarrow {\mathcal {L}}(X)\),
satisfies the following properties:
-
(a)
\(\varPsi _A\) is linear
-
(b)
\( \left\Vert h(A) \right\Vert _{X\leftarrow X} \le \Vert h\Vert _\infty \)
-
(c)
\((g h) (A) = g(A) h(A)\)
-
(d)
For \(x\in {\mathcal {D}}(A)\) it holds \(h(A)x \in {\mathcal {D}}(A)\) and \(A h(A)x = h(A) Ax\)
For the construction of the integrators we make use of filter functions.
Definition 3.13
Let \(\chi \in {\mathcal {C}}_b \left( i {\mathbb {R}}\right) \). We call \(\chi \) a filter of order m, \(m=1,2\), if the following properties are satisfied. There exist \(\vartheta , \varTheta \in {\mathcal {C}}_b \left( i {\mathbb {R}}\right) \) such that for all \(z\in i{\mathbb {R}}\)
In addition, for \(m=2\), \(\chi \) is symmetric, i.e.
Note that (F3) is equivalent to \(\chi (z) = \varphi _1(z) \varTheta (z) \).
By Theorem 3.12 we can define a corresponding class of filter operators that we later use in the averaged schemes.
Theorem 3.14
Let \(\tau >0\) and \(\chi \in {\mathcal {C}}_b \left( i {\mathbb {R}}\right) \) be a filter of order m with \(\vartheta ,\varTheta \) from Definition 3.13. Then we have
Here, \(\pi _i:X\rightarrow X\) denotes the projection onto the i-th component.
Proof
The properties (OF1), (OF3), (OF4) directly follow from the functional calculus and (OF2) is a direct consequence of (OF4). To prove (OF5), we use the fact that we can approximate \(\chi \) uniformly on \(i{\mathbb {R}}\) by even rational functions as \(\lim _{x\rightarrow \pm \infty } \chi (i x)=\,0\), see [27, Sect. 1.6]. Hence, the assertion is true since it is easily verified for functions of the type
with some \(\delta > 0\). \(\square \)
Remark 3.15
-
(a)
An example for \(m=2\) is the short average filter proposed in [8] that we used in (2.11). We note that in this example \(\chi (i x) = \mathop {\mathrm {sinc}}(\tfrac{x}{2})\) holds for all \(x\in {\mathbb {R}}\), which relates our filters to the ones considered in [13, Chapter XIII.].
-
(b)
We further obtain \( \left\Vert \tau A \vartheta (\tau A) \right\Vert _{X\leftarrow X}^2 \le 2 \Vert \vartheta \Vert _\infty \) for \(m=2\) as
$$\begin{aligned} |z \vartheta (z)|^2 = |z^2 \vartheta (z)| \, |\vartheta (z)| \le 2 \Vert \vartheta \Vert _\infty \quad \text{for all} \quad z\in i{\mathbb {R}}\,. \end{aligned}$$
4 Averaged problem
In this section we bound the difference between the solution \({\widetilde{u}}\) of the averaged equation (2.8) and the solution u of (2.2). Note that by Proposition 3.10, a unique classical solution \({\widetilde{u}}\) of (2.8) exists since the assumptions on \(f\) also hold for \({\widetilde{f}}\).
In order to apply (A5a′) we define \(r_X\) via
with \(C_{\rm emb}\) defined in (3.1) and K in (2.1).
Theorem 4.1
Let Assumptions 3.2, to 3.4be valid and consider the averaged nonlinearity \({\widetilde{f}}\) defined in (2.8) with second-order filters. Then there is a \(\tau _0>0\) and a constant \(C_{{\rm av}}>0\) such that for all \(\tau \le \tau _0\)
The constant \(C_\text{av}\) and \(\tau _0\) depend on \(r_X\), \(u_0\), \(t_{\mathrm {end}}\), the finite energy K defined in (2.1), the filter functions, and the embedding constant \(C_{\rm emb}\), but not on \(\tau \).
In particular, \({\widetilde{u}}\) exists on \([0,t_{\mathrm {end}}]\) and is bounded by
Proof
Let \(\widetilde{t^*}>0\) be the maximal existence time of \({\widetilde{u}}\) and define
We first observe that for \(t\le \min \{ t_0, t_{\mathrm {end}}\}\) the variation-of-constants formula yields
with
By Assumption (A5a′) and since \(t\le t_0\), the third term in (4.2) is bounded by
It remains to prove
since these bounds are sufficient to apply a Gronwall lemma which shows the assertion for all \(t\le \min \{ t_0, t_{\mathrm {end}}\}\).
To bound \(I_1\) we use integration-by-parts and (OF3) to obtain
where we used that \( f\big (s,u(s)\big )\) is differentiable in \(X\). By Assumptions (A3′), (A4b′), and the bound (3.4) on \(u'\) we have
This proves (4.3) for \(j=1\).
Using the notation \({\mathbf {u}}(s,\sigma ) = \sigma u(s) + (1-\sigma ) \varPhi u(s) \) and the differentiability (A1′) of \(f\) we get
By (OF3) and integration-by-parts, the first term can be rewritten as
Hence, we have \(\left\Vert I_{2,1}(t) \right\Vert \le C \tau ^2\) by (A2′), (A4a′), and (A4b′).
By assumption (A1′) we also have
Again integration-by-parts and Assumptions (A1′) and (A4b′) yield the desired bound (4.3). Using (4.1) for \(t\le \min \{ t_0, t_{\mathrm {end}}\}\) we obtain for \(\tau \le \tau _0\)
This proves \(t_0 \ge t_{\mathrm {end}}\) and hence (4.1) holds on \([0,t_{\mathrm {end}}]\) for all \(\tau \le \tau _0\). \(\square \)
In the next lemma we show that \({\widetilde{u}}\) inherits the regularity of u uniformly in \(\tau \).
Lemma 4.2
Let Assumptions 3.2to 3.4be valid. Then there is a \(\tau _0>0\) and a constant \({\widehat{C_\text{av}}}>0\) such that for all \(\tau \le \tau _0\)
In particular, \({\widetilde{u}}\) satisfies the generalized finite-energy condition uniformly in \(\tau \le \tau _0\), i.e.,
where \(\tau _0\) and the constants \({\widehat{C_\text{av}}}\) and \({\widetilde{K}}\) depend on \(r_X\), \(u_0\), \(t_{\mathrm {end}}\), the finite energy K defined in (2.1), the filter functions, and the embedding constant \(C_{\rm emb}\), but not on \(\tau \).
Proof
We proceed as in the proof of Theorem 4.1 and define \(t_0\) by
For \(0 \le t \le t_0\), (4.1), (4.2), and (A5b′) imply
With Remark 3.15, similar arguments as before yield \({\mathcal {O}}(\tau )\) bounds for \(\left\Vert A I_1(t) \right\Vert \) and \(\left\Vert A I_2(t) \right\Vert \). By possibly reducing \(\tau _0\) we obtain the result for \(0 \le t \le t_{\mathrm {end}}\). This immediately implies the first bound in (4.5) and the second bound is then obtained from (2.8). \(\square \)
Remark 4.3
Note that Theorem 4.1 and Lemma 4.2 remain true for \(\varPsi = I \) as for this choice \(I_1(t) = 0\) and the proof does not require (F3). This case is of interest for methods (2.5) and (2.6).
5 Abstract assumptions on the methods
In this section we characterize the classes of methods which are covered by our error analysis.
We recall that u denotes the solution of the original problem (2.2) and \({\widetilde{u}}\) the solution of the averaged problem (2.8). Further, we denote the numerical flow by \(S_\tau \) and the defect by \(\delta _n\), i.e., a one-step method is given by
We start with an assumption on the stability of the method.
Assumption 5.1
(Stability) The method applied to (2.8) is stable in the sense that for all \(v \in \mathcal{D}(A),\; w \in X\), \(t \ge 0\),
where \({\mathcal {J}}: {\mathbb {R}}\times \mathcal{D}(A) \times X\rightarrow X\) is bounded by
Next, we consider the consistency.
Assumption 5.2
(Consistency for order one) The method applied to the original equation (2.2) satisfies Assumption 5.1 (with \(\phi =\psi =1\)) and its defect (5.1) satisfies
where \(C>0\) is independent of \(\tau \) and n.
For second-order methods, our analysis requires a particular structure of the defect. Before we state this in an abstract way, we briefly motivate it. Most of the methods we consider are constructed from the variation-of-constants formula
where only the integral term is approximated. Hence, this defect can be expressed as some quadrature error that contains the second derivative in s of
depending on the precise method. Terms of order \(\tau ^3\) do not cause any difficulties. However, terms of lower order exist which, in general, require more careful treatment. From \(f_1\) we obtain the second-order term
and \(f_2\) gives in addition the term
For these terms property (OF4) needs to be used in order to carry over the local convergence order to the global error. Similar terms are obtained for the defect of the splitting scheme (2.4). Together with the integral in (5.4), equations (5.7) and (5.6) give rise to the following general structure of \(\delta _n\).
Assumption 5.3
(Structure of defects for order two) The defect \(\delta _n\) defined in (5.1) of a numerical method applied to the averaged equation (2.8) is of the form
with \(\left\Vert D_n \right\Vert \le C \tau ^3\), where the constant \(C>0\) is independent of \(\tau \) and n. In addition, one of the following sets of conditions is satisfied:
-
(a)
If \(\phi ,\psi \) are filters of order 2, then there exist \(w_n \in X\) and a linear map \(W_n:X\rightarrow {\mathcal {D}}(A)\) which satisfy
$$\begin{aligned}& \left\Vert w_n \right\Vert\le C , &&\Big \Vert \frac{1}{\tau } \big ( w_{n+1} - w_{n} \bigr ) \Big \Vert\le C , \end{aligned}$$(5.9a)$$\begin{aligned} &\left\Vert W_n \right\Vert _{X\leftarrow X} \le C ,&\Big \Vert \frac{1}{\tau } \bigl (W_{n+1} - W_n\bigr ) \Big \Vert\le C , \end{aligned}$$(5.9b)$$\begin{aligned} &\left\Vert A W_n \right\Vert _{X\leftarrow X}\le C , &&&&&&&&&&&&&\end{aligned}$$(5.9c)with a constant C which is independent of \(\tau \) and n such that \(\delta _n^{(i)}\) can be written as
$$\begin{aligned} \delta _n^{(1)} = \tau ^2 \bigl (\tau A \varPsi \bigr ) w_n\,, \qquad \delta _n^{(2)} = \tau ^2 W_n \bigl (\tau A \varPhi \bigr ) A {\widetilde{u}}(t_n) \,, \end{aligned}$$(5.10) -
(b)
If \(\psi =1\) and \(\phi \) is a filter of order 2, then (5.9) and (5.10) hold with \(w_n = 0\) for all n.
Remark 5.4
From (5.9) and (OF2) one can directly derive \(\left\Vert \delta _n \right\Vert \le C \tau ^2\). However, this would only yield a suboptimal first-order bound in the global error.
The following proposition embeds the methods presented in Sect. 2.2 in the abstract framework.
Proposition 5.5
Let Assumptions 3.2 to 3.4 be satisfied.
-
(a)
The Strang splitting methods (2.3a) and (2.3b) applied to the averaged equation (2.8) satisfy Assumptions 5.1, 5.2, and 5.3 (a).
-
(b)
The second-order variant of the Lie splitting (2.4) applied to the averaged equation (2.8) satisfies Assumptions 5.1and 5.3 (a).
-
(c)
The exponential Runge–Kutta method (2.5) applied to the averaged equation (2.8) satisfies Assumptions 5.1, 5.2, and 5.3 (b).
Proof
Assumption 5.1 is easily verified for all schemes. We only prove part (a) for the Strang splitting \(\bigl (A,{\widetilde{f}},A\bigr )\) as the statement for the \(\bigl ({\widetilde{f}},A,{\widetilde{f}}\bigr )\) variant and part (c) can be adapted from this proof. We comment on part (b) below.
Let \(t_{n+\xi } :=t_n + \tau \xi \), \({\widetilde{u}}_{n+\xi } :={\widetilde{u}}(t_{n+\xi } ) \), and \({\widetilde{f}}_{n+\xi } :={\widetilde{f}}(t_{n+\xi } ,{\widetilde{u}}_{n+\xi } ) \). Since we can write the scheme as
the defect is given by
where
\({\widehat{I}}_1\) is the quadrature error of the midpoint rule. It can be written in terms of the Peano kernel \(\kappa _2\) as
with \(\left\Vert {\widehat{D}}^{(1)}_n \right\Vert \le C \tau ^3\). Again using the variation-of-constants formula, (OF2), (A5a′), and (A5b′) we obtain
with \(\big \Vert {\widehat{D}}_n \big \Vert \le C \tau ^3\).
To bound \({\widehat{I}}_2\) recall that \({\widetilde{f}}\) only depends on the first component of \({\widetilde{u}}\). Using (A5a′), the variation-of-constants formula, and \(\pi _1 {\widetilde{f}}_{n+1/2} = 0\), we have
since this is just a quadrature error of the (right) rectangular rule.
The properties (5.9a) to (5.9c) follow directly from Lemma 3.9. Using the first order Peano kernel \(\kappa _1\), Assumption 5.2 is verified by writing
as this yields \(\big \Vert {\widehat{I}}_1 \big \Vert \le C \tau ^2\) for \(\psi = \phi =1 \) by (A1′) and (A3′).
We briefly comment on the scheme (2.4). The defect can be written as
In the first term \({\widehat{I}}_3\) we add and subtract \(\tau e^{\tau /2 A} {\widetilde{f}}_{n+1/2} \) and get the quadrature error of the midpoint rule. The term \({\widehat{I}}_4\) admits a similar structure as in \({\widehat{I}}_1\) and hence Assumption 5.3 can be verified as before. \(\square \)
Remark 5.6
We note that method (2.4) applied to the original equation (2.2) does not satisfy Assumption 5.2.
6 Main result for exponential one-step methods
The following result is the last step towards our main Theorem 6.2. It states the global error of a numerical integrator applied to the averaged equation (2.8) with suitable filters satisfying our assumptions (e.g., all the methods of Sect. 2.2) is second order accurate. As before, u denotes the solution of the original problem (2.2) and \({\widetilde{u}}\) the solution of the averaged problem (2.8).
Theorem 6.1
(Global error of the averaged problem) Let Assumptions 3.2to 3.4be fulfilled. Moreover, let \((u_n)_n\) be the numerical approximations of a scheme applied to the averaged equation (2.8) that satisfies Assumptions 5.1and 5.3. Then there is a \(\tau _0>0\) and a constant \(C_\text{e}>0\) such that for all \(\tau \le \tau _0\)
The constant \(C_\text{e}\) and \(\tau _0\) depend on \(u_0\), \(t_{\mathrm {end}}\), the finite energy K defined in (2.1), the filter functions, and the embedding constant \(C_{\rm emb}\), but are independent of \(\tau \) and n.
Proof
The proof makes use of the error recursion from [12] and adapts techniques from Theorem 5.3 in [2].
Due to definition (5.1) of the defect \(\delta _n\), the global error \({\widetilde{e}}_n = {\widetilde{u}}(t_n) - u_n\) can be written as
By Assumption 5.1, the global error satisfies
The error bound follows from a discrete Gronwall lemma, once we established the bound
with a constant \(C_\delta \) being independent of \(\tau \) and n.
The proof is done by induction on n. For \(n=0\), the statement is obviously true. Hence we assume that for all \(0\le k \le n\) it holds
By Assumption 5.3, the defect is split into three parts, which motivates to write
where
Since \(\left\Vert D_j \right\Vert \le C\tau ^3\) and \(n\tau \le t_{\mathrm {end}}\) we easily see
To bound \({\widetilde{e}}_{n+1}^{(\ell )}\), \(\ell =1,2\), we define
Summation-by-parts, Assumption 5.3, and (OF4) yield
To bound \(E_j(e^{\tau A} - I)\) we exploit a telescopic sum to get
Together with (5.9a) and (OF4) this yields (6.2) for \(\delta _j^{(1)}\) instead of \(\delta _j\).
Next we consider \({\widetilde{e}}_{n+1}^{(2)}\). Again, Assumption 5.3, summation-by-parts, and (OF4) with \(\chi = \varPhi \) yield
Here, we have
The terms can be estimated by (5.9b) and (5.9c)
since \(\left| \varphi _1(z) \right| \le 1\) for \(z \in i{\mathbb {R}}\).
Next we consider \((e^{\tau A} - I) A F_j\) for \(j\le n\). After adding the exact solution we apply the variation-of-constants formula, (A3′), and (4.5), which gives
This yields (6.2) for \(\delta _j^{(2)}\) instead of \(\delta _j\) and together with the results above proves (6.2).
Finally, (5.3), (6.1), (6.2), and \({\widetilde{e}}_0 = 0\) give
A discrete Gronwall Lemma thus yields
for \(\tau \le \tau _0 \le \tfrac{1}{2} \bigl (\tfrac{r_X}{C_\text{e}}\bigr )^{1/2}\) and the induction is closed.\(\square \)
Our main result is the following theorem.
Theorem 6.2
Let Assumptions 3.2to 3.4be fulfilled. Further let \((u_n)_n\) be the numerical approximations of a scheme that satisfies Assumptions 5.1and 5.3.
-
(a)
If the method also satisfies Assumptions 5.2and is applied to the original equation (2.2), then there is a \(\tau _0>0\) and a constant \(C_1>0\) such that for all \(\tau \le \tau _0\)
$$\begin{aligned} \left\Vert u_n - u(t_n) \right\Vert \le C_1 \tau , \qquad 0 \le t_n = n\tau \le t_{\mathrm {end}}\,. \end{aligned}$$ -
(b)
Let \(\phi ,\psi \) such that Assumption 5.3is satisfied. Then there is a \(\tau _0>0\) and a constant \(C_2>0\) such that for all \(\tau \le \tau _0\)
$$\begin{aligned} \left\Vert u_n - u(t_n) \right\Vert \le C_2 \tau ^2, \qquad 0 \le t_n = n\tau \le t_{\mathrm {end}}, \end{aligned}$$if the method is applied to the averaged equation (2.8).
The constants \(C_1,C_2\) and \(\tau _0\) depend on \(u_0\), \(t_{\mathrm {end}}\), the finite energy K defined in (2.1), the filter functions, and the embedding constant \(C_{\rm emb}\), but are independent of \(\tau \) and n.
Proof
Part (a) follows directly from Assumption 5.2 and equation (6.1). For part (b), we simply combine Theorem 4.1 and Theorem 6.1 by the triangle inequality (2.9). \(\square \)
7 Main result for exponential multistep methods
We briefly indicate how to extend the developed theory to the exponential multistep methods of Sect. 2.2.4. The first-order convergence as in part (a) of Theorem 6.2 is easily shown. To get second order, Assumption 5.1 needs to be modified.
For method (2.6), we denote the numerical flow by \(S_\tau (t,v_n,v_{n-1})\) and obtain
where \({\mathcal {J}}_n = {\mathcal {J}}\bigl (t,v_n,v_{n-1},w_n,w_{n-1}\bigr )\) is bounded by
This yields the following convergence result.
Corollary 7.1
Let Assumptions 3.2to 3.4be valid. Consider the numerical approximations \((u_n)_n\) from (2.6) applied to the averaged equation (2.8) with \(\psi = 1\) and a filter \(\phi \) of order 2. Then there is a \(\tau _0>0\) and a constant \(C>0\) such that for all \(\tau \le \tau _0\)
where C and \(\tau _0\) depend on \(u_0\), \(t_{\mathrm {end}}\), the finite energy K defined in (2.1), the filter functions, and the embedding constant \(C_{\rm emb}\), but are independent of \(\tau \) and n.
Proof
We first employ Theorem 4.1 and Lemma 4.2, so again it remains to prove the error in approximating the filtered solution. As in the proof of [19, Thm. 4.3] the defect stems from a quadrature error that yields the dominant terms as in (5.6). Considering the defect
Assumption 5.3 (b) is satisfied and a slight modification of the proof of Theorem 6.1 yields the assertion. \(\square \)
For method (2.7) we have
where \({\mathcal {J}}_n={\mathcal {J}}\left( t,v_n,w_n\right) \) is bounded by
In order to apply the techniques from above we define the modification
and can state the following result.
Corollary 7.2
Let Assumptions 3.2to 3.4be valid and u be the classical solution of (2.2).
Consider the numerical approximations \((u_n)_n\) from (2.7) applied to the averaged equation (2.8) with filters \(\chi _2\psi ,\chi _2 \phi \) where \(\psi ,\phi \) are filters of order 2. Then there is a \(\tau _0>0\) and a constant \(C>0\) such that for all \(\tau \le \tau _0\)
where C and \(\tau _0\) depend on \(u_0\), \(t_{\mathrm {end}}\), the finite energy K defined in (2.1), the filter functions, and the embedding constant \(C_{\rm emb}\), but are independent of \(\tau \) and n.
Proof
Since the method stems from a midpoint rule applied to the variation-of-constants formula the defect is again given with dominant terms similar to (5.6) and (5.7). If we resolve the error recursion, we only obtain every second defect and the propagation is driven by \(e^{2 \tau A}\). As \(e^{z}\) in (F3) is replaced by \(e^{2 z}\), this can be combined to conclude the assertion similar to the proof of Theorem 6.1. \(\square \)
Change history
29 June 2021
The funding information section was missing from this article and should have read Open Access funding enabled and organized by Projekt DEAL
References
Baumstark, S., Faou, E., Schratz, K.: Uniformly accurate exponential-type integrators for Klein-Gordon equations with asymptotic convergence to the classical NLS splitting. Math. Comp. 87(311), 1227–1254 (2018). https://doi.org/10.1090/mcom/3263
Buchholz, S., Gauckler, L., Grimm, V., Hochbruck, M., Jahnke, T.: Closing the gap between trigonometric integrators and splitting methods for highly oscillatory differential equations. IMA J. Numer. Anal. 38(1), 57–74 (2018). https://doi.org/10.1093/imanum/drx007
Buchholz, S.F.: Fehleranalyse von auf trigonometrischen Integratoren basierenden Splittingverfahren für hochoszillatorische, semilineare Probleme. Ph.D. thesis, Karlsruher Institut für Technologie (KIT) (2019), https://doi.org/10.5445/IR/1000088935
Cano, B.: Conservation of invariants by symmetric multistep cosine methods for second-order partial differential equations. BIT 53(1), 29–56 (2013). https://doi.org/10.1007/s10543-012-0393-1
Celledoni, E., Cohen, D., Owren, B.: Symmetric exponential integrators with an application to the cubic Schrödinger equation. Found. Comput. Math. 8(3), 303–317 (2008). https://doi.org/10.1007/s10208-007-9016-7
Cohen, D., Hairer, E., Lubich, C.: Conservation of energy, momentum and actions in numerical discretizations of non-linear wave equations. Numer. Math. 110(2), 113–143 (2008). https://doi.org/10.1007/s00211-008-0163-9
Frisch, U., She, Z.S., Thual, O.: Viscoelastic behaviour of cellular solutions to the Kuramoto–Sivashinsky model. J. Fluid Mech. 168, 221–240 (1986). https://doi.org/10.1017/S0022112086000356
García-Archilla, B., Sanz-Serna, J.M., Skeel, R.D.: Long-time-step methods for oscillatory differential equations. SIAM J. Sci. Comput. 20(3), 930–963 (1999). https://doi.org/10.1137/S1064827596313851
Gauckler, L.: Error analysis of trigonometric integrators for semilinear wave equations. SIAM J. Numer. Anal. 53(2), 1082–1106 (2015). https://doi.org/10.1137/140977217
Gauckler, L., Lu, J., Marzuola, J.L., Rousset, F., Schratz, K.: Trigonometric integrators for quasilinear wave equations. Math. Comp. 88(316), 717–749 (2019). https://doi.org/10.1090/mcom/3339
Goldberg, H., Kampowsky, W., Tröltzsch, F.: On Nemytskij operators in \(L_p\)-spaces of abstract functions. Math. Nachr. 155, 127–140 (1992). https://doi.org/10.1002/mana.19921550110
Grimm, V., Hochbruck, M.: Error analysis of exponential integrators for oscillatory second-order differential equations. J. Phys. A 39(19), 5495–5507 (2006). https://doi.org/10.1088/0305-4470/39/19/S10
Hairer, E., Lubich, C., Wanner, G.: Geometric numerical integration: Structure-preserving algorithms for ordinary differential equations. Springer Series in Computational Mathematics, vol. 31, 2nd edn. Springer-Verlag, Berlin (2006)
Hairer, E., Nørsett, S.P., Wanner, G.: Solving ordinary differential equations I: Nonstiff problems. Springer Series in Computational Mathematics, vol. 8, 2nd edn. Springer-Verlag, Berlin (1993)
Hansen, E., Ostermann, A.: Exponential splitting for unbounded operators. Math. Comp. 78(267), 1485–1496 (2009). https://doi.org/10.1090/S0025-5718-09-02213-3
Hochbruck, M., Leibold, J., Ostermann, A.: On the convergence of Lawson methods for semilinear stiff problems. Numer. Math. 145(3), 553–580 (2020). https://doi.org/10.1007/s00211-020-01120-4
Hochbruck, M., Lubich, C.: A Gautschi-type method for oscillatory second-order differential equations. Numer. Math. 83(3), 403–426 (1999). https://doi.org/10.1007/s002110050456
Hochbruck, M., Ostermann, A.: Explicit exponential Runge–Kutta methods for semilinear parabolic problems. SIAM J. Numer. Anal. 43(3), 1069–1090 (2005). https://doi.org/10.1137/040611434
Hochbruck, M., Ostermann, A.: Exponential multistep methods of Adams-type. BIT 51(4), 889–908 (2011). https://doi.org/10.1007/s10543-011-0332-6
Jahnke, T., Lubich, C.: Error bounds for exponential operator splittings. BIT 40(4), 735–744 (2000). https://doi.org/10.1023/A:10223965196567
Lubich, C.: On splitting methods for Schrödinger–Poisson and cubic nonlinear Schrödinger equations. Math. Comp. 77(264), 2141–2153 (2008). https://doi.org/10.1090/S0025-5718-08-02101-7
McLachlan, R.I., Quispel, G.R.W.: Splitting methods. Acta Numer. 11, 341–434 (2002). https://doi.org/10.1017/S0962492902000053
Pazy, A.: Semigroups of Linear Operators and Applications to Partial Differential Equations. Applied Mathematical Sciences, vol. 44. Springer, New York (1983), https://doi.org/10.1007/978-1-4612-5561-1
Sanz-Serna, J.M.: Mollified impulse methods for highly oscillatory differential equations. SIAM J. Numer. Anal. 46(2), 1040–1059 (2008). https://doi.org/10.1137/070681636
Schmüdgen, K.: Unbounded self-adjoint operators on Hilbert space. Graduate Texts in Mathematics, vol. 265. Springer, Dordrecht (2012), https://doi.org/10.1007/978-94-007-4753-1
Suzuki, M.: Generalized Trotter’s formula and systematic approximants of exponential operators and inner derivations with applications to many-body problems. Comm. Math. Phys. 51(2), 183–190 (1976), https://doi.org/10.1007/BF01609348
Timan, A.F.: Theory of approximation of functions of a real variable. Translated from the Russian by J. Berry. English translation edited and editorial preface by J. Cossar. International Series of Monographs in Pure and Applied Mathematics, Vol. 34. A Pergamon Press Book. The Macmillan Co., New York (1963)
Wang, B., Wu, X.: Global error bounds of one-stage extended RKN integrators for semilinear wave equations. Numer. Algorithms 81(4), 1203–1218 (2019). https://doi.org/10.1007/s11075-018-0585-0
Acknowledgements
The authors thank Ludwig Gauckler for helpful discussions on handling the global error in Theorem 6.1 and Jan Leibold for his careful reading of this manuscript.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Corresponding author
Additional information
This article is part of the topical collection “Waves 2019 invited papers” edited by Manfred Kaltenbacher and Markus Melenk.
Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 258734477 – SFB 1173.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Buchholz, S., Dörich, B. & Hochbruck, M. On averaged exponential integrators for semilinear wave equations with solutions of low-regularity. Partial Differ. Equ. Appl. 2, 23 (2021). https://doi.org/10.1007/s42985-020-00045-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s42985-020-00045-9
Keywords
- Highly oscillatory problems
- Error bounds
- Order-reduction
- Time-integration
- Second-order evolution equations
- Filter functions
- Summation-by-parts formula