[论文解读] Simulated Navier-Stokes trefoil reconnection

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

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

一句话总结

本文在纳维-斯托克斯方程中模拟了受扰动的三叶形涡 knot 的自重连过程，发现在粘度趋近于零时，尽管在有限时间内发生耗散，但其手性（helicity）在较长时间内保持不变。与实验中声称重连过程中手性守恒的观点相反，模拟结果表明手性耗散仅在初始重连阶段之后才开始，暗示实验中对重连 timescale 的理解存在误解。

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.

研究动机与目标

• 在纳维-斯托克斯方程中模拟受扰动的三叶形涡 knot 的自重连过程，复现 Scheeler 等人（2014a）实验中的条件。
• 解决实验观察中声称手性守恒与数值模拟中显示重连期间手性发生有限时间耗散之间的矛盾。
• 探究三叶形拓扑结构是否能导致宏观手性的有限时间耗散，从而为纳维-斯托克斯方程的正则性提供新见解。

提出的方法

• 使用谱方法对三维不可压缩纳维-斯托克斯方程进行数值模拟，采用周期性边界条件。
• 初始涡结构为 (3,2) 三叶形 knot，具有固定环量，并通过四个弱涡环扰动，以确保仅发生单次重连。
• 利用高斯链接积分和涡度场诊断方法计算手性、能量和自 linking 数。
• 通过谱分析、手性分布图以及涡度等值面和涡线的三维可视化，分析手性转移与耗散机制。
• 模拟了更细核心的涡结构，以检验重连 timescale 对初始涡度结构和粘度的依赖性。
• 通过减小粘度趋近于零，研究雷诺数的依赖性，评估 Sobolev 范数和能量耗散的行为。

实验结果

研究问题

• RQ1 纳维-斯托克斯方程对受扰动三叶形涡 knot 的模拟是否能再现氢气泡示踪实验中观察到的重连定性特征？
• RQ2 为何模拟显示有限时间内的手性耗散，而实验却报告重连过程中手性守恒？
• RQ3 重连 timescale 由什么决定——初始涡轨迹、环量，还是核心厚度与涡度分布？
• RQ4 三叶形涡 knot 是否可作为粘性流体中手性有限时间耗散的测试平台，即使能量耗散可忽略不计？
• RQ5 在模拟与实验中，对重连事件的解读有何不同，特别是首次重连阶段的时间判断？

主要发现

• 在显著衰减开始前，手性保持时间比反平行或正交涡重连长 50% 以上。
• 在粘度 ν → 0 时，有限时间内的手性耗散发生在 t_h ≈ 72，表明尽管能量衰减可忽略，仍存在非零耗散率。
• 在大尺度上产生了负的手性，暗示一种机制：大尺度的 writhe-helicity 被转移并持续存在，而小尺度的 twist-helicity 已被耗散。
• 重连 timescale 仅取决于初始轨迹、环路间距和环量 Γ，与初始涡度大小或核心厚度无关。
• 更细核心的模拟中，Sobolev 范数更大，但手性耗散时间与原始模拟一致，证明 timescale 由拓扑结构和环量决定，而非粘度或涡度范数。
• 在中等粘度范围内，ν → 0 的极限并未抑制有限时间耗散，挑战了此类行为需要数学奇点的假设。