Search

We consider the long-time behavior of an explicit tamed exponential Euler scheme applied to a class of parabolic semilinear stochastic partial differential equations driven by additive noise, under a one-sided Lipschitz continuity condition. The setting encompasses nonlinearities with polynomial growth. First, we prove that moment bounds for the numerical scheme hold, with at most polynomial dependence with respect to the time horizon. Second, we apply this result to obtain error estimates, in the weak sense, in terms of the time-step size and of the time horizon, to quantify the error to approximate averages with respect to the invariant distribution of the continuous-time process. We justify the efficiency of using the explicit tamed exponential Euler scheme to approximate the invariant distribution, since the computational cost does not suffer from the at most polynomial growth of the moment bounds. To the best of our knowledge, this is the first result in the literature concerning the approximation of the invariant distribution for SPDEs with non-globally Lipschitz coefficients using an explicit tamed scheme.

We study a class of stochastic semilinear damped wave equations driven by additive Wiener noise. Owing to the damping term, under appropriate conditions on the nonlinearity, the solution admits a unique invariant distribution. We apply semi-discrete and fully-discrete methods in order to approximate this invariant distribution, using a spectral Galerkin method and an exponential Euler integrator for spatial and temporal discretization respectively. We prove that the considered numerical schemes also admit unique invariant distributions, and we prove error estimates between the approximate and exact invariant distributions, with identification of the orders of convergence. To the best of our knowledge this is the first result in the literature concerning numerical approximation of invariant distributions for stochastic damped wave equations.

We construct a positivity-preserving Lie--Trotter splitting scheme with finite difference discretization in space for approximating the solutions to a class of nonlinear stochastic heat equations with multiplicative space-time white noise. We prove that this explicit numerical scheme converges in the mean-square sense, with rate $1/4$ in time and rate $1/2$ in space, under appropriate CFL conditions. Numerical experiments illustrate the superiority of the proposed numerical scheme compared with standard numerical methods which do not preserve positivity.

We introduce a positivity-preserving numerical scheme for a class of nonlinear stochastic heat equations driven by a purely time-dependent Brownian motion. The construction is inspired by a recent preprint by the authors where one-dimensional equations driven by space-time white noise are considered. The objective of this paper is to illustrate the properties of the proposed integrators in a different framework, by numerical experiments and by giving convergence results.

We propose a modification of the standard linear implicit Euler integrator for the weak approximation of parabolic semilinear stochastic PDEs driven by additive space-time white noise. The new method can easily be combined with a finite difference method for the spatial discretization. The proposed method is shown to have improved qualitative properties compared with the standard method. First, for any time-step size, the spatial regularity of the solution is preserved, at all times. Second, the proposed method preserves the Gaussian invariant distribution of the infinite dimensional Ornstein--Uhlenbeck process obtained when the nonlinearity is absent, for any time-step size. The weak order of convergence of the proposed method is shown to be equal to $1/2$ in a general setting, like for the standard Euler scheme. A stronger weak approximation result is obtained when considering the approximation of a Gibbs invariant distribution, when the nonlinearity is a gradient: one obtains an approximation in total variation distance of order $1/2$, which does not hold for the standard method. This is the first result of this type in the literature. A key point in the analysis is the interpretation of the proposed modified Euler scheme as the accelerated exponential Euler scheme applied to a modified stochastic evolution equation. Finally, it is shown that the proposed method can be applied to design an asymptotic preserving scheme for a class of slow-fast multiscale systems, and to construct a Markov Chain Monte Carlo method which is well-defined in infinite dimension. We also revisit the analysis of the standard and the accelerated exponential Euler scheme, and we prove new results with approximation in the total variation distance, which serve to illustrate the behavior of the proposed modified Euler scheme.

Explicit stabilized integrators are an efficient alternative to implicit or semi-implicit methods to avoid the severe timestep restriction faced by standard explicit integrators applied to stiff diffusion problems. In this paper, we provide a fully discrete strong convergence analysis of a family of explicit stabilized methods coupled with finite element methods for a class of parabolic semilinear deterministic and stochastic partial differential equations. Numerical experiments including the semilinear stochastic heat equation with space-time white noise confirm the theoretical findings.

Explicit stabilized integrators are an efficient alternative to implicit or semi-implicit methods to avoid the severe timestep restriction faced by standard explicit integrators applied to stiff diffusion problems. In this paper, we provide a fully discrete strong convergence analysis of a family of explicit stabilized methods coupled with finite element methods for a class of parabolic semilinear deterministic and stochastic partial differential equations. Numerical experiments including the semilinear stochastic heat equation with space-time white noise confirm the theoretical findings.