This talk is a personal how-to (and how-not-to) manual for doing Maths with industry, or indeed with government. The Maths element is essential but lots of other skills and activities are equally necessary. Examples: problem elicitation; understanding the environmental constraints; power analysis; understanding world-views and aligning personal motivations; and finally, understanding the wider systems in which the Maths element will sit. These issues have been discussed for some time in the management science community, where their generic umbrella name is Problem Structuring Methods (PSMs).

I will present a nonlinear version of the open mapping principle which applies to constant-coefficient PDEs which are both homogeneous and weak* stable. An example of such a PDE is the Jacobian equation. I will discuss the consequences of such a result for the Jacobian and its relevance towards an answer to a long-standing problem due to Coifman, Lions, Meyer and Semmes. This is based on joint work with Lukas Koch and Sauli Lindberg.

/

I present a general nonlinear open mapping principle suited to applications to scale-invariant PDEs in regularity regimes where the equations are stable under weak* convergence. As an application I show that, for any $p < \infty$, the set of initial data for which there are dissipative weak solutions in $L^p_t L^2_x$ is meagre in the space of solenoidal L^2 fields. This is based on joint work with A. Guerra (Oxford) and S. Lindberg (Aalto).

We analyze the steady motion of a viscous incompressible fluid in a two- and three-dimensional channel containing an obstacle through the Navier-Stokes equations under different types of boundary conditions. In the 2D case we take constant non-homogeneous Dirichlet boundary data in a (virtual) square containing the obstacle, and emphasize the connection between the appearance of lift and the unique solvability of Navier-Stokes equations. In the 3D case we consider mixed boundary conditions: the inflow is given by a fairly general datum and the flow is assumed to satisfy a constant traction boundary condition on the outlet. In the absence of external forcing, explicit bounds on the inflow velocity guaranteeing existence and uniqueness of such steady motion are provided after estimating some Sobolev embedding constants and constructing a suitable solenoidal extension of the inlet velocity. In the 3D case, this solenoidal extension is built through the Bogovskii operator and explicit bounds on its Dirichlet norm (in terms of the geometric parameters of the obstacle) are found by solving a variational problem involving the infinity-Laplacian.

The talk accounts for results obtained in collaboration with Filippo Gazzola and Ilaria Fragalà (both at Politecnico di Milano).

Solving kinetic or related models with high-dimensional random parameters has been a challenging problem. In this talk, we will discuss how to employ the bi-fidelity stochastic collocation and choose efficient low-fidelity models in order to solve a class of multi-scale kinetic equations with uncertainties, including the Boltzmann equation, linear transport and the Vlasov-Poisson equation. In addition, some error analysis for the bi-fidelity method based on these PDEs will be presented. Finally, several numerical examples are shown to validate the efficiency and accuracy of the proposed method.

In the study of geometric flows it is often important to understand when a flow which converges along a sequence of times going to infinity will, in fact, converge along every such sequence of times to the same limit. While examples of finite dimensional gradient flows that asymptote to a circle of critical points show that this cannot hold in general, a positive result can be obtained in the presence of a so-called Lojasiewicz-Simon inequality. In this talk we will introduce this problem of uniqueness of asymptotic limits and discuss joint work with Melanie Rupflin and Peter M. Topping in which we examined the situation for a geometric flow that is designed to evolve a map describing a closed surface in a given target manifold into a parametrization of a minimal surface.

/

The Landau equation is an important PDE in kinetic theory modelling plasma particles in a gas. It can be derived as a limiting process from the famous Boltzmann equation. From the mathematical point of view, the Landau equation can be very challenging to study; many partial results require, for example, stochastic analysis as well as a delicate combination of kinetic and parabolic theory. The major open question is uniqueness in the physically relevant Coulomb case. I will present joint work with Jose Carrillo, Matias Delgadino, and Laurent Desvillettes where we cast the Landau equation as a generalized gradient flow from the optimal transportation perspective motivated by analogous results on the Boltzmann equation. A direct outcome of this is a numerical scheme for the Landau equation in the spirit of de Giorgi and Jordan, Kinderlehrer, and Otto. An extended area of investigation is to use the powerful gradient flow techniques to resolve some of the open problems and recover known results.

Trading of financial instruments has largely moved away from floor trading and onto electronic exchanges. Orders to buy and sell are queued at these exchanges in a limit-order book. While a full analysis of the dynamics of a limit-order book requires an understanding of strategic play among multiple agents, and is thus extremely complex, so-called zero-intelligence Poisson models have been shown to capture many of the statistical features of limit-order book evolution. These models can be addressed by traditional queueing theory techniques, including Laplace transform analysis. In this work, we demonstrate in a simple setting that another queueing theory technique, approximating the Poisson model by a diffusion model identified as the limit of a sequence of scaled Poisson models, can also be implemented. We identify the diffusion limit, find an embedded semi-Markov model in the limit, and determine the statistics of the embedded semi-Markov model. Along the way, we introduce and study a new type of process, a generalization of skew Brownian motion that we call two-speed Brownian motion.

Abstract: A recent paradigm views deep neural networks as discretizations of certain controlled ordinary differential equations, sometimes called neural ordinary differential equations. We make use of this perspective to link expressiveness of deep networks to the notion of controllability of dynamical systems. Using this connection, we study an expressiveness property that we call universal interpolation, and show that it is generic in a certain sense. The universal interpolation property is slightly weaker than universal approximation, and disentangles supervised learning on finite training sets from generalization properties. We also show that universal interpolation holds for certain deep neural networks even if large numbers of parameters are left untrained, and are instead chosen randomly. This lends theoretical support to the observation that training with random initialization can be successful even when most parameters are largely unchanged through the training. Our results also explore what a minimal amount of trainable parameters in neural ordinary differential equations could be without giving up on expressiveness.

Joint work with Martin Larsson, Josef Teichmann.

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please send email to @email.

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please send email to @email.

I will present a simple model of market microstructure which explains the concavity of price impact. In the proposed model, the local relationship between the order flow and the fundamental price (i.e. the local price impact) is linear, with a constant slope, which makes the model dynamically consistent. Nevertheless, the expected impact on midprice from a large sequence of co-directional trades is nonlinear and asymptotically concave. The main practical conclusion of the model is that, throughout a meta-order, the volumes at the best bid and ask prices change (on average) in favor of the executor. This conclusion, in turn, relies on two more concrete predictions of the model, one of which can be tested using publicly available market data and does not require the (difficult to obtain) information about meta-orders. I will present the theoretical results and will support them with the empirical analysis.