Past

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.