L1

In this talk I will discuss our recent work on two problems. The first problem concerns with capillary rise between rough structures, a fundamental wetting phenomenon that is functionalised in biological organisms and prevalent in geological or man-made materials. Predicting the liquid rise height is more complex than currently considered in the literature because it is necessary to couple two wetting phenomena: capillary rise and hemiwicking. Experiments, simulations and analytic theory demonstrate how this coupling challenges our conventional understanding and intuitions of wetting and roughness. For example, the critical contact angle for hemiwicking becomes separation-dependent so that hemiwicking can vanish for even highly wetting liquids. The rise heights for perfectly wetting liquids can also be different in smooth and rough systems. The second problem concerns with droplets (or condensates) formed via a liquid-liquid phase separation process in biological cells. Despite the widespread importance of surface tension for the interactions between these droplets and other cellular components, there is currently no reliable technique for their measurement in live cells. To address this, we develop a high-throughput flicker spectroscopy technique. Applying it to a class of cellular droplets known as stress granules, we find their interface fluctuations cannot be described by surface tension alone. It is necessary to consider elastic bending deformation and a non-spherical base shape, suggesting that stress granules are viscoelastic droplets with a structured interface, rather than simple Newtonian liquids. Moreover, given the broad distributions of surface tension and bending rigidity observed, different types of stress granules can only be differentiated via large-scale surveys, which was not possible previously and our technique now enables.

It is well known that stochasticity can play a fundamental role in various biochemical processes, such as cell regulatory networks and enzyme cascades. Isothermal, well-mixed systems can be adequately modeled by Markov processes and, for such systems, methods such as Gillespie’s algorithm are typically employed. While such schemes are easy to implement and are exact, the computational cost of simulating such systems can become prohibitive as the frequency of the reaction events increases. This has motivated numerous coarse-grained schemes, where the “fast” reactions are approximated either using Langevin dynamics or deterministically. While such approaches provide a good approximation for systems where all reactants are present in large concentrations, the approximation breaks down when the fast chemical species exist in small concentrations, giving rise to significant errors in the simulation. This is particularly problematic when using such methods to compute statistics of extinction times for chemical species, as well as computing observables of cell cycle models. In this talk, we present a hybrid scheme for simulating well-mixed stochastic kinetics, using Gillespie–type dynamics to simulate the network in regions of low reactant concentration, and chemical Langevin dynamics when the concentrations of all species are large. These two regimes are coupled via an intermediate region in which a “blended” jump-diffusion model is introduced. Examples of gene regulatory networks involving reactions occurring at multiple scales, as well as a cell-cycle model are simulated, using the exact and hybrid scheme, and compared, both in terms of weak error, as well as computational cost. If there is time, we will also discuss the extension of these methods for simulating spatial reaction kinetics models, blending together partial differential equation with compartment based approaches, as well as compartment based approaches with individual particle models.

This is joint work with Andrew Duncan (Imperial), Radek Erban (Oxford), Kit Yates (Bath), Adam George (Bath), Cameron Smith (Bath), Armand Jordana (New York )

Macroscopic Transport in Heterogeneous Porous Materials

Lucy Auton

Solute transport in porous materials is a key physical process in a wide variety of situations, including contaminant transport, filtration, lithium-ion batteries, hydrogeological systems, biofilms, bones and soils. Despite the prevalence of solute transport in porous materials, the effect of microstructure on flow and transport remains poorly understood and improving our understanding of this remains a major challenge.  In this presentation, I consider a two-dimensional microstructure comprising an array of obstacles of smooth but arbitrary shape, the size and spacing of which can vary along the length of the porous medium, allowing for anisotropy.  I use a nontrivial extension to classical homogenisation theory via the method of multiple scales to rigorously upscale the novel problem involving cells of varying area. This results in simple effective continuum equations for macroscale flow and transport where the effect of the microscale geometry on the macroscopic transport and removal is encoded within these simple macroscale equations via effective parameters such as an effective local anisotropic diffusivity and an effective local adsorption rate.  For a simple example geometry I exploit the two degrees of microstructural freedom in this problem, obstacle size and obstacle spacing, to investigate scenarios of uniform porosity but heterogenous microstructure, noting the impact this heterogeneity has on filter efficiency. 

This model constitutes the development of the core framework required to consider other crucial problems such as solute transport within soft porous materials for which there does not currently exist a simple macroscale model where the effective diffusivity and removal depend on the microstructure. Further, via this methodology I will  derive a  bespoke model for fluoride and arsenic removal filters. With this model I will be able to optimise the design of fluoride-removal filters which are being deployed across rural India. The design optimisation will both increase filter lifespan and reduce filter cost, enabling more people to access safe drinking water

Averaged interface conditions: evaporation fronts in porous media

Ellen Luckins

Homogenisation methods are powerful tools for deriving effective PDE models for processes incorporating multiple length-scales. For physical systems in which interface processes are crucial to the overall system, we might ask how the microstructure impacts the effective interface conditions, in addition to the PDEs in the bulk. In this talk we derive an effective model for the motion of an evaporation front through porous media, combining homogenisation and boundary layer analysis to derive averaged interface conditions at the evaporation front. Our analysis results in a new effective parameter in the boundary conditions, which encodes how the shape and speed of the porescale evaporating interfaces impact the overall drying process.

Liouville conformal field theory models two-dimensional gravity with a cosmological constant and conformal matter. In its timelike regime, it reproduces the characteristic negative kinetic term of the conformal factor of the metric in the Einstein-Hilbert action, the sign which infamously makes the gravity path integral ill-defined. In this talk, I will first discuss the perturbative computation of the timelike Liouville partition function around the sphere saddle and propose an all-orders result. I will then turn to the disk and present the bulk 1-point functions of this CFT, and discuss possible interpretations in terms of boundary conditions.

Statics and dynamics of droplets on lubricated surfaces

Zhaohe Dai

The abstract is "Slippery liquid infused porous surfaces are formed by coating surface with a thin layer of oil lubricant. This thin layer prevents other droplets from reaching the solid surface and allows such deposited droplets to move with ultra-low friction, leading to a range of applications. In this talk, we will discuss the static and dynamic behaviour of droplets placed on lubricated surfaces. We will show that the layer thickness and the size of the substrate are key parameters in determining the final equilibrium. However, the evolution towards the equilibrium is extremely slow (on the order of days for typical experimental parameter values). As a result, we suggest that most previous experiments with oil films lubricating smooth substrates are likely to have been in an evolving, albeit slowly evolving, transient state.

Modeling and Design Optimization for Pleated Membrane Filters

Yixuan Sun

Membrane filtration is widely used in many applications, ranging from industrial processes to everyday living activities. With growing interest from both industrial and academic sectors in understanding the various types of filtration processes in use, and in improving filter performance, the past few decades have seen significant research activity in this area. Experimental studies can be very valuable, but are expensive and time-consuming, therefore theoretical studies offer potential as a cost-effective and predictive way to improve on current filter designs. In this work, mathematical models, derived from first principles and simplified using asymptotic analysis, are proposed for pleated membrane filters, where the macroscale flow problem of Darcy flow through a pleated porous medium is coupled to the microscale fouling problem of particle transport and deposition within individual pores of the membrane. Asymptotically-simplified models are used to describe and evaluate the membrane performance numerically and filter design optimization problems are formulated and solved for industrially-relevant scenarios. This study demonstrates the potential of such modeling to guide industrial membrane filter design for a range of applications involving purification and separation.

We consider the four-point correlator of the stress-energy tensor in N=4 SYM, to leading order in inverse powers of the central charge, but including all order corrections in 1/lambda. This corresponds to the AdS version of the Virasoro-Shapiro amplitude to all orders in the small alpha'/low energy expansion. Using dispersion relations in Mellin space, we derive an infinite set of sum rules. These sum rules strongly constrain the form of the amplitude, and determine all coefficients in the low energy expansion in terms of the CFT data for heavy string operators, in principle available from integrability. For the first set of corrections to the flat space amplitude we find a unique solution consistent with the results from integrability and localisation.

Title: Topics on regularity theory for fully non-linear integro-differential equations

Abstract: We will focus on local and non-local Monge Ampere type equations, equations with deforming kernels and convex envelopes of functions with optimal special conditions. We discuss global solutions and their regularity properties.

Title: Quasilinear Conservative Collisional Transport in Kinetic Mean Field models

Abstract: We shall focus the on the interplay of nonlinear analysis  and numerical approximations to mean field models in particle physics where kinetic transport flows in momentum are strongly nonlinearly  modified by macroscopic quantities in classical or spectral density spaces. Two noteworthy models arise: the classical Fokker-Plank Landau dynamics as a low magnetized plasma regimes in the modeling of perturbative non-local high order terms. The other one corresponds to perturbation under strongly magnetized dynamics for fast electrons  in momentum space  give raise to a coupled system of classical kinetic diffusion processes described by the balance equations for electron probability density functions (electron pdf) coupled to the time dynamics on spectral energy waves  (quasi-particles) in a quantum process of their resonant interaction. Both models are rather different, yet there are derived form the Liouville-Maxwell system under different scaling. Analytical tools and some numerical  simulations show a presence of  strong hot tail anisotropy  formation taking the stationary states away from Classical equilibrium solutions stabilization for the iteration in a three dimensional cylindrical model. The semi-discrete schemes preserves the total system mass, momentum and energy, which are enforced by the numerical scheme. Error estimates can be obtained as well.

Work in collaboration with Clark Pennie and Kun Huang

Regularization of linear least squares problems is necessary in a variety of contexts. However, the optimal regularization parameter is usually unknown a priori and is often to be determined in an ad hoc manner, which may involve solving the problem for multiple regularization parameters. In this talk, we will discuss three randomized algorithms, building on the sketch-and-precondition framework in randomized numerical linear algebra (RNLA), to efficiently solve this set of problems. In particular, we consider preconditioners for a set of Tikhonov regularization problems to be solved iteratively. The first algorithm is a Cholesky-based algorithm employing a single sketch for multiple parameters; the second algorithm is SVD-based and improves the computational complexity by requiring a single decomposition of the sketch for multiple parameters. Finally, we introduce an algorithm capable of exploiting low-rank structure (specifically, low statistical dimension), requiring a single sketch and a single decomposition to compute multiple preconditioners with low-rank structure. This algorithm avoids the Gram matrix, resulting in improved stability as compared to related work.

Time-optimal control can be used to improve driving efficiency for autonomous
vehicles and it enables us explore vehicle and driver behaviour in extreme
situations. Due to the computational cost and limited scope of classical
optimal control methods we have seen new interest in applying reinforcement
learning algorithms to autonomous driving tasks.
In this talk we present methods for translating time-optimal vehicle control
problems into reinforcement learning environments. For this translation we
construct a sequence of environments, starting from the closest representation
of our optimisation problem, and gradually improve the environments reward
signal and feature quality. The trained agents we obtain are able to generalise
across different race tracks and obtain near optimal solutions, which can then
be used to speed up the solution of classical time-optimal control problems.

Many problems in engineering can be understood as controlling the bifurcation structure of a given device. For example, one may wish to delay the onset of instability, or bring forward a bifurcation to enable rapid switching between states. In this talk, we will describe a numerical technique for controlling the bifurcation diagram of a nonlinear partial differential equation by varying the shape of the domain or a parameter in the equation. Our aim is to delay or advance a given branch point to a target parameter value. The algorithm consists of solving an optimization problem constrained by an augmented system of equations that characterize the location of the branch points. The flexibility and robustness of the method also allow us to advance or delay a Hopf bifurcation to a target value of the bifurcation parameter, as well as controlling the oscillation frequency. We will apply this technique on systems arising from biology, fluid dynamics, and engineering, such as the FitzHugh-Nagumo model, Navier-Stokes, and hyperelasticity equations.