优化地球三维有限元计算网格的一种方法-兽类学报.pdfVIP

第 卷第 期 中国科学院研究生院学报 ２３ ５ Ｖｏｌ．２３ Ｎｏ．５ 年 月 ２００６ ９ ＪｏｕｒｎａｌｏｆｔｈｅＧｒａｄｕａｔｅＳｃｈｏｏｌｏｆｔｈｅＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ Ｓｅｐｔｅｍｂｅｒ ２００６  文章编号： （ ） １００２１１７５２００６０５０６７６０５ 简 报 优化地球三维有限元计算网格的一种方法 ， １２ １ １ 朱桂芝 张 怀 石耀霖 （中国科学院研究生院计算地球动力学实验室，北京 ； 中国地震局第二监测中心，西安 ） １ １０００４９ ２ ７１００５４ （ 年 月 日收稿； 年 月 日收修改稿） ２００５ ９ ８ ２００５ １２ ６ ， ， ＺｈｕＧＺ ＺｈａｎｇＨ ＳｈｉＹＬ．Ａｍｅｔｈｏｄｏｆｍｅｓｈｇｅｎｅｒａｔｉｏｎｆｏｒ３Ｄｆｉｎｉｔｅｅｌｅｍｅｎｔｅａｒｔｈｃｏｍｐｕｔａｔｉｏｎ．ＪｏｕｒｎａｌｏｆｔｈｅＧｒａｄｕａｔｅ ， ， （）： ＳｃｈｏｏｌｏｆｔｈｅＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ２００６２３５ ６７６～６８０ 摘 要 在研究大尺度的三维地球动力学问题的数值模拟过程中，往往需要对整个球体或部 分球壳生成计算网格，而且网格单元的尺寸对所研究的问题尤为重要 计算网格既满足地球 ． 物理分析的需要，又要避免内部网格过小的浪费 本文提出了一种地球结构化网格的生成和 ． 加密划分方案，可以对整个球体形成较为均匀的网格，并能够对近地表的区域完成网格加密以 达到地球物理分析的需要 生成的网格在简单的试算中，得到了优异的结果 ． ． 关键词 地球，结构化网格单元，有限单元法 中图分类号 Ｐ５４１ １ 引言 在诸如地球发电机、地幔对流的数值模拟和全球地震波传播的谱元法计算等大尺度的三维地球动 力学问题的研究中，都会遇到球面作为边界的网格生成问题 二维有限元地幔对流计算中往往简单地 ． ［］ １ 采用沿弧向和径向均匀剖分 ，如果采用简单的方法按经纬度划分网格，会造成在两极处网格极密而赤 ［］ ２ 道处网格稀疏的状况，这种情况会引起坐标的奇异性问题和近两极网格的收敛性问题 因此，三维地 ． ［，］ ［］ ［］ ３４