Past Industrial and Applied Mathematics Seminar

 I will present a new framework for determining effectively the spectrum and stability of traveling waves on networks with symmetries, such as rings and lattices, by computing master stability curves (MSCs). Unlike traditional methods, MSCs are independent of system size and can be readily used to assess wave destabilization and multi-stability in small and large networks.

We apply techniques of exponential asymptotics to the KdV equation to derive the small-time behaviour for dispersive waves that propagate in one direction.  The results demonstrate how the amplitude, wavelength and speed of these waves depend on the strength and location of complex-plane singularities of the initial condition.  Using matched asymptotic expansions, we show how the small-time dynamics of complex singularities of the time-dependent solution are dictated by a Painlevé II problem with decreasing tritronquée solutions.  We relate these dynamics to the solution on the real line.

In this talk, we will explore three flow configurations that illustrate the behaviour of slow-moving viscous fluids in confined geometries: viscous gravity currents, fracturing of shear-thinning fluids in a Hele-Shaw cell, and rectangular channel flows of non-Newtonian fluids. We will first develop simple mathematical models to describe each setup, and then we will compare the theoretical predictions from these models with laboratory experiments. As is often the case, we will see that even models that are grounded in solid physical principles often fail to accurately predict the real-world flow behaviour. Our aim is to identify the primary physical mechanisms absent from the model using laboratory experiments. We will then refine the mathematical models and see whether better agreement between theory and experiment can be achieved.

Odd systems, characterised by broken time-reversal or parity symmetry, 
exhibit striking transport phenomena due to transverse responses. In this 
talk, I will introduce the concept of odd diffusion, a generalisation of 
diffusion in two-dimensional systems that incorporates antisymmetric tensor 
components. Focusing on systems of interacting particles, I adapt a 
geometric approach to derive effective transport equations and show how 
interactions give rise to unusual transport in odd systems. I present 
effects like enhanced self-diffusion, reversed Hall drift and even absolute 
negative mobility that solely originate in odd diffusion. These results 
reveal how microscopic symmetry-breaking gives rise to emergent, equilibrium 
and non-equilibrium transport, with implications for soft matter, chiral 
active systems, and topological materials.

Structural mechanics plays a crucial role in soft matter physics, mechanobiology, metamaterials, pattern formation, active matter, and soft robotics. What unites these seemingly disparate topics is the natural balance that emerges between elasticity, geometry, and stability. This seminar will serve as a high-level overview of our work on several problems concerning the stability of structures. I will cover three topics: (1) shapeshifting shells; (2) mechanical metamaterials; and (3) elastogranular mechanics.

I will begin by discussing our development of a generalized, stimuli-responsive shell theory. (1) Non-mechanical stimuli including heat, swelling, and growth further complicate the nonlinear mechanics of shells, as simultaneously solving multiple field equations to capture multiphysics phenomena requires significant computational expense. We present a general shell theory to account for non-mechanical stimuli, in which the effects of the stimuli are
generalized into three forms: those that add mass to the shell, those that increase the area of the shell through the natural stretch, and those that change the curvature of the shell through the natural curvature. I will show how this model can capture the morphogenesis of the optic cup, the snapping of the Venus flytrap, leaf growth, and the buckling of electrically active polymer plates. (2) I will then discuss how cutting thin sheets and shells, a process
inspired by the art of kirigami, enables the design of functional mechanical metamaterials. We create linear actuators, artificial muscles, soft robotic grippers, and mechanical logic units by systematically cutting and stretching thin sheets. (3) Finally, if time permits, I will introduce our work on the interactions between elastic and granular matter, which we refer to as elastogranular mechanics. Such interactions occur across all lengths, from morphogenesis, to root growth, to stabilizing soil against erosion. We show how combining rocks and string in the absence of any adhesive we can create large, load bearing structures like columns, beams, and arches. I will finish with a general phase diagram for elastogranular behavior.

Nowadays, the 3D ultrastructure of a fibrous tissue can be reconstructed in order to visualize the complex nanoscale arrangement of collagen fibrils including neighboring proteoglycans even in the stretched loaded state [1]. In particular, experimental data of collagen fibers in human artery layers have shown that the f ibers are not symmetrically dispersed [2]. In addition, it is known that collagen f ibers are cross-linked and the density of cross-links in arterial tissues has a stiffening effect on the associated mechanical response. A first attempt to characterize this effect on the elastic response is presented and the influence of the cross-link density on the mechanical behavior in uniaxial tension is shown [3]. A recently developed extension of the model that accounts for dispersed fibers connected by randomly distributed cross-links is outlined [4]. A simple shear test focusing on the sign of the normal stress perpendicular to the shear planes (Poynting effect) is analyzed. In [5] it was experimentally observed that, in contrast to rubber, semi-flexible biopolymer gels show a tendency to approach the top and bottom faces under simple shear. This so-called negative Poynting effect and its connection with the cross-links as well as the fiber and crosslink dispersion is also examined. 

References 

[1]A. Pukaluk et al.: An ultrastructural 3D reconstruction method for observing the arrangement of collagen fibrils and proteoglycans in the human aortic wall under mechanical load. Acta Biomaterialia, 141:300-314, 2022. 

 [2] G.A. Holzapfel et al.: Modelling non-symmetric collagen fibre dispersion in arterial walls. Journal of the Royal Society Interface, 12:20150188, 2015. 

 [3] G.A. Holzapfel and R.W. Ogden: An arterial constitutive model accounting for collagen content and cross-linking. Journal of the Mechanics and Physics of Solids, 136:103682, 2020. 

 [4] S. Teichtmeister and G.A. Holzapfel: A constitutive model for fibrous tissues with cross-linked collagen fibers including dispersion – with an analysis of the Poynting effect. Journal of the Mechanics and Physics of Solids, 164:104911, 2022. 

[5] P.A. Janmey et al.: Negative normal stress in semiflexible biopolymer gels. Nature Materials, 6:48–51, 2007.

Predictions of complex flows remain a significant challenge for engineering systems. Computationally affordable predictions of turbulent flows generally require Reynolds-Averaged Navier–Stokes (RANS) simulations and Large-Eddy Simulation (LES), the predictive accuracy of which can be insufficient due to non-Boussinesq turbulence and/or unresolved multiphysics that preclude qualitative fidelity in certain regimes. For example, in turbulent combustion, flame–turbulence interactions can lead to inverse-cascade energy transfer, which violates the assumptions of many RANS and LES closures. We present an adjoint-based, solver-embedded data assimilation method to augment the RANS and LES equations using trusted data. This is accomplished using Python-native flow solvers that leverage differentiable programming techniques to construct the adjoint equations needed for optimization. We present applications to shock-tube ignition delay predictions, turbulent premixed jet flames, and shock-dominated nonequilibrium flows and discuss the potential of adjoint-based approaches for future machine learning applications.

A fundamental question in neuroscience is to understand how information is represented in the activity of  tens of thousands of neurons in the brain. Towards this end, low-rank matrix and tensor decompositions are commonly used to identify correlates of behavior in high-dimensional neural data. In this talk I will first present a novel tensor decomposition based on the slice rank which is able to disentangle mixed modes of covarying patterns in data tensors. Second, to compliment this statistical approach, I will present our recent dynamical systems modelling of neural activity over learning. Rather than factorizing data tensors themselves, we instead fit a dynamical system to the data, while constraining the tensor of parameters to be low rank. Together these projects highlight how applications in neural data can inspire new classes of low-rank models.

Determining stationary compact configurations of vorticity described by the 2D Euler equations is a classic problem dating back to the late 19th century. The aim is to find equilibrium distributions of vorticity, in the form  of point vortices, vortex sheets, vortex patches, and hollow vortices. This endeavour has driven the development of mathematical and numerical techniques such as Hamiltonian vortex dynamics and contour dynamics.

In the case of vortex sheets, methods and results are presented for finding rotating equilibria, some in the presence of point vortices. To begin, a numerical approach based on that recently developed by Trefethen, Costa, Baddoo, and others for solving Laplace's equation in the complex plane by series and rational approximation is described. The method successfully reproduces the exact vortex sheet solutions found by O'Neil (2018) and Protas & Sakajo (2020). Some new solutions are found.

The numerical approach suggests an analytical method based on conformal mapping for finding exact closed-form vortex sheet equilibria. Examples are presented.

Finally, new numerical solutions are computed for steady, doubly-connected vortex layers of uniform vorticity surrounding a solid object such that the fluid velocity vanishes on the outer free boundary. While dynamically unrelated, these solutions have mathematical analogy and application to the industrial free boundary problem arising in the dip-coating of objects by a viscous fluid.

The beautiful displays exhibited by fish schools and bird flocks have long fascinated scientists, but the role of their complex behavior remains largely unknown. In particular, the influence of hydrodynamic interactions on schooling and flocking has been the subject of debate in the scientific literature. I will present a model for flapping wings that interact hydrodynamically in an inviscid fluid, wherein each wing is represented as a plate that executes a prescribed time-periodic kinematics. The model generalizes and extends thin-airfoil theory by assuming that the flapping amplitude is small, and permits consideration of multiple wings through the use of conformal maps and multiply-connected function theory. We find that the model predictions agree well with experimental data on freely-translating, flapping wings in a water tank. The results are then used to motivate a reduced-order model for the temporally nonlocal interactions between schooling wings, which consists of a system of nonlinear delay-differential equations. We obtain a PDE as the mean-field limit of these equations, which we find supports traveling wave solutions. Generally, our results indicate how hydrodynamics may mediate schooling and flocking behavior in biological contexts.

I will discuss two classes of effective, macroscopic models in elasticity: (i) 1D models applicable to thin structures, and (ii) homogenized 2D or 3D continua applicable to materials with a periodic microstructure. In both systems, the separation of scales calls for the definition of macroscopic models that slave fine-scale fluctuations to an effective, macroscopic deformation field. I will show how such models can be established in a systematic and rigorous way based on a two-scale expansion that accounts for nonlinear and higher-order (i.e. deformation gradient) effects. I will further demonstrate that the resulting models accurately predict nonlinear effects, finite size effects and localization for a set of examples. Finally, I will discuss two challenges that arise when solving these effective models: (1) missed boundary layer effects and (2) negative stiffness associated with higher-order terms.

We consider two problems in fluid dynamics: the collective locomotion of flying animals and the interaction of vortex rings with fluid interfaces. First, we present a model of formation flight, viewing the group as a material whose properties arise from the flow-mediated interactions among its members. This aerodynamic model explains how flapping flyers produce vortex wakes and how they are influenced by the wakes of others. Long in-line arrays show that the group behaves as a soft, excitable "crystal" with regularly ordered member "atoms" whose positioning is susceptible to deformations and dynamical instabilities. Second, we delve into the phenomenon of vortex ring reflections at water-air interfaces. Experimental observations reveal reflections analogous to total internal reflection of a light beam. We present a vortex-pair--vortex-sheet model to simulate this phenomenon, offering insights into the fundamental interactions of vortex rings with free surfaces.

Biological matter has the fascinating ability to autonomously generate material deformations via intrinsic active forces, where the latter are often present within effectively two-dimensional structures. The dynamics of such “active surfaces” inevitably entails a complex, self-organized interplay between geometry of a surface and its mechanical interactions with the surrounding. The impact of these factors on the self-organization capacity of surfaces made of an active material, and how related effects are exploited in biological systems, is largely unknown.

In this talk, I will first discuss general numerical challenges in analysing self-organising active surfaces and the bifurcation structure of emergent shape spaces. I will then focus on active surfaces with broken up-down symmetry, of which the eukaryotic cell cortex and epithelial tissues are highly abundant biological examples. In such surfaces, a natural interplay arises between active stresses and surface curvature. We demonstrate that this interplay leads to a comprehensive library of spontaneous shape transformations that resemble stereotypical morphogenetic processes. These include cell-division-like invaginations and the autonomous formation of tubular surfaces of arbitrary length, both of which robustly overcome well-known shape instabilities that would arise in analogue passive systems.

Active solids consume energy to allow for actuation and shape change not possible in equilibrium. I will first introduce active solids in comparison with their active fluid counterparts. I will then focus on active solids composed of non-reciprocal springs and show how so-called odd elastic moduli arise in these materials. Odd active solids have counter-intuitive elastic properties and require new design principles for optimal response. For example, in floppy lattices, zero modes couple to microscopic non-reciprocity, which destroys odd moduli entirely in a phenomenon reminiscent of rigidity percolation. Instead, an optimal odd lattice will be sufficiently soft to activate elastic deformations, but not too soft. These results provide a theoretical underpinning for recent experiments and point to the design of novel soft machines.

It is well known that adding even small amounts of  long chain polymers (e.g. few parts per million) to Newtonian solvents can drastically change the flow behaviour by introducing elasticity. In particular,  two decades ago, experiments in curved geometries  demonstrated  that polymer flows can be  chaotic even at vanishingly small Reynolds numbers. The situation in `straight’ flows  such as pressure-driven flow down a channel is less clear  and hence an area of current focus. I will discuss recent progress.

Tension-induced giant actuation in elastic sheets

Dr. Marc Suñé

Buckling is normally associated with a compressive load applied to a slender structure; from railway tracks in extreme heat to microtubules in cytoplasm, axial compression is relieved by out-of-plane buckling. However, recent studies have demonstrated that tension applied to structured thin sheets leads to transverse motion that may be harnessed for novel applications, such as kirigami grippers, multi-stable `groovy-sheets', and elastic ribbed sheets that close into tubes. Qualitatively similar behaviour has also been observed in simulations of thermalized graphene sheets, where clamping along one edge leads to tilting in the transverse direction. I will discuss how this counter-intuitive behaviour is, in fact, generic for thin sheets that have a relatively low stretching modulus compared to the bending modulus, which allows `giant actuation' with moderate strain.

The global COVID19 pandemic has highlighted the lethality and morbidity associated with infectious respiratory diseases. These diseases can lead to devastating syndrome known as the acute respiratory distress syndrome (ARDS) where bacterial/viral infections cause excessive lung inflammation, pulmonary edema, and severe hypoxemia (shortness of breath). Although ARDS patients require artificial mechanical ventilation, the complex biofluid and biomechanical forces generated by the ventilator exacerbates lung injury leading to high mortality. My group has used mathematical and computational modeling to both characterize the complex mechanics of lung injury during ventilation and to identify novel ways to prevent injury at the cellular level. We have used in-vitro and in-vivo studies to validate our mathematical predictions and have used engineering tools to understand the biological consequences of the mechanical forces generated during ventilation. In this talk I will specifically describe how our mathematical/computational approach has led to novel cytoskeletal based therapies and how coupling mathematics and molecular biology has led to the discovery of a gene regulatory mechanisms that can minimize ventilation induced lung injury. I will also describe how we are currently using nanotechnology and gene/drug delivery systems to enhance the lung’s native regulatory responses and thereby prevent lung injury during ARDS.

Slender elastic filaments with intrinsic helical geometry are encountered in a wide range of physical and biological settings, ranging from coil springs in engineering to bacteria flagellar filaments. The equilibrium configurations of helical filaments under a variety of loading types have been well studied in the framework of the Kirchhoff rod equations. These equations are geometrically nonlinear and so can account for large, global displacements of the rod. This geometric nonlinearity also makes a mathematical analysis of the rod equations extremely difficult, so that much is still unknown about the dynamic behaviour of helical rods under external loading.

An important class of simplified models consists of 'equivalent-column' theories. These model the helical filament as a naturally-straight beam (aligned with the helix axis) for which the extensional and torsional deformations are coupled. Such theories have long been used in engineering to describe the free vibrations of helical coil springs, though their validity remains unclear, particularly when distributed forces and moments are present. In this talk, we show how such an effective theory can be derived systematically from the Kirchhoff rod equations using the method of multiple scales. Importantly, our analysis is asymptotically exact in the small-wavelength limit and can account for large, unsteady displacements. We then illustrate our theory with two loading scenarios: (i) a heavy helical rod deforming under its own weight; and (ii) axial rotation (twirling) in viscous fluid, which may be considered as a simple model for a bacteria flagellar filament. More broadly, our analysis provides a framework to develop reduced models of helical rods in a wide variety of physical and biological settings, as well as yielding analytical insight into their tensile instabilities.

Even though most of us tie our shoelaces "wrongly," knots in ropes and filaments have been used as functional structures for millennia, from sailing and climbing to dewing and surgery. However, knowledge of the mechanics of physical knots is largely empirical, and there is much need for physics-based predictive models. Tight knots exhibit highly nonlinear and coupled behavior due to their intricate 3D geometry, large deformations, self-contact, friction, and even elasto-plasticity. Additionally, tight knots do not show separation of the relevant length scales, preventing the use of centerline-based rod models. In this talk, I will present an overview of recent work from our research group, combining precision experiments, Finite Element simulations, and theoretical analyses. First, we study the mechanics of two elastic fibers in frictional contact. Second, we explore several different knotted structures, including the overhand, figure-8, clove-hitch, and bowline knots. These knots serve various functions in practical settings, from shoelaces to climbing and sailing. Lastly, we focus on surgical knots, with a particularly high risk of failure in clinical settings, including complications such as massive bleeding or the unraveling of high-tension closures. Our research reveals a striking and robust power law, with a general exponent, between the mechanical strength of surgical knots, the applied pre-tension, and the number of throws, providing new insights into their operational and safety limits. These findings could have potential applications in the training of surgeons and enhanced control of robotic-assisted surgical devices.

In geophysical and astrophysical flows, we are often interested in understanding the impact of internal waves on the non-wavelike flow. For example, oceanic internal waves generated at the surface and the seafloor transfer energy from the large scale flow to dissipative scales, thereby influencing the global ocean state. A primary challenge in the study of wave-flow interactions is how to separate these processes – since waves and non-wavelike flows can vary on similar spatial and temporal scales in the Eulerian frame. However, in a Lagrangian flow-following frame, temporal filtering offers a convenient way to isolate waves. Here, I will discuss a recently developed method for evolving Lagrangian mean fields alongside the governing equations in a numerical simulation, and extend this theory to allow effective filtering of waves from non-wavelike processes.

Liquid crystal elastomers are rubbery solids containing molecular LC rods that align along a common director. On heating, the alignment is disrupted, leading to a substantial (~50%) contraction along the director.  In recent years, there has been a great deal of interest in fabrication LCE sheets with a bespoke alignment pattern. On heating, these patterns generate  corresponding patterns of contraction that can morph a sheet into a bespoke curved surface such as a cone or face. Moreover, LCEs can also be activated by light, either photothermally or photochemically, leading to similarly large contractions. Stimulation by light also introduces an important new possibility: using spatio-temporal patterns of illumination to morph a single LCE sample into a range of different surfaces. Such stimulation can enable non-reciprocal actuation for viscous swimming or pumping, and control over the whole path taken by the sheet through shape-space rather than just the final destination. In this talk, I will start by with an experimental example of a spatio-temporal pattern of illumination being used to actuate an LCE peristaltic pump. I will then introduce a second set of experiments, in which a monodomain sheet morphs first into a cone, an anti-cone and then an array of cones upon exposure to different patterns of illumination. Finally, I will then discuss the general problem of how to choose a pattern of illumination to morph a director-patterned sheet into an arbitrary surface, first analytically for axisymmetric cases, then numerically for low symmetry cases. This last study exceeds our current experimental capacity, but highlights how, with full spatio-temporal control over the stimulation magnitude, one can choreograph an LCE sheet to undergo almost any pattern of morphing.

applied mathematics.

Mathematical modelling can support decontamination processes in a variety of ways.  In this talk, we focus on the contamination step: understanding how much of a chemical spill has seeped into the Earth or a building material, and how far it has travelled, are essential for making good decisions about how to clean it up.  

We consider an infiltration problem in which a chemical is poured on an initially unsaturated porous medium, and seeps into it via capillary action. Capillarity-driven flow through partially-saturated porous media is often modelled using Richards’ equation, which is a simplification of the Buckingham-Darcy equation in the limit where the infiltrating phase is much more viscous than the receding phase.  In this talk, I will explore the limitations of Richards equation, and discuss some scenarios in which predictions for small-but-finite viscosity ratios are very different to the Richards simplification.

What do fiber polymers and ice sheets have in common? They both flow with a directionally dependent - anisotropic - viscosity. This behaviour occurs in other geophysical flows, such as the Earth's mantle, where a material's microstructure affects its large-scale flow. In ice, the alignment of crystal orientations can cause the viscosity to vary by an order of magnitude, consequently having a strong impact on the flow of ice sheets and glaciers. However, the effect of anisotropy on large-scale flow is not well understood, due to a lack of understanding of a) the best physical approximations to model crystal orientations, and b) how crystal orientations affect rheology. In this work, we aim to address both these questions by linking rheology to crystal orientation predictions, and testing a range of models against observations from the Greenland ice sheet. The results show assuming all grains experience approximately the same stress provides realistic predictions, and we suggest a set of equations and parameters which can be used in large-scale models of ice sheets. 

Traditionally, stochastic processes are modelled one of two ways: a continuum Fokker-Planck approach, where a PDE is solved to determine the time evolution of the probability density, or a Langevin approach, where the SDE describing the system is sampled, and multiple simulations are used to collect statistics. There is also a third way: the functional or path integral. Originally developed by Wiener in the 1920s to model Brownian motion, path integrals were famously applied to quantum mechanics by Feynman in the 1950s. However, they also have much to offer classical stochastic processes (and statistical physics).  

In this talk I will introduce the formalism at a physicist’s level of rigour, and focus on determining the dominant contribution to the path integral when the noise is weak. There exists a remarkable correspondence between the most-probable stochastic paths and Hamiltonian dynamics in an effective potential [1,2,3]. I will then discuss some applications, including reaction pathways conditioned on finite time [2]. We demonstrate that the most probable pathway at a finite time may be very different from the usual minimum energy path used to calculate the average reaction rate. If time permits, I will also discuss the extremely nonlinear crystal dislocation response to applied stress [4].  

[1] Ge, Hao, and Hong Qian. Int. J. Mod. Phys. B 26.24 1230012 (2012)     

[2] Fitzgerald, Steve, et al. J. Chem. Phys. 158.12 (2023).

[3] Honour, Tom and Fitzgerald, Steve. in press J. Phys. A (2024)

[4] Fitzgerald, Steve. Sci. Rep. 6 (1) 39708 (2016)

We develop a comprehensive geometrically-exact theory for an end-loaded elastic rod constrained to deform on a cylindrical surface. By viewing the rod-cylinder system as a special case of an elastic braid, we are able to obtain all forces and moments imparted by the deforming rod to the cylinder as well as all contact reactions. This framework allows us to give a complete treatment of static friction consistent with force and moment balance. In addition to the commonly considered model of hard frictionless contact, we analyse two friction models in which the rod, possibly with intrinsic curvature, experiences either lateral or tangential friction. As applications of the theory we study buckling of the constrained rod under compressive and torsional loads, finding critical loads to depend on Coulomb-like friction parameters, as well as the tendency of the rod to lift off the cylinder under further loading. The cylinder can also have arbitrary orientation relative to the direction of gravity. The cases of a horizontal and vertical cylinder, with gravity having only a lateral or axial component, are amenable to exact analysis, while numerical results map out the transition in buckling mechanism between the two extremes. Weight has a stabilising effect for near-horizontal cylinders, while for near-vertical cylinders it introduces the possibility of buckling purely due to self-weight. Our results are relevant for many engineering and medical applications in which a slender structure winds inside or outside a cylindrical boundary.

The Fokker-Planck equation is one of the major tools of statistical physics in the description of stochastic processes, with numerous applications in physics, chemistry and biology. In the case of heterogeneous diffusions, the formulation of the equation depends on the choice of the discretization of the stochastic integral in the underlying Langevin-equation due to the multiplicative noise. In the Fokker-Planck equation, the choice of discretization then enters as a parameter in the definition of drift and diffusion terms. I show how both short- and long-time limits are affected by this choice. In the long-time limit, the existence of normalizable probability distribution functions is not always guaranteed which can be remedied by invoking elements of infinite ergodic theory. 

[1] S. Giordano, F. Cleri, R. Blossey, Phys Rev E 107, 044111 (2023)

[2] T. Dupont, S. Giordano, F. Cleri, R. Blossey, arXiv:2401.01765 (2024)

Add to calendar

Global scale ocean currents are strongly constrained by the Earth’s rotation, while this effect is generally negligible at small scales. In between, motions with scales from 1-10km are marginally affected by the Earth’s rotation. These intermediate scales, collectively termed the ocean submesoscale, have been hidden from view until recent years. Evidence from field measurements, numerical models, and satellite data have shown that submesoscales play a particularly important role in the upper ocean where they help to control the transport of material between the ocean surface and interior. In this talk I will review some recent work on submesoscale dynamics and their influence on biogeochemistry and accumulation of microplastics in the surface waters.

We study “bottom-up” models for the collective motion of large groups of animals. Similar models can be encoded into (micro)robotic matter, capable of sensing light and processing information. Agents are endowed only with visual sensing and information processing. We study a model in which moving agents reorientate to maximise the path-entropy of their visual environment over paths into the future. There are general arguments that principles like this that are based on retaining freedom in the future may confer fitness in an uncertain world. Alternative “top-down” models are more common in the literature. These typically encode coalignment and/or cohesion directly and are often motivated by models drawn from physics, e.g. describing spin systems. However, such models can usually give little insight into how co-alignment and cohesion emerge because these properties are encoded in the model at the outset, in a top-down manner. We discuss how our model leads to dynamics with striking similarities with animal systems, including the emergence of coalignment, cohesion, a characteristic density scaling anddifferent behavioural phenotypes. The dynamics also supports a very unusual order-disorder transition in which the order (coalignment) initially increases upon the addition of sensory or behavioural noise, before decreasing as the noise becomes larger.

During the last fifteen years, there has been a paradigm shift in the continuum modelling of granular materials; most notably with the development of rheological models, such as the μ(I)-rheology (where μ is the friction and I is the inertial number), but also with significant advances in theories for particle segregation. This talk details theoretical and numerical frameworks (based on OpenFOAM®) which unify these disconnected endeavours. Coupling the segregation with the flow, and vice versa, is not only vital for a complete theory of granular materials, but is also beneficial for developing numerical methods to handle evolving free surfaces. This general approach is based on the partially regularized incompressible μ(I)-rheology, which is coupled to a theory for gravity/shear-driven segregation (Gray & Ancey, J. Fluid Mech., vol. 678, 2011, pp. 353–588). These advection–diffusion–segregation equations describe the evolving concentrations of the constituents, which then couple back to the variable viscosity in the incompressible Navier–Stokes equations. A novel feature of this approach is that any number of differently sized phases may be included, which may have disparate frictional properties. The model is used to simulate the complex particle-size segregation patterns that form in a partially filled triangular rotating drum. There are many other applications of the theory to industrial granular flows, which are the second most common material used after fluids. The same processes also occur in geophysical flows, such as snow avalanches, debris flows and dense pyroclastic flows. Depth-averaged models, that go beyond the μ(I)-rheology, will also be derived to capture spontaneous self-channelization and levee formation, as well as complex segregation-induced flow fingering effects, which enhance the run-out distance of these hazardous flows.

Recent advances in experimental techniques have enabled remarkable discoveries and insight into how the dynamics of thin gas/vapour films can profoundly influence the behaviour of liquid droplets: drops impacting solids can “skate on a film of air” [1], so that they can “bounce off walls” [2,3]; reductions in ambient gas pressure can suppress splashing [4] and initiate the merging of colliding droplets [5]; and evaporating droplets can levitate on their own vapour film [7] (the Leidenfrost effect). Despite these advances, the precise physical mechanisms governing these phenomena remains a topic of debate.  A theoretical approach would shed light on these issues, but due to the strongly multiscale nature of these processes brute force computation is infeasible.  Furthermore, when films reach the scale of the mean free path in the gas (i.e. ~100nm) and below, new nanoscale physics appears that renders the classical Navier-Stokes paradigm inaccurate.

In this talk, I will overview our development of efficient computational models for the aforementioned droplet dynamics in the presence of gas nanofilms into which gas-kinetic, van der Waals and/or evaporative effects can be easily incorporated [8,9].  It will be shown that these models can reproduce experimental observations – for example, the threshold between bouncing and wetting for drop impact on a solid is reproduced to within 5%, whilst a model excluding either gas-kinetic or van der Waals effects is ~170% off!  These models will then be exploited to make new experimentally-verifiable predictions, such as how we expect drops to behave in reduced pressure environments.  Finally, I will conclude with some exciting directions for future wor

[1] JM Kolinski et al, Phys. Rev. Lett.  108 (2012), 074503. [2] JM Kolinski et al, EPL.  108 (2014), 24001. [3] J de Ruiter et al, Nature Phys.  11 (2014), 48. [4] L Xu et al, Phys. Rev. Lett. 94 (2005), 184505. [5] J Qian & CK Law, J. Fluid. Mech. 331 (1997), 59.  [6] KL Pan J. Appl. Phys. 103 (2008), 064901. [7] D Quéré, Ann. Rev. Fluid Mech. 45 (2013), 197. [8] JE Sprittles, Phys. Rev. Lett.  118 (2017), 114502.  [9] MV Chubynsky et al, Phys. Rev. Lett.. 124 (2020), 084501.

Predictive scoring modelling is a common approach to measure financial health and credit worthiness in banking. Whilst the latter is a key factor in making decisions for lending, evaluating financial health helps to identify vulnerable customers trending towards financial hardship, who need support. The current macroeconomic uncertainties amplify the importance of extensive flexibility in modelling data solutions so that they can remain effective and adaptive to a volatile economic environment. This workshop is focused on discussing relevant techniques and mathematical methodologies which can help modernise traditional scoring models and accelerate innovation. In summary, the problem definition in the banking context is how automation and optimality can be achieved in a multiobjective problem where a subset of existing data features should be selected by relevance and uniqueness, assigned scoring weights by importance and how a pool of customers can be categorised accordingly using their individual scores and auto-adjusting thresholds of risk classification scales. The key challenge is imposed by the mutual dependency of the three sub-problems and their objectives. Introducing or removing constraints in any of them can change the feasibility and optimality of the others and the overall solution. It is common for traditional scoring models to be mainly focused on the predictive accuracy and their setup is often defined and revised manually, following ad-hoc exploratory data analysis and business-led decision making. An automated optimisation of the data features’ selection, scoring weights and classification thresholds definition can achieve respectively: ▪ Precise financial health evaluation and book classification under changing economic climate; ▪ Development of innovative data-driven solutions to enhance prevention from financial hardship and bankruptcies.

Capillary forces acting at the surface of a liquid drop can be strong enough to deform small objects and recent studies have provided several examples of elastic instabilities induced by surface tension. Inspired by the windlass mechanism in spider webs, we present a system where a liquid drop sits on a straight fiber and attracts the fiber which thereby coils inside the drop. We then introduce a 2D extension of the mechanism and build a membrane that can extend/contract by a factor of 20.

In order for artificial neural networks to learn a task, one must solve an inverse design problem. What network will produce the desired output? We have harnessed AI approaches to design physical systems to perform functions inspired by biology, such as protein allostery. But artificial neural networks require a computer in order to learn in top-down fashion by the global process of gradient descent on a cost function. By contrast, the brain learns by local rules on its own, with each neuron adjusting itself and its synapses without knowing what all the other neurons are doing, and without the aid of an external computer. But the brain is not the only biological system that learns by local rules; I will argue that the actin cortex and the amnioserosa during the dorsal closure stage of Drosophila development can also be viewed this way.

The key role of physical and mechanical interactions in cancer emerges from a very large variety of data sources and methods - from genomics to bioimaging, from proteomics to clinical records. Thus, learning physics-driven relational information is crucial to characterize its progression at different scales.

In this talk I will discuss how mathematical and computational tools allow for learning  and better understanding of  the mechano-biology of cancer,  thanks to the integration of  patient-specific data and physics-based models. I will present a few applications developed in the last decade in which the development of  digital twins,  empowered by ad-hoc learning tools,  allows us to test new hypotheses,  to assess the model predictions against biological and clinical data, and to aid decision-making in a clinical setting.

Funding from MUR - PRIN 2020, Progetto di Eccellenza 2023-2027 and Regione Lombardia (NEWMED Grant, ID: 117599, POR FESR 2014-2020) is gratefully acknowledged.

We formulate and solve a generalized inverse Navier–Stokes boundary value problem for velocity field reconstruction and simultaneous boundary segmentation of noisy Flow-MRI velocity images. We use a Bayesian framework that combines CFD, Gaussian processes, adjoint methods, and shape optimization in a unified and rigorous manner.
With this framework, we find the velocity field and flow boundaries (i.e. the digital twin) that are most likely to have produced a given noisy image. We also calculate the posterior covariances of the unknown parameters and thereby deduce the uncertainty in the reconstructed flow. First, we verify this method on synthetic noisy images of flows. Then we apply it to experimental phase contrast magnetic resonance (PC-MRI) images of an axisymmetric flow at low and high SNRs. We show that this method successfully reconstructs and segments the low SNR images, producing noiseless velocity fields that match the high SNR images, using 30 times less data.
This framework also provides additional flow information, such as the pressure field and wall shear stress, accurately and with known error bounds. We demonstrate this further on a 3-D in-vitro flow through a 3D-printed aorta and 3-D in-vivo flow through a carotid artery.