The main research activity in this area has been and are performed, by Albert Ruehli1 at IBM Thomas J. Watson Research Center, starting with a publication in 1972. At that time the foundation of the PEEC method was presented, i.e. the calculation of the partial inductances. The PEEC method was extended to more generalized problems, including dielectric material and retardation effect.
The PEEC method is not one of the most common techniques used in EM simulation software or as a research area but it has just been starting to gain recognition and for the first time there is a session at the 2001 IEEE EMC Symposium named after the technique. In the mid-1990s, two researchers from the University of L'Aquila in Italy, Professor Antonio Orlandi and Professor Giulio Antonini, published their first PEEC paper and are now together with Dr. Ruehli considered the top researchers in the area. Starting year 2006, several research projects have been initiated by the faculty of Computer Science and Electrical Engineering of Luleå University of Technology in Sweden in the focus area of PEEC with the emphasis on computer-based solvers for PEEC.
PEEC is widely used for combined electromagnetic and circuit problems in various areas such as power electronics, antenna design, signal integrity analysis, etc. Using PEEC the designed model of a physical structure is transferred from the electromagnetic domain into the circuit domain. Therefore, external electrical components and circuits can be connected to the equivalent circuit which consists of extracted partial elements, in a straightforward manner. Moreover, since the final model consists of circuit elements, various components can easily be excluded from the circuit to simplify the problem while the accuracy is still ensured. For instance, for low-frequency problems, one can safely remove capacitive couplings without degrading the accuracy of the results and hence reduce the problem size and complexity.
The classical PEEC method is derived from the equation for the total electric field at a point2 written as
where E → i {\displaystyle {\vec {E}}^{i}} is an incident electric field, J → {\displaystyle {\vec {J}}} is a current density, A → {\displaystyle {\vec {A}}} is the magnetic vector potential, ϕ {\displaystyle \phi } is the scalar electric potential, and σ {\displaystyle \sigma } the electrical conductivity all at observation point r → {\displaystyle {\vec {r}}} . In the figures on the right, an orthogonal metal strip with 3 nodes and 2 cells, and the corresponding PEEC circuit are shown.
By using the definitions of the scalar and vector potentials, the current- and charge-densities are discretized by defining pulse basis functions for the conductors and dielectric materials. Pulse functions are also used for the weighting functions resulting in a Galerkin type solution. By defining a suitable inner product, a weighted volume integral over the cells, the field equation can be interpreted as Kirchhoff's voltage law over a PEEC cell consisting of partial self inductances between the nodes and partial mutual inductances representing the magnetic field coupling in the equivalent circuit. The partial inductances are defined as
for volume cell α {\displaystyle \alpha } and β {\displaystyle \beta } . Then, the coefficients of potentials are computed as
and a resistive term between the nodes, defined as
The rigorous full-wave version of the PEEC method is called (Lp,P,R,t) PEEC, where Lp is partial inductance, P is the Maxwell potential coefficient (inverse of capacitance), R is resistance, and t is the time-delay. If available, a reduced model of the full-wave version can be used. For example, if the EIP structure is electrically small, the delay term t can be omitted and the model can be reduced to (Lp,P,R) PEEC model. In addition, if the angular frequency w is sufficiently high so that w*Lp >> R, we can omit R term and use approximate (Lp,P) PEEC model. According to various modeling situations, (Lp) and (Lp,R) models are also useful.
Model Order Reduction (MOR) has become an active research topic for circuit models in general and PEEC models in particular. Integration of a PEEC model directly into a circuit simulator is computationally expensive for two main facts. One is that a large number of circuit elements are generated for complex structures at high frequencies, and the other is that the circuit matrices based on modified nodal analysis (MNA) are usually dense due to full inductive and capacitive coupling. In order to model/simulate such problems efficiently, developing compact model representation via model order reduction is desirable for PEEC modeling.
A. E. Ruehli: Equivalent Circuit Models for Three-Dimensional Multiconductor Systems, IEEE Transactions on Microwave Theory and Techniques, Vol. 22 (1974), Nr. 3 ↩
S. Ramo, J. R. Whinnery and T. Van Duzer: Fields and Waves in Communication Electronics, John Wiley and Sons, 1972 ↩