Electromagnetic field solver

Last updated

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.

Contents

Introduction

The extraction of parasitic circuit models is essential for various aspects of physical verification such as timing, signal integrity, substrate coupling, and power grid analysis. As circuit speeds and densities have increased, the need has grown to account accurately for parasitic effects for more extensive and more complicated interconnect structures. In addition, the electromagnetic complexity has grown as well, from resistance and capacitance to inductance, and now even full electromagnetic wave propagation. This increase in complexity has also grown for the analysis of passive devices such as integrated inductors. Electromagnetic behavior is governed by Maxwell's equations, and all parasitic extraction requires solving some form of Maxwell's equations. That form may be a simple analytic parallel plate capacitance equation or may involve a full numerical solution for a complex 3D geometry with wave propagation. In 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 Maxwell's equations is typically solved by one of two classes of methods. The first uses a differential form of the governing equations and requires the discretization (meshing) of the entire domain in which the electromagnetic fields reside. Two of the most common approaches in this first class are the finite difference (FD) and 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

Related Research Articles

<span class="mw-page-title-main">Capacitance</span> Ability of a body to store an electrical charge

Capacitance is the capability of a material object or device to store electric charge. It is measured by the change in charge in response to a difference in electric potential, expressed as the ratio of those quantities. Commonly recognized are two closely related notions of capacitance: self capacitance and mutual capacitance. An object that can be electrically charged exhibits self capacitance, for which the electric potential is measured between the object and ground. Mutual capacitance is measured between two components, and is particularly important in the operation of the capacitor, an elementary linear electronic component designed to add capacitance to an electric circuit.

The boundary element method (BEM) is a numerical computational method of solving linear partial differential equations which have been formulated as integral equations, including fluid mechanics, acoustics, electromagnetics, fracture mechanics, and contact mechanics.

<span class="mw-page-title-main">Finite-difference time-domain method</span>

Finite-difference time-domain (FDTD) or Yee's method is a numerical analysis technique used for modeling computational electrodynamics. Since it is a time-domain method, FDTD solutions can cover a wide frequency range with a single simulation run, and treat nonlinear material properties in a natural way.

The electric circuit extraction or simply circuit extraction, also netlist extraction, is the translation of an integrated circuit layout back into the electrical circuit (netlist) it is intended to represent. This extracted circuit is needed for various purposes including circuit simulation, static timing analysis, signal integrity, power analysis and optimization, and logic to layout comparison. Each of these functions require a slightly different representation of the circuit, resulting in the need for multiple layout extractions. In addition, there may be a postprocessing step of converting the device-level circuit into a purely digital circuit, but this is not considered part of the extraction process.

<span class="mw-page-title-main">Signal integrity</span>

Signal integrity or SI is a set of measures of the quality of an electrical signal. In digital electronics, a stream of binary values is represented by a voltage waveform. However, digital signals are fundamentally analog in nature, and all signals are subject to effects such as noise, distortion, and loss. Over short distances and at low bit rates, a simple conductor can transmit this with sufficient fidelity. At high bit rates and over longer distances or through various mediums, various effects can degrade the electrical signal to the point where errors occur and the system or device fails. Signal integrity engineering is the task of analyzing and mitigating these effects. It is an important activity at all levels of electronics packaging and assembly, from internal connections of an integrated circuit (IC), through the package, the printed circuit board (PCB), the backplane, and inter-system connections. While there are some common themes at these various levels, there are also practical considerations, in particular the interconnect flight time versus the bit period, that cause substantial differences in the approach to signal integrity for on-chip connections versus chip-to-chip connections.

<span class="mw-page-title-main">Computational electromagnetics</span> Branch of physics

Computational electromagnetics (CEM), computational electrodynamics or electromagnetic modeling is the process of modeling the interaction of electromagnetic fields with physical objects and the environment.

<span class="mw-page-title-main">Technology CAD</span>

Technology computer-aided design is a branch of electronic design automation that models semiconductor fabrication and semiconductor device operation. The modeling of the fabrication is termed Process TCAD, while the modeling of the device operation is termed Device TCAD. Included are the modelling of process steps, and modelling of the behavior of the electrical devices based on fundamental physics, such as the doping profiles of the devices. TCAD may also include the creation of compact models, which try to capture the electrical behavior of such devices but do not generally derive them from the underlying physics. SPICE simulator itself is usually considered as part of ECAD rather than TCAD.

<span class="mw-page-title-main">Power network design (IC)</span> Electrical system design based on power distribution

In the design of integrated circuits, power network design is the analysis and design of on-chip conductor networks that distribute electrical power on a chip. As in all engineering, this involves tradeoffs - the network must have adequate performance, be sufficiently reliable, but should not use more resources than required.

E The finite-difference frequency-domain (FDFD) method is a numerical solution method for problems usually in electromagnetism and sometimes in acoustics, based on finite-difference approximations of the derivative operators in the differential equation being solved.

The fast multipole method (FMM) is a numerical technique that was developed to speed up the calculation of long-ranged forces in the n-body problem. It does this by expanding the system Green's function using a multipole expansion, which allows one to group sources that lie close together and treat them as if they are a single source.

<span class="mw-page-title-main">Finite element method</span> Numerical method for solving physical or engineering problems

The finite element method (FEM) is a popular method for numerically solving differential equations arising in engineering and mathematical modeling. Typical problem areas of interest include the traditional fields of structural analysis, heat transfer, fluid flow, mass transport, and electromagnetic potential.

<span class="mw-page-title-main">Partial element equivalent circuit</span>

Partial element equivalent circuit method (PEEC) is partial inductance calculation used for interconnect problems from early 1970s which is used for numerical modeling of electromagnetic (EM) properties. The transition from a design tool to the full-wave method involves the capacitance representation, the inclusion of time retardation and the dielectric formulation. Using the PEEC method, the problem will be transferred from the electromagnetic domain to the circuit domain where conventional SPICE-like circuit solvers can be employed to analyze the equivalent circuit. By having the PEEC model one can easily include any electrical component e.g. passive components, sources, non-linear elements, ground, etc. to the model. Moreover, using the PEEC circuit, it is easy to exclude capacitive, inductive or resistive effects from the model when it is possible, in order to make the model smaller. As an example, in many applications within power electronics, the magnetic field is a dominating factor over the electric field due to the high current in the systems. Therefore, the model can be simplified by just neglecting capacitive couplings in the model which can simply be done by excluding the capacitors from the PEEC model.

<span class="mw-page-title-main">Jacob K. White</span> American electronics engineer and professor

Jacob K. White is the Cecil H. Green Professor of Electrical Engineering and Computer Science at the Massachusetts Institute of Technology. He researches fast numerical algorithms for simulation, particularly the simulation of circuits. His work on the FASTCAP program for three-dimensional capacitance calculation and FASTHENRY, a program for three-dimensional inductance calculations, is highly cited. He has also done extensive work on steady-state simulation of analog and microwave circuits. White was a significant early contributor to the development of Spectre and SpectreRF.

Eigenmode expansion (EME) is a computational electrodynamics modelling technique. It is also referred to as the mode matching technique or the bidirectional eigenmode propagation method. Eigenmode expansion is a linear frequency-domain method.

In electronic design automation, parasitic extraction is the calculation of the parasitic effects in both the designed devices and the required wiring interconnects of an electronic circuit: parasitic capacitances, parasitic resistances and parasitic inductances, commonly called parasitic devices, parasitic components, or simply parasitics.

The multilevel fast multipole method (MLFMM) is used along with method of moments (MoM) a numerical computational method of solving linear partial differential equations which have been formulated as integral equations of large objects almost faster without loss in accuracy. This method is an alternative formulation of the technology behind the MoM and is applicable to much larger structures like radar cross-section (RCS) analysis, antenna integration on large structures, reflector antenna design, finite size antenna arrays, etc., making full-wave current-based solutions of such structures a possibility.

<span class="mw-page-title-main">Frequency selective surface</span>

A frequency-selective surface (FSS) is any thin, repetitive surface designed to reflect, transmit or absorb electromagnetic fields based on the frequency of the field. In this sense, an FSS is a type of optical filter or metal-mesh optical filters in which the filtering is accomplished by virtue of the regular, periodic pattern on the surface of the FSS. Though not explicitly mentioned in the name, FSS's also have properties which vary with incidence angle and polarization as well - these are unavoidable consequences of the way in which FSS's are constructed. Frequency-selective surfaces have been most commonly used in the radio frequency region of the electromagnetic spectrum and find use in applications as diverse as the aforementioned microwave oven, antenna radomes and modern metamaterials. Sometimes frequency selective surfaces are referred to simply as periodic surfaces and are a 2-dimensional analog of the new periodic volumes known as photonic crystals.

<span class="mw-page-title-main">Weng Cho Chew</span> Malaysian-American electrical engineer

Weng Cho Chew is a Malaysian-American electrical engineer and applied physicist known for contributions to wave physics, especially computational electromagnetics. He is a Distinguished Professor of Electrical and Computer Engineering at Purdue University.

Kane Shee-Gong Yee is a Chinese-American electrical engineer and mathematician. He is best known for introducing the finite-difference time-domain method (FDTD) in 1966.

<span class="mw-page-title-main">Method of moments (electromagnetics)</span> Numerical method in computational electromagnetics

The method of moments (MoM), also known as the moment method and method of weighted residuals, is a numerical method in computational electromagnetics. It is used in computer programs that simulate the interaction of electromagnetic fields such as radio waves with matter, for example antenna simulation programs like NEC that calculate the radiation pattern of an antenna. Generally being a frequency-domain method, it involves the projection of an integral equation into a system of linear equations by the application of appropriate boundary conditions. This is done by using discrete meshes as in finite difference and finite element methods, often for the surface. The solutions are represented with the linear combination of pre-defined basis functions; generally, the coefficients of these basis functions are the sought unknowns. Green's functions and Galerkin method play a central role in the method of moments.

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. ↑ Yu, Wenjian; Zhuang, Hao; Zhang, Chao; Hu, Gang; Liu, Zhi (2013). "RWCap: A Floating Random Walk Solver for 3-D Capacitance Extraction of Very-Large-Scale Integration Interconnects". IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems. 32 (3): 353–366. CiteSeerX   10.1.1.719.3986 . doi:10.1109/TCAD.2012.2224346. S2CID   16351864.
  3. ↑ O. M. Ramahi; B. Archambeault (1995). "Adaptive absorbing boundary conditions in finite-difference time-domain applications for EMC simulations". IEEE Trans. Electromagn. Compat. 37 (4): 580–583. doi:10.1109/15.477343.
  4. ↑ J.C. Veihl; R. Mittra (Feb 1996). "An efficient implementation of Berenger's perfectly matched layer (PML) for finite-difference time-domain mesh truncation". IEEE Microwave and Guided Wave Letters. 6 (2): 94. doi:10.1109/75.482000.
  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. ↑ K. Nabors; J. White (November 1991). "Fastcap: A multipole accelerated 3-D capacitance extraction program". IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems. 10 (11): 1447–1459. CiteSeerX   10.1.1.19.9745 . doi:10.1109/43.97624.
  8. ↑ A. Brandt. Multilevel computations of integral transforms and particle interactions with oscillatory kernels. Computer Physics Communications, 65:24-38, 1991.
  9. ↑ J.R. Phillips; J.K. White (October 1997). "A precorrected-FFT method for electrostatic analysis of complicated 3-d structures". IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems. 16 (10): 1059–1072. CiteSeerX   10.1.1.20.791 . doi:10.1109/43.662670.
  10. ↑ S. Kapur; D.E. Long (Oct–Dec 1998). "IES3: Efficient electrostatic and electromagnetic simulation". IEEE Computational Science and Engineering. 5 (4): 60–67. doi:10.1109/99.735896.
  11. ↑ J.M. Song; C.C. Lu; W.C. Chew; S.W. Lee (June 1998). "Fast Illinois Solver Code (FISC)". IEEE Antennas and Propagation Magazine. 40 (3): 27–34. Bibcode:1998IAPM...40...27S. CiteSeerX   10.1.1.7.8263 . doi:10.1109/74.706067.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.