OpenSBLI v3.0: High-fidelity multi-block transonic aerofoil CFD simulations using domain specific languages on GPUs

OpenSBLI is an automatic code-generation framework for compressible Computational Fluid Dynamics (CFD) simulations on heterogeneous computing architectures (previous release: Lusher et al. (2021) [4]). OpenSBLI is coupled to the Oxford Parallel Structured (OPS) Domain Specific Language (DSL), which uses source-to-source translation to enable parallel execution of the code on large-scale supercomputers, including multi-GPU clusters. To date, OpenSBLI has largely been applied to compressible turbulence and shock-wave/boundary-layer interactions on very simple geometries comprised of single mesh blocks with essentially orthogonal grid lines. OpenSBLI has been extended in this new release to target strongly curvilinear cases, including transonic aerofoils using multi-block grids. In addition to multi-block mesh support, more efficient numerical shock-capturing methods and filters have been added to the codebase. Improvements to post-processing, reduced-dimension data output, and coupling to a modal decomposition library are also included. A set of validation cases are presented to showcase the new code features. Furthermore, state-of-the-art wide-span transonic aerofoil simulations on up to N=2.5×109 grid points demonstrate that wider aspect ratios can alter buffet predictions and increase the regularity of the low-frequency shock oscillations by accommodating fully-developed trailing edge flow separation. Spectral Proper Orthogonal Decomposition (SPOD) analysis showed that overly-narrow aerofoil simulations contain additional domain-dependent energy content at a Strouhal number of St≈3 associated with wake modes. PROGRAM SUMMARY Program Title OpenSBLI v3.0 CPC Library link to program files"https://doi.org/10.17632/3sdb6hck2c.2 Developer's repository link:https://github.com/opensbli/opensbli Licensing provisions: GPLv3 Programming language: Python, C/C++, OPS DSL Journal reference of previous version: Comput. Phys. Commun. 267 (2021), 108063 Does the new version supersede the previous version?: Yes. Reasons for the new version: Multi-block mesh support. Summary of revisions: Added multi-block mesh support, cylinder and aerofoil validations, new shock-capturing methods, new filters, improved data output and post-processing. Nature of problem: The compressible 3D Navier-Stokes equations are solved via Implicit Large Eddy Simulation (ILES) or Direct Numerical Simulation (DNS). Solution method: OpenSBLI [1] is a Python-based code-generation system that uses symbol