Method of moments (electromagnetics)

Simulation of negative refraction from a metasurface at 15 GHz for different angles of incidence. The simulations are performed through the method of moments.

The method of moments (MoM), also known as the moment method and method of weighted residuals,Template:Sfn 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,Template:Efn 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.Template:Efn It is one of the most common methods in microwave and antenna engineering.

History

Development of boundary element method and other similar methods for different engineering applications is associated with the advent of digital computing in the 1960s. Prior to this, variational methods were applied to engineering problems at microwave frequencies by the time of World War II. While Julian Schwinger and Nathan Marcuvitz have respectively compiled these works into lecture notes and textbooks, Victor Rumsey has formulated these methods into the "reaction concept" in 1954. The concept was later shown to be equivalent to the Galerkin method. 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.

A scheme and radiation pattern of a log-spiral antenna, designed with a NEC-based modeling software

In the 1960s, early research work on the method was published by Kenneth Mei, Jean van Bladel and Jack Richmond. In the same decade, the systematic theory for the method of moments in electromagnetics was largely formalized by Roger Harrington. While the term "the method of moments" was coined earlier by Leonid Kantorovich and Gleb Akilov for analogous numerical applications, Harrington has adapted the term for the electromagnetic formulation. Harrington published the seminal textbook Field Computation by Moment Methods on the moment method in 1968. 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.Template:Sfn In the 1990s, introduction of fast multipole and multilevel fast multipole methods enabled efficient MoM solutions to problems with millions of unknowns.Template:Sfn

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.Template:Sfn Many non-commercial and public domain codes of different sophistications are also available.Template:Sfn In addition to its use in electrical engineering, the method of moments has been applied to light scattering and plasmonic problems.

Background

Basic concepts

Template:See also An inhomogeneous integral equation can be expressed as: $$L(f)=g$$ where Template:Math denotes a linear operator, Template:Math denotes the known forcing function and Template:Math denotes the unknown function. Template:Math can be approximated by a finite number of basis functions ($f_n$): $$f\approx \sum _n^N{a}_n{f}_n.$$

By linearity, substitution of this expression into the equation yields: $$\sum _n^N{a}_{n}L(f_n)\approx g.$$

We can also define a residual for this expression, which denotes the difference between the actual and the approximate solution: $$R=\sum _n^N{a}_{n}L(f_n)-g$$

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.Template:Sfn After the determination of a suitable inner product for the problem, the expression then becomes: $$\sum _n^N{a}_n⟨w_m,L(f_n)⟩\approx ⟨w_m,g⟩$$

Thus, the expression can be represented in the matrix form: $$\left[{\ell }_{mn}\right]\left[{\alpha }_m\right]=[g_n]$$

The resulting matrix is often referred as the impedance matrix. The coefficients of the basis functions can be obtained through inverting the matrix.Template:Sfn For large matrices with a large number of unknowns, iterative methods such as conjugate gradient method can be used for acceleration.Template:Sfn The actual field distributions can be obtained from the coefficients and the associated integrals.Template:Sfn The interactions between each basis function in MoM is ensured by Green's function of the system.Template:Sfn

Basis and testing functions

Interpolation of function with rooftop basis functions

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.Template:Sfn Choice of Dirac delta function as basis function is known as point-matching or collocation. This corresponds to enforcing the boundary conditions on $N$ discrete points and is often used to obtain approximate solutions when the inner product operation is cumbersome to perform.Template:SfnTemplate:Sfn Other subsectional basis functions include pulse, piecewise triangular, piecewise sinusoidal and rooftop functions.Template:Sfn Triangular patches, introduced by S. Rao, D. Wilton and A. Glisson in 1982, are known as RWG basis functions and are widely used in MoM.Template:Sfn Characteristic basis functions were also introduced to accelerate computation and reduce the matrix equation.

The testing and basis functions are often chosen to be the same; this is known as the Galerkin method.Template:Sfn 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.

Integral equations

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.Template:Sfn For surface problems, common integral equation formulations include electric field integral equation (EFIE), magnetic field integral equation (MFIE)Template:Sfn and mixed-potential integral equation (MPIE).Template:Sfn

Thin-wire equations

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.Template:Sfn Pocklington's equation precedes the computational techniques, having been introduced in 1897 by Henry Cabourn Pocklington. For a linear wire that is centered on the origin and aligned with the z-axis, the equation can be written as: $$\underset{-l/2}{\overset{l/2}{\int }}{I}_z(z^{\prime })[(\frac{d^2}{d{z}^2}+{\beta }^2)G(z,z^{\prime })]d{z}^{\prime }=-j\omega \epsilon E_z^{\text{inc}}(p=a)$$ where $l$ and $a$ denote the total length and thickness, respectively. $G(z,z^{\prime })$ is the Green's function for free space. The equation can be generalized to different excitation schemes, including magnetic frills.Template:Sfn

Hallén integral equation, published by E. Hallén in 1938, can be given as: $$(\frac{d^2}{d{z}^2}+{\beta }^2)\underset{-l/2}{\overset{l/2}{\int }}{I}_z(z^{\prime })G(z,z^{\prime })d{z}^{\prime }=-j\omega \epsilon E_z^{\text{inc}}(p=a)$$

This equation, despite being more well-behaved than the Pocklington's equation,Template:Sfn is generally restricted to the delta-gap voltage excitations at the antenna feed point, which can be represented as an impressed electric field.Template:Sfn

Electric field integral equation (EFIE)

Template:Main The general form of electric field integral equation (EFIE) can be written as: $$\widehat{𝐧}\times {𝐄}^{\text{inc}}(𝐫)=\widehat{𝐧}\times \underset{S}{\int }[\eta jk𝐉({𝐫}^{\prime })G(𝐫,{𝐫}^{\prime })+\frac{\eta }{jk}\{{{\mathrm{\nabla }}_s}^{\prime }\cdot 𝐉({𝐫}^{\prime })\}{\mathrm{\nabla }}^{\prime }G(𝐫,{𝐫}^{\prime })]d{S}^{\prime }$$ where ${𝐄}_{\text{inc}}$ is the incident or impressed electric field. $G(r,r^{\prime })$ is the Green function for Helmholtz equation and $\eta$ represents the wave impedance. The boundary conditions are met at a defined PEC surface. EFIE is a Fredholm integral equation of the first kind.Template:Sfn

Magnetic field integral equation (MFIE)

Another commonly used integral equation in MoM is the magnetic field integral equation (MFIE), which can be written as: $$-\frac{1}{2}𝐉(r)+\widehat{𝐧}\times {\oint }_S𝐉(r^{\prime })\times {\mathrm{\nabla }}^{\prime }G(r,r^{\prime })d{S}^{\prime }=\widehat{𝐧}\times {𝐇}_{\text{inc}}(r)$$

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.Template:Sfn

Other formulations

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.Template:Sfn 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.Template:Sfn Another approach, the volumetric integral equation, necessitates the discretization of the volume elements and is often computationally expensive.Template:Sfn

MoM can also be integrated with physical optics theoryTemplate:Sfn and finite element method.

Green's functions

A microstrip scheme. MoM analysis of such layered structures necessitates the derivation of appropriate Green's functions.

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.Template:Sfn

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.Template:Sfn 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. 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.Template:Sfn Common approaches for evaluating these integrals include tail extrapolation approaches such as weighted-averages method.

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. 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. The associated Green's functions are referred as closed-form Green's functions. The method has also been extended for cylindrically-layered structures.

Rational-function fitting method, as well as its combinations with DCIM, 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. 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.

See also

• Boundary element method
• Characteristic mode analysis
• Discrete dipole approximation
• Fast multipole method
• Finite element method
• Multilevel fast multipole method

References

Bibliography

• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book
• Template:Cite book