Published January 2026 | Published
Journal Article Open

Populating the wall layer, one eddy at a time: Resolvent analysis for Wall-Modelled LES

  • 1. ROR icon Queen's University
  • 2. ROR icon California Institute of Technology
  • 3. ROR icon Stanford University

Abstract

Computational cost precludes direct numerical simulation or wall-resolved large-eddy simulations of non-equilibrium, wall-bounded turbulent flows in realistic conditions. Wall-modelled large-eddy simulations (WMLES) and hybrid RANS/LES methods can be used to analyse these flows at much decreased cost, but require modelling of the near-wall layer and, in particular, a means to address the deficit of turbulent activity, or eddies, in the vicinity of the interface between the outer flow and the wall model. We report a computational framework to populate the wall region with synthetic but realistic eddies and reflect their integrated effect on the flow in the inner layer. Two means of generating spatio-temporal representations for the synthetic eddies are investigated: low-order, resolvent-based representations of the wall layer and a coarse-grained, data-driven spectral proper orthogonal decomposition (SPOD) model, both generated in turbulent channel flow at a friction Reynolds number, Reτ = 1000. The eddy-augmented WMLES models are then tested in Reτ = 5000 and 20,000 channels and compared with experimental and numerical data. The inherent scaling of the resolvent operator can be used to scale the resolvent model to higher Reynolds numbers (and potentially populate new, self-similar eddies as the wall layer grows in inner units), while the SPOD model is energetically optimal for reconstruction of the flow at Reynolds numbers close to that where it is obtained, but degrades as the Reynolds number is increased. The results show that the effect of the introduction of synthetic eddies is twofold: first, a direct contribution to the stress due to the presence of the synthetic eddies and, second, an improved prediction of the normal Reynolds stresses in the inner layer due to an accompanying, coupled reduction in the time- and length-scales of the variation of the URANS-like velocity in the inner layer. Implications and extensions of the method for more complex flows, for example external boundary layers with pressure gradient and separation, are briefly discussed.

Copyright and License

© 2025 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Acknowledgement

This research was supported by the Air Force Office of Scientific Research, United States under award number FA9550-21-1-0443 P00004, monitored by Dr. Gregg Abate. Computational resources were provided by the Digital Research Alliance of Canada.

Files

1-s2.0-S0142727X25003248-main.pdf
Files (3.0 MB)
Name Size Download all
md5:47c42352bb29fbdb578565d1fad47df7
3.0 MB Preview Download

Additional details

Created:
September 26, 2025
Modified:
September 26, 2025