[论文解读] Single to Double Mill Small Noise Transition via Semi-Lagrangian Finite Volume Methods

José A. Carrillo, Axel Klar+1 more|arXiv (Cornell University)|Jul 18, 2014

stochastic dynamics and bifurcation参考文献 42被引用 4

一句话总结

本文提出一种新颖的半拉格朗日有限体积方法，以精确模拟群集动力学中带有随机扰动的动能平均场方程，揭示了在小随机扰动下从单重磨坊到双重磨坊的噪声诱导相变。该方法首次实现了对这种非均匀相变的数值观测，表明在低噪声条件下双重磨坊比单重磨坊更稳定，而在较高噪声水平下则逐渐变得无序。

ABSTRACT

We show that double mills are more stable than single mills under stochastic perturbations in swarming dynamic models with basic attraction-repulsion mechanisms. In order to analyse accurately this fact, we will present a numerical technique for solving kinetic mean field equations for swarming dynamics. Numerical solutions of these equations for different sets of parameters will be presented and compared to microscopic and macroscopic results. As a consequence, we numerically observe a phase transition diagram in term of the stochastic noise going from single to double mill for small stochasticity fading gradually to disordered states when the noise strength gets larger. This bifurcation diagram at the inhomogeneous kinetic level is shown by carefully computing the distribution function in velocity space.

研究动机与目标

• 开发一种高精度数值格式，用于求解带有随机扰动的群集动力学中的动能平均场方程。
• 研究在不同噪声强度下磨坊结构——特别是单重磨坊与双重磨坊——的稳定性。
• 将动能方程的数值解与微观粒子模拟及宏观近似进行比较。
• 识别并表征非均匀群集系统中从有序磨坊状态到无序状态的噪声诱导相变。
• 通过在不同尺度下比较空间密度和速度分布函数，验证动能模型的准确性。

提出的方法

• 采用分裂格式，将空间对流项用半拉格朗日方法处理，速度空间扩散项则采用半隐式有限体积格式。
• 利用算子分裂方法将对流与扩散步骤解耦，从而提高求解福克-普朗克型动能方程时的稳定性和精度。
• 预先计算网格点之间的成对距离，以高效评估空间域中的相互作用力。
• 该格式被扩展以处理均匀与非均匀情形，并通过收敛性分析进行验证。
• 在固定空间点上计算速度分布函数，以追踪噪声逐渐增强时磨坊结构的演化。
• 定义一个衡量与单重磨坊解距离的泛函，并用于量化不同噪声水平下磨坊结构的持久性。

实验结果

研究问题

• RQ1 随着随机噪声的增加，群集动力学中单重磨坊与双重磨坊构型的稳定性如何变化？
• RQ2 动能平均场模型能否准确再现微观模拟中观察到的从有序磨坊到无序状态的转变？
• RQ3 在小噪声条件下，双重磨坊结构在群集系统的相图中起什么作用？
• RQ4 动能方程的数值解与微观粒子模拟及宏观近似在定量上如何比较？
• RQ5 半拉格朗日有限体积格式在多大程度上能够解析速度空间中非均匀的、由噪声驱动的相变？

主要发现

• 在小随机扰动下，双重磨坊比单重磨坊更稳定，可在更广范围的噪声振幅下持续存在。
• 当噪声强度 $ A = 0.123 $ 时，动能模拟与微观模拟结果高度一致，仅在动能解中表现出略高的扩散性。
• 在固定空间点上的速度分布函数随 $ A $ 增大，从单峰（单重磨坊）变为双峰（双重磨坊），最终趋于正态分布。
• 衡量与单重磨坊解距离的泛函表明，与单重磨坊初始条件相比，双重磨坊初始条件在更广范围的 $ A $ 下仍保持与单重磨坊结构的接近性。
• 在 $ A = 0 $ 时，动能方程的径向空间密度与微观和宏观结果吻合良好，但宏观模型因单动量闭合而无法捕捉双重磨坊解。
• 在高噪声（$ A = 3.0 $）下，动能解显示正态分布，与宏观扩散模型一致，证实了向无序状态的相变。