No evidence of vorticity production from initially irrotational turbulent gravitational collapse

技术摘要：初步无旋湍流引力塌缩中无涡度产生证据

问题陈述
引力塌缩被广泛认为会产生湍流，这可能为小尺度发电机驱动磁场放大提供种子，使其在 Jeans 尺度以下增强。然而，小尺度发电机从根本上需要涡旋（旋转）湍流，而引力塌缩主要产生无旋（势流）流。在缺乏磁场或外部驱动力，且采用巴洛特罗皮克状态方程（即压力和密度梯度平行）的情况下，斜压项无法产生涡度。因此，从最初的无旋流中产生涡度的唯一机制是通过与粘性的相互作用。先前的研究表明，塌缩诱导的压缩可能会驱动发电机作用，但目前尚不清楚塌缩流本身是否产生了必要的涡度，或者观察到的发电机效应仅仅是预先存在的初始湍流的衰减。

研究方法
作者使用 Pencil Code 对引力塌缩进行了直接数值模拟（DNS），采用了六阶中心差分和三阶时间步进方案。研究利用了一个周期性立方区域，并采用等熵状态方程（$\gamma = 5/3$）以及经过修正的等温 Bonnor-Ebert 球作为初始密度分布。

关键方法论选择包括：

• 初始条件： 速度场在傅里叶空间中被构建为严格无旋的，具有随机相位和特定的谱斜率。研究未包含磁场，以隔离涡度产生机制。
• 粘性： 作者没有使用冲击粘性（这可能会掩盖结果），而是使用了一个恒定的运动粘度 ($\nu$)。
• 分辨率与参数： 模拟在 $512^3$ 到 $2048^3$ 个网格点的分辨率范围内运行。作者通过改变初始速度振幅 ($u_{ini}$) 和粘度 ($\nu$) 来确定数值可靠的极限。
• 分析： 研究追踪了恩斯特罗菲（enstrophy，$\langle \omega^2 \rangle/2$）的演化，并将其时间导数分解为三个项：发电机项 ($T_{dyn}$)、由粘性和密度梯度驱动的产生项 ($T_{gen}$)，以及耗散项 ($T_{dis}$)。作者还分析了速度、密度和涡度谱，并检查了涡度随初始流振幅的标度关系。

关键结果

1. 粘性依赖性： 发现涡度的产生是一个次级效应，由通过粘性实现的初始流散度的放大所驱动。研究发现，涡度产生与耗散的比率 ($R_{gen}$) 与粘度不成比例；相反，涡度产生在主塌缩发生前的早期阶段最为显著，而不是作为塌缩流本身的后果。
2. 与初始条件的标度关系： 均方根（rms）涡度 ($\omega_{rms}$) 与初始速度振幅的平方 ($u_{ini}^2$) 成比例，而 rms 速度散度与 $u_{ini}$ 成线性比例。这种标度关系在整个模拟过程中持续存在，表明引力塌缩并不直接驱动涡度产生。相反，塌缩放大了初始流的散度，随后粘性将其转化为涡度。
3. 数值可靠性： 作者确定，当网格雷诺数 ($Re_{mesh}$) 超过 0.7 至 2.0 之间的值时（取决于具体运行情况），结果变得不可靠。在可靠范围内，发电机项 ($T_{dyn}$) 保持次要地位，对恩斯特罗菲演化的贡献小于 20%。
4. 谱分析： 涡度谱表现出 $k^4$ 的亚惯性区间，类似于 MHD 中的磁场行为，但作者将其归因于 Jeans 不稳定性驱动的随机速度波动，而非发电机过程。速度谱显示出一个陡峭的惯性区间，且没有出现“瓶颈”效应。
5. 发电机临界涡度： 通过重新审视发电机启动的临界磁涡度雷诺数 ($Pr_M Re_\omega$)，作者确认了在 300 左右存在分叉。这意味着纯无旋流在没有外部驱动或特定初始条件的情况下，无法维持足够的涡度，因此不会作为发电机发挥作用。

意义与主张
论文得出结论：在巴洛特罗皮克状态方程和初始无旋湍流的条件下，引力塌缩并不能有效地产生涡度。观察到的涡度是一个次级效应，是由于粘性将初始流的散度转化为涡度所致，而非塌缩动力学的直接结果。

作者断言，他们的发现挑战了先前关于塌缩诱导小尺度发电机作用的说法。他们认为，此前在塌缩流中观察到的磁场放大现象，很可能是由初始湍流的衰减驱动的，而非由塌缩本身驱动。该研究强调了不使用亚网格尺度建模进行直接数值模拟的必要性，以避免人工产生的涡度。这项工作意味着，在这种塌缩场景下实现自持发电机是非常困难的，且仅靠塌缩流本身不足以产生发电机作用所需的涡旋湍流。