[論文レビュー] 傾き3のMSTACSはCo=0.5で255ステップ目に発散した — 圧縮差分とCourant限界
·6 min read
MSTACS式(18)の圧縮差分を1Dスラブで試しました。再分配なしではCo>0.44で発散します。原因は傾きkとCoの積です。
Kを消したら水が1.86倍圧縮された — 格子を細かくしても減らない界面の誤差
·8 min read
圧縮性VOFの圧縮率因子Kを等圧クロージャから直接導き、Kをゼロにした場合と符号を反転した場合に液体体積がどれだけ削られるかを1次元圧縮で再現します。
Δtを5倍にして1.9倍速くなり、10倍では変わらなかった — 陰的表面張力が開ける窓
·8 min read
毛管時間刻み制約を破っても上限は残る。Janodetら(2025)が密度比1000で測った安定限界と速度向上の最適点を追う。
[論文レビュー] 毛細管の足枷を外したらCFL 0.05が残った — VOF界面移流の本当の上限
·8 min read
表面張力を陰的に解いて毛細管時間刻み制約を破っても、色関数を運ぶ圧縮スキームがCFL 0.05で時間刻みを掴み直します。
液滴シミュレーションがΔtに足を取られる理由 — 表面張力を陰的に解く論文
·5 min read
格子をN=64から256へ細かくするとΔt_σは8分の1。コメントのωΔt_σは常に2.221で、陽的側は0.9倍を超えると発散する。
[論文レビュー] Poissonを使わずに二相流を解く — 一般圧力方程式とMSTACSの完全陽的結合
·7 min read
人工音速を5から15に上げると、安定な時間刻みΔtは1/3になり、体積誤差E_Gは6.57×10^{-8}まで下がります。
界面捕捉法の比較:VOF vs Level Set vs Diffuse Interface
·4 min read
圧縮性の界面接触では、αを保存形に解くと圧力が揺れ出す。5方程式モデルはαの移流だけを非保存のまま一つ残し、平衡を守る。