Simple Quasi Static Moment Method Analysis of Microstrip

Simple Quasi Static Moment Method Analysis of a Microstrip
Computational electromagnetics (CEM) can provide accurate solutions for Z0 and other microstrip properties of interest including plots of `E and `H everywhere in space and `Js and rs on the land or the ground plane. This can be accomplished regardless of the cross-sectional geometry of the microstrip, the thickness of the land or its conductivity. The method of moments (MM) is a very popular CEM technique. It is particularly useful for planar geometries such as microstrip, stripline, conformal antennas, etc. The MM was popularized by R. F. Harrington in 1965 with his book “Field Computations by Moment Methods.” Today, it is one of the most widely used CEM techniques. We’ll illustrate the MM technique with a solution to a quasistatic microstrip immersed in an infinite dielectric as shown:

That is, there is no substrate, per se.
Integral Equation
We’ll imagine that a time harmonic voltage source has been applied across the two conductors:

This causes a charge accumulation as shown. Next, the image method will be employed to create an equivalent problem for the fields in the upper half space ( y ³ 0):

In a previous EM course, you’ve likely learned that the electric potential ?e at a point r in a homogeneous space produced by aline charge density r l (`r) is given by

(For example, see J. Van Bladel, Electromagnetic Fields. New York: Hemisphere Publishing, 1985.) It is very important to realize that this contour C?must include all charge densities in the space, which means we must include both conductors in this integral. To develop an equation from which we can solve for the charge density, we’ll apply the boundary condition

Now, using (2) in (1) and accounting for both the +r l and -rl strips yields

Recall that the unknown in (3) is the line charge density r l . But how do we solve for this function? It varies along the strip so we can’t simply “pull” it out of the integral. Actually, (3) is called an integral equation because the unknown function is located in an integrand. You most likely haven’t encountered such equations before. Integral equations are very difficult to solve analytically. We’ll use a numerical solution method instead.
Basis Function Expansion
In the moment method, we first expand rl in a set of basis functions. For a simple MM solution, here we’ll use pulse basis functions and divide the strips into N uniform sections:


What is not known in (4) are the amplitudes a n of the line charge density expansion. These are just numbers. So, instead of directly solving for the spatial variation of r L in (3), now we’ll just be computing these N numbers, a n. Much simpler! However, we need to allow enough “degrees of freedom” in this basis function expansion (4) so that an accurate solution can be found. This is accomplished by choosing the proper type of expansion functions, a large enough N, etc. The next step in the MM solution is to substitute


and is called the Green’s function. We can interchange the order of integration and summation in (5) since these are linear operators, except perhaps when x = x In this case, the integrand becomes singular. We’ll consider this situation later in this lecture. Then, (5) becomes


Testing the Integral Equation
In (7), we have N unknown coefficients a n to solve for, but only a single equation. We will generate a total of N equations by evaluating (7) at N points along the (top) strip. This process is called “testing” the integral equation. We’ll test (7) at the centers of each of the N segments, xm , giving

This is the final system of equations that we will use to solve for all the coefficients r n .
Matrix Equation
It is helpful to cast (8) into the form of a matrix equation

Vm = V
an = r n

The numerical solution to (9) is accomplished by “filling” or “populating” [V] and [Z], then solving a system of linear, constant coefficient equations. In particular, for
• [V] – choose V = 1 V in (10a), for example.
• [Z] – compute (10c) analytically, if possible, or by numerical integration.
The “filling” of [V] is very simple, while filling [Z] is a bit more difficult. In this quasi-static microstrip example, though, it is possible to evaluate all of the terms analytically since a simple anti-derivative is available. In particular, with the center of the strip located at the origin as shown:

then the electrostatic potential at point r produced by a strip of width Ñ supporting a constant line charge density r l is given by

This integral can be evaluated analytically since the integrand, as it turns out, has a relatively simple anti-derivative. Performing the integration in (11), we find that when r does not lie anywhere on the strip, the potential is

while if x = y = 0 (at the center of the strip), then

Using (12) and (13) in (10c), it can be shown that for m ¹ n

where Ñmn º xm - xn while if m = n, then