τ=0.6时流速只有六分之一 — Guo forcing漏掉半步速度F/2的后果
·4 min read
保留Guo源项却在速度中去掉F/2,有效体积力会缩小为(1−1/2τ)倍。Kolmogorov流中τ=0.6时振幅只有0.166倍。
只改了 τ，静止气体的密度就变了 — Shan–Chen 速度平移 forcing 的 (τ−½)² 应力
·4 min read
给静止的 D1Q3 气柱施加同样的力，只有 Shan–Chen forcing 的密度随 τ 改变。平衡速度的平移留下了一个 (τ−½)²F²/ρ 应力。
老实补上三阶矩后误差反而涨了2.8倍 — 力项真正需要守住的几行
·5 min read
在D2Q9上把力项三阶矩补成动理论的值，误差从0.053%升到0.150%。格子锁住张量一半，该项只在Burnett阶才进动量方程。
把 τ 直接代进去，粘度大了六倍 —— LBM 离散化留下的三处 Δt/2
·4 min read
一次梯形积分和一行变量替换，在粘度、速度、应力三处各留下一个 Δt/2。用标量 ODE 和 D2Q9 把这三处都测了一遍。
把压力直接代入后静止的界面开始颤动 — 非理想LBM forcing的两种形式
·4 min read
静止范德瓦尔斯液滴用格子玻尔兹曼求解时，力项代入压力p还是化学势mu，连续下靠吉布斯关系重合，一到格子界面上便立刻分开。
LBM 里那一半的力去哪了 — 四种 forcing 格式与 1−1/(2τ)
·7 min read
定常Poiseuille通道分不出五种forcing格式，误差到第四位都一样。τ=0.6时只开半格修正，幅值变为1.83倍。