Full discretisation of semi-linear stochastic wave equations driven by multiplicative noise
A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric method for the temporal approximation. This explicit time integrator allows for mean-square error bounds independent of the space discretisation and thus do not suffer from a step size restriction as in the often used Störmer-Verlet-leap-frog scheme. Furthermore, it satisfies an almost trace formula (i. e., a linear drift of the expected value of the energy of the problem). Numerical experiments are presented and confirm the theoretical results.
Key words. Semi-linear stochastic wave equation, Multiplicative noise, Strong convergence, Trace formula, Stochastic trigonometric methods, Geometric numerical integration
AMS subject classifications. 65C20, 60H10, 60H15, 60H35, 65C30
We consider the numerical discretisation of semi-linear stochastic
wave equations of the form
where and , , is a bounded convex domain with polygonal boundary . The “” denotes the time derivative . Assumptions on the smoothness of the nonlinearities and will be given below. The stochastic process is an -valued (possibly cylindrical) -Wiener process with respect to a normal filtration on a filtered probability space . The initial data and are -measurable random variables. We will numerically solve this problem with a linear finite element method in space and a stochastic trigonometric method in time.
We refer to the introductions of [kls10] and [cls13] for the relevant literature on the spatial, respectively temporal, discretisation of stochastic (linear) wave equations. Further, the recent publication [wgt13] presents a full discretisation of the wave equation with additive noise: a spectral Galerkin approximation is used in space and an adapted stochastic trigonometric method, using linear functionals of the noise as in [jk09], is employed in time. Furthermore, the time discretisation of nonlinear stochastic wave equations by stochastic trigonometric methods is analysed in [raey]. Finally, let us mention the recent publication [cqs14] which analyses convergence in of the stochastic trigonometric method applied to the one-dimensional nonlinear stochastic wave equation.
In the present publication, we prove mean-square convergence for the full discretisation to the exact solution to the nonlinear problem (LABEL:swe). Furthermore, using this result, we derive a geometric property of our numerical integrator, namely a trace formula. The trace formula (the linear drift of the expected value of the energy) for the exact solution of (LABEL:swe) as well as for the finite element solution and the completely discrete solution are presented.
Strong approximations of stochastic wave equations are relevant in many real applications. For example, let us consider the motion of a strand of DNA floating in a liquid as presented in [MR2508773] and references therein. The motion of the DNA molecule may be modeled by a wave equation and the impact of the fluid’s molecules may be modeled by a stochastic force acting on the string. When two normally distant parts of the DNA get close enough, biological events, such as release of enzymes, occur. It is thus of interest to consider strong approximation of stochastic wave equations in such a situation.
The paper is organised as follows. We introduce some notations and mention some useful results in the next section. Section LABEL:sect:ms presents a mean-square convergence analysis for our numerical discretisation. A trace formula for the exact and numerical solutions is given in Section LABEL:sect:trace. Finally, numerical experiments illustrating the rates of convergence and the trace formula of the numerical solution are given in the final section.
2 Notations and useful results
Let and be separable Hilbert spaces with norms and respectively. We denote the space of bounded linear operators from to by , and we let be the set of Hilbert-Schmidt operators with norm
where is an arbitrary orthonormal basis of . If , then we write and . Let be a self-adjoint, positive semidefinite operator. We denote the space of Hilbert-Schmidt operators from to by with norm
For the stochastic wave equation (LABEL:swe), we define and denote the -norm by . Further, we set with .
Let be a filtered probability space and the space of -valued square integrable random variables with norm
Next, we define the space , for with norm
where are the eigenpairs of with orthonormal eigenvectors. We also introduce the space
with norm for and . Note that and . In the following we denote the scalar product by and recall the notation for the norm .
Denoting the velocity of the solution to our stochastic partial differential equation by , one can rewrite (LABEL:swe) as
where , , , and . The operator with is the generator of a strongly continuous semigroup of bounded linear operators on , in fact, a unitary group.
Let be a quasi-uniform family of triangulations of the convex polygonal domain with and . Let be the space of piecewise linear continuous functions with respect to which are zero on the boundary of , and let denote the -orthogonal projector and the -orthogonal projector (Ritz projector). Thus,
The discrete Laplace operator is then defined by
We note that . We also define discrete variants of and by
and equipped with the norm . Finally, the finite element approximation of (LABEL:swe) can then be written as
or in the abstract form
where , , and are as before, and with . Note the abuse of notation for the projection and similarly for . This will be used throughout the paper. Again, is the generator of a -semigroup on .
We study the equations (LABEL:swe2) and (LABEL:femswe2) in their mild form
where the semigroups can be expressed as
with , , and .
In order to ensure existence and uniqueness of problem (LABEL:swe) we shall assume that and , with for some regularity parameter , and that the functions and satisfy
for all in the first two inequalities and for all in the last one. Through the text, (or etc.) denotes a generic positive constant that may vary from line to line. We assume that the order of initial regularity so that the discrete initial value is well defined.
Assume that , with and the functions and satisfy (LABEL:assFG) for some . Then there exists a unique solution to the stochastic wave equation (LABEL:swe2) and the finite element equation (LABEL:femswe2) given by the solution of their respective mild equation, i. e., equations (LABEL:exactsol) and (LABEL:exactsolfem).
The proof of this lemma follows from [DaPrato1992, Theorem 7.4], see also the proof of Theorem 2.1 in [raey].
We now collect some results that we will use later on. Sketches of the proofs of these results are collected in the appendix at the end of this paper.
The error estimates for the cosine and sine operators (Corollary 4.2 in [kls10]): Denote and let
Then we have
These will be used to estimate the error contributions from the initial values. In order to deal with the convolution terms in (LABEL:exactsolfem) we single out the following error estimates. Let
Then we have
The temporal Hölder continuity of the sine and cosine operators, see in [cls13]:
together with its continuous version:
The equivalence of and , see the proof of Theorem 4.4 in [kll]: This uses an inverse inequality, hence our assumption about the quasi-uniformity of the mesh family.
The equivalence of the discrete and continuous norm, see in [al13]:
Using the above estimates, one can deduce the following regularity results for the exact solution to our stochastic wave equation (LABEL:swe) and for the exact solution of the finite element approximation (LABEL:femswe1).
Let be the solution to (LABEL:swe), where the initial values satisfy , with , and the functions and satisfy (LABEL:assFG) for some . Then it holds that
and, for ,
The proof of this proposition is very similar to the proof of Proposition LABEL:prop:regFEM given below and is therefore omitted (see also the proofs of Proposition 3.1 and Lemma 3.3 in [raey]).
The next result will be useful in Section LABEL:sect:trace when we will deal with the trace formula of the numerical solution.
Let be the solution to the finite element problem (LABEL:femswe1), where the initial values satisfy , with , and the functions and satisfy (LABEL:assFG) for some . Then it holds that
where we recall that and are the initial position and velocity to the finite element problem.
Proof. Let us start with the first estimate of the norm of and consider the expression
Using the fact that and commute, the boundedness of the cosine operator, together with our assumptions on the initial values for the finite element problem, we get
Similarly, one obtains
To estimate the third term, we use (LABEL:lpl), the assumptions on given in (LABEL:assFG), and the equivalence of the norms stated in (LABEL:equivnorm). First for , we get
because and are bounded. For , we have by (LABEL:lpl)
Finally, Ito’s isometry, equations (LABEL:equivnorm) and (LABEL:lpl), and the assumptions (LABEL:assFG) on give us
All together, for , one thus obtains
and an application of Gronwall’s lemma give the desired bound for .
The proof for the other bound is done in the same way except for a slight difference in the initial values and that in the integrals is replaced by .
We now prove a Hölder regularity property of the finite element solution. We write, for ,
To estimate the first term we use (LABEL:o41cls) to get
for . For we note that and that is bounded in the operator norm. Using a similar argument for the second term, we get the following estimate for the first two terms
for . In order to estimate the third term, we use (LABEL:o41cls), the assumptions on , and the equivalence of the norms given in (LABEL:equivnorm). First for , we obtain
For we have, using (LABEL:o41cls), (LABEL:lpl), (LABEL:equivnorm) and the fact that is bounded in the operator norm
Similarly we get for the fourth term
To estimate terms five and six we use Ito’s isometry, (LABEL:o41cls), (LABEL:lpl), (LABEL:equivnorm) and the assumptions on to get, for ,
For we again use that is bounded in the operator norm.
Collecting the above estimates give us the statement about the regularity of the finite element solution.
3 Mean-square convergence analysis
Recall that the exact solutions to (LABEL:swe2) and (LABEL:femswe2) solve the following equations