Multiscale Modeling and Multirate Time-Integration of Field/Circuit Coupled Problems

This treatise is intended for mathematicians and computational engineers that work on modeling, coupling and simulation of electromagnetic problems. This includes lumped electric networks, magnetoquasistatic field and semiconductor devices. Their coupling yields a multiscale system of partial differential algebraic equations containing device models of any dimension interconnected by the electric network. It is solved in time domain by multirate techniques that efficiently exploit the structure. The central idea is the usage of lumped surrogate models that describe latent model parts sufficiently accurate (e.g. the field model by an inductance) even if other model parts (e.g. the circuit) exhibit highly dynamic behavior. We propose dynamic iteration and a bypassing technique using surrogate Schur complements. A mathematical convergence analysis is given and numerical examples are discussed. They show a clear reduction in the computational costs compared to single rate approaches.

[1]  Rolf Schuhmann,et al.  The Perfect Boundary Approximation Technique Facing the Big Challenge of High Precision Field Computation , 1998 .

[2]  F. Henrotte,et al.  A Circuit Coupling Method Based on a Temporary Linearization of the Energy Balance of the Finite Element Model , 2008, IEEE Transactions on Magnetics.

[3]  Siegfried Selberherr,et al.  Mixed-mode device simulation , 2000 .

[4]  J. Gyselinck,et al.  Modeling of 3-D stranded inductors with the magnetic vector potential formulation and spatially dependent turn voltages of reduced support , 2004, IEEE Transactions on Magnetics.

[5]  Joachim Schöberl,et al.  An algebraic multigrid method for finite element discretizations with edge elements , 2002, Numer. Linear Algebra Appl..

[6]  Ricardo Riaza,et al.  Differential-Algebraic Systems: Analytical Aspects and Circuit Applications , 2008 .

[7]  Andreas Bartel,et al.  A Cosimulation Framework for Multirate Time Integration of Field/Circuit Coupled Problems , 2010, IEEE Transactions on Magnetics.

[8]  Herbert De Gersem,et al.  A finite element model for foil winding simulation , 2000 .

[9]  Janne Roos,et al.  Scientific computing in electrical engineering : SCEE 2008 , 2010 .

[10]  R. Stratton,et al.  Diffusion of Hot and Cold Electrons in Semiconductor Barriers , 1962 .

[11]  Caren Tischendorf,et al.  Structural analysis of electric circuits and consequences for MNA , 2000, Int. J. Circuit Theory Appl..

[12]  Michael Günther,et al.  Preconditioned Dynamic Iteration for Coupled Differential-Algebraic Systems , 2001 .

[13]  A. Bossavit Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism , 1988 .

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

[15]  S. J. Salon,et al.  Finite element analysis of electrical machines , 1995 .

[16]  C. W. Gear,et al.  Multirate linear multistep methods , 1984 .

[17]  Sebastian Schöps,et al.  Structural analysis of electrical circuits including magnetoquasistatic devices , 2011 .

[18]  Wolfgang Hackbusch,et al.  Direct Domain Decomposition using the Hierarchical Matrix Technique , 2002 .

[19]  Diana Estévez Schwarz,et al.  Consistent initialization for index-2 differential algebraic equations and its application to circuit simulation , 2000 .

[20]  Daniel Sundström,et al.  Mathematics in industry , 1985 .

[21]  Caren Tischendorf,et al.  Consistent Initialization for Coupled Circuit-Device Simulation , 2010 .

[22]  Mitiko Miura-Mattausch,et al.  On Local Handling of Inner Equations in Compact Models , 2010 .

[23]  H. De Gersem,et al.  Field-circuit coupling for time-harmonic models discretized by the finite integration technique , 2004, IEEE Transactions on Magnetics.

[24]  Thomas Weiland,et al.  Solution of Maxwell's equations , 1992 .

[25]  Gérard Meunier,et al.  3D line current model of coils and external circuits , 1995 .

[26]  Enzo Tonti,et al.  On the Geometrical Structure of Electromagnetism , 1999 .

[27]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[28]  K. Preis,et al.  On the use of the magnetic vector potential in the finite-element analysis of three-dimensional eddy currents , 1989 .

[29]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[30]  Paola Pietra,et al.  New mixed finite element schemes for current continuity equations , 1990 .

[31]  Ralf Hiptmair,et al.  Current and voltage excitations for the eddy current model , 2005 .

[32]  A. Bossavit Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements , 1997 .

[33]  Prashanth K. Vijalapura,et al.  Fractional step methods for index-1 differential-algebraic equations , 2005 .

[34]  Ernst Hairer,et al.  The numerical solution of differential-algebraic systems by Runge-Kutta methods , 1989 .

[35]  Stefan Kurz,et al.  Field-Circuit Coupling by Means of the Woodbury Formula , 2004 .

[36]  Carlo de Falco,et al.  A Demonstrator Platform for Coupled Multiscale Simulation , 2007 .

[37]  A. Konrad,et al.  Coupled Field-Circuit Problems: Trends and Accomplishments , 1992, Digest of the Fifth Biennial IEEE Conference on Electromagnetic Field Computation.

[38]  Koji Fujiwara,et al.  Improvements of convergence characteristics of Newton-Raphson method for nonlinear magnetic field analysis , 1992 .

[39]  YANQING CHEN,et al.  Algorithm 8 xx : CHOLMOD , supernodal sparse Cholesky factorization and update / downdate ∗ , 2006 .

[40]  E. F. Kaasschieter,et al.  Preconditioned conjugate gradients for solving singular systems , 1988 .

[41]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

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

[43]  Andreas Bartel,et al.  A Concept for Classification of Partial Differential Algebraic Equations in Nanoelectronics , 2008 .

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

[45]  Andreas Bartel,et al.  Parabolic Differential-Algebraic Models in Electrical Network Design , 2005, Multiscale Model. Simul..

[46]  A. Bossavit "Stiff" problems in eddy-current theory and the regularization of Maxwell's equations , 2001 .

[47]  J. R. Brauer,et al.  Anisotropic materials in electromagnetic finite element analysis , 1991 .

[48]  Martin J. Gander,et al.  Optimized waveform relaxation methods for RC type circuits , 2004, IEEE Transactions on Circuits and Systems I: Regular Papers.

[49]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[50]  H. G. Brachtendorf,et al.  Numerical steady state analysis of electronic circuits driven by multi-tone signals , 1996 .

[51]  J. Verwer,et al.  A multirate time stepping strategy for parabolic PDE , 2005 .

[52]  Igor Tsukerman,et al.  Finite Element Differential–Algebraic Systems for Eddy Current Problems , 2002, Numerical Algorithms.

[53]  Markus Brunk,et al.  Positivity preserving discretization of time dependent semiconductor drift-diffusion equations , 2012 .

[54]  T. Steinmetz,et al.  Multiple right-hand side techniques for the numerical simulation of quasistatic electric and magnetic fields , 2008 .

[55]  Thomas Weiland,et al.  A CONSISTENT SUBGRIDDING SCHEME FOR THE FINITE DIFFERENCE TIME DOMAIN METHOD , 1996 .

[56]  Sebastian Schöps,et al.  Multirate Time Integration of Field/Circuit Coupled Problems by Schur Complements , 2012, CSC 2012.

[57]  Hermann A. Haus,et al.  Electromagnetic Fields And Energy , 1989 .

[58]  M. Günther,et al.  The DAE-index in electric circuit simulation , 1995 .

[59]  Andreas Bartel,et al.  Multirate Methods in Electrical Circuit Simulation , 2002 .

[60]  C. W. Gear,et al.  The index of general nonlinear DAEs , 1995 .

[61]  Sebastian Schöps,et al.  A Convergent Iteration Scheme for Semiconductor/Circuit Coupled Problems , 2012 .

[62]  U. Langer,et al.  Boundary and Finite Element Domain Decomposition Methods , 2008 .

[63]  Frank Cameron,et al.  Variable step size time integration methods for transient eddy current problems , 1998 .

[64]  J. Vaananen,et al.  Circuit theoretical approach to couple two-dimensional finite element models with external circuit equations , 1996 .

[65]  Alberto L. Sangiovanni-Vincentelli,et al.  The Waveform Relaxation Method for Time-Domain Analysis of Large Scale Integrated Circuits , 1982, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[66]  Wil H. A. Schilders,et al.  Numerical methods in electromagnetics , 2005 .

[67]  Alberto Valli,et al.  A domain decomposition approach for heterogeneous time-harmonic Maxwell equations , 1997 .

[68]  Thomas Weiland,et al.  Field—circuit coupled models in electromagnetic simulation , 2002 .

[69]  Oszkar Biro,et al.  Various FEM formulations for the calculation of transient 3D eddy currents in nonlinear media , 1995 .

[70]  Falk Ebert,et al.  On Partitioned Simulation of Electrical Circuits using Dynamic Iteration Methods , 2008 .

[71]  Thomas Weiland,et al.  Coupled finite-element, spectral-element discretisation for models with circular inclusions and far-field domains , 2002 .

[72]  J. Maxwell A Dynamical Theory of the Electromagnetic Field , .

[73]  S. Kanerva,et al.  Data transfer methodology between a FEM program and a system simulator , 2001, ICEMS'2001. Proceedings of the Fifth International Conference on Electrical Machines and Systems (IEEE Cat. No.01EX501).

[74]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[75]  Jian-She Wang A nodal analysis approach for 2D and 3D magnetic-circuit coupled problems , 1996 .

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

[77]  Markus Clemens,et al.  Large systems on equations in the discrete electromagnetism : formulations and numerical algorithms , 2005 .

[78]  W. Fu,et al.  A general cosimulation approach for coupled field-circuit problems , 2006, IEEE Transactions on Magnetics.

[79]  Thomas Weiland,et al.  On the Unique Numerical Solution of Maxwellian Eigenvalue Problems in Three-dimensions , 1984 .

[80]  O. Nevanlinna,et al.  Convergence of dynamic iteration methods for initial value problems , 1987 .

[81]  Barry Smith,et al.  Domain Decomposition Methods for Partial Differential Equations , 1997 .

[82]  Barbara I. Wohlmuth,et al.  An iterative substructuring method for Maxwell's equations in two dimensions , 2001, Math. Comput..

[83]  Kay Hameyer,et al.  Methodologies for coupled transient electromagnetic-thermal finite-element modeling of electrical energy transducers , 2002 .

[84]  Inmaculada Higueras,et al.  Differential algebraic equations with properly stated leading terms , 2004 .

[85]  G. Bedrosian,et al.  A new method for coupling finite element field solutions with external circuits and kinematics , 1993 .

[86]  H. V. D. Vorst,et al.  The rate of convergence of Conjugate Gradients , 1986 .

[87]  Michael Guenther,et al.  CAD based electric circuit modeling in industry. Part I: Mathematical structure and index of network equations. Part II: Impact of circuit configurations and parameters , 1999 .

[88]  Thomas Weiland,et al.  Differential Material Matrices for the Finite Integration Technique , 2006 .

[89]  E. Griepentrog,et al.  Differential-algebraic equations and their numerical treatment , 1986 .

[90]  T. Weiland Time Domain Electromagnetic Field Computation with Finite Difference Methods , 1996 .

[91]  H. De Gersem,et al.  Construction of Differential Material Matrices for the Orthogonal Finite-Integration Technique With Nonlinear Materials , 2007, IEEE Transactions on Magnetics.

[92]  Michael Günther,et al.  A charge oriented mixed multirate method for a special class of index-1 network equations in chip design , 2005 .

[93]  A. Konrad Integrodifferential finite element formulation of two-dimensional steady-state skin effect problems , 1982 .

[94]  R. März Differential Algebraic Systems with Properly Stated Leading Term and MNA Equations , 2003 .

[95]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[96]  M. S. Soto,et al.  Numerical analysis of DAEs from coupled circuit and semiconductor simulation , 2005 .

[97]  Caren Tischendorf,et al.  Coupled Systems of Differential Algebraic and Partial Differential Equations in Circuit and Device Simulation , 2003 .

[98]  S. Selberherr Analysis and simulation of semiconductor devices , 1984 .

[99]  Sebastian Schöps,et al.  Coupling and Simulation of Lumped Electric Circuits Refined by 3-D Magnetoquasistatic Conductor Models Using MNA and FIT , 2008 .

[100]  A. Quarteroni Numerical Models for Differential Problems , 2009 .

[101]  Albert E. Ruehli,et al.  WAVEFORM RELAXATION: THEORY AND PRACTICE , 1985 .

[102]  Sebastian Schöps,et al.  DAE-Index and Convergence Analysis of Lumped Electric Circuits Refined by 3-D Magnetoquasistatic Conductor Models , 2010 .

[103]  A. Y. Hannalla,et al.  Numerical analysis of transient field problems in electrical machines , 1976 .

[104]  Thomas Weiland,et al.  Regularization of eddy-current formulations using discrete grad-div operators , 2002 .

[105]  François Henrotte,et al.  A general and natural method to define circuit relations associated with magnetic vector potential formulations , 1999 .

[106]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[107]  Alain Bossavit,et al.  Yee-like schemes on staggered cellular grids: a synthesis between FIT and FEM approaches , 2000 .

[108]  G. Denk Circuit Simulation for Nanoelectronics , 2006 .

[109]  K. Hameyer,et al.  An Efficient Field-Circuit Coupling Based on a Temporary Linearization of FE Electrical Machine Models , 2009, IEEE Transactions on Magnetics.

[110]  U. van Rienen,et al.  Computation of currents induced by ELF electric fields in anisotropic human tissues using the Finite Integration Technique (FIT) , 2005 .

[111]  T. Steinmetz,et al.  Reuse, Recycle, Reduce (3R) -- strategies for the calculation of transient magnetic fields , 2009 .

[112]  Kevin Burrage,et al.  Parallel and sequential methods for ordinary differential equations , 1995, Numerical analysis and scientific computation.

[113]  Stefan Vandewalle,et al.  SOR WAVEFORM , 1997 .

[114]  Bodo Heise Analysis of a fully discrete finite element method for a nonlinear magnetic field problem , 1994 .

[115]  Lauri Kettunen,et al.  Gauging in Whitney spaces , 1999 .

[116]  Kay Hameyer,et al.  An algebraic multigrid method for solving very large electromagnetic systems , 1997 .

[117]  Manfred Kaltenbacher,et al.  Numerical Simulation of Mechatronic Sensors and Actuators , 2004 .

[118]  Leon O. Chua,et al.  Computer-Aided Analysis Of Electronic Circuits , 1975 .

[119]  T. Weiland A discretization model for the solution of Maxwell's equations for six-component fields , 1977 .

[120]  K. Preis,et al.  Finite element analysis of anisotropic and nonlinear magnetic circuits , 1981 .

[121]  Caren Tischendorf Model Design Criteria for Integrated Circuits to Have a Unique Solution and Good Numerical Properties , 2001 .

[122]  M. Arnold,et al.  Errors in the numerical solution ofnonlinear di erential-algebraic systems of index 2 , 1995 .

[123]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

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

[125]  Gene H. Golub,et al.  Matrix computations , 1983 .

[126]  Galina Benderskaya,et al.  Numerical methods for transient field-circuit coupled simulations based on the finite integration technique and a mixed circuit formulation , 2007 .

[127]  Zdzislaw Jackiewicz,et al.  Convergence of Waveform Relaxation Methods for Differential-Algebraic Systems , 1996 .

[128]  Muruhan Rathinam,et al.  Dynamic Iteration Using Reduced Order Models: A Method for Simulation of Large Scale Modular Systems , 2002, SIAM J. Numer. Anal..

[129]  Albert E. Ruehli,et al.  The modified nodal approach to network analysis , 1975 .

[130]  Roswitha März,et al.  Solvability of linear differential algebraic equations with properly stated leading terms , 2004 .

[131]  Frédéric Guyomarc'h,et al.  An Augmented Conjugate Gradient Method for Solving Consecutive Symmetric Positive Definite Linear Systems , 2000, SIAM J. Matrix Anal. Appl..

[132]  Qiushi Chen,et al.  A review of finite element open boundary techniques for static and quasi-static electromagnetic field problems , 1997 .

[133]  S. Vandewalle Parallel multigrid waveform relaxation for parabolic problems , 1993 .