“ 湍流是经典物理学中最重要的未解难题。”
可压缩多相流CFD研究者的工作笔记
START HERE
选择你的第一条路线
为初次访问者保留三个清晰入口：数值基础、可压缩流动或动手实践。
本周热门文章
…
τ=0.6时流速只有六分之一 — Guo forcing漏掉半步速度F/2的后果
保留Guo源项却在速度中去掉F/2,有效体积力会缩小为(1−1/2τ)倍。Kolmogorov流中τ=0.6时振幅只有0.166倍。
混合分数稍有波动，平均温度就降 964 K —— presumed β-PDF 与 Z=0 处的无穷大
在化学计量 Z̄ 处，β-PDF 平均温度为何比 T(Z̄) 低 964 K；当 a<1 时，均匀网格积分为何偏高 115 K，用代码验证
界面只有两格宽时，表面张力少了39% — Korteweg应力的各向异性
平直气液界面上只有切向压力会下凹。下凹的面积就是表面张力，在dz=1的格子上，两格宽的界面会丢掉39%。
[论文评述] 斜率为3的MSTACS在Co=0.5时第255步发散 — 压缩差分与Courant极限
把MSTACS式(18)的压缩差分放到一维平板上测试。没有再分配时Co>0.44即发散，原因是NVD斜率k与Co的乘积。
[论文解读] 把 β 调到 2,界面停在 14 个网格,正弦波却比一阶迎风还差 — Tallois 的压缩型限制器
对体积分数使用 β=2 压缩型限制器,界面在五个周期后仍只有 14 个网格宽;同一限制器作用于正弦波时,误差却大于一阶迎风格式。
只改了 τ，静止气体的密度就变了 — Shan–Chen 速度平移 forcing 的 (τ−½)² 应力
给静止的 D1Q3 气柱施加同样的力，只有 Shan–Chen forcing 的密度随 τ 改变。平衡速度的平移留下了一个 (τ−½)²F²/ρ 应力。
用 Kp 去除逆反应速率常数，2500 K 下 AlO 竟解离了 99.6% — 平衡常数里藏着的 (p°/RT)^Δn
把 Kp 代入 kr = kf/K，Δn≠0 的反应平衡会偏差 (p°/RT)^Δn 倍。用 AlO 解离反应器复现 SI 与 CGS 两种错误
网格加密一倍,误差反而变成三倍 — Jacobi 迭代的增量收敛判据隐藏了什么
用增量判据停止 Jacobi 迭代时,网格越细,残留误差按 1/(1-ρ) 放大。以二维导热平板比较 Jacobi、SOR 与最优 ω
把 v 项和 v² 项相加，两端都对，中间却偏了 42% — Prony 两项阻力律与 Re=60
Bossut 测得阻力正比于 v²，Du Buat 测得正比于 v。Prony 的和 a·v²+b·v 在 Re=60 处交接主导权，并在附近把终端速度高估达 42%。
只给热传导单独的τ,粘性加热却被锁在Pr=1 — 双分布LB中的准平衡g*
Couette通道里,简单双τ无论Pr取何值,θ_max都贴在1.000。加上准平衡g*后,Pr=0.71给出0.710。
[论文解读] 水和空气各占一半,声音只剩24 m/s — Murrone–Guillard 五方程模型中的K∇·u项
含气50%的气泡水声速只有23.6 m/s,比水和空气都慢。去掉五方程模型的K∇·u项,同样的计算会给出745 m/s。
[论文评述] 把一个球顶在前面的四面体没有散开 — Mach 20 碎片簇的钝度指数
在 Mach 20 气流中释放的四个铁球，分离速度随四面体姿态被切开。四万个随机姿态重算了钝度 B 和非对称指数 A。