We present a numerical investigation of stochastic transport for the damped and driven incompressible 2D Euler fluid flows. According to Holm (Proc Roy Soc, 2015) and Cotter et al. (2017), the principles of transformation theory and multi-time homogenisation, respectively, imply a physically meaningful, data-driven approach for decomposing the fluid transport velocity into its drift and stochastic parts, for a certain class of fluid flows. We develop a new methodology to implement this velocity decomposition and then numerically integrate the resulting stochastic partial differential equation using a finite element discretisation. We show our numerical method is consistent.

Numerically, we perform the following analyses on this velocity decomposition. We first perform uncertainty quantification tests on the Lagrangian trajectories by comparing an ensemble of realisations of Lagrangian trajectories driven by the stochastic differential equation, and the Lagrangian trajectory driven by the ordinary differential equation. We then perform uncertainty quantification tests on the resulting stochastic partial differential equation by comparing the coarse-grid realisations of solutions of the stochastic partial differential equation with the ``true solutions'' of the deterministic fluid partial differential equation, computed on a refined grid. In these experiments, we also investigate the effect of varying the ensemble size and the number of prescribed stochastic terms. Further experiments are done to show the uncertainty quantification results "converge" to the truth, as the spatial resolution of the coarse grid is refined, implying our methodology is consistent. The uncertainty quantification tests are supplemented by analysing the L2 distance between the SPDE solution ensemble and the PDE solution. Statistical tests are also done on the distribution of the solutions of the stochastic partial differential equation. The numerical results confirm the suitability of the new methodology for decomposing the fluid transport velocity into its drift and stochastic parts, in the case of damped and driven incompressible 2D Euler fluid flows. This is the first step of a larger data assimilation project which we are embarking on. This is joint work with Colin Cotter, Dan Crisan, Darryl Holm and Igor Shevchenko.

# Stochastic Analysis Seminar

Please note that the list below only shows forthcoming events, which may not include regular events that have not yet been entered for the forthcoming term. Please see the past events page for a list of all seminar series that the department has on offer.

We present a support theorem for subcritical parabolic stochastic partial differential equations (SPDEs) driven by Gaussian noises. In the spirit of the classical theorem by Stroock and Varadhan for ordinary stochastic differential equations, we identify the support of the solution to singular SPDEs with the closure of the union of the support of solutions to approximate and renormalized equations. We implement our approach in the setting of regularity structures and obtain a general result covering a range of singular SPDEs (including $\Phi^4_3$, $\Phi^d_2$, KPZ, PAM (2D+3D), SHE, ...). As a Corollary to our result we obtain the uniqueness of invariant measures for various interesting SPDEs. This is a joint work with Martin Hairer.

Consider dY(t)=f(X(t))dX(t), where X(t) is a pure jump Levy process with finite p-variation norm, 1<= p < 2, and f is a Lipchitz continuous function. Following the geometric solution construction of Levy-driven stochastic differential equations in (Williams 2001), we develop a class of epsilon-strong simulation algorithms that allows us to construct a probability space, supporting both the geometric solution Y and a fully simulatable process Y_epsilon, such that Y_epsilon is within epsilon distance from Y under the uniform metric on compact time intervals with probability 1. Moreover, the users can adaptively choose epsilon’ < epsilon, so that Y_epsilon’ can be constructed conditional on Y_epsilon. This tolerance-enforcement feature allows us to easily combine our algorithm with Multilevel Monte Carlo for efficient estimation of expectations, and adding as a benefit a straightforward analysis of rates of convergence. This is joint with Jose Blanchet, Fei He and Offer Kella.

McKean-Vlasov SDEs are SDEs where the coefficients depend on the law of the solution to the SDE. Their interest is in the links with nonlinear PDEs on one side (the SDE-related Fokker-Planck equation is nonlinear) and with interacting particles on the other side: the McKean-Vlasov SDE be approximated by a system of weakly coupled SDEs. In this talk we consider McKean-Vlasov SDEs with irregular drift: though well-posedness for this SDE is not known, we show a large deviation principle for the corresponding interacting particle system. This implies, in particular, that any limit point of the particle system solves the McKean-Vlasov SDE. The proof combines rough paths techniques and an extended Vanrdhan lemma.

This is a joint work with Thomas Holding.