Particle Level Set + Lattice Boltzmann Method

  LBM方法是对流体的另一种模拟方法，最大的好处就是其计算局部性，通过微观的密度分布函数的传递来表现宏观的密度与速度的变化。但LBM的Free Surface模拟一直没有成为大家的关注点，真正完成这一研究的可能也只有Nils Thuerey，此后出现的LBM Free Surface Flow绝大多数都是采用了Nils的方法，但从未有人对该方法进行进一步的研究。我们在完成了这一方法的重现之后，发现，其最大的缺点就是在内存占用十分严重（由于没有试过其他方法，不知是不是流体模拟算法的通病），32位的机器最多只能模拟150^3左右的空间，而最近在64位机器上尝试了300x150x300的模拟空间，居然达到了10G的内存，当然效果的确也好了很多，见下图：



此后在其上，我们完成了二相流，温度，固液耦合等方面的扩展。开始觉得那些讨论对某一方法进行扩展，而不是深化的论文总是会差那么一些火候，发不了顶级的会议或者期刊。因此考虑着对一些方法进行算法上的改进，或者完成某些方法的混合建模，已达到相互扬长避短的效果，这也是发好文章的一种经典模式。

因此，我考虑着将Particle Level set 引入进来，Nils也曾将LBM和LS相结合，并与其方法做了对比，个人感觉LS+LBM的效果还是不错的，如果用的是PLS，那么效果必将更上一层楼。Particle Level Set 在网上有一份源代码，以及配套的pdf文档，是非常好的学习资料。主要用做界面的追踪，其原理在文档中有着很详细的介绍，在被引入流体模拟之后，成为了一种经典的方法，与NS求解器组合，能够模拟很不错的流体现象。

目前经过试验，完成了一个二维的PLS+LBM耦合，见下图：



其中的速度外插尚未完成，因此效果上还很有问题...但也初具雏形了~