Electromagnetics


There is a joke, or saying, when it comes to designing an oscillator or amplifier:

“amplifiers oscillate, and oscillators amplify”

The statement describes the fact that, when attempting to design one, there are many pitfalls that cause you to have the other happen.  Usually when designing an amplifier oscillations are your enemy, and you tune them out with bypass capacitors and work to ensure the stability criteria is met.  When designing an oscillator, on the other hand, the circuit can readily swing into a stable state if you don’t have a very good design, and end up as just an amplifier.  Granted, with proper design and construction this is avoidable, but when you are building one for the first time or just quickly building one, this can happen.

I propose a similar saying for antenna design, although to make this work, i am going to use “radiator” instead of “antenna”:

“Radiators couple, and couplers radiate”

This only applies, really, to antennas with more than one port, and really can apply to any microwave circuit.  But this way the phrasing works out much better.  But basically says that usually when designing an antenna, you are battling this coupling between ports, trying to get efficient radiation.  Yet on the other hand, when designing a microwave circuit, one of your enemies becomes radiation, and you work to keep that to an absolute minimum, because unwanted radiation is unwanted loss.

One of my favorite things to point out to friends of mine is this fact:

Remember in Back to the Future (I, II, and III) Doc Brown used a “flux capacitor” to travel through time?  Well, what if I told you that every capacitor is a flux capacitor?

Well that’s true.   As a simple example, say I took a capacitor and hooked it up to a battery, as below.  Then I show the capacitor as two metal plates with electric field between them, and a surface “S” parallel to the two plates.

fluxcapacitor.JPG

As this crude analysis points out, the flux is nonzero so long as the voltage across the capacitor – any capacitor  – is nonzero.

Movie fans will point out that this isn’t magnetic flux, to which I have two rebuttals:

1.) the movie never states whether Doc Brown means electric or magnetic flux

2.) fine, you want a magnetic flux capacitor?  I give you, the mystical LC network – the transmission line stub, the antenna, the cavity resonator…etc.

Sorry Robert Zemeckis.

Every now and then I write or modify a matlab function that I need for my research, and will post them to this blog here and there, with the hopes that they will be useful to someone doing antenna/microwave design.  I know, there is the Matlab File Exchange (which is very useful!) , but hey, this is my blog, I want to post my content here.

Below is a function I used to plot smith chart results in matlab.  I used basic plotting code to generate the chart itself, and added a simply plot function to add the impedance locus and constant VSWR circles.  Hope this is helpeful!  Please let me know if anyone finds a bug.  I will note that there is no safeguard in the function for s-parameters that are greater than 1, which if you are plotting the active S-parameters of a multiport device, are possible.  This isn’t a big deal when the S-parameter goes slightly above 1, but if it swings well above 1 you end up with a tiny smith chart and these erratic line segments…it’s a mess.  For passive applications I haven’t come across a bug yet.

smithchart.m


This article from NASA talks about the measurement of magentic portals that connect the earth to the sun every 8 minutes, due to their magnetic fields pushing against each other and forming connections.  It is always interesting to see something happening on a such a large scale and with these huge fields that, as the article points out, many people didn’t believe in only 10 years ago.
I find the interaction of the electromagnetic waves on the scale of the solar system, between immensely large objects (here the earth and the sun) to be fascinating.  I spend so much time dealing with such small scale fields – antennas, microwave circuits, where fields are coerced into existence in these carefully designed structures to behave in a certain way.  Then I read articles on fields stretching 100 million miles between the earth and the sun, or the various effects the sun’s fields have on the earth, such as solar flares – fields covering half our solar system, in unbounded space, whose origins are not always clear or understood.  I mean, the earth’s magnetic field seems such a simple field that we take for granted – yet why it exists in the first place is not fully understood – you can’t exactly cut the earth in half and take a peak at it’s cross-section. It has even been shown to have flipped thousands of times over it’s lifetime – also something that we don’t quite understand the reason for.

I got an email informing me that my masters thesis is now accessible through the University of Massachusetts “scholarworks” site, their digital library of thesis and disseratations.

If you are interested in taking a look at what I spent the last two years working on…you can find it  here.

The thesis was on a project aiming to design GPS antennas that fit into a volume of approximately 1.25X1.25X0.2″.

I have been running some Ansoft HFSS simulations of infinite arrays of antennas, recently.  A peculiar thing I noticed is that some of the impedance data, which is calculated from the S-parameter matrices, had sizable spikes (upwards of 300Ohms) inside a frequency band that, from the S-parameters, looked smooth and well matched.  Clearly something was wrong.  A low VSWR and S11 imply impedance levels close to the source impedance, yet these spikes in the impedance were there.

So I decided that, being able to compute all the data I need from S-parameters alone (as HFSS does), I exported the entire S-matrix to MATLAB and took a look at the impedance values I calculated.

HFSS gives the equation for finding the Z matrix as :

Z = sqrt(Zo*I)*inv(I-S)*(I+S)*sqrt(Zo*I)

where I=identity matrix, and S is the full S-parameter matrix

from this I plotted the exact same thing HFSS was giving for the impedances.  So I started plotting different chunks of the equation (e.g. inv(I-s) alone, (I+S) alone…) and found that the inv(I-S) term was spiking at the exact places my impedance was spiking, while the other terms were well-behaved.  Well, looks like that inversion is messing things up…indicating a singular, or more correctly a nearly-singular matrix (singular matrices have NO inverse defined).  Sure enough, using the function:

cond()

in MATLAB revealed huge condition numbers for the matrix (I-S) where the impedance spikes occurred, 2 orders of magnitude larger than the condition numbers calculated over the rest of the frequency sweep.  So I need to invert a matrix that is nearly-singular….now what??  I need single-value-decomposition (SVD) to take out those values that make the matrix singular and do a quasi-inversion, basically approximate the inversion.  MATLAB includes just such a function:

pinv() — where it’s use is  pinv(matrix, tolerance)

The tolerance can be set to an arbitrary value.  For my purposes, I can stand to lose some accuracy, so I chose a relatively large tolerance (around 0.02 or so) to take out those spikes but preserve the overal shape of the impedance waveforms.  Usually you will want to keep this very small to keep the inversion close to the true value.  But this solves my problem.

But it leads to another conclusion — HFSS, apparently, doesn’t test if its S-parameter matrices are nearly-singular (or at least I haven’t seen a way to have it do so).  Certainly it is using all sorts of inversion approximations in the solver the software uses (you can even do iterative matrix solutions in v11 of HFSS…), I don’t think it is likely that they have skimped on such checking on the data output front end.   Please post comments if you konw of any solutions in HFSS for this.

Here is an interesting Java applet that allows good visualization of the fields of some basic antenna elements.

 http://www.falstad.com/antenna/

You can readily see the field pattern distortions that occur from things like making a dipole or loop longer than half a wavelength in size, the multiple lobes that form, and how the current distribution forms on the structure.  You can for from a triangular current distribution on the dipole, to 10 or 20 sine wave variations along the length.  Really easy way to see what happens to the pattern.  There are also some array examples of end-fire arrays and schelkunoff polynomial alrrays.

My friend Bill Whirity (see link to his website on the side) is making a movie (look for his other movie Broke on netflix!) called The MisInvetntions of Milo Weatherby.  From what I have seen of it, it will be awesome.  It is about these two kids and their, well, misinventions with science. They just finished shooting it, so no idea when it will be out.  Here is a great video from it:

http://billwhirity.com/milo/videos/microwave_implosion.mov

UPDATE:

http://www.billwhirity.com/milo/index.html

Has a trailer and information about the movie.  I have seen it, and it is excellent!

I passed my oral thesis defense today, so I currently have a masters!

In celebration, watch this:

Today! 10AM. I have put a lot of work over the past few months into putting together my thesis, and I am hoping that will pay off.

I found this xkcd comic strangely appropriate — hopefully I am not subconsciously doing something like this during my presentation:

math_paper.png

I guess there aren’t imaginary antennas…

Next Page »