A spectral-based numerical method for Kolmogorov equations in Hilbert spaces

We propose a numerical solution for the solution of the Fokker-Planck-Kolmogorov (FPK) equations associated with stochastic partial differential equations in Hilbert spaces. The method is based on the spectral decomposition of the Ornstein-Uhlenbeck semigroup associated to the Kolmogorov equation. This allows us to write the solution of the Kolmogorov equation as a deterministic version of the Wiener-Chaos Expansion. By using this expansion we reformulate the Kolmogorov equation as a infinite system of ordinary differential equations, and by truncation it we set a linear finite system of differential equations. The solution of such system allow us to build an approximation to the solution of the Kolmogorov equations. We test the numerical method with the Kolmogorov equations associated with a stochastic diffusion equation, a Fisher-KPP stochastic equation and a stochastic Burgers Eq. in dimension 1.Comment: 28 pages, 10 figure

Stability of numerical method for semi-linear stochastic pantograph differential equations

Abstract As a particular expression of stochastic delay differential equations, stochastic pantograph differential equations have been widely used in nonlinear dynamics, quantum mechanics, and electrodynamics. In this paper, we mainly study the stability of analytical solutions and numerical solutions of semi-linear stochastic pantograph differential equations. Some suitable conditions for the mean-square stability of an analytical solution are obtained. Then we proved the general mean-square stability of the exponential Euler method for a numerical solution of semi-linear stochastic pantograph differential equations, that is, if an analytical solution is stable, then the exponential Euler method applied to the system is mean-square stable for arbitrary step-size h > 0 $h>0$ . Numerical examples further illustrate the obtained theoretical results

Numerical Solution of Stochastic Partial Differential Equations with Correlated Noise

In this paper we investigate the numerical solution of stochastic partial differential equations (SPDEs) for a wider class of stochastic equations. We focus on non-diagonal colored noise instead of the usual space-time white noise. By applying a spectral Galerkin method for spatial discretization and a numerical scheme in time introduced by Jentzen $\&$ Kloeden, we obtain the rate of path-wise convergence in the uniform topology. The main assumptions are either uniform bounds on the spectral Galerkin approximation or uniform bounds on the numerical data. Numerical examples illustrate the theoretically predicted convergence rate

On the stochastic Magnus expansion and its application to SPDEs

We derive a stochastic version of the Magnus expansion for the solution of linear systems of It\^o stochastic differential equations (SDEs). The goal of this paper is twofold. First, we prove existence and a representation formula for the logarithm associated to the solution of the matrix-valued SDEs. Second, we propose a new method for the numerical solution of stochastic partial differential equations (SPDEs) based on spatial discretization and application of the stochastic Magnus expansion. A notable feature of the method is that it is fully parallelizable. We also present numerical tests in order to asses the accuracy of the numerical schemes