一种求解瞬时涡流问题的新型有限元算法.pdfVIP

第14卷，第2期 中国传媒大学学报自然科学版 Vo1．14．No．2 2007年6月 JOURNAL OF COMMUNICATION UNIVERSITY OF CHINA SCIENCE AND TECHNOLOGY June，2007 一 种求解瞬时涡流问题的新型有限元算法 李建军 (河北省承德市承德旅游学院基础部，河北067000) 摘 要：本文提出了一种求解瞬时涡流问题的新型有限元算法，并讨论了该算法的误差估计。此方法不仅更真实 地模似物理背景条件，而且其计算格式又是一种解耦的形式，减少了计算量。 关键词：瞬时涡流问题；节点有限元法，中心差分Crank—Nicholson；误差分析。 中图分类号：0241．81，0242．21 文献标识码：A 文章编号：1673—4793(2007)02—0023—07 A New Algorithm of the Nodal Finite Element Method for a Trancient Eddy Current Problem LI Jian-jun (Department of Basic Sciences，Chengde Tourism College，Chengde 067000，China) Abstract：A new decoupled finite element method is suggested to approximate transcient eddy current equations in a three—dimensional polyhedral domain．This method is based on solving a vector and a scalar from the splitting of the electric field by using nodal finite elements．An optimal energy—norm error estimate in finite time is obtained by introducing a projection operator． · Key words：eddy current problem；nodal finite element method；C—N difference scheme；nodal elements， error est m ate 1 引言 涡流问题出自舍人位移电流的Maxwell方程组，即 curlH E+Js， ， J curIE：一l ． (1)1 ( H)：。 这意指针对低频和高导的情况通过忽略电容的影响(或忽略空间电荷)，来提供一种合理的、模拟原问题的 近似模型。与[5]中不同，我们通过引入罚函数，提供了一种采用节点有限元法、基于求解电势与磁势的A 一 法。此方法不仅更真实地模拟物理背景条件，而且其计算格式又是一种解耦的形式。它保证了有限元 解的唯一性，同时又不会造成更多的计算复杂性。本文提出了在时间上基于中心差分Crank—Nicholson的 计算格式，并给出了其误差分析。 收稿日期：2006—12—10 作者简介：李建军(1967一)，男(汉族)，河北承德人，承德旅游学院讲师 中国传媒大学学报自然科学版 第14卷 2 电场分解公式 为了简化数学分析起见，令Y2CR 代表导体所在有界的、连续的多边形区域，其边界为 。我们给出如 下与基本涡流模型的物理背景一致的假设： ( )假设材料参数为线性且不依赖时间，即电导率or