有限单元法的研究全球地幔对流主要影响因素.pdfVIP

有限单元法研究全球地幔对流主要影晌因素 张怀朱桂芝 (中科院研究生院计算地球动力学实验室，北京，100049) 近些年来的地震学研究揭示的全球内部结构，对地幔对流研究提供了一批新的、更加可靠的数据，使人们 可以更有把握的选择物质参数和边界条件，也可以对计算结果提供更强的磐L嫩。地震学观测结果表明，固体地 球内部的确存在对流。而且对流并不是简单的地球化学坚持的分层地幔对流、或地球动力学家所坚持的全地幔 对流模型。因为从实际的地震波层析成像结果来看，确实存在从地球表面一直到达核幔边界的斜插俯冲板片， 但也存在界面对对流的阻滞现象。基于大规模数值分析的结果来看，在某些条件下，确实能够出现以670公啦 相变介面为主的分层地幔对流形式的可能。 全球地幔对流的整体格局究竟是什么样的，整体和局部地幔对流究竟都有哪些具体的形式，以及哪些关键 因素会影响地幔对流等问题，己成为了当今地球科学研究中的一个争论性课题。目前有地球化学和地热学家们 大力支持的分层地幔对流和地震学及地球动力学研究广泛支持的全地幔对流模式。为了协调二者之间的矛盾， 也有一些新模型被提出，如团块模型、拉瓦灯模型和过渡区水过滤模型，随着计算机技术的发展，数值模拟得 al 到长足的进步。(Matyska，Yuen，et 放热钙钛矿相情况下的分层地幔对流，并考虑依赖温度和与辐射热导率有关的依赖深度的变粘滞系数对上升流 和下降流的动力学的影响。Zhong等(2000)在球坐标下研究了变粘度对全地幔对流的结构形态和尺度的重要 影响。本文主要研究大尺度地幔对流模式结构的影响因素，在此主要介绍我们的模型和计算方法。 我们采用基于算子分裂算法的有限元方法，对流体速度压力场和温度场采用耦合算法，使用符合LBB条件 的有限元单元模型，建立一整套关于地幔对流的全新模型，综合各种文献种提出的影响全球地幔对流的关键因 素，对全球地幔对流进行了大规模数值模拟。假设Boussinesq近似成立，地幔中的热对流遵循质量、动量和能 量守恒方程，具体表达如下： 甜，．，=0 一Pj+(qui。j+qujj、)。J+pp96|r=o(2) 可，+“，T，=KT^+H (3) 其中“，是流速，尸是动态压力，印是密度异常，g是重力加速度，77是粘滞系数，71是温度，盯是热 扩散率，H是产热率，X，，表示X对Y的倒数，f和，是空间标号，，是时间，重复的标号表示求和，密度 差异包括温度变化和相变影响因素，户。为参考密度，瓦为表面温度。本文主要针对依赖温度和深度的变粘度 系数、内部产热率和相变对地幔对流模式的影响。 ’ 院知识创新工程信息化建设专项(INFl05一scE一02一12)等资助。 参考文献 vander for Hilst，R．D．，S．Widiyantoro，andE．R．Engdahl，Evidencemantlecirculationfrom deep globaltomography，Nature，386，578—584，1997． influenceof C．Matyska，D．A．Yuen．Thetemperature—dependentandradiativethermal viscosity lowermantle on withthe Planet．SciLett． conductivity dynamics post—perovskitephasechange，Earth inpress(2004)． MT．Roleof andsurfacein ZhongS，Zuber temperature—dep