Skip to content
cfd-lab:~/zh/tags/Diffuse-Interface● online
$ cd ..
grep "#Diffuse-Interface"// 13 篇

#Diffuse-Interface

~58 分钟阅读

界面只有两格宽时,表面张力少了39% — Korteweg应力的各向异性

·4 min read

平直气液界面上只有切向压力会下凹。下凹的面积就是表面张力,在dz=1的格子上,两格宽的界面会丢掉39%。

[论文解读] 把 β 调到 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,把那层厚度重新锁回常数。

方程少了,反而不炸了 — 四方程模型与界面平衡条件

·4 min read

压力本来处处均匀,可守恒型输运一过气体界面就把P晃起来。四方程模型把Γ改成非守恒搬运,这阵振荡也就落到了双精度的舍入误差。

给每一相各自的方程 — Baer–Nunziato 七方程模型

·4 min read

印出的324 bar,是刚性气体两个大项相消后剩下的1.2%。松弛第一步就在推这个差额,液体几乎只卸压力,气体几乎不动。

水面上的呼喊为何传不进水下 — 跨越扩散界面的声传播

·4 min read

空气–水阻抗比约3600,大约99.9%的声音被水面直接弹回。临界角13.4度与弱可压缩插值,都在一维声管里重算过一遍。

[论文综述] 用两个压力先解,再合成一个 — 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

在可压缩的接触界面上,把体积分数α按守恒型来更新,压力就会晃动。五方程模型只把α的输运留成非守恒形式,以此守住压力平衡。