先说下我们qq群:格子玻尔兹曼救星:293267908。每个人都是自己的救星。
铁打的群流水的学生,总有人问的问题,还是得记录一下。
郭照立:
形式上完全相同,是否编程完全相同?
实际上half_way是碰撞之后立马反弹,也就是写入反弹的时候坐标不变。在边界上,直接令 f(k,x,y) = fpost(inv(k),x,y); 其中fpost(inv(k),x,y)是上一步走向固体的分布函数,f(k,x,y) 则是当前由固体走向流体的分布函数,也就是一步就弹回来了,你的xy坐标没有体现固体的位置,反弹前后xy位置都是在流体中的。标准反弹是在流动过程中,分布函数碰到壁面时,壁面处的反弹用上一步的分布函数反弹。也就是这个反弹是在墙壁上进行的,反弹时你的xy坐标就是固体边界的坐标。所以需要两个时间步长。
也就是,半步长反弹过程只用一个时间步,标准反弹格式要用两个时间步。
(持续更新中...)