Gabriel Lord

Lyonell Boulton, Gabriel Lord

For $q>12/11$ the eigenfunctions of the non-linear eigenvalue problem associated to the one-dimensional $q$-Laplacian are known to form a Riesz basis of $L^2(0,1)$. We examine in this paper the approximation properties of this family of functions and its dual, in order to establish non-orthogonal spectral methods for the $p$-Poisson boundary value problem and its corresponding parabolic time evolution initial value problem. The principal objective of our analysis is the determination of optimal values of $q$ for which the best approximation is achieved for a given $p$ problem.

Daniel Stone, Sebastian Geiger, Gabriel Lord

A new class of asynchronous discrete-event simulation schemes for advection-diffusion-reaction equations are introduced, which is based on the principle of allowing quanta of mass to pass through faces of a Cartesian finite volume grid. The timescales of these events are linked to the flux on the the face, and the schemes are self-adaptive, local in time and space. Experiments are performed on realistic physical systems related to porous media flow applications, including a large 3D advection diffusion equation and advection diffusion reaction systems. The results are compared to highly accurate results where the temporal evolution is computed with exponential integrator schemes using the same finite volume discretisation. This allows a reliable estimation of the solution error. Our results indicate a first order convergence of the error as a control parameter is decreased.

Gabriel Lord, Simon J. A. Malham, Anke Wiese

We present numerical schemes for the strong solution of linear stochastic differential equations driven by an arbitrary number of Wiener processes. These schemes are based on the Neumann (stochastic Taylor) and Magnus expansions. Firstly, we consider the case when the governing linear diffusion vector fields commute with each other, but not with the linear drift vector field. We prove that numerical methods based on the Magnus expansion are more accurate in the mean-square sense than corresponding stochastic Taylor integration schemes. Secondly, we derive the maximal rate of convergence for arbitrary multi-dimensional stochastic integrals approximated by their conditional expectations. Consequently, for general nonlinear stochastic differential equations with non-commuting vector fields, we deduce explicit formulae for the relation between error and computational costs for methods of arbitrary order. Thirdly, we consider the consequences in two numerical studies, one of which is an application arising in stochastic linear-quadratic optimal control.

Daniel Stone, Gabriel Lord

A new numerical scheme for conservation equations based on evolution by asynchronous discrete events is presented. During each event of the scheme only two cells of the underlying Cartesian grid are active, and an event is processed as the exact evolution of this subsystem. This naturally leads to and adaptive scheme in space and time. Numerical results are presented which show that the error of the asynchronous scheme decreases to zero as a control parameter is reduced. The construction of the scheme allows it to be expressed as repeated multiplications of matrix exponentials on an initial state vector; thus techniques such as the Goldberg series and the Baker Campbell Hausdorff (BCH) formula can be used to explore the theoretical properties of the scheme. We present the framework of a convergence proof in this manner.

Daniel Stone, Gabriel Lord

Exponential integrators are time stepping schemes which exactly solve the linear part of a semilinear ODE system. This class of schemes requires the approxima- tion of a matrix exponential in every step, and one successful modern method is the Krylov subspace projection method. We investigate the effect of breaking down a single timestep into arbitrary multiple substeps, recycling the Krylov subspace to minimise costs. For these recyling based schemes we analyse the lo- cal error, investigate them numerically and show they can be applied to a large system with 106 unknowns. We also propose a new second order integrator that is found using the extra information from the substeps to form a corrector to increase the overall order of the scheme. This scheme is seen to compare favourably with other order two integrators.