Predicting transport by Lagrangian coherent structures with a high-order method

Recent developments in identifying Lagrangian coherent structures from finite-time velocity data have provided a theoretical basis for understanding chaotic transport in general flows with aperiodic dependence on time. As these theoretical developments are extended and applied to more complex flows, an accurate and general numerical method for computing these structures is needed to exploit these ideas for engineering applications. We present an unstructured high-order hp/spectral-element method for solving the two-dimensional compressible form of the Navier-Stokes equations. A corresponding high-order particle tracking method is also developed for extracting the Lagrangian coherent structures from the numerically computed velocity fields. Two different techniques are used; the first computes the direct Lyapunov exponent from an unstructured initial particle distribution, providing easier resolution of structures located close to physical boundaries, whereas the second advects a small material line initialized close to a Lagrangian saddle point to delineate these structures. We demonstrate our algorithm on simulations of a bluff-body flow at a Reynolds number of Re = 150 and a Mach number of M = 0.2 with and without flow forcing. We show that, in the unforced flow, periodic vortex shedding is predicted by our numerical simulations that is in stark contrast to the aperiodic flow field in the case with forcing. An analysis of the Lagrangian structures reveals a transport barrier that inhibits cross-wake transport in the unforced flow. The transport barrier is broken with forcing, producing enhanced transport properties by chaotic advection and consequently improved mixing of advected scalars within the wake.