5.4 人工粘性与冲击波捕捉技术 在爆炸力学的数值模拟世界里,冲击波不是一条平滑的曲线,而是一道撕裂连续性的刀锋——密度、压力、速度在不足一个网格单元的厚度内完成从静止到超音速的跃迁。这种极端的非线性与不连续性,使得传统的有限差分或有限元方法在面对强激波时,要么产生剧烈振荡(Gibbs现象),要么因数值耗散过大而模糊了波前结构,甚至导致计算崩溃。于是,人工粘性(Artificial Viscosity)应运而生,它不是物理真实的存在,却成为我们驯服爆炸狂澜的“隐形缰绳”。 一、为何需要人工粘性?——从物理现实到数值困境 真实爆炸中,冲击波是熵增的不可逆过程,其内部存在分子尺度的耗散机制。然而,在欧拉或拉格朗日网格的连续介质近似中,我们忽略了微观分子碰撞,仅通过守恒方程描述宏观演化。