## Members

Chris Budd

Matthias Ehrhardt

Melina Freitag

Silvia Gazzola

Ivan Graham

Adrian Hill

Eike MÃ¼ller

Rob Scheichl

Tony Shardlow

Alastair Spence

Euan Spence

# Numerical Analysis and Scientific Computing

Numerical analysis refers to the development and implementation of algorithms for the efficient determination of computational solutions to mathematical problems. Efficiency can, depending on the context, mean algorithmic accuracy, stability, or speed; usually all three of these are important.

Research in numerical analysis at Bath encompasses a wide spectrum of modern topics. A connecting thread which runs through much of our research is the derivation, analysis and implementation of computational techniques for Partial Differential Equations (PDEs) of various types, or their reduction to a system of Differential Algebraic Equations (DAEs).

There is a weekly informal seminar series, followed by a research group lunch, along with regular research workshops organised with industry (for example the Rutherford Appleton Laboratory) and with other academic institutions. There is a weekly informal seminar series, followed by a research group lunch, along with regular research workshops organised with industry (for example the Rutherford Appleton Laboratory) and with other academic institutions. The Numerical Analysis group webpages, contain further details of current events.

Our research in numerical analysis can be arranged under the following broad headings

## Data assimilation

Chris Budd, Melina Freitag, Rob Scheichl

Data assimilation problems arise for example in weather prediction and hydrology where observations of variables are incorporated into a mathematical model. The goal is to improve the accuracy of the model and make optimal predications of the future state of the system. Current data assimilation algorithms such as variational data assimilation and Kalman filters lead to large-scale problems in numerical linear algebra that need to be solved efficiently. Data assimilation belongs to the class of inverse problems, which are typically ill-conditioned and require regularisation methods such as Tikhonov regularisation to stabilise the problem. We collaborate with the UK Met Office on data assimilations in weather prediction.

## Geometric integration

The image shows snapshots of the evolution of a cat shape under a Hamiltonian flow map. The map is symplectic and while the shape of the cat changes its area is preserved. To explore the delicate dynamics of such systems, it is necessary that such geometric features are preserved by the numerical approximation in order to respect the underlying physics. The group in Bath looks at developing general linear methods for Hamiltonian differential equations and adaptive numerical methods for capturing shocks and singularities in PDEs with conservation laws.

- Budd, C.J. and Piggott, M.D. Geometric Integration and its Applications. Handbook of Numerical Analysis Vol.11 (2003). link

## High frequency scattering

Our understanding of acoustic and electromagnetic wave scattering underpins technologies such as radar, sonar, telecommunications, ultrasound, and seismic imaging. Many of these applications require the accurate computation of high frequency waves. However, standard numerical methods become increasingly expensive as the frequency increases, putting very high frequency problems out of range of standard methods. Recent research in Bath has looked at ways to combat this, including the design of 'hybrid numerical-asymptotic methods', and novel preconditioners for use in conjunction with domain decomposition methods. This research currently involves collaboration with Schlumberger, the UK Met Office, BAE Systems, and the Institute for Cancer Research.

- Chandler-Wilde, Simon N. and Graham, Ivan G. and Langdon, Stephen and Spence, Euan A. Numerical-asymptotic boundary integral methods in high-frequency acoustic scattering. Acta Numer. Vol.21 (2012). link

## Moving mesh methods

Numerical methods for PDEs usually require the construction of a mesh on which the PDE is then discretised. However, if a fixed mesh is used during the computation then this can lead to poor resolution of the evolving features of the PDE and also the loss of important qualitative formation such as conservation and balance laws. This is a major problem in many areas of computation of the inherently multi-scale problems that arise in nature, a particular example being meteorology. Budd is actively working in the construction of moving meshes in 3D which progressively move mesh points to resolve emergent features. These have already proved very effective for a wide range of problems and are now used as part of a data assimilation calculation in the UK Met Office Operational Code. Much new work needs to be done to extend these methods to work on the whole sphere and to study the complex PDEs encountered in nature. Moving mesh methods are a lovely mix of differential geometry, dynamical systems theory and numerical analysis with major practical applications.

- Budd, C. and Huang, W. and Russel, R. Adaptivity with moving grids. Acta Numerica Vol.18 (2009). link

## Numerical linear algebra

Melina Freitag, Silvia Gazzola, Alastair Spence

Eigenvalue problems are central to a number of applications in science and industry, such as structural dynamics, quantum mechanics, material science and stability analysis. Often the problems are very large and only a few eigenvalues close to a target value are required. Many methods for solving eigenvalue problems iteratively are available, such as the inverse iteration, the Jacobi-Davidson method, and Krylov-methods such as the Lanczos and Arnoldi method. Linear and nonlinear eigenvalue problems that depend on a scalar parameter arise, for example in aerodynamics, electrical power systems and matrix perturbation theory, and are important in determining bifurcations numerically. It is an important research problem to find critical parameter values for parameter dependent eigenproblems such that, for example, eigenvalues coalesce (and become defective, that is, form a Jordan block), eigenvalues cross the imaginary axis (that is, stable matrices become unstable) or move to another critical region in the complex plane.

## Stochastic differential equations

Ivan Graham, Rob Scheichl, Tony Shardlow

Uncertainty quantification is a fundamental problem in all types of mathematical modelling. Due to the rapid increase in computer power, it is now possible to incorporate uncertainty into many types of differential equations and compute accurate numerical approximations to sample paths or distributions. The group in Bath has interests in stochastic differential equations as an ordinary, delay, or reaction diffusion equation with Brownian forcing, where stochastic calculus is necessary to understand the model. A different type of stochastic differential equation are problems with random data, such as an elliptic PDEs with random porosity coefficient modelling flow through a porous media. The PDE can be understood in the conventional sense for each sample of the coefficients, while advanced techniques from probability yield accurate approximation techniques.

Lord, G. and Powell, C. and Shardlow, T. An introduction to Computational Stochastic PDEs. CUP (2014). link