緩和係数を一つ変えたら収束解が動いた — Rhie–Chow面速度がアンダーリラクゼーションを覚えている場所
·6 min read
SIMPLEでα_uだけ変えると収束した速度場が変わります。Rhie–Chow面速度が前反復値を混ぜる方法が原因で、Majumdar補正で消えます。
粘性項だけ2次にしたら誤差が1.82倍に増えた — semi-implicit Gearスキームが対流項を外挿する理由
·7 min read
BDF2に対流項をN^nのまま組み合わせると1次に落ち、誤差はむしろ増えます。外挿2N^n−N^(n−1)の1行で2次が戻ることをコードで確かめます。
直径を1%読み違えると流量は4%ずれる — 壁の半セルとPoiseuilleの定数
·7 min read
円管の流量は直径の4乗に乗る。壁を半セル内側へ놓는格子も、1838年Poiseuilleのガラス管も、同じ指数4で読める。
格子が知らない壁 — IBMのデルタカーネルとmulti-direct forcing
·8 min read
一度の補正後もスリップは自由流の64%残る。Δs/h=3では隙間の漏れが0.287、0.25では条件数が10^11まで上がる。
圧力が市松模様に振動するとき — SIMPLE と Rhie-Chow の圧力-速度連成
·5 min read
48点の周期格子では、naive の市松模様は 0.80 のまま、Rhie-Chow は 3.1e-16 まで消えます。
圧力を方程式から消した男 — Chorinの射影法とフラクショナルステップ
·6 min read
予測直後の発散はおよそ1e-2、FFT射影の後は1e-13まで落ちる。その圧力Poissonが1ステップの費用の大半を食う。