Applied Mathematics Seminars 

Spring 2017

All seminars (unless otherwise stated) will take place on Tuesdays at 3.00pm in Room Archeology G6 (31-34 Gordon Square). See the map for further details. There will be tea afterwards in Mathematics Room 606 (25 Gordon Street). If you require any more information on the Applied seminars please contact Prof Slava Kurylev e-mail: y.kurylev AT or tel: 020-7679-7896.

17 January 2017

Dr Max Jensen (Sussex University)

Title:  Convergent semi-Lagrangian methods for the Monge-Ampère equation on unstructured grids

In this presentation I will discuss a semi-Lagrangian discretisation of the Monge-Ampère operator on P1 finite element spaces. The wide stencil of the scheme is designed to ensure uniform stability of numerical solutions. Monge-Ampère equations arise for example in the Inverse Reflector Problem where the geometry of a reflecting surface is reconstructed from the illumination pattern on a target screen and the characteristics of the light source.

Monge-Ampère type equations, along with Hamilton-Jacobi-Bellman type equations are two major classes of fully nonlinear second order partial differential equations (PDEs). From the PDE point of view, Monge-Ampère type equations are well understood.  On the other hand, from the numerical point of view, the situation is far from ideal. Very few numerical methods, which can reliably and efficiently approximate viscosity solutions of Monge-Ampère type PDEs on general convex domains. There are two main difficulties which contribute to the situation:

* Firstly, it is well known that the fully nonlinear structure and nonvariational concept of viscosity solutions of the PDEs prevent a direct formulation of any Galerkin-type numerical methods.

* Secondly, the Monge-Ampère operator is not an elliptic operator in generality, instead, it is only elliptic in the set of convex functions and the uniqueness of viscosity solutions only holds in that space.  This convexity constraint, imposed on the admissible space, causes a daunting challenge for constructing convergent numerical methods; it indeed screens out any trivial finite difference and finite element analysis because the set of convex finite element functions is not dense in the set of convex functions

The goal of our work is to develop a new approach for constructing convergent numerical methods for the Monge-Ampère Dirichlet problem, in particular, by focusing on overcoming the second difficulty caused by the convexity constraint. The crux of the approach is to first establish an equivalent (in the viscosity sense) Bellman formulation of the Monge-Ampère equation and then to design monotone numerical methods for the resulting Bellman equation on general triangular grids. An aim in the design of the numerical schemes was to make Howard's algorithm available, which is a globally superlinearly converging semi-smooth Newton solver as this allows us to robustly compute numerical approximations on very fine meshes of non-smooth viscosity solutions. An advantage of the rigorous convergence  analysis of the numerical solutions is the comparison principle for the Bellman operator, which extends to non-convex functions. We deviate from the established Barles-Souganidis framework in the treatment of the boundary conditions to address challenges arising from consistency and comparison. The proposed approach also bridges the gap between advances on numerical methods for these two classes of second order fully nonlinear PDEs.

The contents of the presentation is based on joint work with X. Feng from the University of Tennessee.

24 January 2017

Dr Angela Mihai (Cardiff University)

Title: 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.  This is joint work with Prof. Alain Goriely (University of Oxford).

7 February 2017

Professor Andrew Soward (Newcastle Unviersity)


21  February 2017

Dr Joel Daou (University of Manchester) 

Title: Flame propagation and Taylor dispersion


7 March 2017 

Please see the Departmental Colloquia webpage

21 March 2017

Dr Maria Korotyaeva

Title: The resolvent method for shear waves spectra calculation in 2D phononic crystals

We propose the resolvent method for calculating the horizontally polarized shear waves spectra in the 2D phononic crystal (PC) waveguides: the free PC plate and the PC plate sandwiched between two substrates. Since the propagator over a unit cell approximated by Fourier harmonics in one coordinate can have very large components, we introduce its resolvent as a numerically stable substitute.