公开(公告)号：CN106021828A

公开(公告)日：2016-10-12

申请号：CN201610554629.X

申请日：2016-07-15

Applicant: 华中科技大学

Inventor： 武琛 ,  施保昌 ,  柴振华 ,  赵勇 ,  黄昌盛 ,  汪垒 ,  唐冲

IPC: G06F17/50

Abstract: 本发明公开了一种基于格子‑玻尔兹曼模型的流体模拟方法。所述模拟方法包括，将多孔介质的图像网格化，用fi(x,y,t)表示网格点I(x,y)处，运动速度为ci的粒子所对应的粒子分布；判断粒子运动方向ci是否朝向固壁边界，是则令粒子分布fi(x,y,t)执行反向函数，否则对fi(x,y,t)执行格子‑玻尔兹曼模型的碰撞函数，然后根据演化后的粒子分布fi(x,y,t)获得流体密度ρ’(x,y)和流体速度u’(x,y)；直至满足演化结束条件。本发明通过将多孔介质的图像网格化，并将网格中的流体离散化为运动的粒子，从而根据粒子分布fi(x,y,t)获得流体密度和速度，提高了模拟的运行效率。

公开(公告)号：CN106021828B

公开(公告)日：2017-09-08

申请号：CN201610554629.X

申请日：2016-07-15

Applicant: 华中科技大学

Inventor： 武琛 ,  施保昌 ,  柴振华 ,  赵勇 ,  黄昌盛 ,  汪垒 ,  唐冲

IPC: G06F17/50

Abstract: 本发明公开了一种基于格子‑玻尔兹曼模型的流体模拟方法。所述模拟方法包括，将多孔介质的图像网格化，用fi(x,y,t)表示网格点I(x,y)处，运动速度为ci的粒子所对应的粒子分布；判断粒子运动方向ci是否朝向固壁边界，是则令粒子分布fi(x,y,t)执行反向函数，否则对fi(x,y,t)执行格子‑玻尔兹曼模型的碰撞函数，然后根据演化后的粒子分布fi(x,y,t)获得流体密度ρ’(x,y)和流体速度u’(x,y)；直至满足演化结束条件。本发明通过将多孔介质的图像网格化，并将网格中的流体离散化为运动的粒子，从而根据粒子分布fi(x,y,t)获得流体密度和速度，提高了模拟的运行效率。