{"id":54,"date":"2010-01-18T13:18:25","date_gmt":"2010-01-18T18:18:25","guid":{"rendered":"http:\/\/blogs.vassar.edu\/magnes\/?page_id=54"},"modified":"2012-01-16T11:57:02","modified_gmt":"2012-01-16T16:57:02","slug":"maxwell-fagin","status":"publish","type":"page","link":"https:\/\/pages.vassar.edu\/magnes\/advanced-electromagnetism-phys-341\/maxwell-fagin\/","title":{"rendered":"Modeling Antennas"},"content":{"rendered":"

Final Results:<\/strong><\/h2>\n

The Matlab<\/a> simulation of my final results. \u00a0The PowerPoint Presentation<\/a> associated with the simulation.<\/p>\n

I have produced a simulation of the near and far fields of a two-element Yagi antenna. \u00a0See the previously linked .m and .ppt files for a complete discussion of methodology and modeling.<\/p>\n

Results:<\/h2>\n

\"Final <\/strong><\/p>\n

Figure 1:<\/strong> The calculated electric field strength(in dB) \u00a0in the region surrounding the antenna.\u00a0 Antenna elements are shown to scale.\u00a0 Axis are in meters. \u00a0Note the preferential transmission in one direction. \u00a0An animated version of this plot is available in the linked Matlab file.<\/p>\n

\"Final<\/p>\n

Figure 2:<\/strong> The calculated directivity of the antenna shown in figure 1. \u00a0Note the similarity to the directivity plots show in previous entries for professional models of Yagi antenna.<\/p>\n

——— 04\/21\/10 ————<\/p>\n

Project Update: Yagi Antenna<\/h1>\n

Introduction<\/span><\/p>\n

I am attempting to generate a visual representation of the EM fields produced by a directional Yagi-Uda antenna.\u00a0 I\u2019ve used Yagi-Uda antenna for direction finding before, and I would like to know why arranging antenna elements like this:<\/p>\n

\"yagi\"<\/p>\n

Produces a directional gain pattern like this:<\/p>\n

\"Image<\/p>\n

Yagi-Uda Antenna Background<\/span><\/p>\n

All derivations presented here come from the text \u201cElectromagnetic Waves and Antennas\u201d by Sophocles J. Orfanidis, chapters 14, 21 and 22, or \u201cApplied Electromegnetics\u201d by Fawwaz T. Ulaby, chapter 9.<\/p>\n

A Yagi-Uda antenna contains two types of elements: Driven and parasitic.\u00a0 Driven elements are conductors that are subject to an applied external AC potential source, such as an RF encoder or a function generator.\u00a0 Parasitic elements are left unconnected to any electrical components, and respond only to the incident EM fields of free space.\u00a0 When transmitting a signal, one element of a Yagi-Uda antenna is driven, while the rest remain parasitic.\u00a0 When receiving a signal, all the elements in the antenna remain parasitic, and the signal is measured as the voltage at the midpoint of the second longest element, as shown above.\u00a0 It should be noted that the directivity of a Yagi-Uda antenna is identical in both sending and receiving modes.<\/p>\n

\u201cElectromagnetic Waves and Antennas\u201d presents a computation method for calculating the field produced by any arbitrary antenna configuration, driven by any arbitrary current.\u00a0 It is very long and involved, and I do not yet understand how to apply the method to a multiple element antenna.\u00a0 However, I can apply it to some simpler cases, such as a single element dipole antenna.\u00a0 I hope that by modeling the fields of the component parts of a Yagi-Uda antenna, I will eventually be able to assemble the information into my final answer.<\/p>\n

Theory: Single Element dipole Antenna<\/span><\/p>\n

With most antenna problems, it is easiest to formulate a solution if we are working, not with the fields, but with the retarded potentials, A(r,t) and V(r,t) .\u00a0 Recall that for a given charge distribution \u03c1, and a given current distribution J, the retarded potentials at point r are defined such that:<\/p>\n

\"Equation1\"<\/p>\n

Where V represents all space where \u03c1\u00a0and J are non-zero, and R is the distance between the potential source r\u2019 and the point of interest r.\u00a0 Although most of our work will be done with potentials, recall that at any time, the fields E and B may be calculated from the potentials:<\/p>\n

\"Equation2\"<\/p>\n

If we assume the potentials are driven by some harmonic source (as is the case with most radio applications), then the potentials will assume a harmonic time dependence:<\/p>\n

\"Equation3\"<\/p>\n

Where k is the wave number, \u03c9\/c.\u00a0 Now let us assume a current distribution J for the antenna element of interest.\u00a0 Let us assume a thin dipole antenna on length (L) and radius (a) driven by a time dependant sinusoidal current source located at its midpoint.\u00a0 Assuming the dipole lies along the z-axis, and assuming the current must go to zero at the ends of the conductor, such a situation will produce the current distribution of:<\/p>\n

\"Equation4\"<\/p>\n

Where the 2\u03c0a factor comes from the fact that the current will be distributed around the surface of the dipole. \u00a0Evaluating the magnitude of A for such a current distribution and collecting all constants produces a potential:<\/p>\n

\"Equation5\"<\/p>\n

Where \u03b8\u00a0is the azimuth angle.\u00a0 The key features to note here is that it produces a radial propagating wave\u00a0with a 1\/R falloff and a maximum at an angle broadside to the dipole:<\/p>\n

\"Image<\/p>\n

For a periodic A, the magnitude of the E and B fields are proportional to the magnitude of A, so this plot may also be thought of as a plot of the magnitude of either E or B (though not both, as they differ by a phase constant).<\/p>\n

My next step will be to see what happens to this field in the presence of multiple dipole elements, as well as model their response to an incident field at multiple angles.<\/p>\n

\"Image

Incident EM waves at a 45 degree angle.<\/p><\/div>\n

Note: All plots in this post were generated in Matlab. \u00a0The .m file containing all the tools used to generate these plots is available in the previous post.<\/p>\n

——— 04\/21\/10 ————<\/p>\n

Here<\/a> is a preliminary MATLAB simulation I have put together of my project. \u00a0It contains a graphical display the fields of EM plane waves, as well as the fields generated by a dipole antenna, all time dependent. \u00a0I will be assembling these into the fields of a Yagi-Uda antenna over the next 2 weeks.<\/p>\n

——— 04\/11\/10 ————<\/p>\n

I have found a technique online which allows for the computation of the fields induced inside an array of perfect conductors of an arbitrary shape when subject to a rapidly oscillating EM field. \u00a0The technique makes a great deal of use of what is apparently known as the “Pocklington Integral Equation.”<\/p>\n

I’ve also found the text “Electromagnetic Waves and Antennas” (the full text of which is available for download\u00a0here<\/a>). \u00a0This text devotes two chapters to the application of the Pocklington Integral Equation for antenna of various shapes. \u00a0I’m still working through the material, so I’m not sure if this method will lend itself to easy computation, but this is still the best lead I have so far at how to solve this problem.<\/p>\n

——— 03\/29\/10 ————<\/p>\n

My wave mathematica sheet can be found here<\/a>.<\/p>\n

Based on my experience with this worksheet, I am reconsidering my decision to complete my project with Mathematica. \u00a0Mathematica seems to be great with number crunching, but clumsy and difficult to use when it comes to graphics. \u00a0And since my project is about putting together a user friendly GRAPHICAL interface which can help with VISUALIZING a complicated EM field, I am starting to think that JAVA will produce better results than Mathematica.<\/p>\n

——— 03\/19\/10 ————<\/p>\n

Journal summary: \u00a0Unfortunately, I have not found the journal to be especially useful or helpful.<\/p>\n

If I have questions, I voice them in class or ask them during office hours; if I think of something important, I write it down in my class notes; if I find something interesting, I put it on the blog. \u00a0The journal just seems to be redundant in light of these other media.<\/p>\n

To date, I only have two entries in my journal. \u00a0And they are both of the nature: “I don’t think I quite understand the concept of X. \u00a0Why is this true?” \u00a0And in both cases, the question I had was resolved in a later class.<\/p>\n

——— 03\/03\/10 ————<\/p>\n

Project: Yagi Antenna<\/h1>\n

I spent last summer at Dartmouth working with a team that launched high altitude balloons to study the upper atmosphere. \u00a0The balloons would carry the payloads to ~90,000 feet and then burst, allowing the payloads to return to Earth under parachute, generally within 100 km of the launch site.<\/p>\n

\"Balloon<\/p>\n

The payloads all contained onboard GPS receivers and radios which relayed their position back to us so we could keep track of their location. \u00a0But we found that in the last 1000 feet of the payload’s descent, we would loose radio contact, and only be able to pinpoint the payloads’ landing location to within a few kilometers. \u00a0Since the landing zone was usually in rough mountainous terrain and thick forrest, an uncertainty in the landing position of 1 km was too big to guarantee a recovery of the payload. \u00a0We thus equipped each payload with a short range Emergency Location Transponder (ELT) and used a directional Yagi Antenna to home in onto the payload’s location after it had touched down.<\/p>\n

\"Max_ELT\"<\/p>\n

\"yagi\"<\/p>\n

[Image from RadioElectronic.com.]<\/p>\n

Yagi antenna are a set of three or more parallel conductors, with their optimal dimensions and spacings determined by the wavelength of the radio signal the antenna is meant to receive. \u00a0When arranged properly, the antenna will produce a very powerful response only when oriented towards<\/em> the oncoming radio signal. \u00a0There exist tables and manuals which give the proper dimensions for Yagi antenna at various frequencies, but I have no idea exactly WHY these dimensions should work for a given signal. \u00a0As such, I intend to build a JAVA based program to simulate the response of a Yagi antenna in an oscillating EM field, and perhaps find a way to visualize exactly what the EM field is doing in the vicinity of a Yagi antenna to produce it’s directional behavior.<\/p>\n

——— 02\/16\/10 ————<\/p>\n

Mathematica analysis of refraction at a dielectric threshold can be found here<\/a>.<\/p>\n

——— 02\/09\/10 ————<\/p>\n

My Mathematica solution to Griffiths 4.10 can be found here<\/a>.<\/p>\n

——— 02\/01\/10 ————<\/p>\n

Proposed solution methods to Griffiths 4.1 part 2:<\/p>\n

Method 1:<\/p>\n

1) Assume the Hydrogen atom is composed of an electron held in “orbit” about a proton by electrical attraction, with orbital radius .5 Angstroms.<\/p>\n

2) Calculate the “orbital” energy of the electron as if it were a point mass in a classical system. \u00a0This is the ground state energy of the electron, (Eg).<\/p>\n

3) Calculate the orbital radius associated with the orbital energy Eg+13.6\u00a0eV. \u00a0This radius is the ionization radius (Ri).<\/p>\n

4) As in part 1 of the problem, calculate the applied potential necessary to distort the hydrogen atom to (Ri).<\/p>\n

Method 2: (Conceived by D. Bridgman-Packer)<\/p>\n

1) Make a reasonable assumption for the dimensions of the plates (~ 1 cm^2 for example). \u00a0Energy stored in a parallel plate capacitor (Ecap) can be calculated as:<\/p>\n

Ecap = (C*V^2)\/2<\/p>\n

Where C= (EpsNaught*A\/d)<\/p>\n

2) Assume the stored energy is evenly distributed throughout the space between the plates (VolCap), so energy density as a function of applied voltage is:<\/p>\n

Energy Density = Ecap\/VolCap = (C*V^2)\/(2*A*d) = ((EpsNaught*A\/d)*V^2)\/(2*A*d)<\/p>\n

3) Calculate the applied voltage where the Energy Density exceeds 13.6 eV per cubic angstrom.<\/p>\n

—————————————————-<\/p>\n

A moving charge simulation referring to the first worksheet.<\/p>\nhttp:\/\/www.youtube.com\/watch?v=eWaA9RiLsno<\/a>\n","protected":false},"excerpt":{"rendered":"

Final Results: The Matlab simulation of my final results. \u00a0The PowerPoint Presentation associated with the simulation. I have produced a simulation of the near and far fields of a two-element Yagi antenna. \u00a0See the previously linked .m and .ppt files for a complete discussion of methodology and modeling. Results: Figure 1: The calculated electric field […]<\/p>\n","protected":false},"author":95,"featured_media":0,"parent":9,"menu_order":0,"comment_status":"open","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-54","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/pages\/54","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/users\/95"}],"replies":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/comments?post=54"}],"version-history":[{"count":51,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/pages\/54\/revisions"}],"predecessor-version":[{"id":967,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/pages\/54\/revisions\/967"}],"up":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/pages\/9"}],"wp:attachment":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/media?parent=54"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}