Reference map technique for Lagrangian exploration of coherent structures

Abstract

We explore the application of the reference map technique, originally developed for Eulerian simulation of solid mechanics, in Lagrangian kinematics of turbulent flows. Unlike traditional methods based on explicit particle tracking, the reference map facilitates the calculation of flow maps and gradients without the need for particles. This is achieved through an Eulerian update of the reference map, which records the take-off positions of fluid particles. This approach is found to be mathematically equivalent to the work of Leung (J. Comput. Phys., vol. 230, issue 9, 2011, pp. 3500-3524), who computed the flow map of simple two-dimensional flows using an Eulerian approach. We discuss important modifications necessary for its first application to complex three-dimensional turbulent flows, including the conservative, low-dissipation update of the flow map and the treatment of periodic boundary conditions. We first demonstrate the accuracy of finite-time Lyapunov exponent (FTLE) calculations based on the reference map against the standard particle-based approach in a two-dimensional Taylor-Green vortex. Then we apply it to turbulent channel flow at, where Lagrangian coherent structures identified as ridges of the backward-time FTLE are found to bound vortical regions of flow, consistent with Eulerian coherent structures from the -criterion. The reference map also proves suitable for material surface tracking despite not explicitly tracking particles. This capability can provide valuable insights into the Lagrangian landscape of turbulent momentum transport, complementing Eulerian velocity field analysis. The evolution of initially wall-normal material surfaces in the viscous sublayer, buffer layer and log layer sheds light on the Reynolds stress-generating events from a Lagrangian perspective.