Towards adaptive simulations of turbulent wings at high Reynolds numbers

TL;DR

This paper demonstrates adaptive mesh refinement in high-order spectral-element simulations of turbulent flow over a wing, achieving accurate results with reduced computational cost at high Reynolds numbers.

Contribution

It introduces a non-conformal AMR implementation in Nek5000 for simulating turbulent wings at high Reynolds numbers, reducing computational costs and eliminating steady boundary condition errors.

Findings

Validated flow statistics against experimental data and DNS.

Achieved accurate flow simulations with fewer elements and lower cost.

Enabled study of complex turbulent flows at high Reynolds numbers.

Abstract

Adaptive mesh refinement (AMR) in the high-order spectral-element method code Nek5000 is demonstrated and validated with well-resolved large-eddy simulations (LES) of the flow past a wing profile. In the present work, the flow around a NACA 4412 profile at a chord-based Reynolds number $Re_c=200,000$ is studied at two different angles of attack: 5 and 11 degrees. The mesh is evolved from a very coarse initial mesh by means of volume-weighted spectral error indicators, until a sufficient level of resolution is achieved at the boundary and wake regions. The non-conformal implementation of AMR allows the use of a large domain avoiding the need of a precursor RANS simulation to obtain the boundary conditions (BCs). This eliminates the effect of the steady Dirichlet BCs on the flow, which becomes a relevant source of error at higher angles of attack (specially near the trailing edge and wake regions). Furthermore, over-refinement in the far field and the associated high-aspect ratio elements are avoided, meaning less pressure-iterations of the solver and a reduced number of elements, which leads to a considerable computational cost reduction. Mean flow statistics are validated using experimental data obtained for the same profile in the Minimum Turbulence Level (MTL) wind tunnel at KTH, as well as with a previous DNS simulation, showing excellent agreement. This work constitutes an important step in the direction of studying stronger pressure gradients and higher Reynolds complex flows with the high fidelity that high-order simulations allow to achieve. Eventually, this database can be used for the development and improvement of turbulence models, in particular wall models.