1.利用cfd方法獲取一系列時(shí)刻的流場(chǎng)速度變化，時(shí)間間隔盡量較小

2. 用tecplot將fluent結(jié)果中的速度數(shù)據(jù)批量插值到均勻網(wǎng)格上，可以參考這個(gè)視頻

這里要注意的是就是：

• 如果要計(jì)算backwardFTLE，那么所插值的區(qū)域就必須包括你關(guān)心的區(qū)域以及其前部的區(qū)域

• 如果要計(jì)算backwardFTLE，那么所插值的區(qū)域就必須包括你關(guān)心的區(qū)域以及其后部的區(qū)域

這樣在做粒子推進(jìn)的時(shí)候，粒子就不會(huì)超出你所提供的數(shù)據(jù)區(qū)域

另外插值的網(wǎng)格也可以超出原網(wǎng)格的范圍，因?yàn)樵瓉?lái)的網(wǎng)格不是規(guī)則正方形，超出的部分速度值為0，粒子到那些區(qū)域就不再移動(dòng)。

3. 利用matlab編寫的ftle計(jì)算程序讀取速度數(shù)據(jù)，并進(jìn)行計(jì)算

設(shè)置相關(guān)參數(shù)：

Istart：開始讀取的首個(gè)文件名

N_SNAP:讀取文件個(gè)數(shù)

GxGy 讀取的數(shù)據(jù)網(wǎng)格數(shù)目

f_start要計(jì)算ftle的時(shí)刻的開始

f_end要計(jì)算ftle的時(shí)刻的結(jié)束

f_end要計(jì)算ftle的時(shí)刻的間隔

t_length粒子推進(jìn)的時(shí)間（50個(gè)delta）

tstep時(shí)間間隔

dt這個(gè)是整數(shù)，指每次粒子推進(jìn)幾個(gè)delta

這里注意的是，如果是backwardFTLE 那么f_start一定要大于t_length

如果是forwardFTLE，那么f_end+t_length<N_SNAP

下面還要設(shè)置所關(guān)心的區(qū)域，一般要小于你所提供的數(shù)據(jù)區(qū)間

4. 計(jì)算得到最終有限時(shí)間李雅普諾夫指數(shù)場(chǎng)結(jié)果