计算止回阀流场的二维FTLE(拉格朗日拟序结构)
2023-03-04 17:42 作者:奔影临存森顶犇巽棍轮 | 我要投稿
1.利用cfd方法获取一系列时刻的流场速度变化,时间间隔尽量较小

2. 用tecplot将fluent结果中的速度数据批量插值到均匀网格上,可以参考这个视频

这里要注意的是就是:
如果要计算backwardFTLE,那么所插值的区域就必须包括你关心的区域以及其前部的区域
如果要计算backwardFTLE,那么所插值的区域就必须包括你关心的区域以及其后部的区域
这样在做粒子推进的时候,粒子就不会超出你所提供的数据区域
另外插值的网格也可以超出原网格的范围,因为原来的网格不是规则正方形,超出的部分速度值为0,粒子到那些区域就不再移动。


3. 利用matlab编写的ftle计算程序读取速度数据,并进行计算
设置相关参数:

Istart:开始读取的首个文件名
N_SNAP:读取文件个数
GxGy 读取的数据网格数目
f_start要计算ftle的时刻的开始
f_end要计算ftle的时刻的结束
f_end要计算ftle的时刻的间隔
t_length粒子推进的时间(50个delta)
tstep时间间隔
dt这个是整数,指每次粒子推进几个delta
这里注意的是,如果是backwardFTLE 那么f_start一定要大于t_length
如果是forwardFTLE,那么f_end+t_length<N_SNAP
下面还要设置所关心的区域,一般要小于你所提供的数据区间
4. 计算得到最终有限时间李雅普诺夫指数场结果
