Electromagnetic field solver

Source: Wikipedia, the free encyclopedia.

Electromagnetic field solvers (or sometimes just field solvers) are specialized programs that solve (a subset of) Maxwell's equations directly. They form a part of the field of electronic design automation, or EDA, and are commonly used in the design of integrated circuits and printed circuit boards. They are used when a solution from first principles or the highest accuracy is required.

Introduction

The

layout extraction, analytic formulas for simple or simplified geometry can be used where accuracy is less important than speed. Still, when the geometric configuration is not simple, and accuracy demands do not allow simplification, a numerical solution of the appropriate form of Maxwell's equations
must be employed.

The appropriate form of

finite element (FEM) methods. The resultant linear algebraic system (matrix) that must be solved is large but sparse (contains very few non-zero entries). Sparse linear solution methods, such as sparse factorization, conjugate-gradient, or multigrid methods
can be used to solve these systems, the best of which require CPU time and memory of O(N) time, where N is the number of elements in the discretization. However, most problems in electronic design automation (EDA) are open problems, also called exterior problems, and since the fields decrease slowly towards infinity, these methods can require extremely large N.

The second class of methods are integral equation methods which instead require a discretization of only electromagnetic field sources. Those sources can be physical quantities, such as the surface charge density for the capacitance problem, or mathematical abstractions resulting from applying Green's theorem. When the sources exist only on two-dimensional surfaces for three-dimensional problems, the method is often called method of moments (MoM) or boundary element method (BEM). For open problems, the sources of the field exist in a much smaller domain than the fields themselves, and thus the size of linear systems generated by integral equations methods are much smaller than FD or FEM. Integral equation methods, however, generate dense (all entries are nonzero) linear systems, making such methods preferable to FD or FEM only for small problems. Such systems require O(n2) memory to store and O(n3) to solve via direct Gaussian elimination or, at best, O(n2) if solved iteratively. Increasing circuit speeds and densities require the solution of increasingly complicated interconnect, making dense integral equation approaches unsuitable due to these high growth rates of computational cost with increasing problem size.

In the past two decades, much work has gone into improving both the differential and integral equation approaches, as well as new approaches based on random walk methods.[1][2] Methods of truncating the discretization required by the FD and FEM approaches has greatly reduced the number of elements required.[3][4] Integral equation approaches have become particularly popular for interconnect extraction due to sparsification techniques, also sometimes called matrix compression, acceleration, or matrix-free techniques, which have brought nearly O(n) growth in storage and solution time to integral equation methods.[5][6][7][8][9][10][11]

Sparsified integral equation techniques are typically used in the IC industry to solve capacitance and inductance extraction problems. The random-walk methods have become quite mature for capacitance extraction. For problems requiring the solution of the full Maxwell's equations (full-wave), both differential and integral equation approaches are common.

See also

References

  1. ^ Y. L. Le Coz and R. B. Iverson. A stochastic algorithm for high-speed capacitance extraction in integrated circuits. Solid State Electronics, 35(7):1005-1012, 1992.
  2. S2CID 16351864
    .
  3. .
  4. .
  5. ^ L. Greengard. The Rapid Evaluation of Potential Fields in Particle Systems. M.I.T. Press, Cambridge, Massachusetts, 1988.
  6. ^ V. Rokhlin. Rapid solution of integral equations of classical potential theory. Journal of Computational Physics, 60(2):187-207, September 15, 1985.
  7. .
  8. ^ A. Brandt. Multilevel computations of integral transforms and particle interactions with oscillatory kernels. Computer Physics Communications, 65:24-38, 1991.
  9. .
  10. .
  11. .