界面幅が2セルなら表面張力の39%が消える — Korteweg応力の異方性
平らな気液界面では接線方向の圧力だけが下がります。その面積が表面張力で、dz=1の格子では幅2セルで39%が失われます。
[論文レビュー] βを2にしたら界面は14セルで止まり、正弦波は1次風上より悪化した — Tallois の圧縮型リミタ
体積分率にβ=2の圧縮型リミタを使うと界面は5周後も14セル幅のままですが、同じリミタは正弦波で1次風上より大きな誤差を出します。
σ₁₂に1.00を入れて0.78が返ってきた — 3成分Allen–Cahnで存在しない相が界面に漏れる場所
3成分Allen–CahnでΣc=1を平均投影で合わせると、存在しない相が界面に漏れて表面張力が21.9%削られる。3重罰則項がその穴を塞ぐ。
[論文レビュー] 同じ点の色関数が経路によって0.544ずれた — カール拘束条件とスタッガード格子
400ステップ後、スタッガード格子のカールは1e-12に留まり、コロケート配置では同じ点の色関数が経路で0.544ずれました。
何も起きないはずの問題で圧力が21%跳ねた — Abgrall条件と非保存項
圧力・速度が一様な流れが一様に保たれない原因はフラックスではなく非保存項の離散化です。Re–Abgrall論文のH_u演算子を同じ格子で測り直しました。
速度を上げたら拡散が27%減った — LBM移流拡散モデルの余分なフラックス
LBMの移流拡散モデルはτでDを合わせてもD(1−u²/cs²)でしか広がりません。Chapman–Enskogが落とす余分なフラックスと、それを消す源項。
にじんだ界面を毎ステップ元に戻す — 反拡散シャープニング(IST)を実装する
200ステップ後、数値拡散でヘリウム気泡の縁が6セルから20セルへ膨らむ。IST後処理はε/aだけでその厚さを定数に戻す。
方程式を減らしたら、むしろ壊れなくなった — 4方程式モデルと界面平衡条件
圧力が一様でも保存型の輸送は界面でPを揺らします。4方程式モデルはΓを非保存の形で運び、その振動を丸め誤差まで落とします。
各相にそれぞれの方程式を与える — Baer–Nunziato 7方程式モデル
印字された324 barは、剛性気体の二つの大きな項が1.2%だけ残して打ち消し合った差です。緩和の初手はその差を押します。
水上で叫んだ声が水中に届かない理由 — 拡散界面を越える音響伝播
空気–水のインピーダンス比およそ3600で反射率は99.9%に張り付く。臨界角13.4度と弱圧縮性補間を1D管で再現する。
[論文レビュー] 圧力を二つに分けて解き、一つに戻す — Saurel(2009) 圧縮性多相流の緩和法
5方程式拡散界面モデルの罠を、6方程式の圧力非平衡 + stiff緩和で迂回したSaurel・Petitpas・Berry(2009)のレビュー。
5-Equation Model 実装ガイド：Python で作る多相流ソルバー
体積分率0.5の混合音速は約24 m/sまで落ちます。同じ初期状態なら空気は約374 m/s、水は約1625 m/sです。
界面捕捉法の比較：VOF vs Level Set vs Diffuse Interface
圧縮性の界面接触では、αを保存形に解くと圧力が揺れ出す。5方程式モデルはαの移流だけを非保存のまま一つ残し、平衡を守る。