In micro-fluidics, at small scales where inertial effects become negligible, surface to volume ratios are large and the interfacial processes are extremely important for the overall dynamics. Integral
equation based methods are attractive for the simulations of e.g. droplet-based microfluidics, with tiny water drops dispersed in oil, stabilized by surfactants. In boundary integral formulations for
Stokes flow, jumps in pressure and velocity gradients are naturally taken care of, viscosity ratios enter only in coefficients of the equations, and only the drop surfaces must be discretized and not the volume inside nor in between.

We present numerical methods for drops with insoluble surfactants, both in two and three dimensions. We discretize the integral equations using Nyström methods, and special care is taken in the evaluation of singular and also nearly singular integrals that is needed in the case of close drop interactions. A spectral method is used to solve the advection-diffusion equation on each drop surface that describes the evolution of surfactant concentration. The drop velocity and surfactant concentration couple together through an equation of state for the surface tension coefficient. An adaptive time-stepping strategy is developed for the coupled problem, with the constraint to minimize the number of Stokes solves, since this is the computationally most expensive part.

For high quality discretization of the drops throughout the simulations, a hybrid method is used in two dimensions, offering an arc-length parameterization of the interface. In three dimensions, a
reparameterization procedure is developed to optimize the spherical harmonics representation of the drop, while conserving the drop volume and amount of surfactant.

We present results from some validation tests and illustrate the ability of the numerical methods in different challenging problems.

Decision trees usefully represent the sparse, high dimensional and noisy nature of chemical data from experiments. Having learned a function from this data, we may want to thereafter optimise the function, e.g. for picking the best catalyst for a chemical process. This work studies a mixed-integer non-linear optimisation problem involving: (i) gradient boosted trees modelling catalyst behaviour, (ii) penalty functions mitigating risk, and (iii) penalties enforcing chemical composition constraints. We develop several heuristic methods to find feasible solutions, and an exact, branch and bound algorithm that leverages structural properties of the gradient boost trees and penalty functions. We computationally test our methods on an industrial instance from BASF.
This work was completed in collaboration with Mr Miten Mistry and Dr Dimitris Letsios at Imperial College London and Dr Robert Lee and Dr Gerhard Krennrich from BASF.
 

We present discrete methods for computing low-rank approximations of time-dependent tensors that are the solution of a differential equation. The approximation format can be Tucker, tensor trains, MPS or hierarchical tensors. We will consider two types of discrete integrators: projection methods based on quasi-optimal metric projection, and splitting methods based on inexact solutions of substeps. For both approaches we show numerically and theoretically that their behaviour is superior compared to standard methods applied to the so-called gauged equations. In particular, the error bounds are robust in the presence of small singular values of the tensor’s matricisations. Based on joint work with Emil Kieri, Christian Lubich, and Hanna Walach.

Solving the Stokes equation by an optimal domain decomposition method derived algebraically involves the use of non standard interface conditions whose discretisation is not trivial. For this reason the use of approximation methods such as hybrid discontinuous Galerkin appears as an appropriate strategy: on the one hand they provide the best compromise in terms of the number of degrees of freedom in between standard continuous and discontinuous Galerkin methods, and on the other hand the degrees of freedom used in the non standard interface conditions are naturally defined at the boundary between elements. In this work we introduce the coupling between a well chosen discretisation method (hybrid discontinuous Galerkin) and a novel and efficient domain decomposition method to solve the Stokes system. We present the detailed analysis of the hybrid discontinuous Galerkin method for the Stokes problem with non standard boundary conditions. This analysis is supported by numerical evidence. In addition, the advantage of the new preconditioners over more classical choices is also supported by numerical experiments.

This work was done in collaboration with G. Barrenechea, M. Bosy (Univ. Strathclyde) and F. Nataf, P-H Tournier (Univ of Paris VI)

Symplectic reduction is the natural quotient construction for symplectic manifolds. Given a free and proper action of a Lie group G on a symplectic manifold M, this process produces a new symplectic manifold of dimension dim(M) - 2 dim(G). For non-free actions, however, the result is usually fairly singular. But Sjamaar-Lerman (1991) showed that the singularities can be understood quite precisely: symplectic reductions by non-free actions are partitioned into smooth symplectic manifolds, and these manifolds fit nicely together in the sense that they form a stratification.

Symplectic reduction has an analogue in hyperkähler geometry, which has been a very important tool for constructing new examples of these special manifolds. In this talk, I will explain how Sjamaar-Lerman’s results can be extended to this setting, namely, hyperkähler quotients by non-free actions are stratified
spaces whose strata are hyperkähler.

Oxford Mathematics Public Lectures
Hooke Lecture

Michael Berry - Chasing the dragon: tidal bores in the UK and elsewhere
15 November 2018 - 5.15pm

In some of the world’s rivers, an incoming high tide can arrive as a smooth jump decorated by undulations, or as a breaking wave. The river reverses direction and flows upstream.

Understanding tidal bores involves

· analogies with tsunamis, rainbows, horizons in relativity, and ideas from  quantum physics;

· the concept of a ‘minimal model’ in mathematical explanation;

· different ways in which different cultures describe the same thing;

· the first unification in fundamental physics.

Michael Berry is Emeritus Professor of Physics, H H Wills Physics Laboratory, University of Bristol

5.15pm, Mathematical Institute, Oxford

Please email @email to register.

Watch live:

https://www.facebook.com/OxfordMathematics
https://livestream.com/oxuni/Berry

Oxford Mathematics Public Lectures are generously supported by XTX Markets.