General order conditions for stochastic partitioned Runge-Kutta methods

In this paper stochastic partitioned Runge-Kutta (SPRK) methods are considered. A general order theory for SPRK methods based on stochastic B-series and multicolored, multishaped rooted trees is developed. The theory is applied to prove the order of some known methods, and it is shown how the number of order conditions can be reduced in some special cases, especially that the conditions for preserving quadratic invariants can be used as simplifying assumptions

Symplectic methods for Hamiltonian systems with additive noise

Stochastic systems, phase flows of which have integral invariants, are considered. Hamiltonian systems with additive noise being a wide class of such systems possess the property of preserving symplectic structure. For them, numerical methods preserving the symplectic structure are constructed. A special attention is paid to systems with separable Hamiltonians, to second order differential equations with additive noise, and to Hamiltonian systems with small additive noise

Variational Integrators and Generating Functions for Stochastic Hamiltonian Systems

In this work, the stochastic version of the variational principle is established, important for stochastic symplectic integration, and for structure-preserving algorithms of stochastic dynamical systems. Based on it, the stochastic variational integrators in formulation of stochastic Lagrangian functions are proposed, and some applications to symplectic integrations are given. Three types of generating functions in the cases of one and two noises are discussed for constructing new schemes

Stochastic discrete Hamiltonian variational integrators

Variational integrators are derived for structure-preserving simulation of stochastic Hamiltonian systems with a certain type of multiplicative noise arising in geometric mechanics. The derivation is based on a stochastic discrete Hamiltonian which approximates a type-II stochastic generating function for the stochastic flow of the Hamiltonian system. The generating function is obtained by introducing an appropriate stochastic action functional and its corresponding variational principle. Our approach permits to recast in a unified framework a number of integrators previously studied in the literature, and presents a general methodology to derive new structure-preserving numerical schemes. The resulting integrators are symplectic; they preserve integrals of motion related to Lie group symmetries; and they include stochastic symplectic Runge–Kutta methods as a special case. Several new low-stage stochastic symplectic methods of mean-square order 1.0 derived using this approach are presented and tested numerically to demonstrate their superior long-time numerical stability and energy behavior compared to nonsymplectic methods

Data-driven structure-preserving model reduction for stochastic Hamiltonian systems

In this work we demonstrate that SVD-based model reduction techniques known for ordinary differential equations, such as the proper orthogonal decomposition, can be extended to stochastic differential equations in order to reduce the computational cost arising from both the high dimension of the considered stochastic system and the large number of independent Monte Carlo runs. We also extend the proper symplectic decomposition method to stochastic Hamiltonian systems, both with and without external forcing, and argue that preserving the underlying symplectic or variational structures results in more accurate and stable solutions that conserve energy better than when the non-geometric approach is used. We validate our proposed techniques with numerical experiments for a semi-discretization of the stochastic nonlinear Schr\"odinger equation and the Kubo oscillator.Comment: 43 pages, 16 figure

Accurate stationary densities with partitioned numerical methods for stochastic partial differential equations

We consider the numerical solution, by finite differences, of second-order-in-time stochastic partial differential equations (SPDEs) in one space dimension. New timestepping methods are introduced by generalising recently-introduced methods for second-order-in-time stochastic differential equations to multidimensional systems. These stochastic methods, based on leapfrog and Runge–Kutta methods, are designed to give good approximations to the stationary variances and the correlations in the position and velocity variables. In particular, we introduce the reverse leapfrog method and stochastic Runge–Kutta Leapfrog methods, analyse their performance applied to linear SPDEs and perform numerical experiments to examine their accuracy applied to a type of nonlinear SPDE