Skip to content
cfd-lab:~/zhonline
~/fortune.txt

湍流是经典物理学中最重要的未解难题。

Richard P. Feynman

可压缩多相流CFD研究者的工作笔记

WEEKLYMONCFD方法TUE流体力学WEDCFD方法THU流体力学FRICFD方法SAT论文综述SUN论文综述

本周热门文章

FILTER
#145
SUN
2026.08.30

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

突破毛细时间步约束之后上限依然在。追踪Janodet等(2025)在密度比1000下测得的稳定极限与加速最优点。

5 min read#Surface-Tension #Capillary-Wave #VOF
#144
SAT
2026.08.29

声速CFL为10时预处理器每次sweep把误差放大100倍 — 块预处理器丢掉的那一支箭头

块Gauss-Seidel只用下三角,丢掉的压力↔速度反馈增益恰好等于声速CFL的平方。本文用一维声学系统实测Schur补在那个位置做了什么。

6 min read#Schur-Complement #Preconditioning #Newton-Krylov
#143
FRI
2026.08.28

把网格佩克莱数调到2,解掉到了-0.33 — 可最小化的能量消失的地方

对流项破坏刚度矩阵的对称性,最小总势能原理也随之消失。伽辽金振荡的真正原因,以及稳定化换来的东西。

5 min read#SUPG #Weighted-Residual #FEM
#142
THU
2026.08.27

科里奥利力做的功是零,叶轮却加了1800 J/kg — 旋转坐标系里的能量与转焓

旋转参考系中守恒的不是总焓而是转焓。本文用数值验证科里奥利力为何不做功,以及泵扬程究竟从哪里来。

6 min read#Rothalpy #Rotating-Frame #Coriolis
#141
WED
2026.08.26

把 τ 直接代进去,粘度大了六倍 —— LBM 离散化留下的三处 Δt/2

一次梯形积分和一行变量替换,在粘度、速度、应力三处各留下一个 Δt/2。用标量 ODE 和 D2Q9 把三处都测了一遍。

8 min read#Trapezoidal-Rule #LBM #Viscosity
#140
TUE
2026.08.25

折角的弦没事,切断的弦跳出了9% — 1747年波动方程之争的数值结局

达朗贝尔的行波解与伯努利的正弦级数是同一个答案,只在不连续初始条件下分道扬镳。用代码找出9%过冲的来源。

4 min read#Gibbs-Phenomenon #Wave-Equation #Characteristics
#139
MON
2026.08.24

壁面上是3个,角上是5个 — 数一数LBM边界节点丢失的分布函数

迁移之后空缺的分布函数个数,等于该节点周围固体格子的数量:平壁3个,凹角5个,凸角1个。

4 min read#Bounce-Back #Zou-He #Boundary-Conditions
#138
SUN
2026.08.23

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

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

6 min read#Two-Fluid-Model #Hyperbolicity #Paper-Review
#137
FRI
2026.08.21

把压力直接代入后静止的界面开始颤动 — 非理想LBM forcing的两种形式

用格子玻尔兹曼求解范德瓦尔斯流体时,力项里代入压力p与代入化学势mu在连续下相同,在离散下却分道扬镳。

6 min read#Van-der-Waals #Forcing-Term #Gibbs-Duhem
#136
THU
2026.08.20

让烟先下降1 m,抽力就少了7 Pa — 富兰克林壁炉与浮力预算

烟囱顶端高度没变,一段下降流道却让抽力崩塌。用浮力账本和倒烟阈值把这笔损失算清楚。

6 min read#Stack-Effect #Natural-Convection #Buoyancy
#135
WED
2026.08.19

少取一个积分点,解就发散了 — DG 的求积下限与泰勒基

决定 DG 单元积分最少高斯点数的是被积多项式的次数,不是网格。三点降到两点毫无变化,降到一点则直接发散。

7 min read#Discontinuous-Galerkin #Quadrature #FEM
#134
TUE
2026.08.18

激波停在了起跑线上 —— 守恒形式与原始形式的分岔口

同一个方程写成通量散度还是速度乘梯度,算出的激波速度就不同。网格加密八倍,仍然收敛到错误的值。

6 min read#Conservative-Form #Rankine-Hugoniot #Shock-Capturing
$ 12 / 12 条目 · 共 17.3K 字