Search

Laminar shock-wave/boundary-layer interactions were simulated using OpenSBLI, a Python-based source code generation framework. Shock-capturing was performed by a 5th order finite-difference Weighted Essentially Non-Oscillatory (WENO)-Z scheme applied in characteristic space. Oblique shock conditions were imposed for a shock angle of θ = 32 . 58 ∘ and Mach 2 free-stream, impinging on a laminar flat-plate boundary-layer. Performance of the code was assessed on different architectures for CPU, GPU and Xeon Phi.

This work introduces a new error indicator which can be used to determine areas of insufficient numerical resolution in unfiltered finite difference simulations. The background behind the methodology is that smaller scales (i.e. the flow features with higher wave numbers) are physically characterised by a smaller energy content in comparison with larger scales. This energy should decrease with increasing wavenumber at a minimum rate; if this rate is not attained it likely means that the smaller scales are not being properly resolved on the computational grid of solution points. An approach using spectral techniques is used to formulate two varieties of the error indicator – one integer-valued and one floating point-valued. These values are computed at a finite number of ‘blocks’ which span the domain. The indicator is implemented within the OpenSBLI finite difference-based modelling framework, and evaluated in the context of a three-dimensional Taylor-Green vortex problem and flow past a V2C laminar flow aerofoil.