Applied and Computational Mathematics Research Group

The major research interests of the group are in theoretical and computational fluid mechanics. However, members of the group also undertake research in numerical analysis, structural and solid mechanics, inverse problems and applied analysis.

Our research in the area of applied and computational mathematics is informed by problems at the interface with physical sciences, biological sciences and engineering and there are many productive interdisciplinary collaborations within Cardiff University and further afield.

Our national and international collaborators include research groups at Imperial, Oxford, Cambridge, Warwick, Swansea, Canberra, Curtin (Perth), Perth, Delft, Northwestern, Ljubljana and Montréal.

The group hosts the Cardiff University Student Chapter of the Society for Industrial and Applied Mathematics and the Institute of Mathematics and its Applications (SIAM-IMA Student Chapter) which encompasses postgraduate students and faculty members from across the University who are interested in mathematics or scientific computing and their real-world applications.

The main areas of research within the current group are:

Fluid mechanics

Viscous buckling phenomena; wetting phenomena; free-surface flows; dynamics of liquid films and jets; hydrodynamic stability theory; laminar-turbulent transition mechanisms; boundary-layer and wake flow instabilities; boundary layer flow control; viscoelastic flows and computational rheology; bubble dynamics; kinetic theory models.

Biological fluid mechanics

Adhesion and dessication of biological cells; fluid-structure interaction problems.

Structural and solid mechanics

Numerical analysis of damage; multi-scale modelling, limit states analysis.

Numerical analysis and scientific computing

Finite element methods for elasticity; spectral element methods; immersed boundary method; lattice Boltzmann methods; proper generalized decomposition; dynamic density functional theory; numerical continuation; differential algebraic equation; automatic differentiation.

Applied analysis

Inverse problems in materials modelling; homogenisation and the mechanics of composites.

Applied mathematical modelling

Singular perturbation methods; matched asymptotics; dynamical systems and bifurcation analysis; stochastic processes; physical resolution of singularities in mathematical models.

Head of Group

Tim Phillips photograph

Professor Tim Phillips

Head of School, Mathematics

+44 (0)29 2087 4194

Academic staff

Photograph of Mikhail Cherdantsev

Dr Mikhail Cherdantsev


+44 (0)29 2087 5549
Photograph of Russell Davies

Professor Russell Davies

Honorary Distinguished Professor

+44 (0)29 2087 4827
Photograph of Chris Davies

Dr Chris Davies


+44 (0)29 2087 5531
Dr Usama Kadri photograph

Dr Usama Kadri

Lecturer in Applied Mathematics

029 208 75863
Photograph of Dr Angela Mihai

Dr Angela Mihai

Senior Lecturer in Applied Mathematics

+44 (0)29 2087 5570
Prof John Pryce photograph

Professor John Pryce

Associate Tutor

+44 (0)29 2087 4207
Photograph of Nikos Savva

Dr Nikos Savva

Lecturer in Applied Mathematics

+44 (0)29 2087 5116
Thomas Woolley Profile Picture

Dr Thomas Woolley

Lecturer in Applied Mathematics

02920 870618


All seminars are held at 15:10 in Room M/2.06, Senghennydd Road, Cardiff unless stated otherwise.

Programme organiser and contact: Dr Usama Kadri


4 October 2016

Edd Lewis

Xander Ramage

Waleed Ali

Lattice Boltzmann methods for FENE type fluids

Onset of absolute instability in Stokes layers

Onset of global instabilities in the plane channel flow between compliant walls

11 October 2016

Prof. Victor Shrira (Keele University)

Inertial waves and deep ocean mixing

For the existing pattern of global oceanic circulation to exist there should be sufficiently strong turbulent mixing in the abyssal ocean. It is commonly believed that it is breaking of inertia-gravity internal waves which provides the required mixing. However this belief is not supported by understanding of why internal waves should break so intensely in the abyssal ocean. The specific physical mechanisms causing the breaking have not been identified and investigated.

The talk discusses a very plausible mechanism leading to intense breaking of near inertial waves near the bottom of the ocean. The simultaneous account of both the horizontal component of the Earth rotation and its latitude dependence (the beta-effect) reveals the existence for near inertial waves of wide waveguides attached to the bottom. These waveguides are narrowing in the poleward direction. Near inertial waves propagating poleward get trapped in these waveguides; in the process the waves are focussing more and more in the vertical direction, while simultaneously their group velocity tends to zero and wave induced vertical shear significantly increases. This results in developing of shear instability, and, hence, to wave breaking and local intensification of turbulent mixing in the abyssal ocean. It is showed that similarly to wind wave breaking on a beach the abyssal ocean always represents a “surf zone” for near inertial waves.

18 October 2016

Scott Morgan

Alex McKay

Danny Groves

Stability of oscillatory rotating boundary layers

Viscoelastic lubrication

Contact line dynamics on heterogeneous substrates with mass transfer

25 October 2016

Dr Alex Bespalov (University of Birmingham)

Adaptive algorithms for high-dimensional parametric PDEs

Parametric PDEs are typical in optimisation problems and in mathematical models with inherent uncertainties (e.g., groundwater flow models).Differential operators in such PDEs depend on a large, possibly infinite, number of parameters, and naive application of numerical methods often results in the 'curse of dimensionality'. In this talk, we focus on a specific numerical method for solving such PDEs, namely on the stochastic Galerkin finite element method, for which we present an efficient adaptive algorithm. In this algorithm, we use an adaptive strategy to 'build' a polynomial space over a low-dimensional manifold in the infinitely-dimensional parameter space so that the total discretisation error is reduced most effectively,and thus the 'curse of dimensionality' is avoided. We will discuss the underlying theoretical results and demonstrate the performance of the algorithm in numerical experiments.

1 November 2016

Prof. Stephen Wilson (University of Strathclyde)

Floating plates and evaporating droplets

In this talk I will discuss two different, but not entirely unrelated, problems in which mathematical modelling and analysis can give some insight into real-world problems in fluid mechanics. In the first half of the talk I will discuss a paradigm problem motivated by recent work on fluid-structure interaction in microfluidics, namely a plate floating on the free surface of a viscous fluid [1]. In the second half of the talk I will discuss the lifetime of evaporating droplets [2]-[4]. Hopefully one or both halves of the talk will be interesting to those with an interest in the application of mathematics (and, in particular, asymptotic and numerical methods) to problems in fluid mechanics.

[1] Trinh, P.H., Wilson, S.K., Stone, H.A. A pinned or free-floating rigid plate on a thin viscous film, J. Fluid Mech. 760 407-430 (2014).

[2] Stauber, J.M., Wilson, S.K., Duffy, B.R., Sefiane, K. On the lifetimes of evaporating droplets, J. Fluid Mech. 744 R2 (2014).

[3] Stauber, J.M., Wilson, S.K., Duffy, B.R., Sefiane, K. Evaporation of droplets on strongly hydrophobic substrates, Langmuir 31 (12) 3653-3660 (2015).

[4] Stauber, J.M., Wilson, S.K., Duffy, B.R., Sefiane, K. On the lifetimes of evaporating droplets with related initial and receding contact angles, Phys. Fluids 27 (12) 122101 (2015).

8 November 2016

Dr Usama Kadri (Cardiff)

Acoustic-gravity waves, theory and applications

Acoustic–gravity waves (AGWs) are compression-type waves generated as a response to a sudden change in the water pressure, e.g. due to nonlinear interaction of surface waves, submarine earthquakes, landslides, falling meteorites and objects impacting the sea surface. AGWs can travel at near the speed of sound in water (ca. 1500 m/s), but they can also penetrate through the sea-floor surface doubling their speed, which turns them into excellent precursors. “Acoustic–gravity waves” is an emerging field that is rapidly gaining popularity among the scientific community, as it finds broad utility in physical oceanography, marine biology, geophysics, and water engineering, to name a few.

This talk is an overview on AGWs, with emphasis on four major applications that I will briefly discuss: (1) early detection of tsunami; (2) transportation of water in deep ocean; (3) detection of impacting objects (with a focus on the missing Malaysian airplane MH370); and (4) nonlinear triad resonance theory of AGWs and its possible implications.

15 November 2016

Walid Almalki

Martina Cracco

Alex Safar

Ahmed Alshehri

Using spectral element method for in mantle convection

Falkner-Skan flows of a non-Newtonian fluid

Two-step computer modelling of cellular bodies with intercellular contact

Partial differential equations using extended spectral element methods (XSEM)

22 November 2016

Dr Georgy Kitavtsev (University of Bristol)

Variational approaches to modelling surface energies in thin-film bilayer flows in microfluidics and martensitic twins in nonlinear elasticity

In this talk applications of non-convex variational approaches to two specific examples arising in material science will be presented. As first, it will be shown that the thin-film bilayer flows, when seen as a gradient flow of the total surface energy in the sharp interface limit, lead to a coupled system of lubrication equations equipped with the natural boundary conditions suggested previously by Kriegsmann and Miksis '03. A robust numerical algorithm for the thin-film gradient flow structure and evolution of the bilayer triple junction is then provided and tested on several examples.

In the second example a two-well non-convex Hamiltonian on a 2D atomic lattice describing the square-to-rectangular elastic transformations in shape-memory materials will be proposed as a model problem and subsequently analysed. The two wells (ground states) of the Hamiltonian are prescribed by two rank-one connected martensitic twins, respectively. It turns out that the Hamiltonian allows for a direct control of the discrete second order gradients and for a one-sided comparison with a two-dimensional spin system. Therefore, one can effectively proceed to a continuum limit and describe the explicit structure of the minimisers and their surface energies in the sharp-interface form.

These results is recent joint work with Sebastian Jachalski, Stephan Luckhaus, Dirk Peschka and Angkana Rueland.

6 December 2016

Dr David Sibley (Loughborough University)

Comparisons of a variety of physical and mathematical models for moving contact lines, including motion at the nanoscale

The moving contact line problem occurs when attempting to model the movement of the location where two fluid phases and a solid meet, as occurs when droplets spread (e.g. in inkjet printing), capillaries fill, insects walk on water, or in many other natural or technological instances. The problem exists when using the classical, macroscopic, equations of fluid motion as a singularity occurs in the predicted stresses and thus forces at the contact line, and the velocity is multi-valued. In this talk, we will look at the problem from the macroscale, and consider models that progressively retain more and more nanoscale physical features, culminating in an overview of results from two projects to understand motion at the nanoscale. The first is joint work with associates of the group of Prof Serafim Kalliadasis (Imperial College London) using dynamic density functional theory (DDFT) to explore the effect of the nanoscopic fluid structure on the motion of the contact line, and the second is work at Loughborough to understand droplet spreading where information from the nanoscale is captured by a function that is used in a coarse-grained model. At various points, joint work with Andreas Nold, Nikos Savva, Ben Goddard, Serafim Kalliadasis, Han Yu Yin, Uwe Thiele, and Andrew Archer will be presented.

17 January 2017

Anna Kalogirou (University of East Anglia)

Nonlinear dynamics of water waves and their impact on moving ships

The study of water waves has been an important area of research for years; their significance becomes obvious when looking at ocean and offshore engineering or naval architecture. Local weather and sea conditions can often lead to extreme wave phenomena, e.g. waves with irregular height. Waves with anomalously high amplitudes relative to the ambient waves are called rogue waves and can appear either at the coast or in the open ocean. The aim of this study is to investigate mathematically the generation and interaction of such waves and their impact on wave-energy devices and moving ships. The modelling is demonstrated by analysing variational methods asymptotically and numerically.

A reduced potential flow water-wave model is derived, based on the assumptions of waves with small amplitude and large wavelength. This model consists of a set of modified Benney-Luke equations describing the deviation from the still water surface and the velocity potential at the bottom of the domain. A novel feature in our model is that the dynamics are non-autonomous due to the explicit dependence of the equations on time. Numerical results obtained using a (dis)continuous Galerkin finite element method (DGFEM) are compared to a soliton splash experiment in a long water channel with a contraction at its end, resulting after a sluice gate is removed at a finite time. The removal of the sluice gate is included in the variational principle through a time-dependent gravitational potential.

The Benney-Luke approximation for water waves is also adapted to accommodate nonlinear ship dynamics. The new model consists of the classical water-wave equations, coupled to a set of equations describing the dynamics of the ship. We will first investigate the dynamics of the coupled system linearised around a rest state. For simplicity, we also consider a simple ship structure consisting of V-shaped cross-sections. The model is solved numerically using a DGFEM and the numerical results are compared to observations from experiments in wave tanks that employ geometric wave amplification to create nonlinear rogue-wave effects.

24 January 2017

Jens Eggers (University of Bristol)

Self-similar structure of caustics and shock formation

Caustics are places where the light intensity diverges, and where the wave front has a singularity. We use a self-similar description to derive the detailed spatial structure of a cusp singularity, from where caustic lines originate. We use this insight to study shock formation in the dKP equation, as well as shocks in classical compressible Euler dynamics. The spatial structure of these shocks is that of a caustic, and is described by the same similarity equation.

31 January 2017

Anne Juel (University of Manchester)

Fluid deposition and spreading in POLED applications

Microdroplet deposition is a technology that spans applications from tissue engineering to microelectronics. Our high-speed imaging measurements reveal how sequential linear deposition of overlapping droplets on flat uniform substrates leads to striking non-uniform morphologies for moderate contact angles. We develop a simple physical model, which for the first time captures the post-impact dynamics drop-by-drop: surface-tension drives liquid redistribution, contact-angle hysteresis underlies initial non-uniformity, while viscous effects cause subsequent periodic variations. Motivated by applications to the manufacture of POLED displays, we turn to the spreading of a single droplet within a recessed stadium-shaped pixel. We find that the sloping side wall of the pixel can either locally enhance or hinder spreading depending on whether the topography gradient ahead of the contact line is positive or negative. Locally enhanced spreading occurs via the formation of thin pointed rivulets along the side walls of the pixel through a mechanism similar to capillary rise in sharp corners. We demonstrate that a thin-film model combined with an experimentally measured spreading law, which relates the speed of the contact line to the contact angle, provides excellent predictions of the evolving liquid morphologies. We also show that the spreading can be adequately described by a Cox-Voinov law for the majority of the evolution.

7 February 2017

Peter Duck (University of Manchester)

Three-dimensional boundary states: States beyond the classical form

14 February 2017

Dr Angela Mihai (Cardiff)

Hyperelastic constitutive models for brain tissue

In some soft biological structures, such as brain, liver and fat tissues, strong experimental evidence suggests that the shear modulus increases significantly under increasing compressive strain, but not under tensile strain, while the apparent elastic modulus increases or remains almost constant when compressive strain increases. These tissues also exhibit a predominantly isotropic, incompressible behaviour. Our aim is to capture these seemingly contradictory mechanical behaviours, both qualitatively and quantitatively, within the framework of finite elasticity, by modelling a soft tissue as a homogeneous, isotropic, incompressible, hyperelastic material and comparing our results with available experimental data.

21 February 2017

Susanne Claus


28 February 2017

Kensuke Yokoi, Cardiff School of Engineering

Constrained interpolation profile conservative semi-Lagrangian scheme based on third-order polynomial functions and essentially non-oscillatory (CIP-CSL3ENO) scheme

We propose a fully conservative and less oscillatory multi-moment scheme for the approximation of hyperbolic conservation laws. The proposed scheme (CIP-CSL3ENO) is based on two CIP-CSL3 schemes and the essentially non-oscillatory (ENO) scheme. In this paper, we also propose an ENO indicator for the multi-moment framework, which intentionally selects non-smooth stencil but can efficiently minimize numerical oscillations.

The proposed scheme is validated through various benchmark problems and a comparison with an experiment of two droplets collision/separation. The CIP-CSL3ENO scheme shows approximately fourth-order accuracy for smooth solution, and captures discontinuities and smooth solutions simultaneously without numerical oscillations for solutions which include discontinuities. The numerical results of two droplets collision/separation (3D free surface flow simulation) show that the CIP-CSL3ENO scheme can be applied to various types of fluid problems not only compressible flow problems but also incompressible and 3D free surface flow problems.

7 March 2017

Christopher Davies

Impulsively Excited Disturbances in Non-Uniform Boundary Layers

Results will be reviewed for the linearized disturbance impulse response of non-uniform boundary layers.

Two distinct forms of boundary layer non-uniformity have been studied.
First, we consider the global behaviour of impulsively excited disturbances in temporally steady rotating-disc boundary layers, where there is a spatial inhomogeneity which stems from the radially increasing circumferential velocity.

We then consider the oscillatory Stokes layer that is driven by the time-periodic in-plane motion of a bounding flat plate. This provides a second type of boundary layer non-uniformity, which allows us to address the effects of base-flow unsteadiness upon the global development of disturbances.

14 March 2017

Nigel Peake (Cambridge)

The Aeroacoustic of the Owl

Many species of owl can hunt in acoustic stealth. The question of precisely how the owl actually manages to fly so quietly has remained open, but it has long been appreciated that owls which need to hunt silently possess two unique features, which are not found on any other bird, and indeed are not even found on owls which do not need to hunt silently (e.g small owls which feed on insects, or Fish Owls). First, the microstructure of the feathers on the upper wing surface is exceedingly complex, with an array of hairs and barbs which form a thick canopy just above the nominal wing surface. Second, the wing trailing edge possesses a small flexible and porous fringe which does not seem to have an obvious aerodynamic function.

The research I am going to describe in this talk is part of an ongoing theoretical (at Cambridge, Lehigh University and Florida Atlantic University) and experimental (at Virginia Tech.) program, with the aims of first attempting to understand how the two unique owl features described above actually work to control the noise, and then second of designing an owl-inspired treatment which can be used to significantly reduce aerodynamic noise generation in an engineering context.

The application we have in mind initially is to noise generation by onshore wind turbines, but there are many other contexts in which one wishes to reduce flow-structure noise where these ideas may be useful. In this talk I will give a flavour of the mathematical analysis, the experiments and the engineering applications.

21 March 2017

Pierre Kerfriden (Cardiff School of Engineering)


28 March 2017



4 April 2017

Silvia Gazzola (Bath)

Iterative regularization methods for large-scale inverse problems

Inverse problems are ubiquitous in many areas of Science and Engineering and, once discretised, they lead to ill-conditioned linear systems, often of huge dimensions: regularisation consists in replacing the original system by a nearby problem with better numerical properties, in order to find a meaningful approximation of its solution. After briefly addressing some classical regularisation strategies (such as Tikhonov method) and surveying some standard iterative regularisation methods (such as many Krylov methods), this talk will introduce the recent and promising class of the Krylov-Tikhonov iterative regularisation methods. In particular, strategies for choosing the regularization parameter and the regularization matrix will be emphasised.

25 April 2017

Wade Parsons (Memorial University of Newfoundland, Canada)

A mathematical model for acoustic-gravity waves generated by impulsive surface sources

Impulsive sources at the ocean surface generate propagating compression-type modes known as acoustic–gravity waves (AGWs) that travel in the water column at speeds near the speed of sound in water, i.e. c =1500 m/s leaving a measurable pressure signature. Possible sources include solid objects impacting the water surface, e.g. falling meteorites, landslides, sudden formation of rogue waves, or storm surges. A lot of promising work has been reported on AGWs in the last few years due to sea floor sources, with an emphasise on remote sensing as an early detection of tsunami; and only very recently have surface sources started to attract more attention. Here, we extend some of these studies to the remote sensing of general events generated at the ocean surface.

To this end, we developed an analytical model for AGWs generated by an impulsive source at the free surface (the Green’s function) from which the solution for various sources can be extracted. The results are compared with various solutions in literature whereby the source is located at the sea-floor. For the validation of the model, we carried out experiments in a water tank where neutrally buoyant spheres impacted the water surface, and the generated acoustic modes were recorded from a distance using a hydrophone. The shape of the pressure signature revealed three main regions that are associated with the impact, cavitation, and secondary wave formation. Employing these findings and solving the inverse problem allows remote sensing and prediction of the main source properties.

25 April 2017

Matthew Hunt (Brighton)

Surface waves with external electromagnetic fields

The topic of surface waves has a long history but what has been less known is how electromagnetic fields affect the wave profiles. This has important application in magnetohydrodynamics of the Sun and with electric fields, electrolysis. This talk will contain two parts, the first will describe a weakly 2D wave which is a generalisation of the KP equation called the 2D Benjamin equation. This has beenderived in the case of interfacial waves by Kim et al. We show how this equation naturally arrives in an MHD context. This is the first such nonlinear equation for surface waves which has been derived for 3D MHD. The second part of the talk will describe how one can remove the irrotational aspect of many flows by considering a global constant vorticity and deriving a corresponding free surface equation.

2 May 2017

Danny Groves (School of Mathematics, Cardiff)

Contact line dynamics on heterogeneous substrates with mass transfer

The contact line dynamics of two and three dimensional droplets spreading over chemically heterogeneous substrates are considered. Assuming small slopes and strong surface tension effects, a long wave expansion of the Stokes equations leads to a single equation for the droplet height where a contact line singularity is removed using a slip condition. Under a quasi-static regime we investigate cases where we have mass transfer through the substrate; modelling absorption or perhaps a needle piercing the droplet base. Utilising the method of matched asymptotic expansions, we approximate the solution of the governing partial differential equation by reducing to a set of ordinary differential equations. These ordinary differential equations are contrasted to a full numerical calculation using a pseudospectral collocation method, which, assesses the validity of the approximate method, as well as offer a glimpse into non quasi-static dynamics on spreading.

9 May 2017

Eric Lauga (Cambridge)

The Hydrodynamics of Swimming Bacteria

Many cellular organisms possess flagella, slender whiplike appendages which are actuated in a periodic fashion in fluids and allow the cells to self-propel. In particular, most motile bacteria are equipped with multiple helical rotating flagella which interact through the fluid, synchronise, and can form a tight helical bundle behind a swimming cell. We highlight in this talk two consequences of hydrodynamics for bacterial flagellar filaments. First we show how interactions between flagella mediated by the fluid allow them to repeatedly bundle and unbundle leading to reorientation of the whole cell during so-called `tumble’ events.

We next show how the flagellar flows induced by bacteria which have differentiated to a swarming state are responsible for large-scale fluid circulation at the scale of the whole swarm.

Past events

Applied and Computational Mathematics Seminars 2015-16