Higher Order Half-Explicit Time Integration of Eddy Current Problems Using Domain Substructuring

In this paper domain substructuring is adapted to the nonlinear transient eddy current problem: conductive and nonconductive domains are separately treated for a more efficient time integration. A matrix factorization of the linear (nonconductive) subproblem, e.g., air, is executed on beforehand and used throughout the simulation. For a general 3D problem these non-sparse factors must be replaced by sparse approximations or explicit time integration must be carried out to increase the efficiency of the solution process. This approach is validated using implicit and half-explicit time integration methods. Numerical results underline the feasibility and the possible gain obtained by higher order half-explicit methods.

[1]  K. Preis,et al.  Discontinuous Galerkin Finite Elements in Time Domain Eddy-Current Problems , 2009, IEEE Transactions on Magnetics.

[2]  Linda R. Petzold,et al.  A stabilized explicit Lagrange multiplier based domain decomposition method for parabolic problems , 2008, J. Comput. Phys..

[3]  Markus Clemens,et al.  GPU Accelerated Adams–Bashforth Multirate Discontinuous Galerkin FEM Simulation of High-Frequency Electromagnetic Fields , 2010, IEEE Transactions on Magnetics.

[4]  J. Verwer Explicit Runge-Kutta methods for parabolic partial differential equations , 1996 .

[5]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[6]  Kay Hameyer,et al.  Solution strategies for transient, field-circuit coupled systems , 2000 .

[7]  Rüdiger Weiner,et al.  Half-explicit Runge-Kutta methods for semi-explicit differential-algebraic equations of index 1 , 1993 .

[8]  Luc Giraud,et al.  Iterative versus direct parallel substructuring methods in semiconductor device modelling , 2005, Numer. Linear Algebra Appl..

[9]  Sebastian Schöps,et al.  Dynamic Iteration for Coupled Problems of Electric Circuits and Distributed Devices , 2013, SIAM J. Sci. Comput..

[10]  Herbert De Gersem,et al.  Decomposition and regularization of nonlinear anisotropic curl‐curl DAEs , 2011 .

[11]  Mario Bebendorf,et al.  Hierarchical Matrices: A Means to Efficiently Solve Elliptic Boundary Value Problems , 2008 .

[12]  Peter Deuflhard,et al.  Domain decomposition with subdomain CCG for material jump elliptic problems , 1997 .

[13]  Igor Tsukerman,et al.  A survey of numerical methods for transient eddy current problems , 1993 .

[14]  Theodoros D. Tsiboukis,et al.  A finite difference time domain scheme for transient eddy current problems , 2001 .

[15]  Joachim Schöberl,et al.  Numerical analysis of nonlinear multiharmonic eddy current problems , 2005, Numerische Mathematik.