1 Introduction
Recently, efficient and robust designs of electromechanical energy converters are gaining again in importance because of the transition towards sustainable energy in Europe (‘Energiewende’ in German). Electrical machinery is well understood and developed in industry close to their technical limits, but often without transient analysis or consideration of uncertainties in the design process. Such studies are only carried out late in the development process due to their high computational costs. This may lead to the fact that better or more robust designs are ruled out and not considered further on. One promising way to speed up transient analysis are parallelintime methods.
In contrast to classical timeintegration techniques based on a timestepping approach, i. e., solving sequentially for one time step after the other, parallelintime algorithms allow simultaneous solution across multiple time steps. Starting with the work of Nievergelt Nivergelt1964 , various approaches for parallelintime integration have been explored; a recent review of the extensive literature in this area is Gander2015_Review . The key practical aspect for choosing one of the many timeparallel methods when aiming at adding parallelism to an existing application code is the level of intrusiveness, i. e., the required amount of implementation effort. There are only a few timeparallel methods that are nonintrusive. In this paper, we consider two of these approaches, the Parareal method Lions_etal2001 and the multigridreductionintime (MGRIT) algorithm Falgout_etal_2014 that, in a specific twolevel setting, can be viewed as a Pararealtype algorithm.
2 Eddy current model problem
For an open, bounded domain and , the evolution of electromagnetic fields is governed by Maxwell’s equations on Jackson_1998aa
(1) 
with suitable initial and boundary conditions at time and , respectively. These equations are completed by constitutive relations laws
(2) 
In these equations, is the magnetic field [A/m], the magnetic flux density [T], the electric field [V/m], the electric flux density [C/m], and are the total and source current density [A/m], is the electric charge density [C/m]. All fields are functions of space and time . The material properties , and are the electric conductivity, the electric permittivity and the magnetic permeability, respectively. It is convenient to invert the magnetic material law, i. e., , using the reluctivity , where can be a sufficiently smooth and bounded function of the magnitude , see Heise_1994aa . In the following, we consider only devices where the displacement current can be neglected with respect to the source currents, i. e., . An analysis of this error can be found in Schmidt_2008aa . Assuming yields the socalled magnetoquasistatic approximation or eddy current problem. Eddy currents lead to the skin effect, i. e., currents through a conductor are pushed to the surface if frequency increases (Jackson_1998aa, , Chapter 5.18).
One may introduce the (‘modified’ Emson_1988aa
) magnetic vector potential
such that . Then, inserting the equations into each other yields(3) 
The source is defined as in the inner cable and vanishes elsewhere, denotes the unit vector in direction and the excitation is given by
(4) 
where is the common sawtooth pattern, with teeth and period s KulchytskaRuchka_2018ac . The reluctivity is modeled as vacuum () in and , and is given in by a Spline curve, the conductivity is only nonzero in the tube region ( MS/m).
Finally, the RitzGalerkin approach is employed with lowest order ansatz functions in space. When considering planar 2D problems, edge shape functions only have a component and can be constructed from the nodal shape functions as
(5) 
where refers to the length in direction, . This leads to the equation
(6) 
with the matrices and the righthand side
respectively. The resulting system (6) consists of differentialalgebraic equations of index1 due the vanishing entries of the mass matrix in and , Nicolet_1996aa .
3 Multigrid reduction in time
The multigridreductionintime (MGRIT) algorithm Falgout_etal_2014 is an iterative, parallel method, based on applying multigrid reduction (MGR) MRies_UTrottenberg_1979 ; MRies_etal_1983 principles in time, for solving timestepping problems of the form
(7) 
with initial condition, , at . Note that form (7) can be a system of ODEs, arising, for example, after spatial discretization of a spacetime PDE, or it can be a system of DAEs such as given in Equation (6). Discretizing the time interval on a grid , , with, for notational convenience, constant time step , let be an approximation to for , and let . Then, considering a onestep timeindependent time integration method with timestepping operator, , that takes a solution at time to that at time , along with a timedependent forcing term, , the solution to (7) is defined via timestepping, which can also be represented as a forward solve of the linear system, written in block form as
(8) 
Hence, in the time dimension, this forward solve is completely sequential.
Alternatively, considering the lower block bidiagonal structure, we could apply cyclic reduction, whereby we first solve the Schur complement system,
(9) 
for the value of the solution at every th temporal point, with consistently restricted forcing terms. Then define the solution at the remaining temporal points by local (and parallel) timestepping between those points defined from the Schur complement. MGRIT is based on interpreting this cyclic reduction approach as a twolevel MGR algorithm, enabling parallelism in the solution process (8). Therefore, define a coarse temporal mesh, or (using multigrid terminology) the set of Cpoints, to be those points included in the Schur complement system (9
), with the remaining temporal points as the set of Fpoints. Further define “ideal” interpolation as the map which takes the solution at the Cpoints and yields a zero residual at the Fpoints, with a similar definition for “ideal” restriction. The Schur complement then arises as the standard PetrovGalerkin coarsegrid operator with these definitions of restriction and interpolation. Cyclic reduction can be viewed as a twolevel method with this PetrovGalerkin coarsegrid operator and a block smoother (called Frelaxation) that converges in one iteration. As it is typical in the MGR setting, the MGRIT approach replaces the true Schur complement with a simpler operator (typically of the same form as the original bidiagonal system, but with a time propagator using timestep
), replaces ideal restriction with simple injection, and compensates by adding relaxation. Furthermore, the twolevel method can be extended to multiple levels in a simple recursive manner, and the full approximation storage (FAS) approach ABrandt_1977b can be used to accommodate nonlinear problems.4 Numerical results
We apply classical sequential time stepping and two MGRIT variants to the eddy current model problem (3) with the pulsed excitation (4) on the spacetime domain , with depicted in Fig. 0(b). The spatial domain, , is discretized using degrees of freedom and the backward Euler method is used on a uniform grid with time steps for the time derivative of the spacediscrete timestepping problem (6). The time step on the finest grid, , is chosen to be , and the time step on each coarse grid, , is given by , . Two MGRIT variants are considered: a twolevel Pararealtype method with a coarsening factor of , and a fivelevel method that coarsens uniformly across all grids with a factor of . Thus, the coarsest grid consists of points in time in both cases. On this coarsest temporal grid, time stepping is used. All spatial problems are solved using a direct LU solver.
The MGRIT algorithm was implemented in parallel using Python and Message Passing Interface (MPI). Numerical results were generated on an Intel Xeon Phi cluster consisting of 272 1.4 GHz Intel Xeon Phi processors.
Fig. 3 shows convergence of the two MGRIT variants applied to the eddy current model problem. We see linear convergence for both variants. Comparing the number of spatial timestepping solves required for the two methods to the optimal count of for sequential time stepping, we note that one iteration requires about or spatial solves, respectively, when considering the twolevel Pararealtype method or the fivelevel MGRIT scheme. This large computational overhead is demonstrated in the strong scaling results in Fig. 3. The dotted and solid lines show results for the two and fivelevel methods, respectively, for increasing the number of processors in the temporal dimension only. The dashed line shows the runtime of time stepping on one processor for reference purposes. Results show that the extra work in the MGRIT variants can be effectively parallelized at high processor counts, i. e., more than , with good strong parallel scaling with a speedup of up to a factor of about over sequential time stepping.
5 Conclusions
MGRIT was applied for the first time to the eddy current problem, which yields an index1 DAE after spatial discretization. A speedup of approximately three times could be obtained. A strong scaling investigation shows that the method converges linearly with the number of processors, even for nonstandard, pulsed righthand sides, which has been shown to be problematic for classical Parareal KulchytskaRuchka_2018ac .
Acknowledgements.
The work is supported by the Excellence Initiative of the German Federal and State Governments, the Graduate School of Computational Engineering at TU Darmstadt, and the BMBF in the framework of project PASIROM (grants 05M18RDA and 05M18PXB).Bibliography

(1)
Nievergelt, J.: Parallel methods for integrating ordinary differential equations. Comm. ACM
7, 731–733 (1964)  (2) Gander, M.J.: 50 years of Time Parallel Time Integration. In: Carraro, T., Geiger, M., Körkel, S., Rannacher, R. (eds.) Multiple Shooting and Time Domain Decomposition, pp. 69–113. Springer (2015)
 (3) Lions, J.L., Maday, Y., and Turinici, G.: A “parareal” in time discretization of PDEs, C. R. Acad. Sci. 332, 661–668 (2001)
 (4) Falgout, R.D., Friedhoff, S., Kolev, Tz. V., MacLachlan, S. P., and Schroder, J. B.: Parallel time integration with multigrid. SIAM J. Sci. Comput. 36 (6), C635–C661 (2014)
 (5) Meeker, D.C.: Finite Element Method Magnetics, Version 4.2 (28Feb2018 Build), http://www.femm.info
 (6) Jackson, J.D.: Classical Electrodynamics, 3rd edn. Wiley and Sons, New York (1998)
 (7) Heise, B.: Analysis of a fully discrete finite element method for a nonlinear magnetic field problem. SIAM J. Num. Anal. 31(3), 745–759 (1994).

(8)
Schmidt, K., Sterz, O., Hiptmair, R.: Estimating the eddycurrent modeling error.
IEEE Trans. Magn. 44(6), 686–689 (2008).  (9) Emson, C.R.I., Trowbridge, C.W.: Transient 3d eddy currents using modified magnetic vector potentials and magnetic scalar potentials. IEEE Trans. Magn. 24(1), 86–89 (1988).
 (10) Nicolet, A., Delincé, F.: Implicit RungeKutta methods for transient magnetic field computation. IEEE Trans. Magn. 32(3), 1405–1408 (1996).
 (11) Ries, M. and Trottenberg, U.: MGREin blitzschneller elliptischer Löser. Preprint 277 SFB 72, Universität Bonn (1979)
 (12) Ries, M., Trottenberg, U., and Winter, G.: A note on MGR methods. J. Lin. Alg. Applic. 49, 1–26 (1983)
 (13) Brandt, A.: Multi–level adaptive solutions to boundary–value problems. Math. Comp. 31, 333–390 (1977)
 (14) Gander, M. J., KulchytskaRuchka, I., Niyonzima, I., Schöps, S.: A New Parareal Algorithm for Problems with Discontinuous Sources. Submitted to SISC, arXiv: 1803.05503 (2018)
Comments
There are no comments yet.