Past Differential Equations and Applications Seminar

Many growing plant cells undergo rapid axial elongation with negligible radial expansion. Growth is driven by high internal turgor pressure causing viscous stretching of the cell wall, with embedded cellulose microfibrils providing the wall with strongly anisotropic properties. We present a theoretical model of a growing cell, representing the primary cell wall as a thin axisymmetric fibre-reinforced viscous sheet supported between rigid end plates. Asymptotic reduction of the governing equations, under simple sets of assumptions about the fibre and wall properties, yields variants of the traditional Lockhart equation, which relates the axial cell growth rate to the internal pressure. The model provides insights into the geometric and biomechanical parameters underlying bulk quantities such as wall extensibility and shows how either dynamical changes in wall material properties or passive fibre reorientation may suppress cell elongation. We then investigate how the action of enzymes on the cell wall microstructure can lead to the required dynamic changes in macroscale wall material properties, and thus demonstrate a mechanism by which hormones may regulate plant growth.

Whenever we say the words "fluid flows" or "shape changes" we enter the realm of infinite-dimensional geometric mechanics. Water, for example, flows. In fact, Euler's equations tell us that water flows a particular way. Namely, it flows to get out of its own way as adroitly as possible. The shape of water changes by smooth invertible maps called diffeos (short for diffeomorphisms). The flow responsible for this optimal change of shape follows the path of shortest length, the geodesic, defined by the metric of kinetic energy. Not just the flow of water, but the optimal morphing of any shape into another follows one of these optimal paths.

The lecture will be about the commonalities between fluid dynamics and shape changes and will be discussed in the language most suited to fundamental understanding -- the language of geometric mechanics. A common theme will be the use of momentum maps and geometric control for steering along the optimal paths using emergent singular solutions of the initial value problem for a nonlinear partial differential equation called EPDiff, that governs metamorphosis along the geodesic flow of the diffeos. The main application will be in the registration and comparison of Magnetic Resonance Images for clinical diagnosis and medical procedures.

The new model is that the universal small scale structure of high Reynolds number turbulence is determined by the dynamics of thin evolving shear layers, with thickness of the order of the Taylor micro scale,within which there are the familiar elongated vortices .Local quasi-linear dynamics shows how the shear layers act as barriers to external eddies and a filter for the transfer of energy to their interiors. The model is consistent with direct numerical simulations by Ishihara and Kaneda analysed in terms of conditional statistics relative to the layers and also with recent 4D measurements of lab turbulence by Wirth and Nickels. The model explains how the transport of energy into the layers leads to the observed inertial range spectrum and to the generation of intense structures, on the scale of the Kolmogorov micro-scale.

But the modelling also explains the important discrepancies between data and the Kolmogorov-Richardson cascade concept ,eg larger amplitudes of the smallest scale motions and of the higher moments ,and why the latter are generally less isotropic than lower order moments, eg in thermal convection. Ref JCRHunt , I Eames, P Davidson,J.Westerweel, J Fernando, S Voropayev, M Braza J Hyd Env Res 2010

Wound coils or rolls accumulate essentially flat strip compactly without folding or cutting and typically, strip is wound and unwound a number of times before its end use. The variety of material that is wound into coils or rolls is very extensive and includes magnetic tape, paper, cellophane, plastics, fabric and metals such as aluminium and steel.

Stresses wound into a coil provide its structural integrity via the frictional forces between the wraps. For a coil with inadequate inter-wrap pressure, the wraps may slip or telescope (causing surface scuffing) or the coil may slump and collapse. On the other hand, large internal stresses can cause increased creep and stress relaxation, collapse at the bore, stress wrinkling and rupture of the material in the coil.

Given the range of applications, it is not surprising that the literature on calculating stresses in wound coils is large and has a long history, which goes back at least to the wire winding of gun barrels. However the basic approach of the resulting accretion models, where the residual stress is recalculated each time a layer is added, has remained essentially the same. In this talk, we take a radically different approach in analysing the winding stresses in coils. Instead of the traditional method, we seek to deduce a winding policy that will achieve a target distribution of residual stresses within a coil. In this way, optimising the coiling tension profile is much more straight-forward, by

* Specifying the residue stresses required to avoid operational problems, tight-bore collapses, and other issues such as scuffing, then

* Determining the winding tension profile to produce the required residue stresses.

The question in the title seems to be neglected in the studies of open dynamical systems. It occurred though that the features of dynamics may play a role comparable to the one played by the size of a hole. For instance, the escape through the smaller hole could be faster than through the larger one.

These studies revealed as well a new role of the periodic orbits in the dynamics which could be exactly quantified in some cases. Moreover, this new approach allows to characterize the elements of networks by their dynamical properties (rather than by static ones like centrality, betweenness, etc.)

Bloch Floquet waves are considered in structured media. Such waves are dispersive and the dispersion diagrams contain stop bands. For an example of a harmonic lattice, we discuss dynamic band gap Green’s functions characterised by exponential localisation. This is followed by simple models of exponentially localised defect modes. Asymptotic models involving uniform asymptotic approximations of physical fields in structured media are compared with homogenisation approximations.

We report the numerical realization and demonstration of robustness of certain 2-component structures in Bose-Einstein Condensates in 2 and 3 spatial dimensions with non-trivial topological charge in one of the components. In particular, we identify a stable symbiotic state in which a higher-dimensional bright soliton exists even in a homogeneous setting with defocusing interactions, as a result of the effective potential created by a stable vortex in the other component. The resulting vortex-bright solitary waves, which naturally generalize the recently experimentally observed dark-bright solitons, are examined both in the homogeneous medium and in the presence of parabolic and periodic external confinement and are found to be very robust.

An emerging set of methods enables an experimental dialogue with biological systems composed of many interacting cell types---in particular, with neural circuits in the brain. These methods are sometimes called “optogenetic” because they employ light-responsive proteins (“opto-“) encoded in DNA (“-genetic”). Optogenetic devices can be introduced into tissues or whole organisms by genetic manipulation and be expressed in anatomically or functionally defined groups of cells. Two kinds of devices perform complementary functions: light-driven actuators control electrochemical signals; light-emitting sensors report them. Actuators pose questions by delivering targeted perturbations; sensors (and other measurements) signal answers. These catechisms are beginning to yield previously unattainable insight into the organization of neural circuits, the regulation of their collective dynamics, and the causal relationships between cellular activity patterns and behavior.

In the first part of my presentation, I plan to review several applications modelled by delay differential equations (DDEs) starting from familiar examples such as traffic flow problems to physiology and industrial problems. Although delay differential equations have the reputation to be difficult mathematical problems, there is a renewed interest for both old and new problems modelled by DDEs. In the second part of my talk, I’ll emphasize the need of developing asymptotic tools for DDEs in order to guide our numerical simulations and help our physical understanding. I illustrate these ideas by considering the response of optical optoelectronic oscillators that have been studied both experimentally and numerically.

How does form emerge from cellular processes? Using cell-based mechanical models of growth, we investigated the geometry of leaf vasculature and the cellular arrangements at the shoot apex. These models incorporate turgor pressure, wall mechanical properties and cell division. In connection with experimental data, they allowed us to, on the one hand, account for characteristic geometrical property of vein junctions, and, on the other hand, speculate that growth is locally regulated.

We study dynamics from models for the human tear film in one and two dimensional domains.

The tear film is roughly a few microns thick over a domain on a centimeter scale; this separation of scales makes lubrication models desirable. Results on one-dimensional blinking domains are presented for multiple blink cycles. Results on two-dimensional stationary domains are presented for different boundary conditions. In all cases, the results are sensitive to the boundary conditions; this is intuitively satisfying since the tear film seems to be controlled primarily from the boundary and its motion. Quantitative comparison with in vivo measurement will be given in some cases. Some discussion of tear film properties will also be given, and results for non-Newtonian models will be given as available, as well as some future directions.

If an ideal elastic spring is greatly stretched, it will develop large stresses. However, solid biological tissues are able to grow without developing such large stresses. This is because the cells within such tissues are able to lay down new fibres and remove old ones, fundamentally changing the mechanical structure of the tissue. In many ways, this is analogous to classical plasticity, where materials stretched beyond their yield point begin to flow and the unloaded state of the material changes. Unfortunately, biological tissues are not closed systems and so we are not able to use standard plasticity techniques where we require the flow to be mass conserving and energetically passive.

In this talk, a general framework will be presented for modelling the changing zero stress state of a biological tissue (or any other material). Working from the multiplicative decomposition of the deformation gradient, we show that the rate of 'desired' growth can represented using a tensor that describes both the total rate of growth and any directional biases. This can be used to give an evolution equation for the effective strain (a measure of the difference between the current state and the zero stress state). We conclude by looking at a perhaps surprising application for this theory as a method for deriving the constitutive laws of a viscoelastic fluid.

We study the axisymmetric stretching of a thin sheet of viscous fluid

driven by a centrifugal body force. Time-dependent simulations show that

the sheet radius tends to infinity in finite time. As the critical time is

approached, the sheet becomes partitioned into a very thin central region

and a relatively thick rim. A net momentum and mass balance in the rim leads

to a prediction for the sheet radius near the singularity that agrees with the numerical

simulations. By asymptotically matching the dynamics of the sheet with the

rim, we find that the thickness in the central region is described by a

similarity solution of the second kind. For non-zero surface tension, we

find that the similarity exponent depends on the rotational Bond number B,

and increases to infinity at a critical value B=1/4. For B>1/4, surface

tension defeats the centrifugal force, causing the sheet to retract rather

than stretch, with the limiting behaviour described by a similarity

solution of the first kind.

The FPU lattice is a coupled system of ordinary differential equations in which each atom in a chain is coupled to its nearest neighbour by way of a nonlinear spring.

After summarising the properties of travelling waves (kinks) we use asymptotic analysis to describe more complicate envelope solutions (breathers). The interaction of breathers and kinks will then be analysed. If time permits, the method will be extended to two-dimensional lattices.

Some years ago Hall and Smith in a number of papers developed a theory governing the interaction of vortices and waves in shear flows. In recent years immense interest has been focused on so-called self-sustained processes in turbulent shear flows where the importance of waves interacting with streamwise vortex flows has been elucidated in a number of; see for example the work of Waleffe and colleagues, Kerswell, Gibson, etc. These processes have a striking resemblance to coherent structures observed in turbulent shear flow and for that reason they are often referred to as exact coherent structures. It is shown that the structures associated with the so-called 'lower branch' state, which has been shown to play a crucial role in these self-sustained process, is nothing but a Rayleigh wave vortex interaction with a wave system generating streamwise vortices inside a critical layer. The theory enables the reduction of the 3D Navier Stokes equations to a coupled system for a steady streamwise vortex and an inviscid wave system. The reduced system for the streamwise vortices must be solved with jump conditions in the shear across the critical layer and the position of that layer constitutes a nonlinear pde eigenvalue problem. Remarkable agreement between the asymptotic theory and numerical simulations is found thereby demonstrating the importance of vortex-wave interaction theory in the mathematical description of coherent structures in turbulent shear flows. The theory offers the possibility of drag reduction in turbulent shear flows by controlling the flow to the neighborhood of the lower branch state. The relevance of the work to more general shear flows is also discussed.

Patterns of sources and sinks in the complex Ginzburg-Landau equation Jonathan Sherratt, Heriot-Watt University The complex Ginzburg-Landau equation is a prototype model for self-oscillatory systems such as binary fluid convection, chemical oscillators, and cyclic predator-prey systems. In one space dimension, many boundary conditions that arise naturally in applications generate wavetrain solutions. In some contexts, the wavetrain is unstable as a solution of the original equation, and it proves necessary to distinguish between two different types of instability, which I will

explain: convective and absolute. When the wavetrain is absolutely unstable, the selected wavetrain breaks up into spatiotemporal chaos. But when it is only convectively stable, there is a different behaviour, with bands of wavetrains separated by sharp interfaces known as "sources" and "sinks". These have been studied in great detail as isolated objects, but there has been very little work on patterns of alternating sources and sinks, which is what one typically sees in simulations. I will discuss new results on source-sink patterns, which show that the separation distances between sources and sinks are constrained to a discrete set of possible values, because of a phase-locking condition.

I will present results from numerical simulations that confirm the results, and I will briefly discuss applications and the future challenges. The work that I will describe has been done in collaboration with Matthew Smith (Microsoft Research) and Jens Rademacher (CWI, Amsterdam).

------------------------------

Evolutionary PDEs on stationary and moving surfaces appear in many applications such as the diffusion of surfactants on fluid interfaces, surface pattern formation on growing domains, segmentation on curved surfaces and phase separation on biomembranes and dissolving alloy surfaces.

In this talk I discuss three numerical approaches based on:- (I) Surface Finite Elements and Triangulated Surfaces, (II)Level Set Method and Implicit Surface PDEs and (III) Phase Field Approaches and Diffuse Surfaces.

Modelling phase change in the presence of a flowing thin liquid film

There are numerous physical phenomena that involve a melting solid

surrounded by a thin layer of liquid, or alternatively a solid

forming from a thin liquid layer. This talk will involve two such

problems, namely contact melting and the Leidenfrost phenomenon.

Contact melting occurs, for example, when a solid is placed on a

surface that is maintained at a temperature above the solid melting

temperature. Consequently the solid melts, while the melt layer is

squeezed out from under the solid due to its weight. This process

has applications in metallurgy, geology and nuclear technology, and

also describes a piece of ice melting on a table. Leidenfrost is

similar, but involves a liquid droplet evaporating after being

placed on a hot substrate. This has applications in cooling systems

and combustion of fuel or a drop of water on a hot frying pan.

The talk will begin with a brief introduction into one-dimensional

Stefan problems before moving on to the problem of melting coupled

to flow. Mathematical models will be developed, analysed and

compared with experimental results. Along the way the Heat Balance

Integral Method (HBIM) will be introduced. This is a well-known

method primarily used by engineers to approximate the solution of

thermal problems. However, it has not proved so popular with

mathematicians, due to the arbitrary choice of approximating

function and a lack of accuracy. The method will be demonstrated on

a simple example, then it will be shown how it may be modified to

significantly improve the accuracy. In fact, in the large Stefan

number limit the modified method can be shown to be more accurate

than the asymptotic solution to second order.

Whispering gallery modes in optical resonators have received a lot of attention as a mechanism for constructing small, directional lasers. They are also potentially important as passive optical components in schemes for coupling and filtering signals in optical fibres, in sensing devices and in other applications. In this talk it is argued that the evanescent field outside resonators that are very slightly deformed from circular or spherical is surprising in a couple of respects. First, even very small deformations seem to be capable of leading to highly directional emission patterns. Second, even though the undelying ray families are very regular and hardly differ from the integrable circular or spherical limit inside the resonator, a calculation of the evanescent field outside it is not straightforward.

This is because even very slight nonintegrability has a profound effect on the complexified ray families which guide the external wave to asymptopia. An approach to describing the emitted wave is described which is based on canonical perturbation theory applied to the ray families and extended to comeplx phase space.