Thu, 09 Jun 2016

16:00 - 17:00
L1

IAM Group Meeting

Javier Buldu, Dave Hewett
Abstract

Dave Hewett: Canonical solutions in wave scattering

By a "canonical solution" I have in mind a closed-form exact solution of the scalar wave equation in a simple geometry, for example the exterior of a circular cylinder, or the exterior of an infinite wedge. In this talk I hope to convince you that the study of such problems is (a) interesting; (b) important; and (c) a rich source of (difficult) open problems involving eigenfunction expansions, special functions, the asymptotic evaluation of integrals, and matched asymptotic expansions.

 

Thu, 02 Jun 2016

16:00 - 17:00
L3

The spreading of a surfactant-laden drop down an inclined and pre-wetted substrate - Numerics, Asymptotics and Linear Stability Analysis

Shailesh Naire
(Keele)
Abstract

Surfactants are chemicals that adsorb onto the air-liquid interface and lower the surface tension there. Non-uniformities in surfactant concentration result in surface tension gradients leading to a surface shear stress, known as a Marangoni stress. This stress, if sufficiently large, can influence the flow at the interface.

Surfactants are ubiquitous in many aspects of technology and industry to control the wetting properties of liquids due to  their ability to modify surface tension. They are used in detergents, crop spraying, coating processes and oil recovery. Surfactants also occur naturally, for example in the mammalian lung. They reduce the surface tension within the liquid lining the airways, which assists in preventing the collapse of the smaller airways. In the lungs of premature infants, the quantity of surfactant produced is insufficient as the lungs are under- developed. This leads to a respiratory distress syndrome which is treated by Surfactant Replacement Therapy.

Motivated by this medical application, we theoretically investigate a model problem involving the spreading of a drop laden with an insoluble surfactant down an inclined and pre-wetted substrate.  Our focus is in understanding the mechanisms behind a “fingering” instability observed experimentally during the spreading process. High-resolution numerics reveal a multi-region asymptotic wave-like structure of the spreading droplet. Approximate solutions for each region is then derived using asymptotic analysis. In particular, a quasi-steady similarity solution is obtained for the leading edge of the droplet. A linear stability analysis of this region shows that the base state is linearly unstable to long-wavelength perturbations. The Marangoni effect is shown to be the dominant driving mechanism behind this instability at small wavenumbers. A small wavenumber stability criterion is derived and it's implication on the onset of the fingering instability will be discussed.

Thu, 26 May 2016

16:00 - 17:00
L3

IAM Group Meeting

Mason Porter, Robert Van Gorder
Abstract

A Simple Generative Model of Collective Online Behavior (Mason Porter)

Human activities increasingly take place in online environments, providing novel opportunities for relating individual behaviors to population-level outcomes. In this paper, we introduce a simple generative model for the collective behavior of millions of social networking site users who are deciding between different software applications. Our model incorporates two distinct mechanisms: one is associated with recent decisions of users, and the other reflects the cumulative popularity of each application. Importantly, although various combinations of the two mechanisms yield long-time behav- ior that is consistent with data, the only models that reproduce the observed temporal dynamics are those that strongly emphasize the recent popularity of applications over their cumulative popularity.

This demonstrates --- even when using purely observational data with- out experimental design --- that temporal data-driven modeling can effectively distinguish between competing microscopic mechanisms, allowing us to uncover previously unidentified aspects of collective online behavior.

---

Bubbles, Turing machines, and possible routes to Navier-Stokes blow-up (Robert van Gorder)

Navier-Stokes existence and regularity in three spatial dimensions for an incompressible fluid... is hard. Indeed, while the original equations date back to the 1840's, existence and regularity remains an open problem and is one of the six remaining Millennium Prize Problems in mathematics that were stated by the Clay Mathematics Institute in 2000. Despite the difficulty, a resolution to this problem may say little about real-world fluids, as many real fluid problems do not seem to blow-up, anyway.
In this talk, we shall briefly outline the mathematical problem, although our focus shall be on the negative direction; in particular, we focus on the possibility of blow-up solutions. We show that many existing blow-up solutions require infinite energy initially, which is unreasonable. Therefore, obtaining a blow-up solution that starts out with nice properties such as bounded energy on three dimensional Euclidean space is rather challenging. However, if we modify the problem, there are some results. We survey recent results on averaged Navier-Stokes equations and compressible Navier-Stokes equations, and this will take us anywhere from bubbles to fluid Turing machines. We discuss how such results might give insight into the loss of regularity in the incompressible case (or, insight into how hard it might be to loose regularity of solutions when starting with finite energy in the incompressible case), before philosophizing about whether mathematical blow-up solutions could ever be physically relevant.

Thu, 19 May 2016

16:00 - 17:00
L3

Formulating short-range elastic interactions between dislocations in a continuum framework

Yichao Zhu
(Hong Kong University of Science and Technology)
Abstract

Permanent deformations of crystalline materials are known to be carried out by a large
number of atomistic line defects, i.e. dislocations. For specimens on micron scales or above, it
is more computationally tractable to investigate macroscopic material properties based on the
evolution of underlying dislocation densities. However, classical models of dislocation
continua struggle to resolve short-range elastic interactions of dislocations, which are believed
responsible for the formation of various heterogeneous dislocation substructures in crystals. In
this talk, we start with discussion on formulating the collective behaviour of a row of
dislocation dipoles, which would be considered equivalent to a dislocation-free state in
classical continuum models. It is shown that the underlying discrete dislocation dynamics can
be asymptotically captured by a set of evolution equations for dislocation densities along with
a set of equilibrium equations for variables characterising the self-sustained dislocation
substructures residing on a shorter length scale, and the strength of the dislocation
substructures is associated with the solvability conditions of their governing equilibrium
equations. Under the same strategy, a (continuum) flow stress formula for multi-slip systems
is also derived, and the formula resolves more details from the underlying dynamics than the
ubiquitously adopted Taylor-type formulae.

Thu, 12 May 2016

16:00 - 17:00
L3

Cancelled - Mathematical Problems within the Analysis of Transport Data

Eddie Wilson
(University of Bristol)
Abstract

My main purpose in this talk is try and convey a sense of my enthusiasm for mathematical modelling generally and how I've come to use it in a range of transport applications. For concreteness, I am going to talk in particular about work I have been doing on EPSRC grant EP/K000438/1 (PI: Jillian Anable, Aberdeen) where we are using the UK government's Department for Transport MOT data to estimate mileage totals and study how they are broken down across the population in various different ways. Embedded inside this practical problem is a whole set of miniature mathematical puzzles and challenges which are quite particular to the problem area itself, and one wider question which is rather deeper and more general: whether it is possible (and how) to convert usage data that is low-resolution in time but high-resolution in individuals to knowledge that is high-resolution in time but only expressed at a population level.

Thu, 05 May 2016

16:00 - 17:00
L3

Singular asymptotics of surface-plasmon resonance

Ory Schnitzer
(Imperial College London)
Abstract

Surface plasmons are collective electron-density oscillations at a metal-dielectric interface. In particular, highly localised surface-plasmon modes of nanometallic structures with narrow nonmetallic gaps, which enable a tuneable resonance frequency and a giant near-field enhancement, are at the heart of numerous nanophotonics applications. In this work, we elucidate the singular near-contact asymptotics of the plasmonic eigenvalue problem governing the resonant frequencies and modes of such structures. In the classical regime, valid for gap widths > 1nm, we find a generic scaling describing the redshift of the resonance frequency as the gap width is reduced, and in several prototypical dimer configurations derive explicit expressions for the plasmonic eigenvalues and eigenmodes using matched asymptotic expansions; we also derive expressions describing the resonant excitation of such modes by light based on a weak-dissipation limit. In the subnanometric ``nonlocal’’ regime, we show intuitively and by systematic analysis of the hydrodynamic Drude model that nonlocality manifests itself as a potential discontinuity, and in the near-contact limit equivalently as a widening of the gap. We thereby find the near-contact asymptotics as a renormalisation of the local asymptotics, and in particular a lower bound on plasmon frequency, scaling with the 1/4 power of the Fermi wavelength. Joint work with Vincenzo Giannini, Richard V. Craster and Stefan A. Maier. 

Thu, 28 Apr 2016

16:00 - 17:00
L3

Mathematics and Molecular Biology: The Engineering Approach

Bob Eisenberg
(Rush University)
Abstract

Life is different because it is inherited. All life comes from a blueprint (genes) that can only make proteins. Proteins are studied by more than one hundred thousand scientists and physicians every day because they are so important in health and disease. The function of proteins is on the macroscopic scale, but atomic details control that function, as is shown in a multitude of experiments. The structure of proteins is so important that governments spend billions studying them. Structures are known in exquisite detail determined by crystallographic measurement of more than 105 different proteins. But the forces that govern the movement and function of proteins are not visible in the structure. Mathematics is needed to compute both function and forces so comparison with experiment can be made. Experiments report numbers, typically sets of numbers in the form of graphs. Verbal models, however beautifully written in the biological tradition, do not provide numerical outputs, and so it is difficult to tell which verbal model better fits data.

The mathematics of molecular biology must be multiscale because atomic details control macroscopic function. The device approach of the engineering and English physiological tradition provides the dimensional reduction needed to solve the multiscale problem. Mathematical analysis of hundreds of experiments (reported in some fifty papers) has been successful in showing how some properties of an important class of proteins—ion channels— work. Ion channels are natural nanovalves as important to animals as Field Effect Transistors (FETs) are to computers. I will present the Fermi Poisson approach started by Jinn Liang Liu. The Fermi distribution is used to describe the saturation of space produced by crowded spherical ions. The Poisson equation (and continuity of current) is used to describe long range electrodynamics. Short range correlations are approximated by the Santangelo equation. A fully consistent mathematical description reproduces macroscopic properties of bulk solutions of sodium and calcium chloride solutions. It also describes several different channels (with quite different atomic detailed structures) quite well in a wide range of conditions using a handful of parameters never changed. It is not clear why the model works as well it does, nor is it clear how well the model will work on other channels, transporters or proteins.

Thu, 10 Mar 2016

16:00 - 17:00
L3

Dynamic homogenisation

Richard Craster
(ICL)
Abstract

The aim of this talk is to describe effective media for wave propagation through periodic, or nearly periodic, composites. Homogenisation methods are well-known and developed for quasi-static and low frequency regimes. The aim here is to move to situations of more practical interest where the frequencies are high, in some sense, and to compare the results of the theory with large scale simulations.

Thu, 03 Mar 2016

16:00 - 17:00
L3

Non-linear continuum models for planar extensible beams and pantographic lattices of beams: Heuristic homogenization, experimental and numerical examples of equilibrium in large deformation

Francesco dell'Isola
(Universita di Roma)
Abstract
There are relatively few results in the literature of non-linear beam theory: we recall here the very first classical results by Euler–Bernoulli and the researches stemming from von Kármán for moderately large rotations but small strains. In this paper, we consider a discretized springs model for extensible beams and propose a heuristic homogenization technique of the kind first used by Piola. The homogenized energy obtained has some peculiar features which we start to describe by solving numerically some exemplary deformation problems. Furthermore we consider pantographic structures constituted by the introduced nonlinear beams and study some planar deformation problems. Numerical solutions for these 2D problems are obtained via minimization of energy and are compared via some experimental measurements, in which the importance of elongation phenomena are clearly pointed out. In the conclusions we indicate a list of some mathematical problems which seems worth of consideration. 
 
Indeed Already Piola in 1848 introduces for microscopically discrete systems to be described via a continuum model: i) the micro-macro kinematical map, ii) the identification of micro- macro work functional and iii) the consequent determination of macro-constitutive equations in terms of the micro properties of considered mechanical system.
 
Piola uses, following the standards of his time, a rigorous mathematical deduction process and considers separately one dimensional, two dimensional and three dimensional continua as continua whose reference configuration is a curve, a surface or a regular connected subset of Euclidean three dimensional space. This subdivision of the presented matter is also followed by Cosserat Brothers: how to detect the influence on their works exerted by Piola’s pioneering ones is a historical problem which deserves further in-depth studies.
 
In the present paper we follow the spirit of Piola while looking for Lagrange density functions for a class of non-linear one dimensional continua in planar motion: we focus on modeling phenomena in which both extensional and bending deformations are of relevance.
 
Usually in literature the simultaneous extension and bending deformation of a beam is not considered: however when considering two dimensional continua embedding families of fibers as a model of some specific microstructured mechanical systems (as fiber fabrics or pantographic sheets ) the assumption that the fibers cannot extend while bending is not phenomenologically well-grounded. Therefore, we are led in the second part of the present paper to present some two dimensional continua in which the second gradient of in plane displacement (involving so called geodesic bending) appears in the expression of deformation energy.
 
The modeling assumptions are, in both cases, based on a physically reasonable discrete microstructure of used beams: in engineering literature these microstructures, constituted by extensional and rotational springs and possibly rigid bars, were introduced in order to get discrete Lagrangian approximation of continuum models in linearized regimes.
 
A natural development, involving the study of spatial placements of one dimensional or two dimensional continua or the introduction of three dimensional continua embedding reinforcement fibers will be subject of further investigations.
 
The study of pantographic sheets by means of a micro model based on Cauchy first gradient continuum models involves the choice of relatively small length scale, implying the introduction of numerical models involving finite elements with several millions of degrees of freedom: the computational burden of such models makes their use, at least in the mid term horizon, absolutely inappropriate. The higher gradient reduced order model presented in this paper involves a rather more effective numerical modeling whose performances (as will be shown in a forthcoming paper Giorgio et al. in preparation) are however absolutely comparable.
 
However the problem of formulating intermediate meso modeling, involving a class of Generalised Beam Theories, will be necessarily to be confronted: for instance the deformation of beam sections involving warping, Poisson effects, elastic necking or large shear or twist deformation can definitively be studied via reduced order models not resorting to the most detailed micro Cauchy first gradient models.
 
One should also remark that higher gradient continuum models may require novel integration schemes, more suitable to their intrinsic structure: we expect that isogeometric methods may further increase the effectiveness of the reduced models we present here, especially when completely spatial models will be considered .
Thu, 25 Feb 2016

16:00 - 17:00
L3

Acrobatics of Liquid Ropes

Neil Ribe
(CNRS and Universite Paris-Sud)
Abstract

Honey poured from a sufficient height onto toast undergoes the well-known `liquid rope coiling’ instability.

We have studied this instability using a combination of laboratory experiments, theory, and numerics, with the aim of determining phase diagrams and scaling laws for the different coiling modes. Finite-amplitude coiling has four distinct modes - viscous, gravitational, inertio-gravitational, and inertial - depending on how the viscous forces that resist deformation of the rope are balanced. The inertio-gravitational mode is particularly interesting as it involves resonance between the coiling portion of the rope and its long trailing `tail’. Further experiments using less viscous fluids reveal that the rope can exhibit five different morphologies, of which steady coiling is only one. We determine the detailed phase diagram of these morphologies, which includes a novel `liquid supercoiling’

state in which the coiled cylinder formed by the primary coiling instability undergoes in turn its own complex buckling instability.  We show that the onset of these different patterns is determined by a non-penetrability condition which takes different forms in the viscous, gravitational and inertial limits. To close, we will briefly evoke two additional related phenomena: spiral waves of bubbles generated by coiling, and the `fluid mechanical sewing machine’ in which the fluid falls onto a moving belt.

Subscribe to