[论文解读] Three-dimensional solutions for the geostrophic flow in the Earth's core

Colin M. Hardy, Philip W. Livermore+3 more|arXiv (Cornell University)|Jun 18, 2018

Geomagnetism and Paleomagnetism Studies参考文献 43被引用 3

一句话总结

本文提出了一种广义方法，用于计算地球外核中任意三维泰勒态下的瞬时地转流，修正了泰勒1963年原始公式中因忽略正确边界条件而产生的错误。通过严格施加边界约束，作者推导出一种完全正则的地转流，即使在非轴对称构型下，其在自转轴上也无对数奇点，并首次给出了此类流场的全球三维实例。

ABSTRACT

In his seminal work, Taylor (1963) argued that the geophysically relevant limit for dynamo action within the outer core is one of negligibly small inertia and viscosity in the magnetohydrodynamic equations. Within this limit, he showed the existence of a necessary condition, now well known as Taylor's constraint, which requires that the cylindrically-averaged Lorentz torque must everywhere vanish; magnetic fields that satisfy this condition are termed Taylor states. Taylor further showed that the requirement of this constraint being continuously satisfied through time prescribes the evolution of the geostrophic flow, the cylindrically-averaged azimuthal flow. We show that Taylor's original prescription for the geostrophic flow, as satisfying a given second order ordinary differential equation, is only valid for a small subset of Taylor states. An incomplete treatment of the boundary conditions renders his equation generally incorrect. Here, by taking proper account of the boundaries, we describe a generalisation of Taylor's method that enables correct evaluation of the instantaneous geostrophic flow for any 3D Taylor state. We present the first full-sphere examples of geostrophic flows driven by non-axisymmetric Taylor states. Although in axisymmetry the geostrophic flow admits a mild logarithmic singularity on the rotation axis, in the fully 3D case we show that this is absent and indeed the geostrophic flow appears to be everywhere regular.

研究动机与目标

• 解决泰勒1963年地转流公式长期存在的不准确性问题，该问题错误地处理了非轴对称磁场的边界条件。
• 开发一种计算瞬时地转流的一般方法，使其与完整的三维泰勒态条件一致，包括正确施加边界条件。
• 证明地转流在三维空间中处处正则，即使存在非轴对称磁场，也与轴对称情况下预测的对数奇点相反。
• 提供一种计算高效的框架，通过在磁流体静力学近似下正确从属地转流与磁场，实现长期地磁发电机演化的建模。

提出的方法

• 该方法将地转流重新表述为一种包含完整边界条件的修正二阶常微分方程的解，从而纠正了泰勒原始方法的缺陷。
• 采用伽辽金谱方法，结合球谐函数和径向多项式基函数（雅可比多项式），在无散度、正交归一的基底下表示磁场和速度场。
• 通过模态展开求解磁流体静力学方程来计算非地转流，未知系数通过匹配径向幂次并利用符号计算（如Maple）进行解析求解。
• 通过从总速度中去除柱对称平均的方位分量来提取地转流，确保其唯一性并符合泰勒约束。
• 该方法确保了$L^2$正交归一性与原点处的正则性，并满足外核-地幔边界上的不透性条件。
• 通过将磁场投影到伽辽金基底上，并利用球谐变换和径向积分求解系数，对解进行了验证。

实验结果

研究问题

• RQ1 当考虑非轴对称磁场和正确边界条件时，泰勒1963年原始的二阶常微分方程在地转流中是否仍然有效？
• RQ2 全三维非轴对称泰勒态是否能产生在自转轴上也处处正则的地转流？
• RQ3 为确保洛伦兹力矩在所有半径和时间上恒为零（如泰勒约束所要求），地转流的正确数学公式是什么？
• RQ4 当磁流体静力学方程本身无法约束地转流时，如何在三维空间中唯一确定地转流？
• RQ5 先前模型中边界条件应用错误对地磁发电机长期演化有何影响？

主要发现

• 由于边界条件处理不完整，特别是非轴对称情况下，泰勒1963年原始的地转流方程在一般情况下无效。
• 修正后的方法产生了一种在三维空间中处处正则的地转流，消除了轴对称解中在自转轴上出现的对数奇点。
• 成功计算并验证了首个由非轴对称泰勒态驱动的全球三维地转流实例。
• 通过强制施加泰勒约束的时间连续性，地转流被唯一确定，确保其与演化中的磁场保持一致。
• 该方法通过在磁流体静力学近似下正确从属地转流与磁场，实现了对长期地磁发电机演化的精确、计算高效的建模。
• 使用雅可比多项式作为伽辽金模态，确保了$L^2$正交归一性、原点处的正则性，并能准确投影磁场和速度场。