English

Strong order of convergence of a fully discrete approximation of a linear stochastic Volterra type evolution equation

Numerical Analysis 2014-11-07 v3 Probability

Abstract

In this paper we investigate a discrete approximation in time and in space of a Hilbert space valued stochastic process {u(t)}t[0,T]\{u(t)\}_{t\in [0,T]} satisfying a stochastic linear evolution equation with a positive-type memory term driven by an additive Gaussian noise. The equation can be written in an abstract form as \ddu+(0tb(ts)Au(s)\dds)\ddt=\ddWQ,t(0,T];u(0)=u0H, \dd u + (\int_0^t b(t-s) Au(s) \, \dd s)\, \dd t = \dd W^{_Q}, t\in (0,T]; \quad u(0)=u_0 \in H, where WQW^{_Q} is a QQ-Wiener process on H=L2(D)H=L^2({\mathcal D}) and where the main example of bb we consider is given by b(t)=tβ1/Γ(β),0<β<1. b(t) = t^{\beta-1}/\Gamma(\beta), \quad 0 < \beta <1. We let AA be an unbounded linear self-adjoint positive operator on HH and we further assume that there exist α>0\alpha >0 such that AαA^{-\alpha} has finite trace and that QQ is bounded from HH into D(Aκ)D(A^\kappa) for some real κ\kappa with α1β+1<κα\alpha-\frac{1}{\beta+1}<\kappa \leq \alpha. The discretization is achieved via an implicit Euler scheme and a Laplace transform convolution quadrature in time (parameter Δt=T/n\Delta t =T/n), and a standard continuous finite element method in space (parameter hh). Let un,hu_{n,h} be the discrete solution at T=nΔtT=n\Delta t. We show that (\Eun,hu(T)2)1/2=O(hν+Δtγ), (\E \| u_{n,h} - u(T)\|^2)^{1/2}={\mathcal O}(h^{\nu} + \Delta t^\gamma), for any γ<(1(β+1)(ακ))/2\gamma< (1 - (\beta+1)(\alpha - \kappa))/2 and ν1β+1α+κ\nu \leq \frac{1}{\beta+1}-\alpha+\kappa.

Keywords

Cite

@article{arxiv.1205.5601,
  title  = {Strong order of convergence of a fully discrete approximation of a linear stochastic Volterra type evolution equation},
  author = {Mihály Kovács and Jacques Printems},
  journal= {arXiv preprint arXiv:1205.5601},
  year   = {2014}
}

Comments

To appear in Math. Comp