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

#Multiphase

~130 分钟阅读

只改了 τ,静止气体的密度就变了 — Shan–Chen 速度平移 forcing 的 (τ−½)² 应力

·4 min read

给静止的 D1Q3 气柱施加同样的力,只有 Shan–Chen forcing 的密度随 τ 改变。平衡速度的平移留下了一个 (τ−½)²F²/ρ 应力。

把 v 项和 v² 项相加,两端都对,中间却偏了 42% — Prony 两项阻力律与 Re=60

·5 min read

Bossut 测得阻力正比于 v²,Du Buat 测得正比于 v。Prony 的和 a·v²+b·v 在 Re=60 处交接主导权,并在附近把终端速度高估达 42%。

1 μm 铝颗粒的升温比连续介质公式慢 4.6 倍 — Knudsen 过渡区的 Sherman 调和插值

·6 min read

1500 K、1 atm下1 μm铝颗粒的Kn为0.67。连续介质Nu=2把加热时间低估为1/4.58。Sherman调和平均把它补回。

σ₁₂填了1.00,拿回0.78 — 三组分Allen–Cahn中缺席相渗入界面的位置

·5 min read

三组分Allen–Cahn中用平均投影满足Σc=1,缺席相会渗入界面,使实际表面张力被削去21.9%。三重罚项补上这个缺口。

迁移率放大10倍,界面厚度却没有动 — 守恒型Allen–Cahn的反扩散项

·5 min read

实测反扩散项精确抵消扩散,M从0.02增到0.20时界面厚度仍锁在W附近。该项在格子上进入LBM平衡分布的位置一并写清。

阻力系数削掉 9%,液滴却慢了 35% — 变形与蒸发修正该加在哪里

·5 min read

把球形、TAB 变形、Renksizbulut–Yuen 蒸发修正这三种阻力系数模型用在同一个液滴上,以贯穿距离进行对比。

在本该什么都不发生的算例里,压力跳了 21% — Abgrall 判据与非守恒项

·6 min read

均匀压力与速度的流动之所以保持不住均匀,根源不在通量而在非守恒项的离散。本文在同一网格上复现了 Re–Abgrall 论文的 H_u 算子。

把Δt放大5倍快了1.9倍,放大10倍却原地踏步 — 隐式表面张力打开的窗口

·6 min read

突破毛细约束后上限仍在。密度比1000下,Δt放大5倍快了1.9倍,10倍则原地踏步。Janodet等(2025)测了这个窗。

网格减半后发散来得快了一倍 — 六方程 two-fluid 模型的复数特征值

·5 min read

单压力六方程 two-fluid 模型只要两相有滑移,特征值就会变成复数。界面压力项系数越过1的瞬间,这对慢波回到实轴。

[论文评述] 把对流和压力彻底拆开 — Baer–Nunziato flux 的 TV 分裂

·7 min read

把Baer–Nunziato的flux拆成对流与压力后,声速只留在压力系统里;Noh滞止区在第三步就会把压力打到−0.12。

[论文评述] 液滴几乎把声波全部弹回 — 声阻抗与sharp界面两相求解器

·5 min read

空气-水界面把声能的99.889%反射回去。气相体积分数只要1%,Wood声速就从1500 m/s掉到120.5 m/s。

衣服为什么在5℃也会晾干 — 扩散控制蒸发与热控制沸腾

·5 min read

气温只有5℃时界面饱和质量分数并不为零,所以衣服照样晾干。到了100℃,沸腾速率改由热平衡与Lee模型的松弛系数r来决定。

每一步都把糊掉的界面拉回来 — 亲手实现反扩散锐化(IST)

·5 min read

激波气泡算完200步后,数值扩散让氦气界面从6个网格糊开到整整20个。IST后处理只靠比值ε/a,把那层厚度重新锁回常数。

静止的液滴为何会自己流动 —— 寄生电流与 well-balanced 表面张力

·4 min read

静止液滴的曲率只要偏3%,残差最大就跳到4.8e5;曲率精确时留在4.1e-10。动能大约差出24个数量级,这就是寄生电流。

雨滴为什么不是泪滴形 — Weber数与Taylor类比破碎(TAB)

·4 min read

TAB把液滴当成弹簧振子。变形y一旦超过1就判定破碎;无阻尼时峰值等于2·We/24,临界Weber数因此正好落在12附近。

液滴模拟为何被 Δt 绊住 — 把表面张力隐式求解的论文

·4 min read

网格从N=64收到256,毛细步长Δt_σ缩成八分之一。ωΔt_σ在每个网格上都是2.221,显式积分一过约0.9倍就发散。

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

·4 min read

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

[论文评论] 单级即达二阶精度 — 面向Kapila多相流的GRP求解器

·4 min read

在N=320的光滑正弦波上,单级GRP阶为2.03、L1误差1.55e-5,同网格的Godunov误差是2.31e-2。

水滴为何是圆的,杯沿的水为何会爬升 — 表面张力与接触角

·4 min read

输出里,半径 5 µm 液滴的压力跳变是 28800 Pa,0.1 mm 细管在接触角 30 度时上升 127.1 mm。

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

·4 min read

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

[论文综述] 不解 Poisson 也能算两相流 —— 通用压力方程与 MSTACS 的完全显式耦合

·5 min read

把人工声速 c_s 从 5 提到 15 以后,稳定时间步大约缩到原来的 1/3,体积分数误差 E_G 则降到 6.57×10^{-8}。

20℃的水也会沸腾 — 空化的物理

·4 min read

20℃的饱和蒸气压只有2.34 kPa,60℃会升到19.9 kPa。喉部静压跌破该值时,空化数σ转为负,水便开始沸腾。

SPH 粒子法求解流体 — 从核函数到人工粘性

·5 min read

γ=7的Tait方程会把约1%的密度误差放大成明显的压力摆动。压力梯度若不做粒子对对称化,这份压力就以动量泄漏的形式留下来。

压缩性多相流 Rhie-Chow 插值法详解

·9 min read

在水-气界面上,密度比取 1000 时,算术平均密度大约为 500.5,调和平均则降到 2.0。滤波只作用在驱动压力上。

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

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

可压缩多相流CFD入门:为什么难,为什么重要

·3 min read

在界面上对γ取平均,压力就会晃动。气体分数α接近0.5时,混合物声速降到单相的大约百分之一,这个落差会直接改写时间步长。