scholarly journals Weak convergence rates of splitting schemes for the stochastic Allen–Cahn equation

2019 ◽  
Vol 60 (3) ◽  
pp. 543-582 ◽  
Author(s):  
Charles-Edouard Bréhier ◽  
Ludovic Goudenège
2018 ◽  
Vol 39 (4) ◽  
pp. 2096-2134 ◽  
Author(s):  
Charles-Edouard Bréhier ◽  
Jianbo Cui ◽  
Jialin Hong

Abstract This article analyses an explicit temporal splitting numerical scheme for the stochastic Allen–Cahn equation driven by additive noise in a bounded spatial domain with smooth boundary in dimension $d\leqslant 3$. The splitting strategy is combined with an exponential Euler scheme of an auxiliary problem. When $d=1$ and the driving noise is a space–time white noise we first show some a priori estimates of this splitting scheme. Using the monotonicity of the drift nonlinearity we then prove that under very mild assumptions on the initial data this scheme achieves the optimal strong convergence rate $\mathcal{O}(\delta t^{\frac 14})$. When $d\leqslant 3$ and the driving noise possesses some regularity in space we study exponential integrability properties of the exact and numerical solutions. Finally, in dimension $d=1$, these properties are used to prove that the splitting scheme has a strong convergence rate $\mathcal{O}(\delta t)$.


2013 ◽  
Vol 2013 ◽  
pp. 1-11 ◽  
Author(s):  
Jürgen Geiser

We present some operator splitting methods improved by the use of the Zassenhaus product and designed for applications to multiphysics problems. We treat iterative splitting methods that can be improved by means of the Zassenhaus product formula, which is a sequential splitting scheme. The main idea for reducing the computation time needed by the iterative scheme is to embed fast and cheap Zassenhaus product schemes, since the computation of the commutators involved is very cheap, since we are dealing with nilpotent matrices. We discuss the coupling ideas of iterative and sequential splitting techniques and their convergence. While the iterative splitting schemes converge slowly in their first iterative steps, we improve the initial convergence rates by embedding the Zassenhaus product formula. The applications are to multiphysics problems in fluid dynamics. We consider phase models in computational fluid dynamics and analyse how to obtain higher order operator splitting methods based on the Zassenhaus product. The computational benefits derive from the use of sparse matrices, which arise from the spatial discretisation of the underlying partial differential equations. Since the Zassenhaus formula requires nearly constant CPU time due to its sparse commutators, we have accelerated the iterative splitting schemes.


2014 ◽  
Vol 46 (04) ◽  
pp. 1059-1083 ◽  
Author(s):  
Qifan Song ◽  
Mingqi Wu ◽  
Faming Liang

In this paper we establish the theory of weak convergence (toward a normal distribution) for both single-chain and population stochastic approximation Markov chain Monte Carlo (MCMC) algorithms (SAMCMC algorithms). Based on the theory, we give an explicit ratio of convergence rates for the population SAMCMC algorithm and the single-chain SAMCMC algorithm. Our results provide a theoretic guarantee that the population SAMCMC algorithms are asymptotically more efficient than the single-chain SAMCMC algorithms when the gain factor sequence decreases slower than O(1 / t), where t indexes the number of iterations. This is of interest for practical applications.


2020 ◽  
Vol 54 (6) ◽  
pp. 2199-2227
Author(s):  
Mihály Kovács ◽  
Annika Lang ◽  
Andreas Petersson

The numerical approximation of the mild solution to a semilinear stochastic wave equation driven by additive noise is considered. A standard finite element method is employed for the spatial approximation and a a rational approximation of the exponential function for the temporal approximation. First, strong convergence of this approximation in both positive and negative order norms is proven. With the help of Malliavin calculus techniques this result is then used to deduce weak convergence rates for the class of twice continuously differentiable test functions with polynomially bounded derivatives. Under appropriate assumptions on the parameters of the equation, the weak rate is found to be essentially twice the strong rate. This extends earlier work by one of the authors to the semilinear setting. Numerical simulations illustrate the theoretical results.


Sign in / Sign up

Export Citation Format

Share Document