[论文解读] The baroclinic instability in the context of layered accretion. Self-sustained vortices and their magnetic stability in local compressible unstratified models of protoplanetary disks

Wladimir Lyra, Hubert Klahr|Max Planck Institute for Plasma Physics|Nov 2, 2010

Astrophysics and Star Formation Studies参考文献 26被引用 54

一句话总结

本研究探讨了磁化、可压缩、非层结的原行星盘剪切盒模型中的斜压不稳定（BI），发现当磁场与气体强耦合时，BI生成的涡旋会迅速被磁不稳定性——尤其是磁椭圆不稳定性——破坏。在磁化活跃区，MRI在增长率和饱和幅度上均占主导，使BI失效；但涡旋仅在磁耦合较弱的弱电离死区中可能存活。

ABSTRACT

Turbulence and angular momentum transport in accretion disks remains a topic of debate. With the realization that dead zones are robust features of protoplanetary disks, the search for hydrodynamical sources of turbulence continues. A possible source is the baroclinic instability (BI), which has been shown to exist in unmagnetized non-barotropic disks. We present shearing box simulations of baroclinicly unstable, magnetized, 3D disks, in order to assess the interplay between the BI and other instabilities, namely the magneto-rotational instability (MRI) and the magneto-elliptical instability. We find that the vortices generated and sustained by the baroclinic instability in the purely hydrodynamical regime do not survive when magnetic fields are included. The MRI by far supersedes the BI in growth rate and strength at saturation. The resulting turbulence is virtually identical to an MRI-only scenario. We measured the intrinsic vorticity profile of the vortex, finding little radial variation in the vortex core. Nevertheless, the core is disrupted by an MHD instability, which we identify with the magneto-elliptic instability. This instability has nearly the same range of unstable wavelengths as the MRI, but has higher growth rates. In fact, we identify the MRI as a limiting case of the magneto-elliptic instability, when the vortex aspect ratio tends to infinity (pure shear flow). We conclude that vortex excitation and self-sustenance by the baroclinic instability in protoplanetary disks is viable only in low ionization, i.e., the dead zone. Our results are thus in accordance with the layered accretion paradigm. A baroclinicly unstable dead zone should be characterized by the presence of large-scale vortices whose cores are elliptically unstable, yet sustained by the baroclinic feedback. As magnetic fields destroy the vortices and the MRI outweighs the BI, the active layers are unmodified.

研究动机与目标

• 评估斜压不稳定（BI）是否能在磁化、可压缩、局部的原行星盘模型中生成并维持大尺度涡旋。
• 确定BI、磁旋转不稳定性（MRI）和磁椭圆不稳定性在涡旋稳定性中的相互作用。
• 评估磁场、电阻率和磁场构型（垂直、方位角、净通量为零）在破坏或保持BI驱动涡旋中的作用。
• 检验当BI被视为死区中潜在湍流源时，分层吸积模型是否仍具可行性。
• 隔离在强MRI活动缺失条件下，BI导致自持涡旋的条件。

提出的方法

• 使用Pencil代码在剪切盒近似下执行局部、可压缩、非层结模拟。
• 引入线性熵梯度以驱动斜压不稳定性，通过浮力实现涡旋生成与反馈。
• 采用可压缩流体动力学框架，允许螺旋密度波激发及角动量输运测量。
• 通过改变磁场强度、构型（垂直、方位角、净通量为零）和电阻率，隔离MRI与磁椭圆不稳定性的影响。
• 通过时间解析诊断（包括涡度平方、动能与磁能、以及α参数）追踪涡旋演化、磁场增长与湍流饱和水平。
• 通过分析不稳定波长范围与涡旋核心的通道流形成，识别磁椭圆不稳定性为主要破坏机制。

实验结果

研究问题

• RQ1 在三维磁化、可压缩、非层结的原行星盘模型中，斜压不稳定性能否维持大尺度涡旋？
• RQ2 与MRI相比，磁椭圆不稳定性如何影响BI生成涡旋的稳定性？
• RQ3 磁场强度、构型与电阻率在涡旋存活或破坏中的作用是什么？
• RQ4 MRI是否在磁耦合区域抑制BI，而BI是否仅在弱耦合、低电离区域持续存在？
• RQ5 涡旋核心在径向上是否存在涡度非均匀性，以及这如何影响其对磁不稳定的敏感性？

主要发现

• 包含磁场时，斜压生成的涡旋迅速被破坏，且MRI在增长率和饱和振幅上均占主导。
• 磁椭圆不稳定性与MRI具有几乎相同的不稳定波长范围，但其增长更快，通过形成通道流拉伸并破坏涡旋的空间相干性，从而导致其破坏。
• 当涡旋纵横比趋于无穷大（纯剪切流）时，MRI被识别为磁椭圆不稳定性的一个极限情况。
• 强垂直磁场导致涡旋核心内磁场迅速增长，触发通道流，使涡旋在数个轨道周期内被破坏。
• 当最长盒尺度的磁 Reynolds 数接近1（Re_M ≲ 2）时，电阻率会抑制MRI与磁椭圆不稳定性，使涡旋存活至MRI发展之前。
• 恒定方位角磁场与零净通量磁场同样通过磁椭圆不稳定性导致涡旋破坏，表明该不稳定性在不同磁场构型下均具鲁棒性。