In dynamical systems, it is advantageous to identify regions of flow which can exhibit maximal influence on nearby behaviour. Hyperbolic Lagrangian Coherent Structures have been introduced to obtain two-dimensional surfaces which maximise repulsion or attraction in three-dimensional dynamical systems with arbitrary time-dependence. However, the numerical method to compute them requires obtaining derivatives associated with the system, often performed through the approximation of divided differences, which can lead to significant numerical error and numerical noise. In this paper, we introduce a novel method for the numerical calculation of hyperbolic Lagrangian Coherent Structures using Differential Algebra called DA-LCS. As a form of automatic forward differentiation, it allows direct computation of the Taylor expansion of the flow, its derivatives, and the eigenvectors of the associated strain tensor, with all derivatives obtained algebraically and to machine precision. It does so without a priori information about the system, such as variational equations or explicit derivatives. We demonstrate that this can provide significant improvements in the accuracy of the Lagrangian Coherent Structures identified compared to finite-differencing methods in a series of test cases drawn from the literature. We also show how DA-LCS uncovers additional dynamical behaviour in a real-world example drawn from astrodynamics.
翻译:在动态系统中,确定流动区域对附近行为具有最大影响是有好处的。超双曲拉格朗加的coherent结构已经引入,以获得二维表面,在三维动态系统中,在任意依赖时间的任意依赖下,最大限度地反射或吸引;然而,计算它们的数字方法需要获得与系统相关的衍生物,通常通过分化差异的近似法进行,这可能导致重大数字错误和数字噪音。在本文中,我们采用了一种新颖的方法,用不同方位(DA-LCS)计算双曲拉格朗加的coherent结构的数字。作为一种自动前向差异形式,它允许直接计算Tayll的流动扩张、其衍生物以及相关压力的导体;所有衍生物的代数和机精度都得到。它这样做时没有事先了解系统的信息,例如变异方程式或明确的衍生物。我们证明,这可以大大改进Lagrangian Coherent结构的精确度。作为自动前向差异区分的一种形式,它允许直接计算出Tayl的流、衍生物的衍生物的衍生物、衍生物的衍生物,并直接计算出从DAR-DADA中从一个从真实的动态变化中绘制中绘制的动态动态动态的立法实例。