6次モーメントで15が9に落ちた — D2Q9を決めた三つのエルミート根
·7 min read
LBMコードの c_s²=1/3 と重み 1/36 がどこから来たのかを、ガウス・エルミート積分の次数まで遡って追いかけます。
BGKが発散する低粘性域を生き抜く方法 — MRT格子ボルツマンの衝突
·5 min read
ν=0.004でs_νは1.953になる。BGKはその値を高次モードにも乗せ、コードのMRTはs_e=1.64、s_q=1.7に固定する。
9本の矢に閉じ込めたMaxwell–Boltzmann — LBM平衡分布が多項式である理由
·6 min read
エルミート級数を2次で切ると、九つの速度を持つD2Q9の平衡多項式になり、同じ切断が低速のマッハ0.3の天井を残します。