Applied and Computational Mathematics Research Group

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:

Theoretical and applied fluid mechanics

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, bubble dynamics, constitutive modelling of polymeric liquids.

Mathematics and mechanics of solids

Nonlinear elasticity, contact problems, limit states analysis, constitutive modelling in materials science.

Mathematical biology

The development of mathematical, computational and statistical methods to address biological and medical problems.

Applied analysis

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

Numerical analysis and scientific computing

The development of algorithms for the numerical solution of partial differential applications.

Head of Group

Tim Phillips photograph

Professor Tim Phillips

Head of School, Mathematics

Email:
phillipstn@cardiff.ac.uk
Telephone:
+44 (0)29 2087 4194

Academic staff

Photograph of Mikhail Cherdantsev

Dr Mikhail Cherdantsev

Lecturer

Email:
cherdantsevm@cardiff.ac.uk
Telephone:
+44 (0)29 2087 5549
Photograph of Chris Davies

Dr Chris Davies

Reader

Email:
daviesc9@cardiff.ac.uk
Telephone:
+44 (0)29 2087 5531
Photograph of Russell Davies

Professor Russell Davies

Honorary Distinguished Professor

Email:
daviesr@cardiff.ac.uk
Telephone:
+44 (0)29 2087 4827
Dr Usama Kadri photograph

Dr Usama Kadri

Lecturer in Applied Mathematics

Email:
kadriu@cardiff.ac.uk
Telephone:
029 208 75863
KaterinaPhoto

Dr Katerina Kaouri

Lecturer in Applied Mathematics

Email:
kaourik@cardiff.ac.uk
Telephone:
+44 (0)29 2087 5259
Angela Mihai

Dr Angela Mihai

Senior Lecturer in Applied Mathematics

Email:
mihaila@cardiff.ac.uk
Telephone:
+44 (0)29 2087 5570
Prof John Pryce photograph

Professor John Pryce

Emeritus Professor

Email:
prycejd1@cardiff.ac.uk
Telephone:
+44 (0)29 2087 4207
Photograph of Nikos Savva

Dr Nikos Savva

Lecturer in Applied Mathematics

Email:
savvan@cardiff.ac.uk
Telephone:
+44 (0)29 2087 5116
Thomas Woolley Profile Picture

Dr Thomas Woolley

Lecturer in Applied Mathematics

Email:
woolleyt1@cardiff.ac.uk
Telephone:
02920 870618

Seminars

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

Programme organiser and contact: Dr Usama Kadri

DateSpeakerSeminar

15 May 2018

Scott Morgan (School of Mathematics, Cardiff University)

Stability of Oscillatory Rotating Disk Boundary Layers

The rotating disk boundary layer has long been considered as providing an archetypal model for studying the stability of three-dimensional boundary-layer flows, and the crossflow inflexion point instability mechanism is common to both the rotating disk boundary layer and the flow over a swept wing. Thus the investigation of strategies for controlling the behaviour of disturbances that develop in the rotating disk flow may prove to be helpful for the identification and assessment of aerodynamical technologies that have the potential to maintain laminar flow over swept wings.

We will consider the changes in the stability behaviour that arise when the rotating disk base-flow configuration is altered by imposing a periodic modulation in the rotation rate of the disk surface. Thomas et. al. [Proc. R. Soc. A (2011) 467:2643-2662] have previously demonstrated that Tollmien-Schlichting waves can be stabilised when a similarly induced Stokes layer is conjoined to a plane channel flow.

Current work encompasses three distinct investigatory approaches. Linearised direct numerical simulations have been conducted, using the vorticity-based methods that were first adopted by Davies & Carpenter [J. Comput. Phys (2001) 172:119-165]. These simulations are complemented by a local in time linear stability analysis, that is made possible by imposing an artificial frozen base-flow approximation. This localised analysis is deployed together with a more exact global treatment based upon Floquet theory, which avoids the need for any simplification of the temporal dependency of the base-flow.

1 May 2018

Dr. Katerina Kaouri (School of Mathematics, Cardiff University)

Mathematical modelling for calcium signalling and for water resources management

Calcium signalling is one of the most important mechanisms of information propagation in the body. In embryogenesis the interplay between calcium signalling and mechanical forces is critical to normal embryonic development, but poorly understood. Several types of embryonic cells exhibit calcium-induced contractions and experiments indicate that calcium oscillations and contractions are linked via a two-way feedback mechanism. I will discuss some experiments and appropriate deterministic mathematical models that I have developed with Philip Maini and Jon Chapman (Oxford) and (the experimentalists) Paris Skourides (Univ. of Cyprus) and Neophytos Christodoulou (Cambridge). I will finish the math-bio part of my talk with stochastic models for calcium signalling that I have developed with Ruediger Thul (Nottingham) and Ioannis Lestas (Cambridge).

Furthermore, I will briefly present a water management project I have been working on with collaborators in Cyprus, UK and India. In collaboration with the Cyprus Water Development Department, we have been working on improved recharge protocols for the Germasogeia aquifer supplying water to Limassol, Cyprus. Furthermore, for the case of accidental contamination, we developed a model with which we can predict the extent and speed of contamination in the aquifer for various scenarios of practical interest. This project emerged from the 1st Study Group with Industry in Cyprus.

13 March 2018

Prof. Marco Marletta (School of Mathematics, Cardiff University)

An inverse problem in electromagnetism with partial data

Magnetic induction tomography involves determining the coefficients in the time-harmonic Maxwell equations (the permeability and permittivity) from surface measurements of current and voltage on a body. For practical reasons these measurements are made only on a part of the boundary. In this talk we show that partial boundary measurements uniquely determine the Maxwell system. This work is part of a joint project with Malcolm Brown (COMSC) and Paul Ledger (Swansea School of Engineering).

6 March 2018

Prof. Alain Goriely (Mathematical Institute, University of Oxford)

The mathematics and mechanics of brain morphogenesis

The human brain is an organ of extreme complexity. Its intricate folded shape has fascinated generations of scientists and has, so far, defied a complete description. How does it emerge? How is its shape related to its function? In this talk, I will review our current understanding of brain morphogenesis and its unique place within a general mathematical theory of biological growth. In particular, I will present simple models for basic pattern formation and show how they help us understand brain folding and skull formation.

13 February 2018

Dr. Kostas Soldatos (School of Mathematical Sciences, University of Nottingham)

On the theory of fibre-reinforced materials: past, present and future

The nonlinear theory of fibre-reinforced materials originated at the middle of the 20th century through pioneering attempts of Rivlin and Adkins to model mathematically large elastic deformation of rubber-like articles used commonly in industry, such as pneumatic tyres and fire hose. Today the use of the theory is spread worldwide and, apart from its usefulness in common industrial applications, it also enables modelling and understanding of complicated biological processes, including the behavior and growth of soft and hard biological tissue. This talk is naturally biased by the fact that the theory was essentially nurtured and, for long after the aforementioned pioneering attempts, developed through the research activity and effort of a group of Nottingham-based mathematicians. It aims to outline, briefly and lightly, theoretical features and principal relevant concepts that became known gradually over the years, assisted, and continue to assist the development of the theory. The foundation of the talk is accordingly based on the initial non-polar, and the more recent polar hyperelastic versions of the theory. However, attention is also paid to the the fact that, in the small deformation regime, the non-polar hyperelastic version of the theory reduces naturally to its older linear anisotropic elasticity counterpart, which has been, and is still used extensively in the static and dynamic analysis of advanced fibre-reinforced structural composites. If time allows, relevant modelling developments that concern plastic behaviour of fibre-reinforced solids and fluid-like behavior of fibre-reinforced resins will also be referred to.

6 February 2018

Dr. Rosemary Dyson (School of Mathematics, University of Birmingham)

Fibre-reinfoced fluids: from plants to extracellular matrix and beyond

Many biological systems depend on an underlying mechanical anisotropy to give the system required functional properties. This anisotropy is often created via fibres embedded within a ground matrix. For example cellulose microfibres within plant cell walls which enable directional pressure driven expansion and collagen fibres within extracellular matrix which guide cell behaviour. Similar ideas can be exploited within a synthetic biology context to investigate the properties of biological molecules via spectroscopy. We employ a common mathematical framework to study these diverse problems, which we discuss here.

23 January 2018

Dr. Georges Limbert (Engineering and Environment, University of Southampton)

A mechanistic approach to skin biophysics learning from mathematical and computational models

Besides the brain, no other organ of the human body plays such a central role in our everyday biological and social life than the skin. After all, this interface is the first line of defence of our body against the external environment and acts as a physical interface. It controls many types of exchanges between our inner and outside worlds which take the form of mechanical, thermal, biological, chemical and electromagnetic processes. Moreover, the skin tells a story about our health, age, past traumas, emotions, ethnicity, and our social and physical environments. Considering the place of the skin in our life and its multiple physiological functions, understanding its complex physiology and biophysics in health, disease and trauma has become, particularly in the last two decades, a broad multidisciplinary research arena.

To unravel some of the secrets of such a complex organ new experimental, imaging and computational techniques are needed and novel mechanistic theories explaining particular mechanobiological processes need to be formulated and put to the test. Developing and exploiting such an integrated framework underpin many aspects of our research which aims to understand the interplay between the microstructural and material properties of the skin, particularly as they evolve over the life course. As mounting evidence suggests, the skin microstructure can play a critical role in how macroscopic deformations are modulated at the microscopic level. These structural mechanisms are also at the heart of skin tribology by being part of, and conditioning mechanical load transmission and the nature of surface physics interactions. Skin biophysics is therefore fundamental to many industrial sectors from biomedical devices, personal care and cosmetic products to vehicle safety, textile, sport equipment, wearable electronics and tactile surfaces.

In this talk, I will present some of the latest modelling approaches we develop to gain a mechanistic understanding of the interplay between the material and structural properties of the skin, and ultimately, to exploit this knowledge for a variety of clinical and industrial applications. Examples will include computational contact homogenisation procedures to study skin friction, constitutive modelling of skin ageing and analysis of skin surface instabilities to understand mechanisms of wrinkle formation

16 January 2018

Prof. Frederic Dias (School of Mathematics and Statistics, University College Dublin)

What makes ocean waves go rogue in the real world?

The study of extreme ocean waves is a rapidly expanding area of research worldwide. Although much work in this area is based on modeling and experiments in controlled wave tanks, the starting point of all studies is wave observation in the natural world. During this talk, we will provide some evidence of extreme wave events, describe the main mechanisms for their generation and conclude with what we believe makes ocean waves go rogue in the real world.

19 December 2017

Prof. Till Bretschneider (Warwick Systems Biology Center, University of Warwick)

Image-based modelling of cell dynamics

Modern live-cell fluorescence microscopy enables us to visualise dynamic cellular processes in unprecedented detail. I will present ongoing research projects which are concerned with bringing together i) image analysis methods for tracking cells and their movements as well as quantifying spatio-temporal patterns of fluorescently labelled cellular constituents, and ii) mathematical models to investigate regulatory mechanisms of cellular biochemistry and mechanics.

12 December 2017

Dr. Andrey Melnik (School of Mathematics and Statistics, University of Glasgow)

Constitutive modelling of myocardium and other fibre-reinforced soft tissues

Nonlinear solid mechanics is used to model normal function and pathological conditions in soft tissues. In many cases, the mechanical behaviour of tissues can be adequately represented using an idealised elastic material. For instance, myocardium can be regarded as a passive (non-linear) hyperelastic solid with pronounced anisotropic properties due to its complex microstructure.

In the first part of the talk we examine basic ideas utilised in structural and semi-structural constitutive models, discuss the multiplicative decomposition framework for growth and remodeling (G&R), and highlight some challenges in using these elements.

The second part of the talk is dedicated to the Generalised Structure Tensor (GST) approach, which is used to formulate constitutive models for anisotropic fibre-reinforced composited with fibre distribution or dispersion [Gasser et al. JRS’06]. The GST approach has been so far successfully applied to models based on invariants I4 and I5, which capture the effect of deformation on each fibre family in isolation. We extend the GST approach to models based on the invariant I8, which couples two fibre families. Using the Holzapfel-Ogden model for myocardium, we demonstrate that accounting for fibre dispersion in the I8 term can have a significant effect on the predicted material response and may also reduce material symmetry.

27 November 2017

Dr. Eldad Avital (School of Engineering and Materials Sciences, Queen Mary University of London)

Challenges in Fluid-Structure Interaction from Renewable Energy Application to Bio-Fluids

In the first part of the talk we will briefly describe a recently developed computational methodology coupling fluid-flow simulation using finite volume with structural dynamics simulation using a combined finite-discrete element method. Applications in water flow as of sediments and hydro-kinetic turbines will be discussed, followed by looking at recent designs for hydrogen fusion power and the challenges. Bio-fluid applications of red blood cell flow and upper renal system will also be illustrated. In the second part of the talk we will discuss our aerodynamic method for blade design calling for continuous surface curvature. It will be computationally and experimentally shown to increase aerodynamic efficiency, particularly at high incidence where stall delay can be achieved. Reduction of tonal noise for low Reynolds number blade sections will also be addressed.

21 November 2017

Prof. Alexander Korobkin (School of Mathematics, University of East Anglia)

Diffraction of hydroelastic waves by a circular cylinder

Linear problem of wave diffraction is studied for a circular cylinder mounted at the sea bed and piercing the fluid surface which is covered by ice plate of infinite extent. The water depth is constant. The ice plate is modeled by a thin elastic plate of constant thickness clamped to the surface of the cylinder. One-dimensional incident hydroelastic wave of small amplitude propagates towards the cylinder. and is diffracted on the cylinder.  Deflection of the ice plate and the bending stresses in it are determined by two methods: (a) using the integral Weber transform in radial direction, (b) using the vertical modes for the fluid of constant depth with the rigid bottom and elastic upper boundary. The solution by the second method is straightforward but we cannot prove that the solution is complete because the properties of the vertical modes are not known yet.  The solution by the Weber transform is more complicated but this solution is unique. In this talk we will show that these two solutions are identical. This result justifies the method of the vertical modes in the hydroelastic wave diffraction problems.

14 November 2017

Dr. John A. Mackenzie (Department of Mathematics and Statistics, University of Strathclyde)

An Adaptive Moving Mesh Method for Geometric Evolutions Laws and Bulk-Surface PDEs: Application to a Model of Cell Migration and Chemotaxis

In this talk I will consider the adaptive numerical solution of curve-shortening flow with a driving force. An adaptive moving mesh approach is used to distribute the mesh points in the tangential direction. This ensures that the resulting meshes evolve smoothly in time and are well adjusted to resolve areas of high curvature. Experiments will be presented to highlight the improvement in accuracy obtained using the new method in comparison with uniform arc-length mesh distributions. We will also discuss the use of the evolving adaptive curve mesh in the adaptive generation of bulk meshes for the solution of bulk-surface PDEs in time dependent domains.

The main motivation for developing these computational tools is the modelling of single cell migration and chemotaxis. Chemoattractant gradients are usually considered in terms of sources and sinks that are independent of the chemotactic cell. However, recent interest has focused on “self-generated” gradients, in which cell populations create their own local gradients as they move. Here we consider the interplay between chemoattractants and single cells. To achieve this we model the breakdown of extracellular attractants by membrane-bound enzymes. Model equations are parameterised using published estimates from Dictyostelium cells chemotaxing towards cyclic AMP. We find that individual cells can substantially modulate their local attractant field under physiologically appropriate conditions of attractant and enzymes. This means the attractant concentration perceived by receptors can be a small fraction of the ambient concentration. This allows efficient chemotaxis in chemoattractant concentrations that would be saturating without local breakdown.

7 November 2017

Dr. Richard Hewitt (School of Mathematics, University of Manchester)

Localised streaks in a Blasius boundary layer

Streaks are common feature of perturbed boundary-layer flows. They play a central role in transient growth mechanisms and are a building block of exact coherent structures. Most theoretical work has focused on streaks that are periodic in the spanwise direction, but in this work we consider a single spatially localised streak embedded into a Blasius boundary layer. For small streak amplitudes, we show the perturbation can be described in terms of a set of eigenmodes that correspond to an isolated streak/roll structure. These modes are new, and arise from a bi-global eigenvalue calculation; they decay algebraically downstream and may be viewed as the natural three-dimensional extension of the two-dimensional Libby \& Fox (1963, JFM vol. 17) solutions. Despite their bi-global nature, we show that a subset of these eigenmodes is fundamentally related to both the Libby & Fox solutions, and those presented by Luchini (1996, JFM vol. 327), as derived for (spanwise) periodic disturbances at small spanwise wavenumber. This surprising connection is made by an analysis of the far-field decay of the bi-global state. We also address the downstream development of nonlinear streaks, confirming that the aforementioned eigenmodes are recovered as the streak/roll decays downstream. Some comparisons are made with available experimental data.

24 October 2017

Dr. Václav Klika (Department of Mathematics, Czech Technical University in Prague)

On modelling self-organisation in real systems

Nowadays models for self-organisation are being used in systems with a great degree of complexity and across disciplines. We show that the widely used Turing model is sensitive to inputs, type of domain growth, but also to the precision of model formulation itself. Hence a great care is needed when applying Turing's model for self-organisation to real problems. For this purpose we consider derivation of evolution equations within non-equilibrium thermodynamic to identify physically relevant formulations. Only then we subject these models to a detailed mathematical analysis. We offer possible extensions of the concept of self-organisation to more general situations and discuss its physical interpretation.

The essence and importance of these ideas is illustrated on the reaction-diffusion-advection system, where we indicate that such a system should be preferred from both physical and mathematical viewpoint. Further we point to the indispensable role of physical viewpoint during relevant model formulations. Using the non-equilibrium thermodynamic framework physically consistent extensions of Turing model are revealed as well as functional constraints for present parameters.

18 October 2017

Scott Morgan Scott Morgan (School of Mathematics, Cardiff University)

Stability of Oscillatory Rotating Disk Boundary Layers

The rotating disk boundary layer has long been considered as providing an archetypal model for studying the stability of three-dimensional boundary-layer flows, and the crossflow inflexion point instability mechanism is common to both the rotating disk boundary layer and the flow over a swept wing. Thus the investigation of strategies for controlling the behaviour of disturbances that develop in the rotating disk flow may prove to be helpful for the identification and assessment of aerodynamical technologies that have the potential to maintain laminar flow over swept wings.

We will consider the changes in the stability behaviour that arise when the rotating disk base-flow configuration is altered by imposing a periodic modulation in the rotation rate of the disk surface. Thomas et. al. [Proc. R. Soc. A (2011) 467:2643-2662] have previously demonstrated that Tollmien-Schlichting waves can be stabilised when a similarly induced Stokes layer is conjoined to a plane channel flow.

Current work encompasses three distinct investigatory approaches. Linearised direct numerical simulations have been conducted, using the vorticity-based methods that were first adopted by Davies & Carpenter [J. Comput. Phys (2001) 172:119-165]. These simulations are complemented by a local in time linear stability analysis, that is made possible by imposing an artificial frozen base-flow approximation. This localised analysis is deployed together with a more exact global treatment based upon Floquet theory, which avoids the need for any simplification of the temporal dependency of the base-flow.

10 October 2017

Dr. Thomas Woolley (School of Mathematics, Cardiff University)

Patterns, cellular movement and brain tumors

I present three pieces of work that illustrate the power of mathematics as a tool for understanding biology. Although the applications appear to be disparate the underlying mathematics is very similar.

I begin by looking at theoretical and experimental pattern formation, with emphasis on whisker formation in mice. Here, reaction-diffusion equations are used to provide insights into how the wavelength of the whiskers are controlled.

Next, I consider the phenomena of blebbing cells. Initially, I use a diffusion equation to understand the motion of muscle stem cells and illustrate how old cells fundamentally move differently to old cells. This is then extended to include solid mechanics, which allows us to link the structural properties of the cell to their motion.

Finally, reaction-diffusion equations are used to understand the formation of brain tumours. Critically, the cells move at different speeds in white and grey matter, including this information can lead to very different migration patterns of the tumours.

3 October 2017

Dr. Davide Crivelli (School of Engineering, Cardiff University)

Detecting damage with waves: an overview of structural health monitoring research in Engineering

Acoustic emission is a promising technique for monitoring damage on structures. It provides a non-invasive, real-time and passive system for detecting cracks and defects in a range of structures as they happen. Applications range from civil structures, such as bridges and railways; rotating machinery; aerospace, such as impact and damage detection in advanced composite structures.

We face multiple challenges however when the amount of data from an Acoustic Emission system has to be interpreted in a reliable way. Wave propagation, source characterization, interaction with damage and sensor characteristics make traditional signal interpretation and source localization a challenging task.

This talk will give a general overview of acoustic emission for structural health monitoring. It will present a series of case studies, describing how Cardiff research has tackled the main challenges of signal localization and interpretation, and how the research has helped bringing the technology closer to commercial applications. It will also present a collaboration with Maths stemmed from acoustic emission research and applied to low frequency acoustic gravity waves in the search of MH370.

Past events

Applied and Computational Mathematics Seminars 2016-17

Applied and Computational Mathematics Seminars 2015-16