我的指导老师建议我看一下格子玻尔兹曼方法作为科学计算的介绍,并可能是本科荣誉论文的主题。我编写了一些 Java 代码(有点非传统,但我现在只知道 Java 和 Python),并且我有一个简单的 2D 模拟。在测试 2D 通道流时,我发现该行为有很多问题。
节点处的速度不稳定。一般来说,它们都接近于零。通过第 10,000 个时间步,它们的数量级为
尽管人们会期望抛物线速度分布,但速度的 x 分量的大小沿 x 轴在正和负之间切换。

我使用的参考资料清楚地解释了一般算法,但我还没有找到关于无滑移边界条件实现的清晰描述。这就是我认为我的问题所在。我正在尝试实现并网边界条件(第 3-4 页)并模拟文章中描述的泊肃叶流。
无滑移边界条件究竟是如何初始化的?对,现在我假设您根据初始密度和速度计算平衡函数。
如果这不是问题,我的代码可能会有更大的问题。如果有人好奇,源代码位于此处。请记住,这对我来说是实现 LBM 的最直观的方法,但不一定是最有效的。