-
地球的自转轴在空间相对黄轴做周期性的类圆锥形运动,Baron[1]通过整理天文学数据得到了25 960 a的地轴进动周期。这一发现被广泛应用于古气候学和天文学研究中,但在地球外核的数值仿真模拟研究中,由于先前的计算机算力不足,相关软件支持较少,人们通常都忽略了地轴进动对地球磁场仿真模拟的影响。目前,国内外对引入地轴进动的发电机数值模拟研究仍较少。自文献[2]公布初代发电机模型以来,国际上陆续发表了KS95、KB97等改进模型[3-4],但绝大多数模型均忽略地轴进动产生的庞加莱力。德国亥姆霍兹德累斯顿罗森多夫研究中心正在进行一项由进动驱动的下一代地球外核实验模拟项目[5];文献[6]建立了该项目的数值模拟模型,但最终结果没有出现运动发电机现象;文献[7]对地球内核与磁场进行了研究。本文借助已有的地球发电机基准模型,在改进部分参数的基础上引入地轴进动参数,通过数值计算研究其对地球主磁场变化的影响。
-
地球发电机模型的求解方程为:
式中,
为密度; 为速度场; 为时间; 为哈密顿算子; 为修正后的包含离心力的压力; 为真空下磁导率; 为磁场; 为重力加速度; 为自转角速度; 为进动角速度; 为应变率张量; 为相对半径矢量。应变率张量的计算如下: 式中,
为外核粘滞系数; 为狄拉克函数; 、 表示 在不同坐标轴上的分量; 表示不同坐标轴的单位向量。 地轴进动是日、月对地球隆起部分的万有引力产生的旋转力矩导致地轴围绕黄轴在圆锥面上的匀速旋转。庞加莱力用于描述行星自转速度发生变化时所受的额外力,在本文中为日、月引力[8]。旋转坐标系中观察到的流体加速度为:
式中,
表示单位质量的科里奥利力,总是垂直于速度方向,不做功,因此不改变动能的大小,只改变流动的方向; 表示单位质量的离心力; 表示单位质量的庞加莱力,在物体角速度的振幅或方向随时间变化时产生,在地球发电机模型中认为其是由地轴的进动产生的。在地球发电机模型的数值模拟中,庞加莱力由无量纲的庞加莱数表示,计算如下: -
1)无量纲参数
本文基于文献[9]提出的发电机基准模型,将普朗特数设置为
,磁普朗特数设置为 ,在此基础上分别设置了埃克曼数 、瑞利数 和 、 两种参数组合。 2)热扩散率
Christensen模型中热扩散率设置为恒定的
,本文通过Gastine的地球流体内核热扩散率理论公式[10]来模拟地球内部的热扩散率,计算如下: 式中,
为相对半径; 为单位密度; 为地球外核半径的相对值,其计算如下: 式中,
为内核半径; 为外核半径。 3)内核参数
文献[11]利用地震波传播数据计算得到内核超速旋转的速率约为每百万年0.1°~1°,本文取该值为1°。由于其是无量纲的旋转速率参数,
时,超速旋转速率应按时间尺度大小转换为 ;在 时,超速旋转速率为 。此外,内核导电性与外核保持一致[12]。 4)进动参数
本文设置地轴进动周期为25 960 a,与黄轴夹角为23°26′。
-
为了便于描述,定义模型1为E=1.0 ×10-3、Ra=2.0×105的无进动模型;模型2定义为E=1.0 ×10-3、Ra=2.0×105,且以25 960 a为周期的进动模型;模型3定义为E=8.0×10-4、Ra=2.5×105的无进动模型;模型4定义为E=8.0×10-4、Ra=2.5×105,且以25 960 a为周期的进动模型。每个模型均运行6个粘性时间尺度,粘性时间尺度的计算如下:
式中,
为地球外核厚度。 基于上述对文献[9]发电机基准模型的改进,通过比较各模型计算得到的外核动能、磁能,以及发电机数值计算的磁场变化,研究引入地轴进动后对发电机模型的影响。极向-环向分解后的动能
和磁能 分别计算如下: 式中,u为速度;V为单位网格体积;
分别为极向-环向分解的级数、阶数; 为密度扰动; 、 分别为动力场极、环向位能; 、 分别为磁场极、环向位能; 为半径。 核-幔边界处的磁场强度通常使用埃尔泽塞尔数表示,其大小等于无量纲磁场强度的平方,计算如下:
式中,
为磁场的方均根(root mean square,RMS)值; 为内核边界处磁扩散率; 为外核边界处密度。 磁雷诺数为磁流体内部的速度
和磁场 的RMS值之比,计算如下: 式中,
为速度场的RMS值; 为外核厚度。 -
图 1给出了模型1、2在粘性时间尺度内的外核极向动能;图 2给出了模型3、4在粘性时间尺度内的外核极向动能。由图1、2可知,对于极向动能,随着地轴进动的引入,极向上的动能趋向稳定,模型2相对模型1的平均极向动能增强了1.5%,模型4相对模型3的平均极向动能增强了2.7%,说明地轴进动对极向动能的影响是非常微弱的。
图 1 模型1、2在6个粘性时间尺度内的外核极向动能
Figure 1. Poloidal Kinetic Energy of the Outer Core Calculated by Models 1 and 2 in 6 Viscous Time Scales
图 2 模型3、4在6个粘性时间尺度内的外核极向动能
Figure 2. Poloidal Kinetic Energy of the Outer Core Calculated by Models 3 and 4 in 6 Viscous Time Scales
图 3给出了模型1、2在粘性时间尺度内的外核环向动能;图 4给出了模型3、4在粘性时间尺度内的外核环向动能。由图3、4可知,对于环向动能,引入地轴进动将产生较大影响。模型2相比模型1平均环向动能增加了22.7%,模型4相比模型3平均环向动能增加了13.9%,这表明进动项增强了磁流体的环向运动,使得环向动能增大。
图 3 模型1、2在6个粘性时间尺度内的外核环向动能
Figure 3. Toroidal Kinetic Energy of the Outer Core Calculated by Models 1 and 2 in 6 Viscous Time Scales
图 4 模型3、4在6个粘性时间尺度内的外核环向动能
Figure 4. Toroidal Kinetic Energy of the Outer Core Calculated by Models 3 and 4 in 6 Viscous Time Scales
由图1~4可以看出,当模型处于无进动的状态时,较大的瑞利数引起模型内热对流加强使模型在极向与环向上的动能产生更频繁且无序的波动。
-
图 5给出了模型1、2在粘性时间尺度内的外核极向磁能;图 6给出了模型3、4在粘性时间尺度内的外核极向磁能。图 7给出了模型1、2在粘性时间尺度内的外核环向磁能;图 8给出了模型3、4在粘性时间尺度内的外核环向磁能。由图5~8可知,由于地轴进动的引入,外核极向、环向磁能稳定性将明显强于原始模型。对于外核极向磁能,模型2相比模型1的平均外核极向磁能减小了16.2%,模型4相比模型3则减少了14.0%。而对于环向磁能,模型2相比模型1的平均环向磁能仅增加了4.7%,而模型4相比模型3的平均环向磁能增加了7.4%。
图 5 模型1、2在6个粘性时间尺度内的外核极向磁能
Figure 5. Poloidal Magnetic Energy of the Outer Core Calculated by Models 1 and 2 in 6 Viscous Time Scales
图 6 模型3、4在6个粘性时间尺度内的外核极向磁能
Figure 6. Poloidal Magnetic Energy of the Outer Core Calculated by Models 3 and 4 in 6 Viscous Time Scales
-
图 9给出了模型1、2在核-幔边界处的磁场强度,图 10给出了模型3、4在核-幔边界处的磁场强度。由图9、10可知,引入进动后磁场强度的变化幅度明显降低。通过数值计算,发现模型2较模型1的平均磁场强度降低了9.1%,而模型4较模型3的平均磁场强度降低了8.6%。
-
用偶极子(axial dipole,
)强度与偶极子和非偶极子(non-axial-dipole, )强度之和的比表示偶极性[3]。图 11为模型1、2在6个粘性时间尺度内的核-幔边界处的偶极子分量占比,图 12为模型3、4在6个粘性时间尺度内的核-幔边界处的偶极子分量占比。 图 11 6个粘性时间尺度内模型1、2在核-幔边界处的偶极子分量占比
Figure 11. Dipolarity at Core-Mantle Boundary of Models 1 and 2 in 6 Viscous Time Scales
图 12 6个粘性时间尺度内模型3、4在核-幔边界处的偶极子分量占比
Figure 12. Dipolarity at Core-Mantle Boundary of Models 3 and 4 in 6 Viscous Time Scales
由图11、12可知,由于进动项的引入,磁场的偶极性减弱了,但其减少程度还不足以使该发电机变成多极子为主导的情况,偶极性数值仍满足类地发电机成立的条件[3]。模型2较模型1的平均偶极性减少了9.5%,模型4较模型3减少了6.2%。
-
图 13给出了模型1、2在6个粘性时间尺度内的磁雷诺数,图 14给出了模型3、4在6个粘性时间尺度内的磁雷诺数。
图 13 模型1、2在6个粘性时间尺度内的磁雷诺数
Figure 13. Magnetic Reynolds Number of Models 1 and 2 in 6 Viscous Time Scales
图 14 模型3、4在6个粘性时间尺度内的磁雷诺数
Figure 14. Magnetic Reynolds Number of Models 3 and 4 in 6 Viscous Time Scales
由图13、14可以看出,进动项的引入明显增大了磁雷诺数的数值。通过数值计算发现,模型2较模型1的平均磁雷诺数增加了7.3%,模型4较模型3增加了5.2%。这一结论与§3.1、§3.2中得到的增强环向动能并减弱极向磁能的结果相一致。
-
Christensen基准模型存在主磁场西向漂移的现象,这种现象随着地轴进动的引入变得更加剧烈和显著。本文计算了在6个粘性时间尺度下4个模型的磁场西向漂移圈数,结果如表 1所示。由表 1可知,自东向西旋转的地轴进动驱动了发电机模型中磁场的西向漂移,地磁场西漂速度相比无进动模型明显加快,是地轴进动模型的磁场表现出的一个显著特点。环向动能在发电机模型中更多表现在磁场的西向漂移上,而进动的引入使得环向动能较无进动模型有了明显的增大,进而导致了磁场的西向漂移速度增大。
表 1 模型1~4在6个粘性时间尺度的西向漂移圈数
Table 1. Circles of West Drifts of 4 Models in 6 Viscous Time Scales
模型 西向漂移圈数 1 32±3 2 49±3 3 35±3 4 44±3 -
本文借助已有的地球发电机基准模型,在改进部分参数的基础上引入地轴进动参数,通过数值计算研究其对地球主磁场变化的影响。地轴进动对Christensen基准模型的磁流体运动与磁场具有稳定作用,并使得环向动能稳定在一个较高的数值水平,而对于极向磁能,则是稳定在一个较低的数值水平。在此基础上,联合核-幔边界处的磁场强度和磁雷诺数的比较,表明进动项使得磁流体的能量更倾向于向动能转换,特别是环向动能。而环向动能在发电机模型中更多地体现在磁场的西向漂移速率上,通过对西向漂移圈数的统计也证实了向西漂移速率的明显加快。此外,本文还分析了核-幔边界处的偶极性,发现进动项的引入会降低模型的偶极性,但仍属于类地发电机模型。
本文研究的模型选用的埃克曼数量级较大,因此无法评估其他参数发电机模型中进动项的影响,但通过模型的测试表明,进动项的加入很可能会导致埃克曼数量级较小的类地发电机模型变为多极子主导的发电机模型。
-
摘要: 在地球外核磁流体发电机数值模拟研究中,与地轴进动相关的庞加莱力通常视其影响极小而在求解发电机方程组中将其忽略,但实际上地轴进动的周期相对地球磁场的倒转周期来说仍是一个值得考虑的项。基于地球发电机数值模拟基准模型,通过引入周期为25 960 a的地轴进动速度,比较分析不同埃克曼数和瑞利数参数组合模型中地轴进动对发电机数值模型的作用。结果表明,地轴进动使得球壳磁流体的动能及磁能波动稳定在一个较原始模型更小的范围内,但增大球壳磁流体的环向动能10%~20%以上,这导致磁场西向漂移速率明显加快;通过比较各模型间的核-幔边界磁场强度和磁雷诺数,发现地轴进动项的引入使得磁流体的能量更倾向于向动能转换;通过比较核-幔边界的偶极性,发现地轴进动项的引入会降低磁流体发电机数值模型的偶极性,但对于所选的原始模型,偶极性的减小程度还不足以使其成为非类地发电机模型。通过数值模拟研究,发现地轴进动项的加入很可能会导致埃克曼数量级较小的类地发电机模型变为多极子主导的发电机模型,因此,在更为准确的磁流体发电机数值模拟研究中,必须考虑地轴进动项。Abstract:
Objectives The Poincaré force related to the precession of Earth's axis is usually omitted in the solution of the dynamo equations as its effects are minimal. But in fact, the period of Earth's axial precession is still a term that is worth considering compared to the period of Earth's magnetic field reversal. Methods This paper compares and analyzes different dynamo models with Ekman number and Rayleigh number by introducing Earth's axis precession velocity with a period of 25 960 years. Results It is found that Earth's axial precession stabilizes the kinetic energy and magnetic energy fluctuation of spherical shell magnetic fluid in a smaller range than the benchmark dynamo, and increases the toroidal kinetic energy of spherical shell magnetic fluid by more than 10%-20%, which leads to significantly accelerated westward drift rate of the magnetic field. Based on the magnetic field strength and magnetic Reynolds number at the core-mantle boundary of those dynamos, it is found that the introduction of the precession item of Earth's axis makes magnetic energy more tend to kinetic energy conversion. And by comparing the core-mantle boundary dipolarity, it is found that the introduction of the precession item will reduce the dipolarity. But for the dynamos in this article, the influence is not strong enough to make it become less Earth-like. Conclusions From the comparison and analysis of the dynamo models, the introduction of Earth's axis precession would possibly change the dynamo model with a small Ekman number to a multipole dominant dynamo model. It can be inferred that the term of Earth's axis precession should be considered in more precise research of the numerical dynamo simulation for more. -
Key words:
- geomagnetic field /
- geodynamo /
- numerical simulation /
- the Earth's axis precessio
-
表 1 模型1~4在6个粘性时间尺度的西向漂移圈数
Table 1. Circles of West Drifts of 4 Models in 6 Viscous Time Scales
模型 西向漂移圈数 1 32±3 2 49±3 3 35±3 4 44±3 -
[1] Baron G C. A Discourse on the Revolutions of the Surface of the Globe, and the Changes Thereby Produced in the Animal Kingdom[M]. Philadelphia: Lea&Blanchard, 1831 [2] Glatzmaier G A, Roberts P H. A Three-Dimensional Self-Consistent Computer Simulation of a Geomagnetic Field Reversal[J]. Nature, 1995, 377(6546): 203-209 doi: 10.1038/377203a0 [3] Christensen U R, Aubert J, Hulot G. Conditions for Earth-Like Geodynamo Models[J]. Earth and Planetary Science Letters, 2010, 296(3): 486-496 [4] Kageyama A, Sato T. Computer Simulation of a Magnetohydrodynamic DynamoⅡ[J]. Physics of Plasmas, 1995, 2(5): 1421-1431 doi: 10.1063/1.871485 [5] Kuang W, Bloxham J. An Earth-Like Numerical Dynamo Model[J]. Nature, 1997, 389: 371-374 doi: 10.1038/38712 [6] Stefani F, Albrecht T, Gerbeth G, et al. Towards a Precession Driven Dynamo Experiment[J]. Magnetohydrodynamics, 2014, 51(2): 275-283 [7] 邓洪涛. 地球内核与磁场[J]. 武汉大学学报·信息科学版, 2010, 35(7): 854-856 http://ch.whu.edu.cn/article/id/989 Deng Hongtao. The Inner Core and Geomagnetic Field[J]. Geomatics and Information Science of Wuhan University, 2010, 35(7): 854-856 http://ch.whu.edu.cn/article/id/989 [8] Glatzmaier G A. Introduction to Modelling Convection in Planets and Stars[M]. Princeton: Princeton University Press, 2013 [9] Christensen U R, Aubert J, Cardin P, et al. A Numerical Dynamo Benchmark[J]. Physics of the Earth and Planetary Interiors, 2001, 128(1): 24-34 [10] Thomas G, Johannes W, Ankit B, et al. Welcome—Magic 5.6 Documentation[EB/OL]. [2019-3-29]. https://magic-sph.github.io/ [11] Lauren W, Jessica I. Reconciling the Hemispherical Structure of Earth's Inner Core with Its Super-Rotation[J]. Nature Geoscience, 2011, 4(4): 264-267 doi: 10.1038/ngeo1083 [12] Wicht J. Inner-Core Conductivity in Numerical Dynamo Simulations[J]. Physics of the Earth and Planetary Interiors, 2002, 132(4): 281-302 doi: 10.1016/S0031-9201(02)00078-X -