[論文レビュー] 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.

研究の動機と目的

• 閉じ込められた極性アクティブ流体におけるアクティブな乱流に至る動的遷移を調査すること。この分野では、このような経路がまだ十分に理解されていない。
• 任意の境界形状を有する複雑な幾何形状および大規模な流体力学系のシミュレーションにおける計算的課題を克服すること。
• 秩序ある渦と完全に発達した乱流の間の中間状態（周期的、混沌的、準周期的）を特定し、それらを特徴づけること。
• 任意の閉じ込め幾何形状に対して高分解能で効率的にToner-Tu-Swift-Hohenberg方程式を解くことが可能なGPUベースの数値スキームの開発および検証すること。

提案手法

• 任意の境界条件を想定した極性アクティブ流体のToner-Tu-Swift-Hohenberg（TTSH）方程式を高効率に統合するため、GPUアクセラレートされたスペクトル解法を開発した。
• スペクトル法における aliasing を回避するために、新しいマスク処理アルゴリズムを実装した。この手法では、バイナリーマスクの低域透過フィルタ処理を施したDFTを計算し、逆DFTを実行した後、二乗することで非負性を保証した。
• 空間マスクを定義することで、流体領域と境界領域を区別する統一フレームワークを採用し、円形およびダムベル型の閉じ込めを含む任意の幾何形状に対応した。
• GPUクラスタ（NVIDIA A6000、Kugui、nekoya/ai）を用いた高分解能シミュレーションを実施し、大規模格子においてCPUと比較して約60倍の高速化を達成した。
• 混沌たるダイナミクスを特徴づけるためにリャプノフ指数を計算した。GPU並列化のおかげで、膨大な計算量を要するこの解析が可能になった。
• 細菌乱流の実験的観察に基づいて境界条件を設定し、物理的妥当性を確保した。

実験結果

リサーチクエスチョン

• RQ1 閉じ込め領域のサイズを増加させた場合、極性アクティブ流体においてどのような動的遷移が生じるか？
• RQ2 閉じ込められた極性流体におけるアクティブな乱流への経路は、ナビエ-ストークス方程式やアクティブネマチック系における既知の経路とどのように異なるか？
• RQ3 周期的状態および混沌的状態は、秩序ある渦から完全に発達した乱流への遷移において果たす役割は何か？
• RQ4 特に円形とダムベル型の閉じ込め形状において、幾何学的形状が遷移の性質や位置にどのように影響を与えるか？
• RQ5 統一されたGPUベースの数値フレームワークは、任意の境界形状を有する大規模なアクティブ流体ダイナミクスを正確にシミュレートできるか？

主な発見

• 単一の静止渦から振動する渦対への最初の遷移はヒステリシス的であり、亜臨界定常分岐（subcritical Hopf bifurcation）に一致する。
• 閉じ込め半径が増加するに従い、系は以下の順序で遷移を経る：周期的振動 → 混沌的ダイナミクス → 無理数比周波数を有する準周期的状態 → アクティブな乱流。
• アクティブな乱流への遷移は、周期的状態と混沌的状態の間を繰り返す再入的遷移を含む、複雑で単調でない経路を経由する。
• ダムベル型の閉じ込めでは、フェルマー磁性体-反強磁性体転移に類似した渦秩序転移が観察されたが、その幾何依存性は、以前の実験報告とは定性的に異なる。
• GPUベースの解法は、大規模シミュレーションにおいてCPUと比較して約60倍の高速化を達成し、一時的および混沌的ダイナミクスの高分解能解析を可能にした。
• リャプノフ指数解析により、混沌たるダイナミクスの発生が確認された。正の指数は、乱流領域において初期条件に敏感であることを示した。