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.

Thu, 18 Feb 2016

16:00 - 17:00
L3

Interactions of noise and discontinuities: transitions and qualitative changes

Rachel Kuske
(University of British Colombia)
Abstract

While there have been recent advances for analyzing the complex deterministic
behavior of systems with discontinuous dynamics, there are many open questions about
understanding and predicting noise-driven and noise-sensitive phenomena in the
non-smooth context.  Stochastic effects can often change the picture dramatically,
particularly if multiple time scales are present.  We demonstrate novel approaches
for exploring and explaining surprising phenomena driven by the interplay of
nonlinearities, delays, randomness, in specific applications with piecewise smooth
dynamics - nonlinear models of balance,  relay control, and impacting dynamics.
Effective techniques typically depend on the combination of mathematical techniques,
multiple scales techniques, and phenomenological intuition from seemingly unrelated
canonical models of biophysics, mechanics, and chemical dynamics.  The appropriate
strategy is not always immediately obvious from the area of application or model
type. This gap may follow from the limited attention that stochastic models with
discontinuous dynamics have received in the past, or it may be the reason for this
limited attention.  Combining the geometrical perspective with asymptotic approaches
in physical and phase space appears to be a critical part of developing effective
approaches.

Thu, 11 Feb 2016

16:00 - 17:00
L3

Wave-particle coupling in fluid mechanics: bouncing droplets and flapping swimmers

Anand Oza
Abstract
Roughly a decade ago, Yves Couder and coworkers demonstrated that droplets walking on a vibrating fluid bath exhibit several features previously thought to be peculiar to the microscopic quantum realm, including single-particle diffraction, tunneling, quantized orbits, and wave-like statistics in a corral. We here develop an integro-differential trajectory equation for these walking droplets with a view to gaining insight into their subtle dynamics. We then rationalize the emergence of orbital quantization in a rotating frame by assessing the stability of the orbital solutions. In the limit of large vibrational forcing, the chaotic walker dynamics gives rise to a coherent statistical behavior with wave-like features.
 
I will then describe recent efforts to model the dynamics of interacting flapping swimmers. Our study is motivated by recent experiments using a one-dimensional array of wings in a water tank, in which the system adopts “schooling modes” characterized by specific spatial phase relationships between swimmers. We develop a discrete dynamical system that models the swimmers as airfoils shedding point vortices, and study the existence and stability of steady solutions. We expect that our model may be used to understand how schooling behavior is influenced by hydrodynamics in more general contexts.
 
Thu, 04 Feb 2016

16:00 - 17:00
L3

Group Meeting

Barbara Mahler, Thomas Woolley, Julian A. Garcia Grajales
Abstract

Barbara Mahler: 15+5 min

Thomas Woolley: 15+5 min

Julian A. Garcia Grajales: 15+5 min
 

Thu, 28 Jan 2016

16:00 - 17:00
L3

Predictive simulations for optimisation of inhaled drug delivery

Laura Nicolaou
(ICL)
Abstract

Respiratory illnesses, such as asthma and chronic obstructive pulmonary disease, account for one in five deaths worldwide and cost the UK over £6 billion a year. The main form of treatment is via inhaled drug delivery. Typically, however, a low fraction of the inhaled dose reaches the target areas in the lung. Predictive numerical capabilities have the potential for significant impact in the optimisation of pulmonary drug delivery. However, accurate and efficient prediction is challenging due to the complexity of the airway geometries and of the flow in the airways. In addition, geometric variation of the airways across subjects has a pronounced effect on the aerosol deposition. Therefore, an accurate model of respiratory deposition remains a challenge.

High-fidelity simulations of the flow field and prediction of the deposition patterns motivate the use of direct numerical simulations (DNS) in order to resolve the flow. Due to the high grid resolution requirements, it is desirable to adopt an efficient computational strategy. We employ a robust immersed boundary method developed for curvilinear coordinates, which allows the use of structured grids to model the complex patient-specific airways, and can accommodate the inter-subject geometric variations on the same grid. The proposed approach reduces the errors at the boundary and retains the stability guarantees of the original flow solver.

A Lagrangian particle tracking scheme is adopted to model the transport of aerosol particles. In order to characterise deposition, we propose the use of an instantaneous Stokes number based on the local properties of the flow field. The effective Stokes number is then defined as the time-average of the instantaneous value. This effective Stokes number thus encapsulates the flow history and geometric variability. Our results demonstrate that the effective Stokes number can deviate significantly from the reference value based solely on a characteristic flow velocity and length scale. In addition, the effective Stokes number shows a clear correlation with deposition efficiency.

Fri, 19 Jun 2015
11:30
L5

iceCAM project with G's-Fresh

Alasdair Craighead
(G's-Fresh)
Abstract

G’s Growers supply salad and vegetable crops throughout the UK and Europe; primarily as a direct supplier to supermarkets. We are currently working on a project to improve the availability of Iceberg Lettuce throughout the year as this has historically been a very volatile crop. It is also by far the highest volume crop that we produce with typical weekly sales in the summer season being about 3m heads per week.

In order to continue to grow our business we must maintain continuous supply to the supermarkets. Our current method for achieving this is to grow more crop than we will actually harvest. We then aim to use the wholesale markets to sell the extra crop that is grown rather than ploughing it back in and then we reduce availability to these markets when the availability is tight.

We currently use a relatively simple computer Heat Unit model to help predict availability however we know that this is not the full picture. In order to try to help improve our position we have started the IceCAM project (Iceberg Crop Adaptive Model) which has 3 aims.

  1. Forecast crop availability spikes and troughs and use this to have better planting programmes from the start of the season.
  2. Identify the growth stages of Iceberg to measure more accurately whether crop is ahead or behind expectation when it is physically examined in the field.
  3. The final utopian aim would be to match the market so that in times of general shortage when price are high we have sufficient crop to meet all of our supermarket customer requirements and still have spare to sell onto the markets to benefit from the higher prices. Equally when there is a general surplus we would only look to have sufficient to supply the primary customer base.

We believe that statistical mathematics can help us to solve these problems!!

Tue, 23 Jun 2015

17:00 - 18:00
L6

Almost small absolute Galois groups

Arno Fehm
(Konstanz)
Abstract

Already Serre's "Cohomologie Galoisienne" contains an exercise regarding the following condition on a field F: For every finite field extension E of F and every n, the index of the n-th powers (E*)^n in the multiplicative group E* is finite. Model theorists recently got interested in this condition, as it is satisfied by every superrosy field and also by every strongly2 dependent field, and occurs in a conjecture of Shelah-Hasson on NIP fields. I will explain how it relates to the better known condition that F is bounded (i.e. F has only finitely many extensions of degree n, for any n - in other words, the absolute Galois group of F is a small profinite group) and why it is not preserved under elementary equivalence. Joint work with Franziska Jahnke.

*** Note unusual day and time ***

Subscribe to