[论文解读] Route to turbulence via oscillatory states in polar active fluid under confinement

Sora Shiratani, Kazumasa A. Takeuchi+1 more|arXiv (Cornell University)|Apr 6, 2023

Micro and Nano Robotics被引用 4

一句话总结

本研究通过GPU加速模拟Toner-Tu-Swift-Hohenberg模型，揭示了受限条件下极性活性流体中主动湍流的新路径。随着受限尺寸的增加，系统经历了一条滞后性转变路径，从静止涡旋态演变为振荡态，随后进入混沌与准周期动力学，最终达到主动湍流——这一路径与传统的湍流形成机制截然不同。

ABSTRACT

We report a novel route to active turbulence, observed in numerical simulations of a polar active fluid model under confinement. To deal with large-scale computations with arbitrary geometries, we developed a GPU-based scheme that can be used for any boundary shape in a unified manner. For the circular confinement, as the radius was increased, we found a series of transitions first from a single stationary vortex to an oscillating pair of vortices, then through reentrant transitions between oscillatory and chaotic dynamics before finally reaching the active turbulence. The first transition turned out to be hysteretic, with the emergence of the oscillatory state consistent with the subcritical Hopf bifurcation. In dumbbell-shaped boundaries composed of two overlapping circles, we observed a transition comparable to the ferromagnetic-antiferromagnetic vortex-order transition reported in previous experiments, but the transition point turned out to show a qualitatively different geometry dependence.

研究动机与目标

• 研究受限极性活性流体中导致主动湍流的动力学转变过程，该领域内此类转变路径尚不清晰。
• 克服在复杂几何形状与大规模流体动力学系统中模拟的计算挑战，尤其针对任意边界形状的系统。
• 识别并表征有序涡旋与完全发展湍流之间的中间态——振荡态、混沌态与准周期态。
• 开发并验证一种基于GPU的数值方法，能够高效求解任意受限几何形状下的Toner-Tu-Swift-Hohenberg方程，实现高分辨率模拟。

提出的方法

• 开发了一种基于GPU加速的谱求解器，用于高效积分极性活性流体在任意边界条件下的Toner-Tu-Swift-Hohenberg（TTSH）方程。
• 实现了一种新型掩膜处理算法，以防止谱方法中的混叠现象，方法为：对二值掩膜进行低通滤波离散傅里叶变换（DFT），经逆DFT后平方，确保结果非负。
• 通过定义区分流体区域与边界的时空掩膜，该方法构建了适用于任意几何形状（包括圆形与哑铃形）的统一框架。
• 在GPU集群（NVIDIA A6000，Kugui，nekoya/ai）上执行高分辨率模拟，对大规模网格实现约60倍于CPU的加速。
• 通过GPU并行化实现大规模计算，计算了李雅普诺夫指数以表征混沌动力学。
• 边界条件基于细菌湍流的实验观测结果设定，确保物理合理性。

实验结果

研究问题

• RQ1 在极性活性流体中，随着受限区域尺寸的增大，会发生哪些动力学转变？
• RQ2 受限极性流体中主动湍流的形成路径与Navier-Stokes系统或活性向列系统中的已知路径有何不同？
• RQ3 振荡态与混沌态在有序涡旋向完全发展湍流转变过程中扮演何种角色？
• RQ4 受限几何形状（特别是圆形与哑铃形）如何改变转变的性质与位置？
• RQ5 基于GPU的统一数值框架能否准确模拟具有任意边界形状的大规模活性流体动力学？

主要发现

• 从单个静止涡旋首次转变为振荡涡旋对的过程具有滞后性，与亚临界Hopf分岔一致。
• 随着受限半径增大，系统经历一系列转变：周期性振荡 → 混沌动力学 → 准周期态（具有无理频率比） → 主动湍流。
• 主动湍流的形成通过一条复杂且非单调的路径实现，其间振荡态与混沌态发生再入式转变。
• 在哑铃形受限条件下，观察到类似铁磁-反铁磁转变的涡旋有序转变，但其几何依赖性与以往实验报告存在定性差异。
• 基于GPU的求解器在大规模模拟中实现约60倍于CPU的加速，使高分辨率分析瞬态与混沌动力学成为可能。
• 李雅普诺夫指数分析证实了混沌的出现，正值表明在湍流区域中对初始条件具有敏感依赖性。