Computing Lagrangian coherent structures from their variational theory

Using the recently developed variational theory of hyperbolic Lagrangian coherent structures (LCSs), we introduce a computational approach that renders attracting and repelling LCSs as smooth, parametrized curves in two-dimensional flows. The curves are obtained as trajectories of an autonomous ordinary differential equation for the tensor lines of the Cauchy-Green strain tensor. This approach eliminates false positives and negatives in LCS detection by separating true exponential stretching from shear in a frame-independent fashion. Having an explicitly parametrized form for hyperbolic LCSs also allows for their further in-depth analysis and accurate advection as material lines. We illustrate these results on a kinematic model flow and on a direct numerical simulation of two-dimensional turbulence.