This article may be too technical for most readers to understand.(September 2021) |
The method of moments (MoM), also known as the moment method and method of weighted residuals, [1] 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, [lower-alpha 1] 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.
For many applications, the method of moments is identical to the boundary element method. [lower-alpha 2] It is one of the most common methods in microwave and antenna engineering.
Development of boundary element method and other similar methods for different engineering applications is associated with the advent of digital computing in the 1960s. [6] Prior to this, variational methods were applied to engineering problems at microwave frequencies by the time of World War II. [7] While Julian Schwinger and Nathan Marcuvitz have respectively compiled these works into lecture notes and textbooks, [8] [9] Victor Rumsey has formulated these methods into the "reaction concept" in 1954. [10] The concept was later shown to be equivalent to the Galerkin method. [7] In the late 1950s, an early version of the method of moments was introduced by Yuen Lo at a course on mathematical methods in electromagnetic theory at University of Illinois. [11]
In the 1960s, early research work on the method was published by Kenneth Mei, Jean van Bladel [12] and Jack Richmond. [13] In the same decade, the systematic theory for the method of moments in electromagnetics was largely formalized by Roger Harrington. [14] While the term "the method of moments" was coined earlier by Leonid Kantorovich and Gleb Akilov for analogous numerical applications, [15] Harrington has adapted the term for the electromagnetic formulation. [7] Harrington published the seminal textbook Field Computation by Moment Methods on the moment method in 1968. [14] The development of the method and its indications in radar and antenna engineering attracted interest; MoM research was subsequently supported United States government. The method was further popularized by the introduction of generalized antenna modeling codes such as Numerical Electromagnetics Code, which was released into public domain by the United States government in the late 1980s. [16] [17] In the 1990s, introduction of fast multipole and multilevel fast multipole methods enabled efficient MoM solutions to problems with millions of unknowns. [18] [19] [20]
Being one of the most common simulation techniques in RF and microwave engineering, the method of moments forms the basis of many commercial design software such as FEKO. [21] Many non-commercial and public domain codes of different sophistications are also available. [22] In addition to its use in electrical engineering, the method of moments has been applied to light scattering [23] and plasmonic problems. [24] [25] [26]
An inhomogeneous integral equation can be expressed as: where L denotes a linear operator, g denotes the known forcing function and f denotes the unknown function. f can be approximated by a finite number of basis functions ():
By linearity, substitution of this expression into the equation yields:
We can also define a residual for this expression, which denotes the difference between the actual and the approximate solution:
The aim of the method of moments is to minimize this residual, which can be done by using appropriate weighting or testing functions, hence the name method of weighted residuals. [27] After the determination of a suitable inner product for the problem, the expression then becomes:
Thus, the expression can be represented in the matrix form:
The resulting matrix is often referred as the impedance matrix. [28] The coefficients of the basis functions can be obtained through inverting the matrix. [29] For large matrices with a large number of unknowns, iterative methods such as conjugate gradient method can be used for acceleration. [30] The actual field distributions can be obtained from the coefficients and the associated integrals. [31] The interactions between each basis function in MoM is ensured by Green's function of the system. [32]
Different basis functions can be chosen to model the expected behavior of the unknown function in the domain; these functions can either be subsectional or global. [33] Choice of Dirac delta function as basis function is known as point-matching or collocation. This corresponds to enforcing the boundary conditions on discrete points and is often used to obtain approximate solutions when the inner product operation is cumbersome to perform. [34] [35] Other subsectional basis functions include pulse, piecewise triangular, piecewise sinusoidal and rooftop functions. [33] Triangular patches, introduced by S. Rao, D. Wilton and A. Glisson in 1982, [36] are known as RWG basis functions and are widely used in MoM. [37] Characteristic basis functions were also introduced to accelerate computation and reduce the matrix equation. [38] [39]
The testing and basis functions are often chosen to be the same; this is known as the Galerkin method. [29] Depending on the application and studied structure, the testing and basis functions should be chosen appropriately to ensure convergence and accuracy, as well as to prevent possible high order algebraic singularities. [40]
Depending on the application and sought variables, different integral or integro-differential equations are used in MoM. Radiation and scattering by thin wire structures, such as many types of antennas, can be modeled by specialized equations. [41] For surface problems, common integral equation formulations include electric field integral equation (EFIE), magnetic field integral equation (MFIE) [42] and mixed-potential integral equation (MPIE). [43]
As many antenna structures can be approximated as wires, thin wire equations are of interest in MoM applications. Two commonly used thin-wire equations are Pocklington and Hallén integro-differential equations. [44] Pocklington's equation precedes the computational techniques, having been introduced in 1897 by Henry Cabourn Pocklington. [45] For a linear wire that is centered on the origin and aligned with the z-axis, the equation can be written as: where and denote the total length and thickness, respectively. is the Green's function for free space. The equation can be generalized to different excitation schemes, including magnetic frills. [46]
Hallén integral equation, published by E. Hallén in 1938, [47] can be given as:
This equation, despite being more well-behaved than the Pocklington's equation, [48] is generally restricted to the delta-gap voltage excitations at the antenna feed point, which can be represented as an impressed electric field. [46]
The general form of electric field integral equation (EFIE) can be written as: where is the incident or impressed electric field. is the Green function for Helmholtz equation and represents the wave impedance. The boundary conditions are met at a defined PEC surface. EFIE is a Fredholm integral equation of the first kind. [42]
Another commonly used integral equation in MoM is the magnetic field integral equation (MFIE), which can be written as:
MFIE is often formulated to be a Fredholm integral equation of the second kind and is generally well-posed. Nevertheless, the formulation necessitates the use of closed surfaces, which limits its applications. [42]
Many different surface and volume integral formulations for MoM exist. In many cases, EFIEs are converted to mixed potential integral equations (MFIE) through the use of Lorenz gauge condition; this aims to reduce the orders of singularities through the use of magnetic vector and scalar electric potentials. [49] [50] In order to bypass the internal resonance problem in dielectric scattering calculations, combined-field integral equation (CFIE) and Poggio—Miller—Chang—Harrington—Wu—Tsai (PMCHWT) formulations are also used. [51] Another approach, the volumetric integral equation, necessitates the discretization of the volume elements and is often computationally expensive. [52]
MoM can also be integrated with physical optics theory [53] and finite element method. [54]
Appropriate Green's function for the studied structure must be known to formulate MoM matrices: automatic incorporation of the radiation condition into the Green's function makes MoM particularly useful for radiation and scattering problems. Even though the Green function can be derived in closed form for very simple cases, more complex structures necessitate numerical derivation of these functions. [55]
Full wave analysis of planarly-stratified structures in particular, such as microstrips or patch antennas, necessitate the derivation of Green's functions that are peculiar to these geometries. [50] [56] This can be achieved in two different methods. In the first method, known as spectral-domain approach, the inner products and convolution operation for MoM matrix entries are evaluated in the Fourier space with analytically-derived spectral-domain Green's functions through Parseval's theorem. [57] [58] [59] The other approach is based on the use of spatial-domain Green's functions. This involves the inverse Hankel transform of the spectral-domain Green's function, which is defined on the Sommerfeld integration path. Nevertheless, this integral cannot be evaluated analytically, and its numerical evaluation is often computationally expensive due to the oscillatory kernels and slowly-converging nature of the integral. [60] Common approaches for evaluating these integrals include tail extrapolation approaches such as weighted-averages method. [61]
Other approaches include the approximation of the integral kernel. Following the extraction of quasi-static and surface pole components, these integrals can be approximated as closed-form complex exponentials through Prony's method or generalized pencil-of-function method; thus, the spatial Green's functions can be derived through the use of appropriate identities such as Sommerfeld identity. [62] [63] [64] This method is known in the computational electromagnetics literature as the discrete complex image method (DCIM), since the Green's function is effectively approximated with a discrete number of image dipoles that are located within a complex distance from the origin. [65] The associated Green's functions are referred as closed-form Green's functions. [63] [64] The method has also been extended for cylindrically-layered structures. [66]
Rational-function fitting method, [67] [68] as well as its combinations with DCIM, [64] can also be used to approximate closed-form Green's functions. Alternatively, the closed-form Green's function can be evaluated through method of steepest descent. [69] For the periodic structures such as phased arrays and frequency selective surfaces, series acceleration methods such as Kummer's transformation and Ewald summation is often used to accelerate the computation of the periodic Green's function. [70] [71]
A metamaterial is a type of material engineered to have a property, typically rarely observed in naturally occurring materials, that is derived not from the properties of the base materials but from their newly designed structures. Metamaterials are usually fashioned from multiple materials, such as metals and plastics, and are usually arranged in repeating patterns, at scales that are smaller than the wavelengths of the phenomena they influence. Their precise shape, geometry, size, orientation, and arrangement give them their "smart" properties of manipulating electromagnetic, acoustic, or even seismic waves: by blocking, absorbing, enhancing, or bending waves, to achieve benefits that go beyond what is possible with conventional materials.
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.
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-field integral equation is a relationship that allows the calculation of an electric field generated by an electric current distribution.
In classical electromagnetism, reciprocity refers to a variety of related theorems involving the interchange of time-harmonic electric current densities (sources) and the resulting electromagnetic fields in Maxwell's equations for time-invariant linear media under certain constraints. Reciprocity is closely related to the concept of symmetric operators from linear algebra, applied to electromagnetism.
Computational electromagnetics (CEM), computational electrodynamics or electromagnetic modeling is the process of modeling the interaction of electromagnetic fields with physical objects and the environment using computers.
Electromagnetic field solvers are specialized programs that solve 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.
Characteristic modes (CM) form a set of functions which, under specific boundary conditions, diagonalizes operator relating field and induced sources. Under certain conditions, the set of the CM is unique and complete (at least theoretically) and thereby capable of describing the behavior of a studied object in full.
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.
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.
Generalized pencil-of-function method (GPOF), also known as matrix pencil method, is a signal processing technique for estimating a signal or extracting information with complex exponentials. Being similar to Prony and original pencil-of-function methods, it is generally preferred to those for its robustness and computational efficiency.
Roger Fuller Harrington is an American electrical engineer and professor emeritus at Syracuse University. He is best known for his contributions to computational electromagnetics with his development of method of moments (MoM). Harrington's 1968 book, Field Computation by Moment Methods, is regarded as a pivotal textbook on the subject.
In physics, the Weyl expansion, also known as the Weyl identity or angular spectrum expansion, expresses an outgoing spherical wave as a linear combination of plane waves. In a Cartesian coordinate system, it can be denoted as
Yuen Tze Lo was a Chinese American electrical engineer and academician. He was a professor emeritus at the Department of Electrical and Computer Engineering at University of Illinois at Urbana–Champaign. He is best known for his contributions to the theory and design of antennas. He is the editor of the textbook series, Antenna Handbook.
In electromagnetism, surface equivalence principle or surface equivalence theorem relates an arbitrary current distribution within an imaginary closed surface with an equivalent source on the surface. It is also known as field equivalence principle, Huygens' equivalence principle or simply as the equivalence principle. Being a more rigorous reformulation of the Huygens–Fresnel principle, it is often used to simplify the analysis of radiating structures such as antennas.
Tapan Kumar Sarkar was an Indian-American electrical engineer and Professor Emeritus at the Department of Electrical Engineering and Computer Science at Syracuse University. He was best known for his contributions to computational electromagnetics and antenna theory.
Propagation graphs are a mathematical modelling method for radio propagation channels. A propagation graph is a signal flow graph in which vertices represent transmitters, receivers or scatterers. Edges in the graph model propagation conditions between vertices. Propagation graph models were initially developed by Troels Pedersen, et al. for multipath propagation in scenarios with multiple scattering, such as indoor radio propagation. It has later been applied in many other scenarios.
The charge-based formulation of the boundary element method (BEM) is a dimensionality reduction numerical technique that is used to model quasistatic electromagnetic phenomena in highly complex conducting media with a very large number of unknowns. The charge-based BEM solves an integral equation of the potential theory written in terms of the induced surface charge density. This formulation is naturally combined with fast multipole method (FMM) acceleration, and the entire method is known as charge-based BEM-FMM. The combination of BEM and FMM is a common technique in different areas of computational electromagnetics and, in the context of bioelectromagnetism, it provides improvements over the finite element method.
Donald Robert Wilton is an American electrical engineer who is a professor emeritus at Department of Electrical and Computer Engineering at University of Houston. He is best known for his contributions to the field of computational electromagnetics.
Kenneth Kwan Hsiang Mei was a Chinese-American electrical engineer and academic, who was a professor emeritus at Department of Electrical Engineering and Computer Science at University of California, Berkeley. From 1994 until his death, he was also a professor of electrical engineering at City University of Hong Kong. He is best known for his contributions to computational electromagnetics.