Eulerian Methods for Visualizing Continuous Dynamical Systems using Lyapunov Exponents

TL;DR

This paper introduces an improved Eulerian numerical method for constructing forward flow maps in continuous dynamical systems, enabling efficient computation of Lyapunov exponents and coherent structures.

Contribution

The authors develop a forward-in-time Eulerian approach that simplifies implementation and enhances efficiency in computing Lyapunov exponents and flow structures, especially for ISLE calculations.

Findings

Efficient computation of FTLE, FSLE, and ISLE using the new Eulerian method.

The method allows on-the-fly determination of flow maps, improving over previous approaches.

Theoretical link established between FTLE and ISLE fields.

Abstract

We propose a new Eulerian numerical approach for constructing the forward flow maps in continuous dynamical systems. The new algorithm improves the original formulation developed in [23, 24] so that the associated partial differential equations (PDEs) are solved forward in time and, therefore, the forward flow map can now be determined on the fly. Due to the simplicity in the implementations, we are now able to efficiently compute the unstable coherent structures in the flow based on quantities like the finite time Lyapunov exponent (FTLE), the finite size Lyapunov exponent (FSLE) and also a related infinitesimal size Lyapunov exponent (ISLE). When applied to the ISLE computations, the Eulerian method is in particularly computationally efficient. For each separation factor r in the definition of the ISLE, typical Lagrangian methods require to shoot and monitor an individual set of ray trajectories. If the scale factor in the definition is changed, these methods have to restart the whole computations all over again. The proposed Eulerian method, however, requires to extract only an isosurface of a volumetric data for an individual value of r which can be easily done using any well-developed efficient interpolation method or simply an isosurface extraction algorithm. Moreover, we provide a theoretical link between the FTLE and the ISLE fields which explains the similarity in these solutions observed in various applications.