We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

Frequency domain Bernstein-Bezier finite element solver for modelling short waves in elastodynamics

Benatia, N. and El Kacimi, A. and Laghrouche, O. and El Alaoui, M. and Trevelyan, J. (2022) 'Frequency domain Bernstein-Bezier finite element solver for modelling short waves in elastodynamics.', Applied Mathematical Modelling, 102 . pp. 115-136.


This work presents a high-order Bernstein-Bézier finite element (FE) discretisation to accurately solve time harmonic elastic wave problems on unstructured triangular mesh grids. Although high-order FEs possess many advantages over standard FEs, the computational cost of matrix assembly is a major issue in high-order computations. A key ingredient to address this drawback is to resort to low complexity procedures in building the local high order FE matrices. This is achieved in this work by exploiting the tensorial property of Bernstein polynomials and applying the sum factorisation method for curved elements. An efficient implementation of the analytical rules for affine elements is also proposed. Furthermore, element-level static condensation of the interior degrees of freedom is performed to reduce the memory requirements. Additionally, the applicability of the method with a variable polynomial order, based on a simple a priori indicator, is investigated. The computational complexities of sum factorisation, analytical rules and standard quadrature are first evaluated, in terms of the CPU time against the polynomial order. The analysis shows that the achieved numerical complexities compare favourably to those expected theoretically. A significant runtime saving is also obtained by using analytical rules and sum factorisation. The performance of the Bernstein-Bézier FEs is then assessed on various benchmark tests, over a wide range of frequencies. Results from the elastic wave scattering problem demonstrate the effectiveness of this method in coping with the pollution error, and its accuracy in resolving high order evanescent wave modes. Additionally, a wave transmission problem with high wave speeds contrast and a curved interior interface is considered, where a simple a priori indicator is proposed to assign the variable polynomial order. The study provides evidence of the great benefit of a non uniform p-refinement in reducing the computational cost and enhancing accuracy.

Item Type:Article
Full text:Publisher-imposed embargo until 01 October 2022.
(AM) Accepted Manuscript
Available under License - Creative Commons Attribution Non-commercial No Derivatives 4.0.
File format - PDF
Publisher Web site:
Publisher statement:© 2021 This manuscript version is made available under the CC-BY-NC-ND 4.0 license
Date accepted:23 September 2021
Date deposited:05 October 2021
Date of first online publication:01 October 2021
Date first made open access:01 October 2022

Save or Share this output

Look up in GoogleScholar