Next generation finite element methods for wave problems

Lead Research Organisation: University of Reading
Department Name: Mathematics and Statistics

Abstract

Efficient and accurate simulation of wave phenomena is a key enabling technology across science and engineering. Applications span diverse areas, and include the whole of acoustics and noise control, non-destructive testing and ultrasonic and microwave technologies for medical imaging, problems of seismic and radar propagation and imaging, and even quantum scale simulations. But even though the underlying partial differential equations are usually linear and well understood, wave phenomena are complex and hard to simulate whenever the wavelength is small compared to the diameter of the region to be simulated.A main and standard computational tool for simulations of wave problems is the so-called finite element method. The idea of the method is to break up the computational domain into small elements and to approximate the solution on each of them in a simple way, e.g. as a linear variation. However, this gives accurate solutions only if the diameter of each element is small compared to the wavelength. Thus the number of elements needed and the associated computational cost and storage is infeasible if the diameter of the region to be simulated is very large compared to the wavelength, as it is for very many complex problems of wave propagation and scattering, e.g. seismic wave propagation for hydrocarbon exploration.Recently, there has been strong international interest in novel finite element formulations that try to solve this problem by representing the wave field on each element by functions that are themselves waves. This allows much bigger element sizes and so a significant reduction of the computational cost. However, these novel finite element methods are still in their infancy and it is poorly understood how to implement them in an optimal way. For example, one key open problem is the question of which wave functions to use. Another open question is how to achieve numerical stability, i.e. an algorithm whose results are not garbled by effects resulting from the limited accuracy that computers have. These and other questions are particularly unclear for three dimensional problems, although most practical applications are three dimensional.The fellowship addresses this wide open research area. Building upon novel ideas about how to locally model wave phenomena in a stable way it combines fundamental research in diverse areas of applied and computational mathematics in order to develop the next generation of finite element methods for wave problems. These new methods have the potential to be orders of magnitude faster than current methods allowing for numerical simulations of phenomena that are currently out of reach. In close collaboration with partners in science and industry the new methods will be applied to exciting research problems in science and engineering. In particular, a major part of the hydrocarbon exploration business is enabled through the modelling and inversion of large scale 3D seismic and electromagnetic data sets, and Schlumberger Cambridge Research will be a key project partner. Throughout the fellowship annual international workshops on next generation finite element methods for wave problems will be organised, at Reading and Schlumberger. These will bring together leading researchers in the area of numerical wave simulations from academia and industry and will drive this research area forward by intensifying collaborations and developing and exploring application areas for these methods.Numerical wave simulations are an essential technology in science and engineering. Innovations in many areas depend upon the ability to simulate complex wave phenomena. The UK is one of the leading countries for wave-related research. This fellowship will enhance this role by building up an internationally outstanding research group on novel finite element methods for wave problems that will have a strong impact on wave-related research and applications long after the duration of the fellowship.

Publications

10 25 50

publication icon
Nigam N (2011) Numerical integration for high order pyramidal finite elements in ESAIM: Mathematical Modelling and Numerical Analysis

Related Projects

Project Reference Relationship Related To Start End Award Value
EP/H004009/1 01/10/2009 01/02/2011 £795,615
EP/H004009/2 Transfer EP/H004009/1 01/02/2011 30/09/2014 £608,157
 
Description Grant was transferred to UCL. See EP/H004009/2 for report.
Exploitation Route Grant was transferred to UCL. See EP/H004009/2 for report.
Sectors Aerospace, Defence and Marine

 
Description Grant was transferred to UCL. See EP/H004009/2 for report.