High weak order methods for stochastic equations based on modified equations

Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (mean-square stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

High order numerical approximation of the invariant measure of ergodic SDEs

Assyr Abdulle, Gilles Vilmart

We introduce new sufficient conditions for a numerical method to approximate with high order of accuracy the invariant measure of an ergodic system of stochastic differential equations, independently of the weak order of accuracy of the method. We then present a systematic procedure based on the framework of modified differential equations for the construction of stochastic integrators that capture the invariant measure of a wide class of ergodic SDEs (Brownian and Langevin dynamics) with an accuracy independent of the weak order of the underlying method. Numerical experiments confirm our theoretical findings

Society for Industrial and Applied Mathematics2014

Numerical methods for multiscale problems

Assyr Abdulle

The subject of this workshop was numerical methods that preserve geometric properties of the flow of an ordinary or partial differential equation: symplectic and multisymplectic integrators for Hamiltonian systems, symmetric integrators for reversible systems, methods preserving first integrals and numerical methods on manifolds, inclusive of Lie group methods and integrators for constrained mechanical systems, and methods for problems with highly oscillatory solutions. The unifying theme was structure preservation: not just the how?' but also why?', where?' and what for?'. The motivation for developing structure-preserving algorithms for special classes of problems originates independently in such diverse areas of research as astronomy, molecular dynamics, mechanics, control theory, theoretical physics and numerical analysis, with important contributions from other areas of both applied and pure mathematics. Moreover, it turns out that preservation of geometric properties of the flow not only produces an improved qualitative behaviour, but also allows for a significantly more accurate long-time integration than with general-purpose methods. In addition to the construction of geometric integrators, an important aspect of geometric integration is the light it sheds on the relationship between geometric properties of a numerical method and favourable error propagation in long-time integration. A very successful organising principle is backward error analysis, whereby the numerical one-step map is interpreted as (almost) the flow of a modified differential equation. In this way, geometric properties of the numerical integrator translate seamlessly into structure preservation on the level of the modified equation. Much insight and rigourous error estimates over long time intervals can then be obtained by combining backward error analysis with the KAM theory and related perturbation theories for Hamiltonian and reversible systems. While this approach has been very successful for ordinary differential equations, much less is currently known about highly oscillatory systems and geometric integrators for partial differential equations. Geometric numerical integration has been an active interdisciplinary research area since the last decade. Although the subject is in a lively phase of intensive development, the results so far are substantive and they impact on a wide range of application areas and on our understanding of core issues in computational mathematics. This is evidenced by the monographs \cite{HLW:GNI2002,LR:SMH2004}.

EMS2006

Analysis and numerical simulation of viscoelastic flows

Andrea Bonito

Mathematical and numerical aspects of viscoelastic flows are investigated here. Two simplified mathematical models are considered. They are motivated by a splitting algorithm for solving viscoelastic flows with free surfaces. The first model is a simplified Oldroyd-B model. Existence on a fixed time interval is proved in several Banach spaces provided the data are small enough. Short time existence is also proved for arbitrarily large data in Hölder spaces for the time variable. These results are based on the maximal regularity property of the Stokes operator and on the analycity behavior of the corresponding semi-group. A finite element discretization in space is then proposed. Existence of the numerical solution is proved for small data, as well as a priori error estimates, using an implicit function theorem framework. Then, the extension of these results to a stochastic simplified Hookean dumbbells model is discussed. Because of the presence of the Brownian motion, existence in a fixed time interval, provided the data are small enough, is proved only in some of the Banach spaces considered previously. The dumbbells' elongation is split in two parts, one satisfying a standart stochastic differential equation, the other satisfying a partial differential equation with a stochastic source term. A finite element discretization in space is also proposed. Existence of the numerical solution is proved for small data, as well as a priori error estimates. A numerical algorithm for solving viscoelastic flows with free surfaces is also described. This algorithm is based on a splitting method in time and two different meshes are used for the space discretization. Convergence of the numerical model is checked for the pure extensional flow and the filling of a pipe. Then, numerical results are reported for the stretching of a filament and for jet buckling.

EPFL2006

Numerical methods for multiscale problems

Assyr Abdulle

EMS2006