Basic Principles of The Method of Moments - Revision history

Kazem Sabet at 21:19, 23 July 2018

2018-07-23T21:19:38Z

Kazem Sabet: /* Numerical Solution of Integral Equations */

2018-07-15T23:15:52Z

‎Numerical Solution of Integral Equations

Kazem Sabet: /* The Rectangular Mesh Advantage */

2017-06-11T16:07:19Z

‎The Rectangular Mesh Advantage

Kazem Sabet: /* The Rectangular Mesh Advantage */

2017-06-11T16:06:52Z

‎The Rectangular Mesh Advantage

Kazem Sabet: /* Discretization of Electric & Magnetic Currents */

2017-06-11T16:05:06Z

‎Discretization of Electric & Magnetic Currents

Kazem Sabet: /* Discretization Of Electric & Magnetic Currents */

2017-06-11T16:03:10Z

‎Discretization Of Electric & Magnetic Currents

Kazem Sabet: /* Multilayer Green’s Functions */

2017-06-11T16:01:20Z

‎Multilayer Green’s Functions

Kazem Sabet: /* Multilayer Green’s Functions */

2017-06-11T16:00:52Z

‎Multilayer Green’s Functions

Kazem Sabet: /* Discretization Of Wire Structures */

2017-06-11T15:59:38Z

‎Discretization Of Wire Structures

Kazem Sabet: /* Discretization Of Wire Structures */

2017-06-11T15:58:46Z

‎Discretization Of Wire Structures

@@ Line 5: / Line 5: @@
 <tr>
 </table>
 == An Overview of the Method of Moments ==

@@ Line 171: / Line 171: @@
 == Numerical Solution of Integral Equations ==
-The planar integral equations derived earlier can be solved numerically by discretizing the unknown currents using a proper meshing scheme. The original functional equations are reduced to discretized linear algebraic equations over elementary cells. The unknown quantities are found by solving this system of linear equations, and many other [[parameters]] can be computed thereafter. This method of numerical solution of integral equations is known as the Method of Moments (MoM). In this method, the unknown electric and magnetic currents are represented by expansions of basis functions as follows:
+The planar integral equations derived earlier can be solved numerically by discretizing the unknown currents using a proper meshing scheme. The original functional equations are reduced to discretized linear algebraic equations over elementary cells. The unknown quantities are found by solving this system of linear equations, and many other parameters can be computed thereafter. This method of numerical solution of integral equations is known as the Method of Moments (MoM). In this method, the unknown electric and magnetic currents are represented by expansions of basis functions as follows:
 :<math>J(r) = \sum_{n=1}^N I_n^{(J)} f_n^{(J)} (r)</math>

@@ Line 245: / Line 245: @@
 <table>
 <tr>
-<td> [[File:PMOM25.png|thumb|480px|Pairs of rooftop basis functions that have identical MoM interactions.]] </td>
+<td> [[File:PMOM25.png|thumb|360px|Pairs of rooftop basis functions that have identical MoM interactions.]] </td>
 </tr>
 </table>

@@ Line 224: / Line 224: @@
 :<math> Z_{ij}^{(\mu \nu)} = \iiint_{V_i} d\nu f_i^{(\mu)}(r) \cdot \iiint_{V_j}d\nu ' \overline{\overline{G}}_{\mu \nu}(r|r') \cdot f_j^{(v)}(r') </math>
 where the spatial-domain dyadic Green's functions are a function of the observation and source coordinates, '''r'''and '''r' '''. The MoM matrix elements can indeed be interpreted as interactions between two elementary basis functions '''f<sub>i</sub>(r)''' and '''f<sub>j</sub>(r')''' on that particular background structure. The spatial-domain dyadic Green's functions can themselves be expressed in terms of the spectral-domain dyadic Green's functions as follows:
 :<math> \overline{\overline{G}}_{\mu \nu}(r|r') = \frac{1}{(2\pi)^2} \int\limits_{-\infty}^{\infty} \int\limits_{-\infty}^{\infty}  \tilde{\overline{\overline{G}}}_{\mu \nu} (k_p, z|z') e^{-j[k_x(x-x')+k_y(y-y')]} \, dk_x \, dk_y , \quad {k_p}^2 = {k_x}^2 + {k_y}^2 </math>
 where the doubly infinite integration is performed with respect to the spectral variables k<sub>x</sub> and k<sub>y</sub>. As can be seen from the above expression, the spatial-domain dyadic Green's functions are functions of z, z', as well as (x-x') and (y-y'). The MoM matrix elements can now be transformed into the spectral domain as
 :<math> Z_{ij}^{(\mu \nu)} = \dfrac{1}{(2\pi)^2} \int\limits_{-\infty}^{\infty} \int\limits_{-\infty}^{\infty} \tilde{f}_i^{(\mu)} (k_x, k_y) \cdot \tilde{\overline{\overline{G}}}_{\mu \nu} (k_{\rho}, z|z') \cdot \tilde{f}_j^{(\nu)} (k_x, k_y) \, dk_x \, dk_y </math>
 where the tilde symbol signifies the Fourier transform of a function defined as
 :<math> \tilde{f}(k_x, k_y) = \dfrac{1}{(2\pi)^2} \int\limits_{-\infty}^{\infty} \int\limits_{-\infty}^{\infty} f(x,y) e^{j(k_x x + k_y y)} \, dx \, dy </math>
 Rectangular cells have simple Fourier transforms. The rooftop basis functions are triangular functions in the direction of current flow and constant in the perpendicular direction. This means that their Fourier transform is a product of a sinc-squared function along one spectral direction and a sinc function along the other. You can see from the figure below that if one deals with a rectangular mesh of identical cells (all equal and parallel), then the interactions among the rooftop basis functions become a functions of the index differences and not the absolute indices:
 :<math> Z_{(i,k)|(j,l)} = Z \Big\langle f_{i,k}(x,y)| f_{j,l}(x', y') \Big\rangle = Z_{(i-j)|(k-l)} </math>
 In the above equation, the vectorial rooftop basis functions have explicit, double indices: i and k along the local X and Y directions, respectively, for the test (observation) basis function, and j and l along the local X and Y directions, respectively, for the expansion (source) basis function. Thus, uniform rectangular cells, i.e. structured rectangular cells of identical size aligned in the same direction, can speed up the planar MoM simulation significantly due to these symmetry and the invariance properties. For example, all the self-interactions are identical regardless of the location of a rooftop basis function. This reduces the matrix fill process for a total of N rooftop basis functions from an N2 process to one of order N.
-[[File:PMOM25.png]]
+<table>
+<tr>
-Pairs of rooftop basis functions that have identical MoM interactions.
+<td> [[File:PMOM25.png|thumb|480px|Pairs of rooftop basis functions that have identical MoM interactions.]] </td>
 == Computing The Near Fields in Planar MoM ==

@@ Line 212: / Line 212: @@
 The volume polarization currents in 2.5-D MoM have a vertical direction along the Z-axis. These are discretized using prismatic basis functions that have either a rectangular or triangular base with a constant profile along the Z-axis.
 <table>
 <tr>
-<td> [[File:image065_tn.png|thumb|480px|Prismatic basis functions built over single triangular and rectangular cells.]] </td>
+<td> [[File:image065_tn.png|thumb|480px|Prismatic basis function with triangular cross section.]] </td>
 </tr>
 </table>
 == The Rectangular Mesh Advantage ==

@@ Line 198: / Line 198: @@
 Similar expressions can be derived for the T<sup>(EM)</sup>, U<sup>(HJ)</sup> and Y<sup>(HM)</sup>elements of the MoM matrix.
-== Discretization Of Electric &amp; Magnetic Currents ==
+== Discretization of Electric &amp; Magnetic Currents ==
 The right choice of the basis functions to represent the elementary currents is very important. It will determine the accuracy and computational efficiency of the resulting numerical solution. Rooftop basis functions are one of the most popular types of basis functions used in a variety of MoM formulations. The surface currents (whether electric or magnetic) are discretized using 2D rooftop basis functions shown in the figure below:
-[[File:image055_tn.png]]
+<table>
+<tr>
-Rooftop or RWG basis functions built over two rectangular, triangular or mixed cells.
+<td> [[File:image055_tn.png|thumb|480px|Rooftop or RWG basis functions built over two rectangular, triangular or mixed cells.]] </td>
 The rooftop basis functions are defined over two adjacent cells with a common edge of length. If the two cells are triangular, then the so-called RWG functions are obtained. It is also possible to define rooftop functions over two adjacent rectangular cells or two adjacent rectangular and triangular cells with a common edge. On a rectangular cell, the function is defined as having a (descending or ascending) linear profile in one direction and a constant profile in the other perpendicular direction.
@@ Line 211: / Line 213: @@
 [[File:image065_tn.png]][[File:image066_tn.png]]
 Prismatic basis functions built over single triangular and rectangular cells.

@@ Line 129: / Line 129: @@
 <table>
 <tr>
-<td> [[Image:PMOM14.png|thumb|480px|A typical planar layered structure.]] </td>
+<td> [[Image:PMOM14.png|thumb|360px|A typical planar layered structure.]] </td>
 </tr>
 </table>

@@ Line 125: / Line 125: @@
 == Multilayer Green’s Functions ==
 The Green’s functions are the solutions of boundary value problems when they are excited by an elementary source. This is usually assumed to be an infinitesimally small vectorial point source. In order for Green’s functions to be computationally useful, they must have analytical closed forms like a mathematical expression, or one should be able to compute them using a recursive process. It turns out that only very few boundary value problems have closed-form Green’s functions. Planar layered structures with laterally infinite extents are one of those few cases, which can be represented by recursive dyadic Green's functions.
 In general, a structure may support both electric ('''J''') and magnetic ('''M''') currents. The total electric ('''E''') and magnetic ('''H''') fields can be expressed in terms of the electric and magnetic currents in the following way:
@@ Line 133: / Line 138: @@
 :<math>H = H^{inc} + \iiint\limits_V \overline{\overline{G}}_{HJ}(r|r') \cdot J(r') \, dv' + \iiint\limits_V \overline{\overline{G}}_{HM}(r|r') \cdot M(r') \, dv'</math>
 where '''G<sub>EJ</sub>''', '''G<sub>EM</sub>''', '''G<sub>HJ</sub>''', '''GH<sub>M</sub>''' are the dyadic Green’s functions for the electric and magnetic currents due to electric and magnetic current source, respectively, and '''E<sup>i</sup>''' and '''H<sup>i</sup>''' are the incident or impressed electric and magnetic fields, respectively. In these equations, '''r''' is the position vector of the observation point and '''r'''' is the position vector of the source point. V is the volume that contains all the sources and the volume integration is performed with respect to the primed coordinates. The incident or impressed fields provide the excitation of the structure. They may come from an incident plane wave or a gap source on a microstrip line, a short dipole, etc. The complexity of the Green’s functions depends on what is considered as the background structure. If you remove all the unknown currents from the structure, you are left with the background structure.

@@ Line 113: / Line 113: @@
 <table>
 <tr>
-<td> [[Image:19_meshing_tn.gif|thumb|480px|Meshing a wire with rooftop (triangular) basis functions.]] </td>
+<td> [[Image:19_meshing_tn.gif|thumb|360px|Meshing a wire with rooftop (triangular) basis functions.]] </td>
 </tr>
 </table>

@@ Line 103: / Line 103: @@
 The right choice of the basis functions that are used to represent the elementary currents is very important. It will determine the accuracy and computational efficiency of the resulting numerical solution. Rooftop basis functions are one of the more popular types of basis functions used in a variety of MoM formulations. The simplest rooftop function is the one-dimensional triangular functions defined as in the figure below:
-[[File:18_meshing_tn.gif]]
+<table>
 This function provides a linear interpolation of the unknown currents or fields in one dimension. Note that the function vanishes at it two ends. This is a desirable feature for basis functions that represent electric currents on metallic wires as the current must vanish at the two ends of a wire. The total current on the wire can be approximated in a linear fashion by a set of one-dimensional rooftop functions as shown in the figure below:
-[[File:19_meshing_tn.gif]]
+<table>
 This can be written as
 :<math> I(l) = \sum_{n=1}^N a_n f_n(l) \mathbf{\hat{s}_n} </math>
 where l is the length coordinate along the wire with l=0 at its start point. <math>f_n(l)</math> is the scaled and translated version of the linear basis function <math>f(l)</math> shown in the previous figure. <math>\mathbf{\hat{s}_n}</math> is the unit vector along wire.