逆反応定数を Kp で割ったら 2500 K で AlO が 99.6% 解離した — 平衡定数に潜む (p°/RT)^Δn
·6 min read
kr = kf/K に Kp を入れると、Δn≠0 の反応だけ平衡が (p°/RT)^Δn 倍ずれます。AlO 解離反応器で SI・CGS の誤りを再現します
セル番号を混ぜただけでLU-SGSが2.8倍遅くなった — RCM再順序付けと開始セルが決めるスイープ方向
·8 min read
非構造格子のセル番号がLU-SGSの収束を変える理由と、RCM再順序付け・開始セル選びの効果を576セルの実験で確かめます
距離計算を574倍減らしたら接触104組が消えた — DEMのセクタ寸法を決めるのは1個の粒子
·7 min read
粒子2400個でセクタを1/4にすると距離計算は574倍減り、接触は104組消えます。安全条件は shell·s ≥ d_max だけです。
refine を stream の一行下へ動かしたら誤差が7.2倍になった — LBM AMR の再帰的時間進行順序
·6 min read
同じ格子で refine を stream の直後に移すと誤差が7.2倍になる。0.657%が4.704%へ増え、ずれは継ぎ目に残る。
クアッドツリー AMR — 格子を必要な場所だけ育てる
·7 min read
レベル7の一様格子は128セルだ。ガウスパルス一つなら葉は41枚、32%まで落ちる。四分木と傾き閾値τが、どこを割るかを決める。
Lax–Wendroffが段差で震えた — Flux Limiterで2次TVD移流スキームを作る
·5 min read
CFL 0.5・500段でLax–Wendroffのオーバーシュートは0.181。minmodとvan Leer、MCは0。
CFD線形ソルバーの選び方：JacobiからBiCGSTABまでの収束速度比較
·5 min read
N=60の1次元ポアソンでは、Jacobiが2000回後も残差2.11e-2、CGは60回で1.87e-13に届きます。
壁距離を O(N log N) で解く方法：Eikonal 方程式と Fast Marching
·7 min read
1000万セルと50万壁面なら距離計算は5兆回になる。Fast Marchingは壁距離をO(N log N)のEikonal解として出す。