Explicit Time Integration of Transient Eddy Current Problems

For time integration of transient eddy current problems commonly implicit time integration methods are used, where in every time step one or several nonlinear systems of equations have to be linearized with the Newton-Raphson method due to ferromagnetic materials involved. In this paper, a generalized Schur-complement is applied to the magnetic vector potential formulation, which converts a differential-algebraic equation system of index 1 into a system of ordinary differential equations (ODE) with reduced stiffness. For the time integration of this ODE system of equations, the explicit Euler method is applied. The Courant-Friedrich-Levy (CFL) stability criterion of explicit time integration methods may result in small time steps. Applying a pseudo-inverse of the discrete curl-curl operator in nonconducting regions of the problem is required in every time step. For the computation of the pseudo-inverse, the preconditioned conjugate gradient (PCG) method is used. The cascaded Subspace Extrapolation method (CSPE) is presented to produce suitable start vectors for these PCG iterations. The resulting scheme is validated using the TEAM 10 benchmark problem.

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

[2]  S. Schops,et al.  Higher Order Half-Explicit Time Integration of Eddy Current Problems Using Domain Substructuring , 2011, IEEE Transactions on Magnetics.

[3]  André Nicolet,et al.  Implicit Runge-Kutta methods for transient magnetic field computation , 1996 .

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

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

[6]  Akihisa Kameari,et al.  Calculation of transient 3D eddy current using edge-elements , 1990 .

[7]  T. Weiland,et al.  Subspace projection extrapolation scheme for transient field simulations , 2004, IEEE Transactions on Magnetics.

[8]  Zhuoxiang Ren Influence of the RHS on the convergence behaviour of the curl-curl equation , 1996 .

[9]  Traianos V. Yioultsis Explicit finite-difference time-domain method with perfectly matched layers for transient low-frequency eddy-current analysis , 2002 .

[10]  Thomas Weiland,et al.  Transient eddy-current calculation with the FI-method , 1999 .

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

[12]  Oszkar Biro,et al.  Discontinuous Galerkin formulation for eddy‐current problems , 2008 .

[13]  Koji Fujiwara,et al.  RESULTS FOR BENCHMARK PROBLEM 10 (STEEL PLATES AROUND A COIL) , 1990 .