## Method of Moments Solver for Metal and Dielectric Structures

Method of Moments computation technique for metal and dielectric antennas.

Antennas using dielectric substrate consists of a metal part and a dielectric part. The first step in the computational solution of electromagnetic problems is to discretize Maxwell's equations. The process results in this matrix-vector system:

`$V=ZI$`
• V — Applied voltage vector. This signal can be voltage or power applied to the antenna or an incident signal falling on the antenna.

• I — Current vector that represents current on the antenna surface.

• Z — Interaction matrix or impedance matrix that relates V to I. For calculating the interaction matrix, the effect of metal and dielectric parts in an antenna are taken separately.

Antenna Toolbox™ uses method of moments (MoM) to calculate the interaction matrix and solve system equations.

### MoM Formulation

The MoM formulation is split into three parts.

#### Discretization of Dielectrics

Discretization enables the formulation from the continuous domain to the discrete domain. This step is called meshing in antenna literature. In the MoM formulation, the metal surface of the antenna is meshed into triangles and the dielectric volume is meshed into tetrahedrons.

#### Basis Functions

Basis functions are used to represent unknown quantities. In the case of antennas using dielectrics, the unknown quantities are the surface current on the metal structure and flux density due to dielectric volume. Antenna Toolbox uses Rao-Wilton-Glisson (RWG) [2] basis functions. For basis functions for a metal structure in an antenna refer, Method of Moments Solver for Metal Structures.

For the dielectric volume of the antenna, Antenna Toolbox uses a zeroth order edge basis function to model the flux density.

The figure shows an edge-based basis function. The vector variation is perpendicular to the base edge AB (or $\overline{l}$). The vector of the edge CD (or $\overline{p}$) defines the basis function. Within a tetrahedron, the basis function is a constant field given by

`$\stackrel{\to }{f}=c\stackrel{\to }{p}$`

• c – normalization coefficient.

• p – vector of the edge defining the basis function.

#### Interaction Matrix

The interaction matrix is a complex dense symmetric matrix. For a metal-dielectric antenna, there are two sets of basis functions and four interactions. To fill out the interaction matrix, calculate the free-space Green's function between all the basis functions on the antenna surface. The final interaction matrix equations are:

• `ZMM` – metal to metal interaction. For a pure metal structure, you only calculate this symmetric square matrix.

`${Z}_{mn}^{MM}=\left(\frac{j\omega \mu }{4\pi }\right)\underset{S}{\int }\underset{S}{\int }{\stackrel{\to }{f}}^{M}{}_{m}\left(\stackrel{\to }{r}\right).{\stackrel{\to }{f}}^{M}{}_{n}\left(\stackrel{\to }{{r}^{\prime }}\right)gd\stackrel{\to }{{r}^{\prime }}d\stackrel{\to }{r}-\left(\frac{j}{4\pi \omega \epsilon }\right)\underset{S}{\int }\underset{S}{\int }\left(\nabla .{\stackrel{\to }{f}}^{M}{}_{m}\right)\left(\nabla .{\stackrel{\to }{f}}^{M}{}_{n}\right)gd\stackrel{\to }{{r}^{\prime }}d\stackrel{\to }{r}$`
• `ZDD` – dielectric to dielectric interaction. For pure dielectric structures, you only calculate this symmetric square matrix.

`$\begin{array}{l}{\stackrel{^}{Z}}_{mn}^{DD}=\sum _{p-1}^{P}\sum _{{p}^{\prime }-1}^{{P}^{\prime }}\frac{{K}_{p}}{{\stackrel{^}{\epsilon }}_{p}}\underset{{V}_{D}}{\int }{\stackrel{\to }{f}}_{mp}\left(\stackrel{\to }{r}\right)\cdot \text{ }\text{ }{\stackrel{\to }{f}}_{n{p}^{\prime }}\left(\stackrel{\to }{r}\right)d\stackrel{\to }{r}\\ \text{\hspace{0.17em}}\text{ }\text{\hspace{0.17em}}\text{ }-\frac{{\omega }^{2}{\mu }_{0}}{4\pi }\sum _{p-1}^{P}\sum _{{p}^{\prime }-1}^{{P}^{\prime }}{K}_{p}{K}_{{p}^{\prime }}\underset{{V}_{D}}{\int }\underset{{V}_{{D}^{\prime }}}{\int }g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right){\stackrel{\to }{f}}_{mp}\left(\stackrel{\to }{r}\right)\cdot {\stackrel{\to }{f}}_{n{p}^{\prime }}\left(\stackrel{\to }{{r}^{\prime }}\right)d\stackrel{\to }{r}d\stackrel{\to }{{r}^{\prime }}\\ \text{ }\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{4\pi {\epsilon }_{0}}\sum _{q-1}^{Q}\sum _{{q}^{\prime }-1}^{{Q}^{\prime }}{\stackrel{^}{K}}_{q}{\stackrel{^}{K}}_{{q}^{\prime }}\underset{{\Omega }_{q}}{\int }\underset{{\Omega }_{{q}^{\prime }}}{\int }g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right){f}_{\perp mq}\left(\stackrel{\to }{r}\right){f}_{\perp n{q}^{\prime }}\left(\stackrel{\to }{{r}^{\prime }}\right){d}_{s}{d}_{{s}^{\prime }}\text{\hspace{0.17em}}\text{ }\text{ }\text{ }m,n=1,....,N\end{array}$`
• `ZMD` and `ZDM` – These matrices calculate the interaction between metal and dielectric. This matrix is not a symmetrical square matrix.

`$\begin{array}{l}{Z}_{mn}^{MD}=-\frac{{\omega }^{2}{\mu }_{0}}{4\pi }\sum _{p-1}^{2}\sum _{{p}^{\prime }-1}^{{P}^{\prime }}{K}_{p}\underset{t}{\int }\underset{{V}_{{D}^{\prime }}}{\int }{\stackrel{\to }{f}}_{n}^{M}\left(\stackrel{\to }{r}\right)\cdot \text{ }\text{ }{\stackrel{\to }{f}}_{m{p}^{\prime }}\left(\stackrel{\to }{{r}^{\prime }}\right)g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right)d\stackrel{\to }{{r}^{\prime }}{d}_{s}\\ \text{\hspace{0.17em}}\text{ }\text{\hspace{0.17em}}\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{4\pi {\epsilon }_{0}}\sum _{p-1}^{2}\sum _{q-1}^{Q}{\stackrel{^}{K}}_{q}\underset{{t}_{D}}{\int }\underset{{\Omega }_{{q}^{\prime }}}{\int }\left({\nabla }_{s}\cdot {\stackrel{\to }{f}}_{n}^{M}\left(\stackrel{\to }{r}\right)\right){f}_{\perp mq}\left(\stackrel{\to }{{r}^{\prime }}\right)g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right){d}_{{\Omega }^{\prime }}{d}_{s}\text{\hspace{0.17em}}\text{ }\text{ }\\ \text{ }\text{ }\text{\hspace{0.17em}}\text{ }m=1,....,{N}_{D};n=1,...,{N}_{M}\end{array}$`
`$\begin{array}{l}{Z}_{mn}^{DM}=-\frac{j\omega {\mu }_{0}}{4\pi }\sum _{p-1}^{2}\sum _{{p}^{\prime }-1}^{{P}^{\prime }}{K}_{{p}^{\prime }}\underset{{V}_{D}}{\int }\underset{{S}_{{D}^{\prime }}}{\int }{\stackrel{\to }{f}}_{np}\left(\stackrel{\to }{r}\right)\cdot {\stackrel{\to }{f}}_{m}^{M}\left(\stackrel{\to }{{r}^{\prime }}\right)g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right){d}_{{s}^{\prime }}d\stackrel{\to }{r}\\ \text{\hspace{0.17em}}\text{ }\text{\hspace{0.17em}}\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{4\pi {\epsilon }_{0}\omega }\sum _{p-1}^{2}\sum _{q-1}^{Q}{\stackrel{^}{K}}_{q}\underset{{\Omega }_{q}}{\int }\underset{{S}_{{D}^{\prime }}}{\int }{f}_{\perp nq}\left(\stackrel{\to }{r}\right)\cdot \left({\nabla }_{s}\cdot {\stackrel{\to }{f}}_{m}^{M}\left(\stackrel{\to }{{r}^{\prime }}\right)\right)g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right){d}_{{s}^{\prime }}{d}_{\Omega }\text{ }\text{ }\\ \text{ }\text{ }\text{\hspace{0.17em}}\text{ }m=1,....,{N}_{D};n=1,...,{N}_{M}\end{array}$`

where

• $g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right)=\frac{\mathrm{exp}\left(-jkR\right)}{R},R=|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|$ is the free space Green's function.

• $K=\frac{{\stackrel{^}{\epsilon }}^{±}-{\epsilon }_{0}}{{\stackrel{^}{\epsilon }}^{±}}$ is the complex dielectric constant within every tetrahedron.

• ${\stackrel{^}{K}}_{q}={K}_{+}-{K}_{-}$ is the differential contrast on every face of the tetrahedron.

For a composite metal structure, you must calculate all four matrices.

### Neighbor Region

The figure shows a typical interaction matrix for a metal structure `ZMM` with 256 basis functions.

From the interaction matrix plot, you observe that the matrix is diagonally dominant. The dielectric interaction matrix is also diagonally dominant. As you move further away from the diagonal, the magnitude of the terms decreases. This behavior is same as the Green's function behavior. The Green's function decreases as the distance between r and r' increases. Therefore, it is important to calculate the region on the diagonal and close to the diagonal accurately.

This region on and around the diagonal is called neighbor region. For a metal-dielectric antenna, the neighborhood region is based on the average size of the tetrahedron.

For neighboring region details for metal antennas refer, Method of Moments Solver for Metal Structures.

### Singularity Extraction

Along the diagonal, r and r' are identical and the defined Green's function becomes singular. To remove the singularity, extraction is performed on these terms. The equations for the singularity extraction of the `ZMM` matrix are:

`$\begin{array}{l}\underset{{t}_{p}}{\int }\underset{{t}_{q}}{\int }\left({\stackrel{\to }{\rho }}_{i}.{{\stackrel{\to }{\rho }}^{\prime }}_{j}\right)g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right)ds\text{'}ds=\underset{{t}_{p}}{\int }\underset{{t}_{q}}{\int }\frac{\left({\stackrel{\to }{\rho }}_{i}.{{\stackrel{\to }{\rho }}^{\prime }}_{j}\right)}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}ds\text{'}ds+\underset{{t}_{p}}{\int }\underset{{t}_{q}}{\int }\frac{\left(\mathrm{exp}\left(-jk|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|\right)-1\right)\left({\stackrel{\to }{\rho }}_{i}.{{\stackrel{\to }{\rho }}^{\prime }}_{j}\right)}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}ds\text{'}ds\\ \underset{{t}_{p}}{\int }\underset{{t}_{q}}{\int }g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right)ds\text{'}ds=\underset{{t}_{p}}{\int }\underset{{t}_{q}}{\int }\frac{1}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}ds\text{'}ds+\underset{{t}_{p}}{\int }\underset{{t}_{q}}{\int }\frac{\left(\mathrm{exp}\left(-jk|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|\right)-1\right)}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}ds\text{'}ds\end{array}$`

The two integrals on the right-hand side of the equations, called potential or static integrals are found using analytical results [3].

The equations for the singularity extraction of the `ZDD` matrix are:

`$\begin{array}{l}\underset{{V}_{D}}{\int }\underset{{V}_{{D}^{\prime }}}{\int }g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right)d\stackrel{\to }{r}d\stackrel{\to }{{r}^{\prime }}=\underset{{V}_{D}}{\int }\underset{{V}_{q}}{\int }\frac{1}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}d\stackrel{\to }{r}d\stackrel{\to }{{r}^{\prime }}+\underset{{V}_{D}}{\int }\underset{{V}_{q}}{\int }\frac{\left(\mathrm{exp}\left(-jk|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|\right)-1\right)}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}d\stackrel{\to }{r}d\stackrel{\to }{{r}^{\prime }}\\ \underset{{\Omega }_{q}}{\int }\underset{{\Omega }_{{q}^{\prime }}}{\int }g\left(\stackrel{\to }{r},\stackrel{\to }{{r}^{\prime }}\right)d\Omega d\Omega \text{'}=\underset{{\Omega }_{D}}{\int }\underset{{\Omega }_{q}}{\int }\frac{1}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}d\Omega d\Omega \text{'}+\underset{{S}_{D}}{\int }\underset{{S}_{q}}{\int }\frac{\left(\mathrm{exp}\left(-jk|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|\right)-1\right)}{|\stackrel{\to }{r}-\stackrel{\to }{{r}^{\prime }}|}d\Omega d\Omega \end{array}$`

### Finite Arrays

The MoM formulation for finite arrays is the same as for a single antenna element. The main difference is the number of excitations (feeds). For finite arrays, the voltage vector is now a voltage matrix. The number of columns are equal to the number of elements in the array.

For example, the voltage vector matrix for a `2x2` array of rectangular patch antenna (with and without dielectric substrate) has four columns as each antenna can be excited separately.

### Infinite Array

To model an infinite array, you change the MoM to account for the infinite behavior. To do so, you replace the free-space Green's functions with periodic Green's functions. The periodic Green's function is an infinite double summation.

Green's FunctionPeriodic Green's Function

`$\begin{array}{l}g=\frac{{e}^{-jkR}}{R}\\ R=|\stackrel{\to }{r}-{\stackrel{\to }{r}}^{\prime }|\end{array}$`

`$\begin{array}{l}{g}_{\text{periodic}}=\sum _{m=-\infty }^{\infty }\sum _{n=-\infty }^{\infty }{e}^{j{\varphi }_{mn}}\frac{{e}^{-jk{R}_{mn}}}{{R}_{mn}}\\ {R}_{mn}=\sqrt{{\left(x-{x}^{\prime }-{x}_{m}\right)}^{2}+{\left(y-{y}^{\prime }-{y}_{n}\right)}^{2}+{\left(z-{z}^{\prime }\right)}^{2}}\\ {\varphi }_{mn}\text{ }\text{ }\text{\hspace{0.17em}}\text{ }=-k\left({x}_{m}\mathrm{sin}\theta \mathrm{cos}\phi +{y}_{n}\mathrm{sin}\theta \mathrm{sin}\phi \right)\\ {x}_{m}\text{ }\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ }=m\cdot {d}_{x},\text{}{y}_{n}=n\cdot {d}_{y}\end{array}$`

dx and dy are the ground plane dimensions that define the x and y dimensions of the unit cell. θ and Φ are the scan angles.

Comparing the two Green's functions, you observe an additional exponential term that is added to the infinite sum. The Φmn accounts for the scanning of the infinite array. The periodic Green's function also accounts for the effect of mutual coupling.

## References

[1] Harringhton, R. F. Field Computation by Moment Methods. New York: Macmillan, 1968.

[2] Rao, S. M., D. R. Wilton, and A. W. Glisson. “Electromagnetic scattering by surfaces of arbitrary shape.” IEEE. Trans. Antennas and Propagation, Vol. AP-30, No. 3, May 1982, pp. 409–418.

[3] Wilton, D. R., S. M. Rao, A. W. Glisson, D. H. Schaubert, O. M. Al-Bundak. and C. M. Butler. “Potential Integrals for uniform and linear source distribution on polygonal and polyhedral domains.” IEEE. Trans. Antennas and Propagation. Vol. AP-30, No. 3, May 1984, pp. 276–281.

[4] Balanis, C.A. Antenna Theory. Analysis and Design. 3rd Ed. New York: John Wiley & Sons, 2005.