[論文レビュー] Simulated Navier-Stokes trefoil reconnection

Robert M. Kerr|arXiv (Cornell University)|Sep 10, 2015

Fluid Dynamics and Turbulent Flows参考文献 22被引用数 4

ひとこと要約

本稿は、ナビエ＝ストークス方程式における摂動を加えたトライフォイル・ボルテックス・ノットの自己再結合をシミュレーションし、粘性が0に近づく際、有限時間内にヘリシティが散逸するが、その前には長期間にわたりヘリシティが保存されることを示した。実験的主張とは対照的に、ヘリシティの保存が再結合の初期段階で成立せず、後続の段階でヘリシティの減少が生じることから、実験における再結合時間スケールの解釈が誤られていた可能性がある。

ABSTRACT

The evolution and self-reconnection of a perturbed trefoil vortex knot is simulated, then compared to recent experimental measurements (Scheeler et al. 2014a). Qualitative comparisons using three-dimensional vorticity isosurfaces and lines, then quantitative comparisons using the helicity. To have a single initial reconnection, as in the experiments, the trefoil is perturbed by 4 weak vortex rings. Initially there is a long period with deformations similar to the experiment during which the energy, continuum helicity and topological self-linking number are all preserved. In the next period, once reconnection has clearly begun, a Reynolds number independent fraction of the initial helicity is dissipated in a finite time. In contrast, the experimental analysis finds that the helicity inferred from the trajectories of hydrogen bubbles is preserved during reconnection. Since vortices reconnect gradually in a classical fluid, it is suggested that the essential difference is in the interpretation of the reconnection timescales associated with the observed events. Both the time when reconnection begins, and when it ends. Supporting evidence for the strong numerical helicity depletion is provided by spectra, a profile and visualisations of the helicity that show the formation of negative helicity on the periphery of the trefoil. A single case with the same trajectory and circulation, but a thinner core, replicates this helicity depletion despite larger Sobolev norms, showing that the reconnection timescale is determined by the initial trajectory and circulation of the trefoil, not the initial vorticity. This case also shows that the very small viscosity, $ν ightarrow 0$ mathematical restrictions upon finite-time dissipative behavior do not apply to this range of modest viscosities.

研究の動機と目的

• Scheelerら（2014a）の実験条件を再現するように、摂動を加えたトライフォイル・ボルテックス・ノットのナビエ＝ストークス方程式における自己再結合をシミュレートすること。
• 実験観察ではヘリシティ保存が報告されている一方で、数値シミュレーションでは再結合中に有限時間内にヘリシティが散逸することを示す、両者の矛盾を解明すること。
• トライフォイルの特異な位相的性質が、マクロなヘリシティの有限時間内散逸を可能にするかどうかを調査し、ナビエ＝ストークス方程式の正則性に関する知見を得ること。

提案手法

• 周期的境界条件を用いたスペクトル法を用いて、3次元非圧縮性ナビエ＝ストークス方程式の数値的シミュレーションを実施。
• 固定された循環を有する(3,2)トライフォイル・ノットとして初期ボルテックス配置を生成し、4つの弱いボルテックス・リングによる摂動を加えて単一再結合を確保。
• ガウスのリンク積分および渦度場の診断を用いて、ヘリシティ、エネルギー、自己リンク数を計算。
• スペクトル、ヘリシティプロファイル、および渦度等値面・等値線の3次元可視化を用いて、ヘリシティの移転および散逸メカニズムを分析。
• 再結合時間スケールが初期渦度構造および粘性に依存するかを検証するため、より細いコアを持つボルテックスケースをシミュレーション。
• レイノルズ数の依存性を検証するため、粘性を0に近づけ、ソボレフノルムおよびエネルギー散逸の挙動を評価。

実験結果

リサーチクエスチョン

• RQ1 摂動を加えたトライフォイル・ボルテックス・ノットのナビエ＝ストークスシミュレーションは、水素バブルトレーサーを用いた実験で観察された再結合の定性的な特徴を再現するか？
• RQ2 なぜ実験では再結合中にヘリシティ保存が報告されているのに対し、シミュレーションでは有限時間内にヘリシティ散逸が生じるのか？
• RQ3 再結合時間スケールは、初期ボルテックス軌道、循環、またはコア厚さおよび渦度分布に依存するのか？
• RQ4 エネルギー散逸が無視できる範囲でも、トライフォイル・ボルテックス・ノットが有限時間内にヘリシティを散逸させるテストベッドとして機能可能か？
• RQ5 特に最初の再結合段階のタイミングに関して、シミュレーションと実験における再結合イベントの解釈はどのように異なるのか？

主な発見

• ヘリシティは、反平行または直交するボルテックス再結合よりも50%以上長期間にわたり保存され、顕著な減衰が生じるまでに至る。
• 有限時間内にヘリシティ散逸が発生し、粘性ν → 0のときt_h ≈ 72で観測される。エネルギーの減少が微小であっても、非ゼロの散逸率が存在することを示唆している。
• 大スケールで負のヘリシティが生成され、これは大スケールのライム・ヘリシティが、小スケールのねじれヘリシティが散逸した後も持続するメカニズムを示唆している。
• 再結合時間スケールは、初期軌道、ループ間隔、循環Γにのみ依存し、初期渦度の大きさやコア厚さには依存しない。
• より大きなソボレフノルムを有する細いコアのシミュレーションでも、同じヘリシティの枯渇が再現され、時間スケールが粘性や渦度ノルムではなく、位相的性質と循環に支配されることを証明した。
• この中程度の粘性範囲では、ν → 0の極限でも有限時間内散逸が抑制されず、数学的特異性がこのような挙動に必要であるという仮説に反する。