4.3.1.1 自驱动机制 4.3.1.1 自驱动机制:当粒子“拒绝躺平”——一个被低估的数值耗散陷阱与它的实时修复方案 你有没有调试过这样一个场景: 在活性粒子模拟中,你精心设置了自驱动速度 $v0 = 0.8$,粒子质量 $m = 1$,推进力方向由内部取向角 $\thetai(t)$ 决定,驱动力写成 $\mathbf{F}i^{\text{drive}} = v0 \, (\cos\thetai,\, \sin\thetai)$;你用了标准的Verlet积分,加了Langevin热噪声,也耦合了短程排斥势——一切看起来都“教科书般正确”。可运行5000步后,粒子集体减速:平均速率从0.79跌到0.62,再跑1万步,竟滑落到0.41。