Fast Solution of 3-D Eddy-Current Problems in Multiply Connected Domains by a, v-φ and t-φ Formulations With Multigrid-Based Algorithm for Cohomology Generation

被引:0
作者
Moro, Federico [1 ]
Napov, Artem [2 ]
Pellikka, Matti [3 ]
Smajic, Jasmin [4 ]
Codecasa, Lorenzo [5 ]
机构
[1] Univ Padua, Dipartimento Ingn Ind, I-35131 Padua, Italy
[2] Univ Libre Bruxelles, Serv Metrol Nucl, B-1050 Brussels, Belgium
[3] Grundium Ltd, Tampere 33720, Finland
[4] Swiss Fed Inst Technol, Swiss Fed Inst Technol, Inst Electromagnet Fields, CH-8092 Zurich, Switzerland
[5] Politecn Milan, Dipartimento Elettron Informaz & Bioingn, I-20133 Milan, Italy
关键词
Eddy currents; AC problem; finite element method; multiply connected; electromagnetic; multigrid; cohomology; OMEGA FORMULATION; CELL METHOD; COMPUTATION; HOMOLOGY; CUTS;
D O I
10.1109/ACCESS.2022.3216876
中图分类号
TP [自动化技术、计算机技术];
学科分类号
0812 ;
摘要
The fast solution of three-dimensional eddy current problems is still an open problem, especially when real-size finite element models with millions of degrees of freedom are considered. In order to lower the number of degrees of freedom a magnetic scalar potential can be used in the insulating parts of the model. This may become difficult when the model geometry presents some conductive parts which are multiply connected. In this work a multigrid-based algoritm is proposed that allows for a calculation in linear-time of cohomology, which is needed to introduce the scalar potential without cuts. This algorithm relies on an algebraic multigrid solver for curl-curl field problems, which ensures optimal computational complexity. Numerical results show that the novel algorithm outperforms state-of-the-art methods for cohomology generation based on homological algebra. In addition, based on this algoritm, novel a, v-phi and t-phi formulations to analyze three-dimensional eddy current problems in multiply connected domains are proposed. Both formulations, after discretization by the cell method, lead to a complex symmetric system of linear equations amenable to fast iterative solution by Krylov-subspace solvers. These formulations are able to provide very accurate numerical results, with a minimum amount of degrees of freedoms to represent the eddy current model. In this way the computational performance is improved compared to the classical A, V-A formulation typically implemented in finite element software for electromagnetic design.
引用
收藏
页码:112416 / 112432
页数:17
相关论文
共 33 条
  • [1] Edge element formulations of eddy current problems[J]. Bíró, O. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1999(3-4)
  • [2] Explicit, consistent, and conditionally stable extension of FD-TD to tetrahedral grids by FIT[J]. Codecasa, Lorenzo;Politi, Marco. IEEE TRANSACTIONS ON MAGNETICS, 2008(06)
  • [3] Cuts for the magnetic scalar potential in knotted geometries and force-free magnetic fields[J]. Crager, JC;Kotiuga, PR. IEEE TRANSACTIONS ON MAGNETICS, 2002(02)
  • [4] Douglas N., 2018, FINITE ELEMENT EXTER
  • [5] Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities[J]. Geuzaine, Christophe;Remacle, Jean-Francois. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2009(11)
  • [6] An algorithm to construct the discrete cohomology basis functions required for magnetic scalar potential formulations without cuts[J]. Henrotte, F;Hameyer, K. IEEE TRANSACTIONS ON MAGNETICS, 2003(03)
  • [7] Current and voltage excitations for the eddy current model[J]. Hiptmair, R;Sterz, O. INTERNATIONAL JOURNAL OF NUMERICAL MODELLING-ELECTRONIC NETWORKS DEVICES AND FIELDS, 2005(01)
  • [8] Auxiliary space preconditioning for edge elements[J]. Hiptmair, Ralf;Xu, Jinchao. IEEE TRANSACTIONS ON MAGNETICS, 2008(06)
  • [9] International Compumag Society, PROBL 3 BATH PLAT 2
  • [10] Kameari A., 1989, Applied Electromagnetics in Materials, P225