τ=0.6で流速が6分の1になった — Guo forcingで半ステップ速度F/2を忘れると
·5 min read
Guoソース項を入れたまま速度からF/2を外すと、実効的な力が(1−1/2τ)倍に縮みます。Kolmogorov流ではτ=0.6で0.166倍です。
τだけ変えたら静止した気体の密度が変わった — Shan–Chen速度シフトforcingの(τ−½)²応力
·6 min read
静止したD1Q3気体柱に同じ力をかけると、Shan–Chen forcingだけが τ で密度を変える。平衡速度のシフトが (τ−½)²F²/ρ の応力を残します。
3次モーメントを正直に埋めたら誤差が2.8倍になった — 力項が守るべき行
·7 min read
LBMの力項の3次モーメントを運動論の値で埋めると誤差が逆に増える。格子が固定する成分と、その項が届く次数を一緒に読む。
τをそのまま入れたら粘性が6倍になった — LBMの離散化が残したΔt/2の三箇所
·6 min read
台形則による積分一回と変数変換一行が、粘性・速度・応力の三箇所にΔt/2を残します。スカラーODEとD2Q9で三箇所すべてを測りました。
圧力をそのまま入れたら静止した界面が震えた — 非理想LBM forcingの二つの形
·5 min read
ファンデルワールス流体を格子ボルツマンで解くとき、力の項に圧力pを入れるのと化学ポテンシャルmuを入れるのは連続では同じでも離散で分かれる。
LBM で力の半分はどこへ消えたのか — forcing スキーム四種と 1−1/(2τ)
·9 min read
定常Poiseuilleでは五つのforcing方式の誤差が一致する。τ=0.6で半シフトだけを入れると振幅は1.83倍になる。