Eulerian Methods for Visualizing Continuous Dynamical Systems using Lyapunov Exponents (1603.06446v2)

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.