Mihály Kovács

Mihály Kovács, Stig Larsson, Fredrik Lindgren

We consider the stochastic Allen--Cahn equation perturbed by smooth additive Gaussian noise in a spatial domain with smooth boundary in dimension $d\le 3$, and study the semidiscretisation in time of the equation by an Euler type split-step method. We show that the method converges strongly with a rate $O(Δt^{\frac12}) $. By means of a perturbation argument, we also establish the strong convergence of the standard backward Euler scheme with the same rate.

Daisuke Furihata, Mihály Kovács, Stig Larsson et al.

We consider the stochastic Cahn-Hilliard equation driven by additive Gaussian noise in a convex domain with polygonal boundary in dimension $d\le 3$. We discretize the equation using a standard finite element method in space and a fully implicit backward Euler method in time. By proving optimal error estimates on subsets of the probability space with arbitrarily large probability and uniform-in-time moment bounds we show that the numerical solution converges strongly to the solution as the discretization parameters tend to zero.

David Bolin, Kristin Kirchner, Mihály Kovács

The numerical approximation of solutions to stochastic partial differential equations with additive spatial white noise on bounded domains in $\mathbb{R}^d$ is considered. The differential operator is given by the fractional power $L^β$, $β\in(0,1)$, of an integer order elliptic differential operator $L$ and is therefore non-local. Its inverse $L^{-β}$ is represented by a Bochner integral from the Dunford-Taylor functional calculus. By applying a quadrature formula to this integral representation, the inverse fractional power operator $L^{-β}$ is approximated by a weighted sum of non-fractional resolvents $(I + t_j^2 L)^{-1}$ at certain quadrature nodes $t_j>0$. The resolvents are then discretized in space by a standard finite element method. This approach is combined with an approximation of the white noise, which is based only on the mass matrix of the finite element discretization. In this way, an efficient numerical algorithm for computing samples of the approximate solution is obtained. For the resulting approximation, the strong mean-square error is analyzed and an explicit rate of convergence is derived. Numerical experiments for $L=κ^2-Δ$, $κ> 0$, with homogeneous Dirichlet boundary conditions on the unit cube $(0,1)^d$ in $d=1,2,3$ spatial dimensions for varying $β\in(0,1)$ attest the theoretical results.

Boris Baeumer, Mihály Kovács, Harish Sankaranarayanan

We give stability and consistency results for higher order Grünwald-type formulae used in the approximation of solutions to fractional-in-space partial differential equations. We use a new Carlson-type inequality for periodic Fourier multipliers to gain regularity and stability results. We then generalise the theory to the case where the first derivative operator is replaced by the generator of a bounded group on an arbitrary Banach space.

Adam Andersson, Mihály Kovács, Stig Larsson

We prove a weak error estimate for the approximation in space and time of a semilinear stochastic Volterra integro-differential equation driven by additive space-time Gaussian noise. We treat this equation in an abstract framework, in which parabolic stochastic partial differential equations are also included as a special case. The approximation in space is performed by a standard finite element method and in time by an implicit Euler method combined with a convolution quadrature. The weak rate of convergence is proved to be twice the strong rate, as expected. Our convergence result concerns not only functionals of the solution at a fixed time but also more complicated functionals of the entire path and includes convergence of covariances and higher order statistics. The proof does not rely on a Kolmogorov equation. Instead it is based on a duality argument from Malliavin calculus.

Mihály Kovács, Felix Lindner, René L. Schilling

We present an abstract framework to study weak convergence of numerical approximations of linear stochastic partial differential equations driven by additive Lévy noise. We first derive a representation formula for the error which we then apply to study space-time discretizations of the stochastic heat equation, a Volterra-type integro-differential equation, and the wave equation as examples. For twice continuously differentiable test functions with bounded second derivative (with an additional condition on the second derivative for the wave equation) the weak rate of convergence is found to be twice the strong rate. The results extend earlier work by two of the authors as we consider general square-integrable infinite-dimensional Lévy processes and do not require boundedness of the test functions and their first derivative. Furthermore, the present framework is applicable to both hyperbolic and parabolic equations, and even to stochastic Volterra integro-differential equations.

Monika Eisenmann, Mihály Kovács, Raphael Kruse et al.

In this paper we introduce a randomized version of the backward Euler method, that is applicable to stiff ordinary differential equations and nonlinear evolution equations with time-irregular coefficients. In the finite-dimensional case, we consider Carathéodory type functions satisfying a one-sided Lipschitz condition. After investigating the well-posedness and the stability properties of the randomized scheme, we prove the convergence to the exact solution with a rate of $0.5$ in the root-mean-square norm assuming only that the coefficient function is square integrable with respect to the temporal parameter. These results are then extended to the numerical solution of infinite-dimensional evolution equations under monotonicity and Lipschitz conditions. Here we consider a combination of the randomized backward Euler scheme with a Galerkin finite element method. We obtain error estimates that correspond to the regularity of the exact solution. The practicability of the randomized scheme is also illustrated through several numerical experiments.

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}-α+κ$.

David Bolin, Kristin Kirchner, Mihály Kovács

The numerical approximation of the solution to a stochastic partial differential equation with additive spatial white noise on a bounded domain is considered. The differential operator is assumed to be a fractional power of an integer order elliptic differential operator. The solution is approximated by means of a finite element discretization in space and a quadrature approximation of an integral representation of the fractional inverse from the Dunford-Taylor calculus. For the resulting approximation, a concise analysis of the weak error is performed. Specifically, for the class of twice continuously Fréchet differentiable functionals with second derivatives of polynomial growth, an explicit rate of weak convergence is derived, and it is shown that the component of the convergence rate stemming from the stochasticity is doubled compared to the corresponding strong rate. Numerical experiments for different functionals validate the theoretical results.

Jianbo Cui, Mihály Kovács, Derui Sheng

We study the numerical approximation of a class of degenerate parabolic stochastic partial differential equations on non-compact metric graphs, which naturally arise in the asymptotic analysis of Hamiltonian flows under small noise perturbations. The numerical discretization of these equations faces several challenges, including the non-compactness of the graph, the degeneracy of the differential operator near vertices, and the non-symmetry of the associated bilinear form. To address these issues, we propose a multi-step numerical strategy combining graph truncation, localized coefficient regularization, and finite element spatial discretization. By incorporating localization techniques, tightness arguments, and resolvent estimates, we establish the strong convergence of the proposed scheme in a weighted $L^2$-space. Our results provide a systematic methodology that is potentially extensible to more general non-compact graphs and degenerate operators.

Mihály Kovács, Gyula Molnár, Máté András Száraz

We consider Gaussian Random Fields on metric graphs defined implicitly as the stationary solution to a fractional SPDE driven by Gaussian white noise. Sampling from the finite element approximation requires the Cholesky factorization of the mass matrix, causing non-linear execution time explosions and massive memory fill-in on large graphs. Hence, we combine Neumann-Neumann graph decomposition with mass matrix lumping and demonstrate empirically, that our approach preserves exact theoretical convergence rates established in [8] while achieving multi-order speedups and massive memory reductions.

Boris Baeumer, Mihály Kovács, Mark M. Meerschaert et al.

This paper derives physically meaningful boundary conditions for fractional diffusion equations, using a mass balance approach. Numerical solutions are presented, and theoretical properties are reviewed, including well-posedness and steady state solutions. Absorbing and reflecting boundary conditions are considered, and illustrated through several examples. Reflecting boundary conditions involve fractional derivatives. The Caputo fractional derivative is shown to be unsuitable for modeling fractional diffusion, since the resulting boundary value problem is not positivity preserving.

Boris Baeumer, Mihály Kovács, Harish Sankaranarayanan

We identify the stochastic processes associated with one-sided fractional partial differential equations on a bounded domain with various boundary conditions. This is essential for modelling using spatial fractional derivatives. We show well-posedness of the associated Cauchy problems in $C_0(Ω)$ and $L_1(Ω)$. In order to do so we develop a new method of embedding finite state Markov processes into Feller processes and then show convergence of the respective Feller processes. This also gives a numerical approximation of the solution. The proof of well-posedness closes a gap in many numerical algorithm articles approximating solutions to fractional differential equations that use the Lax-Richtmyer Equivalence Theorem to prove convergence without checking well-posedness.