欢迎光临散文网 会员登陆 & 注册

有限时间李亚普诺夫指数识别LCS的灵感来源

2023-02-15 13:50 作者:奔影临存森顶犇巽棍轮  | 我要投稿

考虑一个动力学系统

一段时间内点的映射为

映射特性

假设不动点为xc,则v(xc)=0, 随时间增加轨迹渐进于不动点的为稳定流形,随时间减小到渐进于不动点为不稳定流形。

Example:钟摆

钟摆示意图


动力学方程
状态空间形式
单摆的相空间

分析这个系统,具有(0,0)(π,0) (-π,0)等不动点。

(π,0)和(-π,0)都是双曲不动点,他们的稳定和不稳定流形形成了一个分隔线,图中的蓝色曲线。

在这个蓝色曲线以内出发的轨迹,最终单摆会来回摆动,轨迹形成周期振荡,如下图:

在蓝色轨迹以外出发的,不再来回振荡,而是不断旋转,相位不断增大了

(π,0)是一个双曲不动点,它的稳定流形和不稳定流形起到了分界线的作用,将相平面图以(π,0)为中心扩展一下:

在稳定流形两侧布置点,随时间它们会不断远离同理,在不稳定流形附近布置点,如果反向时间积分的话,它们也会不断远离

稳定流形两侧的粒子轨迹随时间增加不断分开

所以,在流体问题中,可以通过定量确定分离度来确定稳定流形和不稳定流形。

FTLE计算方法

考虑一个t0-t0+T 的流动映射

t0时刻x附近有一个点y

T0+T时刻两者分离的距离

它的幅值:L2范数

分离距离范数
Cauchy-Green deformation tensor


Cauchy-Green deformation tensor的最大特征向量和特征值对应了最大的拉伸方向和该方向长度,为

最大特征值

定义有限时间李雅普诺夫指数

有限时间李雅普诺夫指数

那么有限时间李雅普诺夫指数最大的位置,就是在那里附近放两个粒子,一段时间会分离的最大的地方,那么该位置就可能是稳定流形存在的位置。

流体力学FTLE具体计算方法

步骤1:在t0时刻在网格点上布置粒子,获得粒子的位置坐标,并根据流场计算数据插值出各个粒子所在位置的速度。

步骤2:利用时间推进方法(欧拉法,四阶龙格库塔法等),获得下一时刻(t0+Δt)这些粒子的位置。通过计算获得(t0+Δt)时刻的流场信息,获得该时刻粒子的速度。按照该方法,每次前进Δt,直到完成时间T。那么t0时刻布撒的粒子,在(t0+T)时刻的位置可以获得

步骤3:根据以上数据和差分法,计算出每一个网格点的变形梯度:

步骤4:计算柯西格林张量,并计算其特征值,最大特征值则为FTLE。

有限时间李亚普诺夫指数识别LCS的灵感来源的评论 (共 条)

分享到微博请遵守国家法律