×

使用边缘元素的三维静磁问题的高效预处理迭代线性解算器。 (英语) Zbl 1488.65070号

摘要:对于三维(3-D)大规模静磁问题的数值计算,迭代求解器是首选的,因为使用稀疏直接求解器需要大量内存。本文采用最近提出的用于静磁问题的库仑磁矢势(MVP)公式,利用边缘单元进行有限元离散,其中所得到的线性系统是对称的,但条件不好。为了有效地求解此类线性系统,我们利用迭代Krylov子空间解算器,构造了三个新的块预条件器,它们是从传统的块Jacobi、Gauss-Seidel和约束预条件器导出的。文中还讨论了所提出的预条件子的谱特性和实际实现细节。然后,给出了实际仿真的数值例子,以说明所提方法的效率和准确性。

MSC公司:

65层10 线性系统的迭代数值方法
65F08个 迭代方法的前置条件
65N30型 偏微分方程边值问题的有限元、Rayleigh-Riz和Galerkin方法
78A30 静电和磁力静力学
78M10个 有限元、伽辽金及相关方法在光学和电磁理论问题中的应用
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 赵永明和傅文华,一种新的库仑规矢量磁位公式,用于使用边缘元素的三维涡流场分析,IEEE Trans。马格纳。,53(6) (2017), 9400704.
[2] 赵永明,傅文华,利用带边缘元素的矢量磁势在时域中实现鲁棒全波麦克斯韦解算器,IET Sci。测量。技术。,11(6)(2017),第746-752页。
[3] 赵彦和唐振中,三维瞬态全波麦克斯韦问题的对称场路耦合公式,IEEE Trans。马格纳。,55(6) (2019), 7201804.
[4] 何顺良、赵义勇、傅文南和周培平,边元在考虑电感和电容效应的三维电磁场分析中的应用,IEEE Trans。马格纳。,52(3) (2016), 7400504.
[5] 赵永明和傅文南,一种新的适用于所有频率的稳定全波麦克斯韦解算器,IEEE Trans。马格纳。,53(6) (2017), 7200704.
[6] 赵云霞、何世浩、傅文武,一种新型快速无网格变形方法及其在电磁器件优化设计中的应用,IEEE Trans。马格纳。,50(11)(2014),7200904。
[7] 赵永明和傅文南,利用边缘元素求解三维静磁问题的库仑规范新公式,IEEE Trans。马格纳。,53(6) (2017), 9100104.
[8] Y.ZHAO和Z.TANG,基于对偶A和H公式的开放边界静磁问题的改进平衡误差估计,IEEE Trans。马格纳。,55(6) (2019), 7201105.
[9] 赵勇和唐振中,使用无需源场计算的对偶公式精确提取绕组电感,IEEE Trans。马格纳。,55(6) (2019), 7201504.
[10] C.GAIER和H.HAAS,使用H公式对静磁场和瞬态涡流场进行边缘元素分析,IEEE Trans。马格纳。,31(3)(1995),第1460-1463页。
[11] A.BERMüDEZ、R.RODRíGUEZ和P.SALGADO,用标量势解决静磁问题的有限元方法,SIAM J.Numer。分析。,46(3)(2008),第1338-1363页·兹比尔1171.78010
[12] 李玉立,孙松琪。I.DAI和W.C.CHEW,用广义库仑规范求解静磁问题的双旋度方程的矢量解,IEEE Trans。马格纳。,51(8) (2015), 7002306.
[13] H.BATELAAN和A.TONOMURA,《阿哈罗诺夫·博姆效应:微妙主题的变化》,《物理学》。今天,62(9)(2009),第38-43页。
[14] O.BíRóAND K.PREIS,《磁场矢量势在三维涡流有限元分析中的应用》,IEEE Trans。马格纳。,25(4)(1989),第3145-3159页。
[15] K.PREIS、I.BARDI、O.BíR O.、C.MAGELE、G.VRISK和K.R.RICHTER,三维静磁场的不同有限元公式,IEEE Trans。马格纳。,28(2)(1992),第1056-1059页。
[16] O.BíR O.,K.PREIS和K.R.RICHTER,《关于在三维静磁问题的节点和边缘有限元分析中使用磁矢量势》,IEEE Trans。马格纳。,32(3)(1996),第651-654页。
[17] T.A.DAVIS、S.RAJAMANICKAM和W.M.SID-LAKHDAR,稀疏线性系统直接方法综述,《数值学报》。,25(2016),第383-566页·Zbl 1346.65011号
[18] P.R.AMESTOY,T.A.DAVIS和I.S.DUFF,算法837:AMD,近似最小度排序算法,ACM Trans。数学。软质。,30(3)(2004),第381-388页·Zbl 1070.65534号
[19] P.DI BARBA、I.PERUGIA和A.SAVINI,磁场二维模拟混合有限元的最新经验,COMPEL,17(5)(1998),第674-681页·Zbl 0932.78011号
[20] 陈振英、徐永和邹建平,一种自适应边缘元方法及其在磁静力学鞍点问题中的收敛性,数值。方法。偏微分方程,28(5)(2012),第1643-1666页·Zbl 1259.78039号
[21] Y.SAAD,《稀疏线性系统的迭代方法》,第2版,SIAM,宾夕法尼亚州费城,2003年·Zbl 1031.65046号
[22] X.-M.GU,T.-Z.HUANG,L.LI,H.-B.LI,T.SOGABE和M.CLEMENS,电磁模拟中复杂对称线性系统的COCG和COCR方法的准最小残差变体,IEEE Trans。微型。《理论技术》,62(12)(2014),第2859-2867页。
[23] X.-M.GU,M.CLEMENS,T.-Z.HUANG和L.LI,复杂对称线性系统的SCBiCG算法类及其在几个电磁模型问题中的应用,计算。物理学。社区。,191(2015),第52-64页·兹比尔1344.65041
[24] I.PERUGIA、V.SIMONCINI和M.ARIOLI,静磁问题混合近似中的线性代数方法,SIAM J.Sci。计算。,21(3)(1999),第1085-1101页·兹比尔1040.78012
[25] C.C.PAIGE和M.A.SAUNDERS,稀疏线性方程组的求解,SIAM J.Numer。分析。,12(4)(1975年),第617-629页·Zbl 0319.65025号
[26] R.W.FREUND和N.M.NACHTIGAL,简化Lanczos和QMR算法软件,应用。数字。数学。,19(3)(1995),第319-341页·Zbl 0853.65041号
[27] P.DI BARBA、A.SAVINI和I.PERUGIA,电磁装置中场和力模拟的混合有限元,IEEE Trans。马格纳。,34(5)(1998),第3572-3575页。
[28] P.ALOTTO和I.PERUGIA,《静磁学中的树-树隐式凝聚》,IEEE Trans。马格纳。,36(4)(2000),第1523-1526页。
[29] I.PERUGIA,二维静磁问题的基于现场的混合公式,SIAM J.Numer。分析。,34(6)(1997),第2382-2391页·Zbl 0896.65087号
[30] 曹振华,非对称鞍点问题的约束Schur补预条件,应用。数字。数学。,59(1)(2009),第151-169页·Zbl 1161.65037号
[31] I.PERUGIA和V.SIMONCINI,混合有限元公式的块对角和不定对称预条件,数值。线性代数应用。,7(7-8)(2000),第585-616页·Zbl 1051.65038号
[32] J.CHEN,Z.CHEN,T.CUI和L.ZHANG,带电路/场耦合的涡流模型的自适应有限元方法,SIAM J.Sci。计算。,32(2)(2010),第1020-1042页·Zbl 1211.78029号
[33] E.CHOW和Y.SAAD,不定矩阵ILU预条件的实验研究,J.Compute。申请。数学。,86(2)(1997),第387-414页·Zbl 0891.65028号
[34] T.A.MANTEUFFEL,正定线性系统的不完全因式分解技术,数学。计算。,34(150)(1980),第473-497页·Zbl 0422.65018号
[35] D.K.SALKUYEH,M.MASOUDI和D.HEZARI,关于鞍点问题的广义移位分裂预dictioner,Appl。数学。莱特。,48(2015),第55-61页·Zbl 1325.65048号
[36] F.P.ALI BEIK,M.BENZI和S.-H.AZIZI CHAPARPORDI,稳定鞍点问题的块对角和块三角迭代格式及预条件,J.Compute。申请。数学。,326(2017),第15-30页·Zbl 1370.65015号
[37] C.-J.LIN和J.J.MORé,记忆有限的不完全Cholesky因子分解,SIAM J.Sci。计算。,21(1)(1999),第24-45页·兹伯利0941.65033
[38] 张彦,黄忠忠,金义峰,李丽丽,用多参数的柔性不完全Cholesky因子分解来控制预条件中非零元素的个数,Numer。线性代数应用。,19(3)(2012),第555-569页·兹比尔1274.65088
[39] Y.SAAD和M.H.SCHULTZ,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7(3)(1986),第856-869页·Zbl 0599.65018号
[40] S.C.EISENSTAT、H.C.ELMAN和M.H.SCHULTZ,非对称线性方程组的变分迭代方法,SIAM J.Numer。分析。,20(2)(1983年),第345-357页·Zbl 0524.65019号
[41] P.JIRáNEK和M.ROZLOZNíK,更简单GMRES的自适应版本,Numer。《算法》,53(1)(2010),第93-112页·Zbl 1188.65033号
[42] X.-M.GU,T.-Z.HUANG,B.CARPENTIERI,L.LI和C.WEN,求解非厄米线性系统的BiCORSTAB和GPBiCOR方法的混合迭代算法,Com-put。数学。申请。,70(12)(2015),第3019-3031页·Zbl 1443.65042号
[43] E.丹尼森,有限螺线管的轴向磁场,2005年。在线提供地址:https://tiggerntatie.github.io/emagnet/solenoids/solenoidonaxis.htm。
[44] G.AKOUN和J.-P.YONNET,两个立方磁体之间作用力的三维解析计算,IEEE Trans。马格纳。,20(5)(1984),第1962-1964页。
[45] V.SIMONCINI和D.B.SZYLD,不精确Krylov子空间方法理论及其在科学计算中的应用,SIAM J.Sci。计算。,25(2)(2003),第454-477页·兹比尔1048.65032
[46] Y.XI,R.LI和Y.SAAD,稀疏对称矩阵的低秩修正代数多层预条件器,SIAM。《矩阵分析杂志》。申请。,37(1)(2016),第235-259页·Zbl 1376.65036号
[47] P.GATTO和J.S.HESTHAVEN,通过分层低阶近似对hp-FEM矩阵进行有效预处理,科学杂志。计算。,72(1)(2017),第49-80页·Zbl 1371.65028号
[48] M.BEBENDORF,《层次矩阵:有效解决椭圆边值问题的方法》,《计算科学与工程讲义》,第63卷,施普林格-弗拉格出版社,柏林-海德堡,2008年·Zbl 1151.65090号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。