UCL Discovery
UCL home » Library Services » Electronic resources » UCL Discovery

Gradient jump penalty stabilisation of spectral/hp element discretisation for under-resolved turbulence simulations

Moura, RC; Cassinelli, A; Da Silva, AFC; Burman, E; Sherwin, SJ; (2022) Gradient jump penalty stabilisation of spectral/hp element discretisation for under-resolved turbulence simulations. Computer Methods in Applied Mechanics and Engineering , 388 , Article 114200. 10.1016/j.cma.2021.114200.

[thumbnail of Gradient_Jump_Penalty-Oct2021.pdf] Text
Gradient_Jump_Penalty-Oct2021.pdf - Accepted version
Access restricted to UCL open access staff until 26 October 2022.

Download (9MB)

Abstract

One of the strengths of the discontinuous Galerkin (DG) method has been its balance between accuracy and robustness, which stems from DG’s intrinsic (upwind) dissipation being biased towards high frequencies/wavenumbers. This is particularly useful in high Reynolds-number flow simulations where limitations on mesh resolution typically lead to potentially unstable under-resolved scales. In continuous Galerkin (CG) discretisations, similar properties are achievable through the addition of artificial diffusion, such as spectral vanishing viscosity (SVV). The latter, although recognised as very useful in CG-based high-fidelity turbulence simulations, has been observed to be sub-optimal when compared to DG at intermediate polynomials orders (P ≈ 3). In this paper we explore an alternative stabilisation approach by the introduction of a continuous interior penalty on the gradient discontinuity at elemental boundaries, which we refer to as a gradient jump penalisation (GJP). Analogous to DG methods, this introduces a penalisation at the elemental interfaces as opposed to the interior element stabilisation of SVV. Detailed eigenanalysis of the GJP approach shows its potential as equivalent (sometimes superior) to DG dissipation and hence superior to previous SVV approaches. Through eigenanalysis, a judicious choice of GJP’s P-dependent scaling parameter is made and found to be consistent with previous apriori error analysis. The favourable properties of the GJP stabilisation approach are also supported by turbulent flow simulations of the incompressible Navier-Stokes equation, as we achieve high-quality flow solutions at P = 3 using GJP, whereas SVV performs marginally worse at P = 5 with twice as many degrees of freedom in total.

Type: Article
Title: Gradient jump penalty stabilisation of spectral/hp element discretisation for under-resolved turbulence simulations
DOI: 10.1016/j.cma.2021.114200
Publisher version: https://doi.org/10.1016/j.cma.2021.114200
Language: English
Additional information: This version is the author accepted manuscript. For information on re-use, please refer to the publisher’s terms and conditions.
Keywords: High-order discretisation, Spectral/hp element method, Stabilised methods, Velocity splitting scheme, dispersion–diffusion analysis, eigenanalysis
UCL classification: UCL
UCL > Provost and Vice Provost Offices > UCL BEAMS
UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences
UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences > Dept of Mathematics
URI: https://discovery.ucl.ac.uk/id/eprint/10137213
Downloads since deposit
1Download
Download activity - last month
Download activity - last 12 months
Downloads by country - last 12 months

Archive Staff Only

View Item View Item