No abstract

No abstract

No abstract

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Highlights• Discretization and residual error refinement indicators for DG methods are analysed.• p-adaptive simulations of steady and unsteady configurations are performed.• Performance is compared based on convergence history and refinement distribution.• Computational cost and possible implementation issues are discussed.

We assess the performance of a high-order discontinuous Galerkin (DG) approach to simulate wall-bounded turbulent flows for a range of Reynolds numbers. First, a plane channel flow configuration with turbulent heat transfer at Re τ = 640 is considered. Second, a detached flow configuration is investigated. This configuration represents the flow over periodically arranged hills. Four bulk Reynolds numbers Re b are considered: 2 800, 10 595, 19 000, and 37 000. For Re b = 2 800 direct numerical simulation (DNS) is used. Large-eddy simulations (LES) based on the WALE subgrid-scale model are carried out for the higher Re b . The simulation results are compared to reference data from CFD and experiment. hp-convergence analyses are performed which demonstrate the superior performance of increasing the polynomial order as compared to refining the mesh. It appears from this work that the use of a subgrid modelling approach together with local hp-adaptation could greatly improve the accuracy of the solution without penalising the computational cost of the simulation.