[论文解读] Geodynamo simulations with vigorous convection and low viscosity

Nathanaël Schaeffer, D. Jault+2 more|arXiv (Cornell University)|May 31, 2017

Geomagnetism and Paleomagnetism Studies被引用 4

一句话总结

本研究通过超低粘度（E = 1e-7）和高磁 Reynolds 数（Rm > 500）的地球发电机模拟，揭示了强磁场主导的偶极场，其磁能比动能高出一个数量级。模拟显示，切线圆柱体内存在扭曲的极地涡旋、稳定层结以及频繁的扭转阿尔芬波，而大尺度、准地转的 m=1 涡旋在无边界强迫下自发形成，其动力学受科里奥利-浮力平衡控制，且磁场所对齐的流动占主导地位。

ABSTRACT

We analyze a suite of three convection-driven dynamo simulations in a rapidly rotating spherical shell. These three simulations form a sequence along which the viscosity is decreased, while maintaining a large magnetic Reynolds number (Rm > 500). Our most extreme case is characterized by an Ekman number E = 1e-7 and magnetic Prandtl number Pm = 0.1. In this strong-field, dipole-dominated dynamo, the magnetic energy is one order of magnitude larger than the kinetic energy. The spatial distribution of magnetic intensity is highly heterogeneous, and a stark dynamical contrast exists between the interior and the exterior of the tangent cylinder (the cylinder parallel to the axis of rotation that circumscribes the inner core). In the interior, the magnetic field is strongest, and is associated with a vigorous twisted polar vortex, whose dynamics may occasionally lead to the formation of a reverse polar flux patch at the surface of the shell. Furthermore, the strong magnetic field also allows accumulation of light material within the tangent cylinder, leading to stable stratification there. Torsional Alfven waves are frequently triggered in the vicinity of the tangent cylinder and propagate towards the equator. Outside the tangent cylinder, the magnetic field inhibits the growth of zonal winds and the kinetic energy is mostly non-zonal. Spatio-temporal analysis indicates that the low-frequency, non-zonal flow is quite geostrophic (columnar) and predominantly large-scale: an m=1 eddy spontaneously emerges in our most extreme simulations, without any heterogeneous boundary forcing. Our spatio-temporal analysis further reveals that (i) the low-frequency, large-scale flow is governed by a balance between Coriolis and buoyancy forces – magnetic field and flow tend to align, minimizing the Lorentz force; (ii) the high-frequency flow obeys a balance between magnetic and Coriolis forces; (iii) the convective plumes mostly live at an intermediate scale, whose dynamics is driven by a 3-term MAC balance – involving Coriolis, Lorentz and buoyancy forces. However, small-scale (E^{1/3}) quasi-geostrophic convection is still observed in the regions of low magnetic intensity.

研究动机与目标

• 研究在极低粘度和高磁 Reynolds 数极端条件下对流驱动地球发电机的行为。
• 理解强磁场如何影响流动组织化，特别是切线圆柱体区域内的影响。
• 研究在无非均质边界强迫条件下，大尺度非纬向流动的形成及其动力学平衡机制。
• 表征磁能与动能在多个时间尺度上的时空动力学特征。

提出的方法

• 在快速旋转的球形壳中进行三次对流驱动发电机模拟，逐步降低粘度，同时保持 Rm > 500。
• 在最极端情况下采用 Ekman 数 E = 1e-7 和磁 Prandtl 数 Pm = 0.1，以实现高强度磁场。
• 分析磁能与动能的空间分布，重点关注切线圆柱体内外的异质性及动力学差异。
• 进行时空分解，识别由不同力平衡控制的流动模态：科里奥利-浮力平衡、科里奥利-磁力平衡以及 MAC（科里奥利-洛伦兹-浮力）平衡。
• 追踪扭转阿尔芬波的传播及其与切线圆柱体和磁场结构的关联。
• 识别在低频流动分量中自发出现的大尺度 m=1 涡旋，且无外部强迫作用。

实验结果

研究问题

• RQ1 超低粘度（E = 1e-7）如何影响地球发电机中磁能与动能的组织化？
• RQ2 在无边界强迫条件下，哪些动力学机制控制大尺度非纬向流动（如 m=1 涡旋）的形成与持续？
• RQ3 磁场如何影响切线圆柱体内的流动结构与稳定性，特别是在层结与波传播方面的关联？
• RQ4 低频大尺度流动由何种力平衡控制？其与高频动力学有何差异？
• RQ5 在磁强度较低的区域，对流射流与小尺度运动在多大程度上仍持续存在？其行为由何种平衡控制？

主要发现

• 最极端模拟中，磁能比动能高出一个数量级，表明存在强偶极主导的发电机作用。
• 磁强度高度异质，最强磁场集中于切线圆柱体内，与扭曲的极地涡旋及偶发的反向极磁通量斑块相关。
• 切线圆柱体区域因轻质物质积聚而呈现稳定层结，由强磁场维持。
• 扭转阿尔芬波频繁在切线圆柱体附近激发，并向赤道方向传播，表明存在活跃的磁耦合。
• 切线圆柱体外部，磁场抑制了纬向风的增强，动能主要表现为非纬向流动，流动以大尺度、准地转、柱状运动为主。
• 低频大尺度流动由科里奥利力与浮力的平衡主导，磁场与流动方向对齐以最小化洛伦兹力；高频动力学遵循科里奥利-磁力平衡；中等尺度射流则服从三元 MAC 平衡。