Jacques Printems

Arnaud Debussche, Jacques Printems

In this paper we study the approximation of the distribution of $X_t$ Hilbert--valued stochastic process solution of a linear parabolic stochastic partial differential equation written in an abstract form as $$ dX_t+AX_t dt = Q^{1/2} d W_t, \quad X_0=x \in H, \quad t\in[0,T], $$ driven by a Gaussian space time noise whose covariance operator $Q$ is given. We assume that $A^{-α}$ is a finite trace operator for some $α>0$ and that $Q$ is bounded from $H$ into $D(A^β)$ for some $β\geq 0$. It is not required to be nuclear or to commute with $A$. The discretization is achieved thanks to finite element methods in space (parameter $h>0$) and implicit Euler schemes in time (parameter $Δt=T/N$). We define a discrete solution $X^n_h$ and for suitable functions $ϕ$ defined on $H$, we show that $$ |\E ϕ(X^N_h) - \E ϕ(X_T) | = O(h^{2γ} + Δt^γ) $$ \noindent where $γ<1- α+ β$. Let us note that as in the finite dimensional case the rate of convergence is twice the one for pathwise approximations.

Mihály Kovács, Jacques Printems

In this paper we investigate a discrete approximation in time and in space of a Hilbert space valued stochastic process $\{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 $$ \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 $W^{_Q}$ is a $Q$-Wiener process on $H=L^2({\mathcal D})$ and where the main example of $b$ we consider is given by $$ b(t) = t^{β-1}/Γ(β), \quad 0 < β<1. $$ We let $A$ be an unbounded linear self-adjoint positive operator on $H$ and we further assume that there exist $α>0$ such that $A^{-α}$ has finite trace and that $Q$ is bounded from $H$ into $D(A^κ)$ for some real $κ$ with $α-\frac{1}{β+1}<κ\leq α$. The discretization is achieved via an implicit Euler scheme and a Laplace transform convolution quadrature in time (parameter $Δt =T/n$), and a standard continuous finite element method in space (parameter $h$). Let $u_{n,h}$ be the discrete solution at $T=nΔt$. We show that $$ (\E \| u_{n,h} - u(T)\|^2)^{1/2}={\mathcal O}(h^ν + Δt^γ), $$ for any $γ< (1 - (β+1)(α- κ))/2 $ and $ν\leq \frac{1}{β+1}-α+κ$.