Past Differential Equations and Applications Seminar

In this talk I shall discuss special polynomials associated with rational solutions of the Painlevé equations and of the soliton equations which are solvable by the inverse scattering method, including the Korteweg-de Vries, Boussinesq and nonlinear Schrodinger equations. Further I shall illustrate applications of these polynomials to vortex dynamics and rogue waves.

The Painlevé equations are six nonlinear ordinary differential equations that have been the subject of much interest in the past thirty years, and have arisen in a variety of physical applications. Further the Painlevé equations may be thought of as nonlinear special functions. Rational solutions of the Painlevé equations are expressible in terms of the logarithmic derivative of certain special polynomials. For the fourth Painlevé equation these polynomials are known as the generalized Hermite polynomials and generalized Okamoto polynomials. The locations of the roots of these polynomials have a highly symmetric (and intriguing) structure in the complex plane.

It is well known that soliton equations have symmetry reductions which reduce them to the Painlevé equations, e.g. scaling reductions of the Boussinesq and nonlinear Schrödinger equations are expressible in terms of the fourth Painlevé equation. Hence rational solutions of these equations can be expressed in terms of the generalized Hermite and generalized Okamoto polynomials.

I will also discuss the relationship between vortex dynamics and properties of polynomials with roots at the vortex positions. Classical polynomials such as the Hermite and Laguerre polynomials have roots which describe vortex equilibria. Stationary vortex configurations with vortices of the same strength and positive or negative configurations are located at the roots of the Adler-Moser polynomials, which are associated with rational solutions of the Kortweg-de Vries equation.

Further, I shall also describe some additional rational solutions of the Boussinesq equation and rational-oscillatory solutions of the focusing nonlinear Schrödinger equation which have applications to rogue waves.

Three-wave interactions form the basis of our understanding of many

nonlinear pattern forming systems because they encapsulate the most basic

nonlinear interactions. In problems with two comparable length scales, such

as the Faraday wave experiment with multi-frequency forcing, consideration

of three-wave interactions can explain the presence of the spatio-temporal

chaos found in some experiments, enabling some previously unexplained

results to be interpreted in a new light. The predictions are illustrated

with numerical simulations of a model partial differential equation.

Isostatic mounts are used in applications like telescopes and robotics to move and hold part of a structure in a desired pose relative to the rest, by driving some controls rather than driving the subsystem directly. To achieve this successfully requires an understanding of the coupled space of configurations and controls, and of the singularities of the mapping from the coupled space to the space of controls. It is crucial to avoid such singularities because generically they lead to large constraint forces and internal stresses which can cause distortion. In this paper we outline design principles for isostatic mount systems for dynamic structures, with particular emphasis on robots.

Systems with delays frequently appear in engineering. The presence of delays makes system analysis and control design very complicated. In this talk, the standard H-infinity control problem of time-delay systems will be discussed. The emphasis will be on systems having an input or output delay. The problem is solved in the frequency domain via reduction to a one-block problem and then further to an extended Nehari problem using a simple and intuitive method. After solving the extended Nehari problem, the original problem is solved. The solvability of the extended Nehari problem (or the one-block problem) is equivalent to the nonsingularity of a delay-dependent matrix and the solvability conditions of the standard H-infinity control problem with a delay are then formulated in terms of the existence of solutions to two delay-independent algebraic Riccati equations and a delay-dependent nonsingular matrix.

15th Biennial OXFORD / CAMBRIDGE MEETING

PROGRAMME FOR THE

‘WOOLLY OWL TROPHY’

Invited Judges

John Harper

(Victoria University of Wellington, NZ)

Arash Yavari

(Georgia Tech, Atlanta, USA)

Sharon Stephen

(University of Birmingham, UK)

10:45 Morning Coffee The Maths Inst Common Room

Solving partial differential equations (PDEs) on curved surfaces is

important in many areas of science. The Closest Point Method is a new

technique for computing numerical solutions to PDEs on curves,

surfaces, and more general domains. For example, it can be used to

solve a pattern-formation PDE on the surface of a rabbit.

A benefit of the Closest Point Method is its simplicity: it is easy to

understand and straightforward to implement on a wide variety of PDEs

and surfaces. In this presentation, I will introduce the Closest

Point Method and highlight some of the research in this area. Example

computations (including the in-surface heat equation,

reaction-diffusion on surfaces, level set equations, high-order

interface motion, and Laplace--Beltrami eigenmodes) on a variety of

surfaces will demonstrate the effectiveness of the method.

Voltammetry is a powerful method for interrogating electrochemical systems. A voltage is applied to an electrode and the resulting current response analysed to determine features of the system under investigation, such as concentrations, diffusion coefficients, rate constants and thermodynamic potentials. Here we will focus on ac voltammetry, where the voltage signal consists of a high frequency sine-wave superimposed on a linear ramp. Using multiple scales analysis, we find analytical solutions for the harmonics of the current response and show how they can be used to determine the system parameters. We also include the effects of capacitance due to the double-layer at the electrode surface and show that even in the presence of large capacitance, the harmonics of the current response can still be isolated using the FFT and the Hanning window.

Flows involving immiscible liquids are encountered in a variety of industrial and natural processes. Recent applications in micro- and nano-fluidics have led to a significant scientific effort whose aim (among other aspects) is to enable theoretical predictions of the spatiotemporal dynamics of the interface(s) separating different flowing liquids. In such applications the scale of the system is small, and forces such as surface tension or externally imposed electrostatic forces compete and can, in many cases, surpass those of gravity and inertia. This talk will begin with a brief survey of applications where electrohydrodynamics have been used experimentally in micro-lithography, and experiments will be presented that demonstrate the use of electric fields in producing controlled encapsulated droplet formation in microchannels.

The main thrust of the talk will be theoretical and will mostly focus on the paradigm problem of the dynamics of electrified falling liquid films over topographically structured substrates.

Evolution equations will be developed asymptotically and their solutions will be compared to direct simulations in order to identify their practicality. The equations are rich mathematically and yield novel examples of dissipative evolutionary systems with additional effects (typically these are pseudo-differential operators) due to dispersion and external fields.

The models will be analysed (we have rigorous results concerning global existence of solutions, the existence of dissipative dynamics and an absorbing set, and analyticity), and accurate numerical solutions will be presented to describe the large time dynamics. It is found that electric fields and topography can be used to control the flow.Time permitting, I will present some recent results on transitions between convective to absolute instabilities for film flows over periodic topography.

The notion of critical mass in research is one that has been around for a long time without proper definition. It has been described as some kind of threshold group size above which research standards significantly improve. However no evidence for such a threshold has been found and critical mass has never been measured -- until now.

We present a new, simple, sociophysical model which explains how research quality depends on research-group structure and in particular on size. Our model predicts that there are, in fact, two critical masses in research, the values of which are discipline dependent. Research quality tends to be linearly dependent on group size, but only up to a limit termed the 'upper critical mass'. The upper critical mass is interpreted as the average maximum number of colleagues with whom a given individual in a research group can meaningfully interact. Once the group exceeds this size, it tends to fragment into sub-groups and research quality no longer improves significantly with increasing size. There is also a

lower critical mass, which small research groups should strive to achieve for stability.

Our theory is tested using empirical data from RAE 2008 on the quantity and quality of research groups, for which critical masses are determined. For pure and applied mathematics, the lower critical mass is about 2 and 6, respectively, while for statistics and physics it is 9 and 13. The upper critical mass, beyond which research quality does not significantly improve with increasing group size, is about twice the lower value.

We develop a theory of impact of viscoelastic spheres with adhesive

interactions. We assume that the collision velocities are not large to

avoid the fracture and plastic deformation of particles material and

microscopic relaxation time is much smaller than the collision duration.

The adhesive interactions are described with the use of Johnson, Kendall

and Roberts (JKR) theory, while dissipation is attributed to the

viscoelastic behavior of the material. For small impact velocities we

apply the condition of a quasi-static collision and obtain the

inter-particle force. We show that this force is a sum of four

components, having in addition to common elastic, viscous and adhesive

force, the visco-adhesive cross term. Using the derived force we compute

the coefficient of normal restitution and consider the application of our

theory to the collisions of macro and nano-particles.

Human behaviour can show surprising properties when looked at from a collective point of view. Data on collective behaviour can be gleaned from a number of sources, and mobile phone data are increasingly becoming used. A major challenge is combining behavioural data with health data. In this talk I will describe our approach to understanding behaviour change related to change in health status at a collective level.

A mathematical model of Olufsen [1,2] has been extended to study periodic pulse propagation in both the systemic arteries and the pulmonary arterial and venous trees. The systemic and pulmonary circulations are treated as separate, bifurcating trees of compliant and tapering vessels. Each model is divided into two coupled parts: the larger and smaller vessels. Blood flow and pressure in the larger arteries and veins are predicted from a nonlinear 1D cross-sectional area-averaged model for a Newtonian fluid in an elastic tube. The initial cardiac output is obtained from magnetic resonance measurements.

The smaller blood vessels are modelled as asymmetric structured trees with specified area and asymmetry ratios between the parent and daughter arteries. For the systemic arteries, the smaller vessels are placed into a number of separate trees representing different vascular beds corresponding to major organs and limbs. Womersley's theory gives the wave equation in the frequency domain for the 1D flow in these smaller vessels, resulting in a linear system. The impedances of the smallest vessels are set to a constant and then back-calculation gives the required outflow boundary condition for the Navier--Stokes equations in the larger vessels. The flow and pressure in the large vessels are then used to calculate the flow and pressure in the small vessels. This gives the first theoretical calculations of the pressure pulse in the small `resistance' arteries which control the haemodynamic pressure drop.

I will discuss the effects, on both the forward-propagating and the reflected components of the pressure pulse waveform, of the number of generations of blood vessels, the compliance of the arterial wall, and of vascular rarefaction (the loss of small systemic arterioles) which is associated with type II diabetes. We discuss the possibilities for developing clinical indicators for the early detection of vascular disease.

References:

1. M.S. Olufsen et al., Ann Biomed Eng. 28, 1281-99 (2000)

2. M.S. Olufsen, Am J Physiol. 276, H257--68 (1999)

The mechanics of thin elastic or viscous objects has applications in e.g. the buckling of engineering structures, the spinning of polymer fibers, or the crumpling of plates and shells. During the past decade the mathematics, mechanics and physics communities have witnessed an upsurge of interest in those issues. A general question is to how patterns are formed in thin structures. In this talk I consider two illustrative problems: the shapes of an elastic knot, and the stitching patterns laid down by a viscous thread falling on a moving belt. These intriguing phenomena can be understood by using a combination of approaches, ranging from numerical to analytical, and based on exact equations or low-dimensional models.

"Most drivers will recognize the scenario: you are making steady progress along the motorway when suddenly you come to a sudden halt at the tail end of a lengthy queue of traffic. When you move off again you look for the cause of the jam, but there isn't one. No accident damaged cars, no breakdown, no dead animal, and no debris strewn on the road. So what caused everyone to stop?" RAC news release (2005)

The (by now well-known) answer is that such "phantom traffic jams" exist as waves that propagate upstream (opposite to the driving direction) - so that the vast majority of individuals do not observe the instant at which the jam was created - yet what exactly goes on at that instant is still a matter of debate. In this talk I'll give an overview of empirical data and models to describe such spatiotemporal patterns. The key property we need is instability: and using the framework of car-following (CF) models, I'll show how different sorts of linear (convective and absolute) and nonlinear instability can be used to explain empirical patterns.

Rubbers and biological soft tissues undergo large isochoric motions in service, and can thus be modelled as nonlinear, incompressible elastic solids. It is easy to enforce incompressibility in the finite (exact) theory of nonlinear elasticity, but not so simple in the weakly nonlinear formulation, where the stress is expanded in successive powers of the strain. In linear and second-order elasticity, incompressibility means that Poisson's ratio is 1/2. Here we show how third- and fourth-order elastic constants behave in the incompressible limit. For applications, we turn to the propagation of elastic waves in soft incompressible solids, a topic of crucial importance in medical imaging (joint work with Ray Ogden, University of Aberdeen).

Predicting the dynamics of foams requires input from geometry and both Newtonian and non-Newtonian fluid mechanics, among many other fields. I will attempt to give a flavour of this richness by discussing the static structure of a foam and how it allows the derivation of dynamic properties, at least to leading order. The latter includes coarsening due to gas diffusion, liquid drainage under gravity, and the flow of the bubbles themselves.

Several stochastic simulation algorithms (SSAs) have been recently proposed for modelling reaction-diffusion processes in cellular and molecular biology. In this talk, two commonly used SSAs will be studied. The first SSA is an on-lattice model described by the reaction-diffusion master equation. The second SSA is an off-lattice model based on the simulation of Brownian motion of individual molecules and their reactive collisions. The connections between SSAs and the deterministic models (based on reaction-diffusion PDEs) will be presented. I will consider chemical reactions both at a surface and in the bulk. I will show how the "microscopic" parameters should be chosen to achieve the correct "macroscopic" reaction rate. This choice is found to depend on which SSA is used. I will also present multiscale algorithms which use models with a different level of detail in different parts of the computational domain.

This talk will focus on a family of Bayesian inference algorithms built around Gaussian processes. We firstly introduce an iterative Gaussian process for multi-sensor inference problems. Extensions to our algorithm allow us to tackle some of the decision problems faced in sensor networks, including observation scheduling. Along these lines, we also propose a general method of global optimisation, Gaussian process global optimisation (GPGO). This paradigm is extended to the Bayesian decision problem of sequential multi-scale observation selection. We show how the hyperparameters of our system can be marginalised by use of Bayesian quadrature and frame the selection of the positions of the hyperparameter samples required by Bayesian quadrature as a sequential decision problem, with the aim of minimising the uncertainty we possess about the values of the integrals we are approximating.

I will discuss a few problems  that involve randomness , chosen randomly  (?) from the following : (i) the probability of a coin landing on a side  (ii) optimal strategies for throwing accurately, (iii)  the statistical mechanics of a ribbon, (iv) the intermittent dynamics of a growing polymeric assembly (v) fat tails from feedback.



When modeling biochemical reactions within cells, it is vitally important to take into account the effect of intrinsic noise in the system, due to the small copy numbers of some of the chemical species. Deterministic systems can give vastly different types of behaviour for the same parameter sets of reaction rates as their stochastic analogues, giving us an incorrect view of the bifurcation diagram.

Stochastic Simulation Algorithms (SSAs) exist which draw exact trajectories from the Chemical Master Equation (CME). However, these methods can be very computationally expensive, particularly where there is a separation of time scales of the evolution of some of the chemical species. Some of the species may react many times on a time scale for which others are highly unlikely to react at all. Simulating all of these reactions of the fast species is a waste of computational effort, and many different methods exist for reducing the system to one which only contains the slow variables.

In this talk we will introduce the conditional Gillespie algorithm, a method for sampling directly from the conditional distribution on the fast variables, given a static value for the slow variables. Using this, we will go on to describe the constrained Gillespie approach, which uses simulations of the CG algorithm to estimate the drift and diffusion terms of the effective dynamics of the slow variables.

If there is time at the end, I will briefly describe my work on another project, which involves full sampling of the posterior distributions in various problems in data assimilation using Monte Carlo Markov Chain (MCMC) methods.

In my talk I will introduce the concept of spectral discrete solitons

(SDSs): solutions of nonlinear Schroedinger type equations, which are localized on a regular grid in frequency space. In time domain such solitons correspond to periodic trains of pulses. SDSs play important role in cascaded four-wave-mixing processes (frequency comb generation) in optical fibres, where initial excitation by a two-frequency pump leads to the generation of multiple side-bands. When free space diffraction is taken into consideration, a non-trivial generalization of 1D SDSs will be discussed, in which every individual harmonic is an optical vortex with its own topological charge. Such excitations correspond to spatio-temporal helical beams.

A kinetic theory for swarming systems of interacting individuals will be described with and without noise. Starting from the the particle model \cite{DCBC}, one can construct solutions to a kinetic equation for the single particle probability distribution function using distances between measures \cite{dobru}. Analogously, we will discuss the mean-field limit for these problems with noise.

We will also present and analys the asymptotic behavior of solutions of the continuous kinetic version of flocking by Cucker and Smale The large-time behavior of the distribution in phase space is subsequently studied by means of particle approximations and a stability property in distances between measures. It will be shown that the solutions concentrate exponentially fast their velocity to their mean while in space they will converge towards a translational flocking solution.

The Mathematical Institute invites you to attend the Inaugural Lecture of Professor Gui-Qiang G. Chen. Professor in the Analysis of Partial Differential Equations. Examination Schools, 75-81 High Street, Oxford, OX 4BG.

There is no charge to attend but registration is required. Please register your attendance by sending an email to @email specifying the number of people in your party. Admission will only be allowed with prior registration.

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

ABSTRACT

While calculus is a mathematical theory concerned with change, differential equations are the mathematician's foremost aid for describing change. In the simplest case, a process depends on one variable alone, for example time. More complex phenomena depend on several variables – perhaps time and, in addition, one, two or three space variables. Such processes require the use of partial differential equations. The behaviour of every material object in nature, with timescales ranging from picoseconds to millennia and length scales ranging from sub-atomic to astronomical, can be modelled by nonlinear partial differential equations or by equations with similar features. The roles of partial differential equations within mathematics and in the other sciences become increasingly significant. The mathematical theory of partial differential equations has a long history. In the recent decades, the subject has experienced a vigorous growth, and research is marching on at a brisk pace.

In this lecture, Professor Gui-Qiang G. Chen will present several examples to illustrate the origins, developments, and roles of partial differential equations in our changing world.

The presence of additives, which may or may not be surface-active, can have a dramatic influence on interfacial flows. The presence of surfactants alters the interfacial tension and drives Marangoni flow that leads to fingering instabilities in drops spreading on ultra-thin films. Surfactants also play a major role in coating flows, foam drainage, jet breakup and may be responsible for the so-called ``super-spreading" of drops on hydrophobic substrates. The addition of surface-inactive nano-particles to thin films and drops also influences the interfacial dynamics and has recently been shown to accelerate spreading and to modify the boiling characteristics of nanofluids. These findings have been attributed to the structural component of the disjoining pressure resulting from the ordered layering of nanoparticles in the region near the contact line. In this talk, we present a collection of results which demonstrate that the above-mentioned effects of surfactants and nano-particles can be captured using long-wave models.

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.