5.2.2.1 欧拉网格求解 5.2.2.1 欧拉网格求解:当压力泊松方程在非均匀网格上“拒不解耦”——一个被低估的离散一致性陷阱与三行半代码的救赎 你有没有在凌晨三点盯着屏幕上那团永不收敛的压力残差曲线发呆? 它像一条垂死的蛇,在 $10^{-3}$ 附近反复痉挛,忽而跳到 $10^{-2}$,忽而又滑回 $10^{-4}$,却始终不肯沉入 $10^{-6}$ 的寂静深渊。 你检查了时间步长——够小;验证了边界条件——无懈可击;重编译了稀疏求解器——OpenFOAM、SU2、自研矩阵库全试过;甚至把初始场从零场换成解析解扰动场……它依然在笑。 不是数值噪声,不是硬件浮点误差,不是并行通信失配——是离散算子在非均匀欧拉网格上悄然背叛了连续世界的守恒契约。 这不是调试失败,是建模失贞。