[论文解读] 把 β 调到 2,界面停在 14 个网格,正弦波却比一阶迎风还差 — Tallois 的压缩型限制器
·4 min read
对体积分数使用 β=2 压缩型限制器,界面在五个周期后仍只有 14 个网格宽;同一限制器作用于正弦波时,误差却大于一阶迎风格式。
σ₁₂填了1.00,拿回0.78 — 三组分Allen–Cahn中缺席相渗入界面的位置
·5 min read
三组分Allen–Cahn中用平均投影满足Σc=1,缺席相会渗入界面,使实际表面张力被削去21.9%。三重罚项补上这个缺口。
[论文评述] 同一点的色函数沿两条路径相差0.544 — 旋度约束与交错网格
·6 min read
同一点的色函数沿两条路径差了0.544。交错网格把旋度按在约1e-12,同位布置在同一剪切层上跑400步后,旋转噪声已经盖过界面。
在本该什么都不发生的算例里,压力跳了 21% — Abgrall 判据与非守恒项
·6 min read
均匀压力与速度的流动之所以保持不住均匀,根源不在通量而在非守恒项的离散。本文在同一网格上复现了 Re–Abgrall 论文的 H_u 算子。
速度调高后扩散少了27% — LBM对流扩散模型中多出来的通量
·5 min read
LBM对流扩散模型即便用τ标定了D,实际仍以D(1−u²/cs²)扩散。Chapman–Enskog漏下的多余通量,以及消除它的源项。
每一步都把糊掉的界面拉回来 — 亲手实现反扩散锐化(IST)
·5 min read
激波气泡算完200步后,数值扩散让氦气界面从6个网格糊开到整整20个。IST后处理只靠比值ε/a,把那层厚度重新锁回常数。
给每一相各自的方程 — Baer–Nunziato 七方程模型
·4 min read
印出的324 bar,是刚性气体两个大项相消后剩下的1.2%。松弛第一步就在推这个差额,液体几乎只卸压力,气体几乎不动。
[论文综述] 用两个压力先解,再合成一个 — Saurel(2009) 可压缩多相流松弛法
·4 min read
用六方程压力非平衡 + 刚性松弛绕开五方程扩散界面模型陷阱的 Saurel·Petitpas·Berry(2009) 综述。
5-Equation Model 实现指南:用 Python 编写多相流求解器
·3 min read
体积分数取 0.5 时,混合声速大约只有 24 m/s。同一初始状态下,空气约 374 m/s,水约 1625 m/s。
界面捕捉技术对比:VOF vs Level Set vs Diffuse Interface
·4 min read
在可压缩的接触界面上,把体积分数α按守恒型来更新,压力就会晃动。五方程模型只把α的输运留成非守恒形式,以此守住压力平衡。