并行自适应有限元软件平台phg 及其应用 - lsec.pdf

并并行行自自适适应应有有限限元元软软件件平平台台 PHG 及及其其应应用用 张林波, 郑伟英, 卢本卓, 崔涛, 冷伟 and 林灯 Citation: 中国科 : 信息科 46, 1442 (2016 ); doi: 10.1360/N112016-00066 View online: http://engine.s i hina. om/doi/10.1360/N112016-00066 View Table of Contents:http://engine.s i hina. om/publisher/s p/journal/SSI/46/10 Published by the 《中国科 》杂志社 Articles you may be interested in 第一原理实空间并行 自适应计算程序设计原理 中国科 : 信息科 46, 1421 (2016); 纳米集成电路互连线建模和光刻仿真中的大规模并行计算方法 中国科 : 信息科 46, 1372 (2016); 原子S hrǒdinger方程直接求解的相关函数-势谐-广义Laguerre函数方法 中国科 B辑: 化 26, 249 (1996); 一种以软件体系结构为中心的网构软件开发方法 中国科 F辑: 信息科 36, 1100 (2006); 虚拟川滇—基于千万网格并行有限元计算的区域强震演化过程数值模型设计和构建 中国科 D辑: 地球科 39, 260 (2009); 中国科学 信息科学 年 第 卷 第 期 SCIENTIA SINICA Informationis 高性能科学计算若干前沿问题研究专刊 并行自适应有限元软件平台 及其应用 * 张林波 郑伟英 卢本卓 崔涛 冷伟 林灯 中国科学院数学与系统科学研究院科学与工程计算国家重点实验室(LSEC), 北京 100190 * 通信作者. E-mail: zlb@ 收稿日期: 2016–03–29; 接受日期: 2016–07–06; 网络出版日期: 2016–10–25 国家自然科学基金(批准号:、国家重点基础研究发展计划 (973 计划) (批 准号: 2011CB309703)、国家高技术研究发展计划 (863 计划) (批准号: 2012AA01A30901)、国家磁约束聚变能发展研究专项 (批 准号: 2015GB110003) 和中国科学院国家数学与交叉科学研究中心 (NCMIS) 资助项目 摘要 是一个由中国科学院 科学与工程计算国家重点实验室 研制 的开源并行自适应有限元程序开发平台 与国际上其他同类型软件相比 其主要特征有 提供基 于最新顶点单元二分的非结构协调四面体网格并行局部加密 支持 自适应有限元计算 支 持大规模并行和动态负载平衡 提供大型稀疏线性方程组直接求解和灵活的预处理迭代求解 本 文介绍 平台的主要模块和核心算法 以及基于 平台开发的一些并行有限元应用程序 它 们包括 大型变压器铁损模拟、离子通道中的离子输运过程模拟、集成电路互连线寄生参数提取、冰 盖模拟、弹性波 方程谱元计算等 关键词 有限元 自适应网格 并行计算 离子通道 电磁涡流 寄生参数提取 冰盖模拟 地震 波模拟 引言 有限元方法是求解偏微分方程 (组) 的重要数值方法, 它通过将求解区域划分成由简单几何形状 的单元构成的网格, 并在每个单元中用简单函数(通常是多项式) 构成的基函数对解进行近似, 从而将 偏微分方程 (组) 转化为代数方程组来实现数值求解. 在有限元方法中, 影响计算结果精度的主要因素 [1, 2] 是网格单元的分布和基函数. 自适应有限元方法最初由Babuˇska 等提出 , 它提供了一套系统的方 法, 在求解过程中自动根据数值解的性质反复对网格单元的分布及基函数进行调整, 从而逐步形成适 合于所计算的问题的拟最优的网格. 以定常问题为例, 自适应有限元方法从一个初始网格出发,