On Capturing Laminar/Turbulent Regions Over a Wing Using WMLES

TL;DR

This paper investigates grid resolution requirements in wall-modeled large-eddy simulations of flow over a wing, focusing on accurately capturing both laminar and turbulent regions and proposing a variable grid approach based on boundary layer thickness.

Contribution

It introduces a grid generation method based on boundary layer thickness from RANS simulations to improve laminar and turbulent region predictions in WMLES.

Findings

Refined near-wall grid captures laminar skin friction accurately.

Grid based on boundary layer thickness improves laminar region prediction.

Transition location can be delayed with boundary layer-informed grid adaptation.

Abstract

Wall-modeled large-eddy simulation (WMLES) is performed for flow over a wing with a focus on documenting grid resolution requirements to predict both the laminar and turbulent regions accurately. Flow over a spanwise extruded NACA0012 airfoil at 0-degree angle of attack and freestream chord-based Reynolds numbers of 3 million is simulated using an unstructured-grid finite-volume solver. An equilibrium wall function with the first off-wall grid point as the exchange location is used. Two scenarios are simulated wherein either the entire airfoil surface, or only a portion of it, is assumed turbulent based on linear stability calculations. For the latter scenario, the regular no-slip wall boundary condition is imposed in laminar regions. Using the same grid that was employed in the fully turbulent case, WMLES captured the skin friction in the turbulent region close to the RANS results. However, the skin friction is not well captured in the laminar region, due to far few grid points inside the boundary layer. When the near-wall grid was refined in the wall normal direction, the laminar region was captured accurately, but the turbulent region was not resolved well because the first off-wall point moved below the buffer-layer. To reconcile differing resolution requirements in the laminar and turbulent regions, a grid was generated based on the varying boundary layer thickness estimated from a precursor RANS simulation with a specified transition location. This grid produced improved results in the laminar region but resulted in a delayed transition location. Unsteady disturbances with the most amplified frequencies were introduced upstream of the neutral point. The transition process in the laminar region and the skin friction in the turbulent region were then captured satisfactorily.