# Amar G. Bose: 6.312 Lecture 15

DR. AMAR G. BOSE: Today we have only one topic, radiation impedance. And we'll use-- oh, by the way, I just remembered that driving home last Thursday, somehow once in awhile you review what you've done. And I recall that we got into a little side discussion of viscosity. And it all of a sudden flashed before me that I wrote something on the board. I hope I said what's correct, but what I wrote on board was definitely not correct.

We talked about any situation which had velocity shear, so to speak, that this is the velocity as you go up on the y-axis and the wave is going this way. The amount of change of the velocity as you go this way is what gives rise to a pressure here to move through a viscous medium, if the medium has viscosity. And the expression is something like this. This is not particularly important for the subject, but I recalled that I think I wrote a different thing down. du dy. And I think I wrote dy dx or something.

In other words, the amount of so-called velocity shear is what gives rise to the pressure. And this thing is called ultimately coefficient of viscosity, or just plain viscosity. The bigger the viscosity obviously, the bigger the pressure for a given average velocity.

Now, so radiation impedance. You know what it is to calculate the impedance into two terminals of an electrical network or a mechanical network. You know what it is to calculate the impedance looking into a tube, open end, or closed, or whatever is hanging on the other end. Radiation impedance is no mystery, it is simply looking out o something at the rest of the world.

In other words, the tube that you're able to compute the input impedance, z. If you looked at the impedance at this point here, this would be a radiation impedance looking out at the rest of the world. That's all it is.

If it's a specific acoustic impedance, the ratio of pressure, complex amplitudes of course, for impedance. Pressure to velocity at this point is the radiation impedance. It's called radiation only because this is launching the wave into the rest of space. It's called another impedance, input impedance, for this tube over here that's all.

Now, what the radiation impedance looks like in practical devices has a significant role to play on the equipment that you will use to drive the devices.

And we will see soon that the radiation impedance looking out of, for example, something like this, or something like a piston, or a loudspeaker in an infinite baffle. This will tell us a lot about the problems of how do you launch, how do you get some sound out here? It's harder to get it out at some frequencies than it is at other frequencies. And that's because looking out here, this fellow sees a different impedance at different frequencies.

Now, you can go ahead and solve the wave equations here to get this radiation impedance. One way, for example, would be to-- if this was an infinite baffle, to do what we had learned before. Think of this as an array of point sources here, and come out here and get an expression for the pressure as a function of the space variables. In other words, assume that this is driven by some u sub s. That gives you enough to calculate the pressure out of here by super position of a lot of little point sources in integral, which we developed before.

Then when you get that pressure, look at the expression as you go down very close to the surface. That gives you the pressure at each point, which won't necessarily be the same at each point on this-- let's say it's a circular piston.

From that you can add up and get the total force on the piston. And now you have enough to see what the impedance is. If you have the total force on the piston and you know its velocity, you have its mechanical impedance from which you can get all the other impedances just by area factors. So that's one way to get it.

Another way would be to solve the wave equation out here in space subject to the boundary condition that the normal velocity 0 along this infinite plane, except here and it's use of s from here to here. You get the ratio of the pressure to the velocity, or pressure to volume velocity, or force to velocity, whichever impedance you want. And that is the radiation impedance.

Now, funny things happen when you do this. This is going ahead, but we're going to see it in a couple of minutes. When you do this, you don't get something like for impedance of a plane wave going along on the x-axis or a wave travelling down a tube, you see a real impedance if the tube is infinitely long. You don't see that when you look out from practical devices.

If you did, we would have a great time in terms of antennas, in terms of loudspeakers. It would be much easier to radiate certain frequencies than it now is. And there are funny things going on near here. Meaning within a few wavelengths of whatever frequency you're driving that cause a lot of energy to be exchanged. It goes out here and it comes back. It goes out and it comes back. And just like we saw on the power lines how the power company doesn't like that because it has to generate a lot of extra current to give you the same power if you wish.

Well, you're going to see why a piston doesn't like that either. And it's much nicer if all the action here were in phase. If the velocity and the pressure were in phase, and then it would be great because all the motion would contribute to power that pressure-- a real part of p times u, and that's power that's going and it doesn't come back. So in that case, it's radiated power.

So what happens in this area nearby here is called, well, let's say the area itself is called the near field of a radiating source. You have near fields for electromagnetic antennas. In the same, very same sense you pump into a dipole or whatever you want. You pump current in here and you develop a voltage across here.

Well, voltage isn't in phase with the current at some frequencies. For any kind of an antenna. And so you have energy exchanging between-- if the voltage isn't in phase here, we know that means stored energy. And so the fields around that antenna store energy. They exchange it back to the source. And that creates a problem, which you have to get around.

In many antenna problems where you're transmitting a radio signal, it's not a big problem because it's very narrow band and you can treat it as the modulation of let's say, plus and minus 15 kilohertz on an FM station that's 100 megahertz. For all practical purposes, the 100 megahertz is only changing by plus and minus 15. So you can tune the antenna so that the impedance is real looking in here at that frequency, and then you're all set.

But when you need something broadband, you can't do it. And if you try in the audio region, 24 decades, 20 to 20 kilohertz, to move a piston there and tune, you can't tune it to be having-- looking at a real impedance out here, no matter what kind of apparatus you put outside.

So it turns out that it's really easy to look at all these effects in terms of a very simple radiator, which embodies the basics of the problems that you find in all the complicated ones, but the math is very simple. And we get the insight from that, and then you can just carry it over to the situations that I talked about here. And that is the pulsating sphere.

Basically, we have a sphere of some radius r0, and what we're going to do is look at the world the way that pulsating sphere does. Let's say he's expanding. Don't ask how. Imagine it's a balloon or something and it's elastic. Its radius is uniformly pulsating in and out with some velocity, and that generates a pressure. And you try to push the air away. The ratio of the pressure to the velocity is a specific acoustic impedance, and we just want to find that.

There's no more mathematics in this than you've already done on your homework. We know this is now radiating out into space, so there's no objects out there. So there's no p minus coming back. So the complex amplitude of the pressure p of r and omega, we can express as p plus over r. And I'm out here somewhere now for the moment. My r is out here. p plus over r e to the minus j omega t minus-- sorry. Omega r over c.

Now, I put a p plus there just to remind us that the whole p is constituted from p plus. There is no p minus. Now, we want to find the ratio between pressure and velocity. Well, we don't have velocity, but we know that the velocity, of course, satisfies the same kind of wave equation and it's going to have a complex amplitude u of r and omega. Since I don't need any time functions here, I'm looking for impedances. So we can just bring in that u of r and omega and relate it to the p by Newton's law. And Newton's law in time said partial of pressure with respect to x. Let's see, partial pressure with respect to t is equal to-- oh, sorry. Is equal to minus rho 0 partial of u with respect to t. u is the velocity. These are instantaneous quantities.

Now, this whole thing, if I let little p be cap P e to the j omega t. Little u be cap U e to the j omega t. Then, of course, this time domain expression goes over to in complex amplitudes, partial of p of r and omega with respect to r is equal to minus rho 0.

Well, when I take the derivative with respect to t, that of course in the complex, brings down j omega rho 0 and u of r and omega.

In other words, u of t is u of r and omega, e to the j omega t. And you put that in and drop out the j omega. We've done this many times that's why I did it just by putting it down almost here.

OK, so if we have this, we now have a relationship between this and the complex amplitude of the velocity. And the ratio of this to this is what we're looking for, the impedance. So all's we have to do is use this expression in here and that will bring us to the relation.

Partial of p of r and omega with respect to r. OK, I take this and I get a minus 1 over r times p-- I'm going to take a little shortcut now because I've done this, I think, in class once before. And I know you've done it in a homework problem. But if there are questions, ask.

p of r and omega. In other words. If I differentiate with respect to this r, it's a product. I get minus 1 over r squared p plus times all this. But this much of it is p of omega. So it's minus 1 over r times p of r and omega. Plus-- Or minus now the next term. I differentiate that with respect to r. I get j omega over c times the same thing, p of r and omega equals.

Minus. This minus sign should've been here from here. Minus j omega rho 0 u of r and omega. OK.

All minus signs, they'll go out. So p of r and omega, let's see if we can do this all at once.

OK, over u of r and omega.

Now, by the way, when I say p of r and omega and u of r and omega, I had my sphere here, which was radius r0, and I was out here somewhere at radius r. I'm computing the impedance looking out here at the moment until I specialize and say that r in here is equal to r0, which I will do. But I'm just looking at the impedance at any point out here at the moment.

And this is z of r and omega. I won't call it the radiation impedance yet. It's just the impedance seen by the wave looking out from here.

OK, p divided by u is-- all the minus signs are out. That is 1 over r plus j omega over c downstairs. And j omega rho 0 upstairs.

If I multiply that to get it into standard form, multiply top and bottom by c. I will get j omega plus c over r, j omega rho 0 c, which I can express since rho zero c is z0. z0 times j omega over j omega plus c over r.

Now, it's interesting to look at this expression. By the way, I can-- well, I'll do it later.

Look what happens to this expression when omega is very large compared to c over r. And by the way, what does that mean? Omega very, very large compared to c over r means omega over c very, very large compared to 1 over r. And that means that omega over c is 2 pi over lambda. 2 pi over lambda much, much greater than 1 over r. Or r much, much greater than lambda over 2 pi. Again, the lambda over 2 pi pops up everywhere.

So if you're out here much, much greater than sixth of a wavelength. Namely, a few wavelengths at least. What happens is this term dominates this. And look what happens then. For omega much, much greater-- well, either way you want to put it. For omega much, much greater than c over r, z of r and omega goes to z0.

In other words, when we get far enough away, this spherical wave that's traveling out thinks it's a plane wave. And it thinks that because what really happens if you look at it, let's say a spherical wave going out. This little part here, if I'm computing what the impedance is seen by him. It very much depends on what his neighbors are doing.

If his neighbors aren't doing anything. In other words, if I have a tube here and I'm looking at the impedance here and nothing is moving out here, then I see a totally different impedance than if I had another tube here with the same wave going down it, and it was pushing forward too. So it constrains this thing not I have divergence.

If I had a surface like this and then we're looking at the impedance, what would happen? It notices that this surface next to it is doing that. And he's radiating like that. That's compression.

And so if I had that surface, I would have an opposite sign. Not looking like a mass, but looking like a compression.

If I have a convex surface, the next surfaces around it are radiating out in that direction, and there's divergence. This fellow now sees a mass out there. He sees a compliance if there's convergence, he sees a mass if there's divergence. And if everything is going straight, he sees z0 real. So what that little element sees depends upon this.

Now, if the sphere is very small, there's a lot of divergence. Those arrows are going way out, and he sees a lot of mass.

If if were the other way around, he would see a lot of compression. If the wave was going like so. If it was straight, it's straight. So you can almost look at a surface and you can see what kind of impedance that surface is going to see.

If with respect to the wavelength, the thing looks sort of flat, it's real. If it's looking the other way, it's a compression. If it's looking from a convex surface, it's a mass.

Now, let's take a look at this at the other end. Namely, if r gets very, very small. When r gets very, very small, we'll say-- you know what it's going to be small compared to lambda over 2 pi. But I think you can work that out.

If it gets very small, this term gets very large. This one gets neglected. And look what happens. The whole impedance becomes imaginary. So for r-- well, I said omega here. So I'll say omega much, much less.

Yeah, the way I looked at this one was varying frequencies for a fixed r. So let me look at this one the same way. You can choose anything you want, anyway you wish. You can choose a fixed frequency and different r, or a fixed r and different frequencies. Since I started this way, I'll do it again this way.

For omega much, much less than c over r, which means that r was much, much less than lambda over 2 if you want to look at it that way. z, this impedance, z of r and omega, goes to z0 times j omega times r over c. j omega r over c times z0. And you could write this-- since z0 is rho 0 c, you could write this j omega rho 0 r. It's totally imaginary. This z is a ratio of pressure to velocity. And in that analogy, this thing is an inductor, which we know as a mass for acoustics. So the thing behaves like a mass.

Now, if it behaves like a mass, remember what the piston did-- I mean, the open tube did. We were able to compute the impedance for a short open tube. And we were then able to go back and compute it another way, which as if this was a block of air in here that had no compression, and it just moved like that. And we got the same answer. So a mass, when we say something behaves like a mass, it represents a block of air that's just going back and forth. And that's what happens when you're looking from a surface here. It's divergent. Namely, there's air moving in and moving out, moving in and moving out. And you've got to move that darn piston, or whatever is driving this thing, to shuffle all that air back and forth. And as we're going to see in a few minutes, that's the dominant motion that you have on the piston. Moving it back and forth.

Well, it turns out if what you move the piston back and forth with, you put some electrical energy in over here. And if this process from here to the motion is inefficient as it's going to be spectacularly so in normal transducers, then you're putting a huge amount of energy in here just to move that thing back and forth, that mass back and forth, and very little energy is getting radiated. So there are shocking figures.

I mean, for example, as we'll see when we get to the loudspeaker model, I shouldn't give this away, but it'll be as spectacular when we get to it. You have 100 watt amplifier in your hi-fi system. About 99 watts of that thing go to heating the voice coil, which is exactly what you don't want. Because if it heats, the adhesives may age, and the coil may unwind, and then you've got a mess and no sound. So you've bought this big amplifier to heat your room and the loudspeaker. And 1% of it gets out the other end. So you don't like to see things like that with mass on the end. But that's exactly what's happening. You know if the impedance looks like a mass, there is going to be air just flopping back and forth here that has nothing to do with radiation.

OK, so let's take a look at this expression a little closer. This is a specific omega. It's the omega-axis. If I put that in form of an s, let j omega equals s, then I can explore this on the real-- anywhere in the complex plane for the moment. So let me write this as an expression.

z of r and-- now I'm going to specialize for the moment to r0. Let's talk about the let r equal r0. Let's talk about the surface of this sphere, whatever r it is. Let r equal to r0. r0 and s equals z0 from right here, s over s plus c over r.

Let me make a little-- this is the s plane. Here are the poles and zeroes of this function. Quick. Quick.

AUDIENCE: [INAUDIBLE].

DR. AMAR G. BOSE: Great, OK. So the magnitude of something like this is if I'm at this point on the omega-axis here, I'd get the magnitude by the magnitude of this vector divided by the magnitude of this vector. Pole zero patterns never tell you the constant that's out front. It's normalized on that. So I get something that the s becomes-- when the s is equal to j omega, I get this vector divided by the vector, the pole vector from j omega to this point. So I can right away see that the magnitude of the impedance starts out at 0. And as I get way, way up here, the vector from the 0 to omega and the vector from the pole of omega become equal. So it gets to 1 or z0. So this thing comes up in magnitude as frequency goes along. This is magnitude of z.

I can now if you want, call this radiation impedance because we specialized to r0. So you might want to just remind yourself that that's also called radiation impedance by putting an r on there. So that marches along here to z0, the radiation impedance as we go on.

So very small impedance at the origin. Remember now, we made the assumption when we looked at the open tube that p was equal to 0 outside here. And I told you that it wasn't quite 0, but it's very small compared to z0 inside. And of course, this whole thing, when omega is small enough, is very small. This impedance here is z0 multiplied by something that's very small. And so if it's very small compared to z0, the standing waves in here, gamma is almost 1. The standing waves in here don't know the difference between it being z0-- being p equals 0 and p equals some tiny, tiny, tiny little quantity. That's the justification for what we did earlier when we said an open tube has a 0 pressure outside. Which you know, of course, couldn't be true if you had any kind of pressure in here. And you know that an open tube t standing wave at a given frequency, the pressure is 0. We made the assumption it's 0 here. Rises slowly and the velocity is the constant that goes down.

OK, now let's look-- we'll come back to this expression. But let's first look at the real and imaginary parts of the impedance as a function of omega. Impedance as a function of omega is over here.

That's amazing. OK.

Now, z radiation r and omega is equal to z0 j omega over j omega plus c over r0. OK. Let's get that in terms of its real and imaginary parts. And you know how to do that with complex numbers now.

You'd multiply top and bottom by the conjugate of downstairs. In which case you get omega squared plus c over r0 quantity squared. And upstairs, you would get conjugate. That would be minus j omega times j omega. That is minus j omega times j omega is plus omega squared for this times that. Plus j omega c over r0.

OK, now it's easy to get the real part and the imaginary part. Real part is just this. And the imaginary part is this. I could write this as a sum of this plus this.

So taking the real and the imaginary parts of an impedance like this give me the equivalent circuit of something-- two things in series. Some of impedances add. So if I wanted to make a model of this, I'd have something which was real represented by that, plus in series with it something that was imaginary modeled by that.

So let's just look at these two terms. If I did this on a linear scale, omega, the real part which start out as omega squared. And it would go to a constant when omega got greater than-- much, much greater then c over r0, it goes to 1 or z0. I'll just normalize it and say I'm just looking at-- I don't care what the z0 is. So I'll just say it goes to 1.

If you want, this is z real part-- well, real part of z of r over z0. Goes over to 1.

Now, the imaginary part starts-- that's this part over here in the red. That imaginary part starts as omega, not as omega squared.

Look, as you go to high frequencies, it goes to 0 as 1 over omega. So if you actually plot it out, you'll get something that looks like this. And this would be imaginary part of zr over z0. So at the very low frequencies, what happens is the whole impedance that you see-- because these two are in series. In other words, this model here, when I wrote the sum up, I can now write zr of r0 and omega and say that that's looking into a box here, which is real part plus imaginary part. Because the sum of impedances mean that they're in series.

So it says the real part is negligible compared to the imaginary part. As you go down in frequency for a given size ball, the world looks imaginary. And that's exactly what we saw before because as you go down in frequency relative to the wavelength that ball looks very small. It has an enormous amount of curvature. Has enormous amount of divergence. And so the world looks imaginary to it. And it's awful hard to radiate power into an imaginary world. You can't do it.

If your impedance is all imaginary, try very hard. Put all the voltage you want across it. No power.

Instantaneous, yes. It all comes back a little bit later. So no average power, which is the radiated power. Average power-- when you're talking about radiation, the average power that flows into a radiation impedance is the power that's radiated. Normally up until now, you would say the average power flowing into a network, into a box is the power dissipated. When we looked at energy, the real consideration only was for average power that it didn't come back. When we looked at power going through a line, if it didn't return at any point in that cycle, we called it dissipated because that's all we knew at that time. It was getting absorbed somewhere.

Average power, when you're looking at a radiation situation is the power you want. That's the communication signal.

OK, if you looked at these things like you'll quite often see them in textbooks, you'll see this thing on a logarithmic scale. And of course, this thing goes as omega squared. The real part goes down as omega squared. So that 60 b per octave and it goes like that. And then the imaginary part goes as omega, which is-- sorry. This is 12 db. 12 db per octave. Because it's 1 over omega squared, and then it's constant. this is 60 b per octave, and it'll go down also at 60 b per octave.

The imaginary part, where is it? Imaginary part. It goes as 1 over omega. So 6 and 6. This is the way you would see the radiation impedances listed.

Now, very interesting when I first came across it and saw what the results were, I was amazed. That is for a sphere. Very neat, very clean. No approximations. You can calculate it.

If you actually went back and took the piston in an infinite plane, what you find out-- there's a little bit of a change in the numbers scale here. Not too much for a radius similar. The two sizes similar. But what you basically find out is that it looks just like this except here after you've crossed over here there are-- this, by the way, is about c over r somewhere in here. Very close to this point for omega. You find that there's some ripples in here that die out like this. And there's corresponding ripples in the real part that die out like that.

But the shapes are very, very similar. They both go down like this together. You can take a piston in a tube, in a very long tube, and look at its radiation impedance. It'll look like this.

Take a piston in infinite baffle, it'll look like this. And of course, this sphere does look like that.

And you can see these. I want you, in this case, to look-- because it's a good summary-- at the Baranick's chapter five. He has the graphs for these. Just look at them and think about them a little bit. There is a problem in the book as I remember. The graphs for a piston in an infinite baffle and for a piston in a tube are reversed. In other words, the legends are wrong. So you just have to reverse them. But they're so similar in their appearance that you would never know this. Unless you happened to bump into a situation where you had to calculate it and you wanted to check your calculations with what was in the book. Then you'd find it. Otherwise you wouldn't know it. So it is interesting that these more complicated situations look similar.

Now let's take a look at the world of radiation impedances in terms of sort of representing the world as a box with two terminals on it, and see what we would put in there.

In other words, the world now is-- pressure, let's say is the cross variable, velocity is the through variable. And inside this box is the radiation impedance that represents the rest of the world.

To find out what is inside of that box, we only need to go back to the radiation impedance that we calculated and make believe that all of that is contained in a box. And look at the impedance and see what would be in there. So we already computed the radiation impedance. What would you like to put in the box? And why?

Looking at the top of this section right here, there's the radiation impedance. Now, you cannot answer the question I just gave you without answering or asking, at least another question. What's your analogy? You're going from one discipline to another. I'm now going over, if I didn't say it, to electrical where I have to make an analogy.

So the analogy that I want at the moment is voltage analogous to pressure. Which means of course, current analogous to velocity.

So this is the impedance of a network where impedance is defined as voltage over current. What's the network? Somebody's got to know the answer to this.

AUDIENCE: [INAUDIBLE].

DR. AMAR G. BOSE: What's that?

AUDIENCE: Conductor and series resistor.

DR. AMAR G. BOSE: Conductor and series with a resistor? Well, let's try it and see what happens. That's the nice thing about this. You can always just think of something and then check it and see if it checks out or it doesn't. Everybody agree that that checks out?

AUDIENCE: No.

DR. AMAR G. BOSE: Where doesn't it check? What made you say no?

AUDIENCE: They should be parallel, 0, parallel. Parallel resistors.

AUDIENCE: [INAUDIBLE].

DR. AMAR G. BOSE: 0 and omega 0. This thing has 0 impedance at omega 0? No way, so that can't be it. To have a 0 at impedance at negative 0, what do we have to have besides a short circuit?

AUDIENCE: [UNINTELLIGIBLE PHRASE]

DR. AMAR G. BOSE: [UNINTELLIGIBLE PHRASE]. All right, what else?

AUDIENCE: Resistor [INAUDIBLE].

DR. AMAR G. BOSE: Resistor? Hard to do anything, but put it in parallel.

Now, two elements in parallel, two different kind of elements, be they an l and a c, an l and an r, a c and an r. They have a pole at the point where the two impedances are equal and opposite. So those two impedances are equal and opposite when s is equal to minus c over r. That's the poles. The poles means the whole thing goes infinite, the denominator goes zero. So that happens at s equal minus c over r as can see from here. Here's the pole. It's how we got it in the first place.

So this thing has a pole then at that frequency. It has a 0 that's obvious at the origin. Let's see if we can figure out what the elements must be.

What's this element?

AUDIENCE: Z naught.

DR. AMAR G. BOSE: Z naught. Sure. Why? Because when you go to high enough frequencies, you look at this network here. When you go to high enough frequency, this impedance becomes millions of times larger than this one. All the current goes down here, so this fellow sees a resistance. He can't tell that he's here. He doesn't draw any more parallel. So looking at the network, you know that this network goes t this value, whatever the r is at high frequencies. And you know that we figured out before that the radiation impedance at high frequencies, s gets very, very large, goes to z0. So the resistor then, which this network goes to at high frequencies, must be z0.

OK, and so the inductor now? Well, if you have one, you know what the pole frequency is. You can all figure out what the conductor is. You can figure it out from the pole frequency. You can figure it out another way if you wish. You can go to very low frequencies and at very, very low frequencies, this thing is out the circuit because-- I mean, out of the system. He's very small compared to this. And so the whole impedance goes to j omega r0 over c. j omega r0 over c.

AUDIENCE: Z naught.

DR. AMAR G. BOSE: Times z0. z0 j omega r0 over c, or j omega if I put in for z0 rho 0 c, j omega rho 0 r0. So this fellow here then can be expressed as rho 0 r0. And that's the network that matches that radiation impedance.

Any questions up to now? Because this is one trivial almost example algebraically, but it has the heart of radiation impedance in it.

If we understand just the things we've talked about, here we will have a very good feeling and not have to be afraid of radiation impedances in any field. It's just a matter of how complex it is to calculate them. And I'll talk about calculating them and making models of them in a minute. But here's one in which we are very fortunate in the respect that A, it's simple to calculate. And B, which isn't not obvious to you yet, it has a model that's exact, with resistors and, in this case, inductors.

By the way, suppose I want a model now of this thing in which-- and we'll need that, in which velocity is the across variable and pressure is the through variable. Would you please give me the model? This is the box.

AUDIENCE: [INAUDIBLE].

DR. AMAR G. BOSE: Value?

AUDIENCE: 1 over z0.

DR. AMAR G. BOSE: Z0. Value?

AUDIENCE: [UNINTELLIGIBLE PHRASE].

DR. AMAR G. BOSE: [UNINTELLIGIBLE] Done. Now turns out, why will you need this? Well, you'll need this because the most convenient models-- for example, if you want to talk about what drives that piston over there, you have an electromagnetic transducer driving it. The natural, and we'll see were forced into it, across variable in that will be-- well, not say we're forced into it. The natural one will be one with velocity as the across variable. And you want to get the whole circuit together, acoustics and the mechanical and the electrical. And you'll need a circuit of this form to hook it on the end. There you are. Thus the use of the dual networks that we were talking about.

OK, now I just told you I was going to do something. Oh, yeah. First, because it's just so important, at least I feel it is, review what we've done here. And then let me extend this to the much more interesting and general cases.

What we did here is we had a sphere. We computed. In essence, I went though step-- I concluded the impedance out here. But what I was after is the ratio of pressure to velocity on the surface as a function of omega. It's now fixed at r0. Got the pressure to velocity ratio out here, which is the impedance seen this sphere as a function of frequency. Took that function and said, oh, now imagine that it's a network. Looked at the network. Constructed a model. Analyzed this in terms of its real and imaginary parts, so that we would understand when we need it why it's hard to drive something at low frequencies when there's no real part. Without a real part to the impedance, you don't radiate any power.

Now, by the way, this is just like the motor thing we talked about, the industries-- you can go into a big machine shop and go into where the power comes and find a whole bank of capacitors sitting there. Because they want to have a real-- they want to have a unity power factor for their build. And they would put capacitors across here.

You could tune such a device with another acoustical device, it turns out. I'll let you think about that, how you might take a piston in an infinite baffle, for example. What could you build out here that could tune this thing so that that piston saw a real impedance? This is a very interesting problem. I don't think you'll find it in textbooks or anything. It's not difficult. It's within the realm of things you have done in the class.

What structure could you build out here, such that at a given frequency-- you pick the frequency-- that this piston would see a real whirl. And that would be very interesting because I told you the efficiency will go way up then.

Now, since it's only at one frequency, I hope you're in the business of building sirens or tone generators because that's all you're going to get. But it is a very nice problem. It'll exercise your creativity. Find a structure that would do to your impedance what this machine shops do to their impedance, so that they see a real one and have to pay less.

For example, do you ever think of why a digital watch makes that unpleasant ring? It could do other things. It could even talk to you and tell you the time. You can make an extremely effective resonator, highly efficient when you want to only broadcast one frequency. The minute you want bandwidth, totally different thing in terms of power. So your watch makes just a beep. The model is an LC circuit, high Q LC circuit. It doesn't need much power to drive it and there isn't much in there.

So whenever you have a job to design something that's very narrow in frequency, you can do it with efficiency. Which, if it was broadband, you have another problem.

OK, now I keep going on side tracks and I forgot what the heck I want to do.

I want to talk first about a little peculiarity if you wish, of-- I hope I didn't erase it. Let's see what's under here.

Well, the pressure. There's a peculiarity of the thing, it doesn't lose us any generality in our thinking, but the sphere that we analyzed has a situation that you don't encounter otherwise. Namely, it has no near field, if you want to say it, in pressure. The complex amplitude as you go out in pressure just gets smaller. That's all. So in time, if you look at that, if you radiated a wave that was a triangle here, and you neglect the fact that this end of the triangle was a little further away than the other end. That's a small thing compared to the big R. That triangle just goes out and it gets smaller and smaller. So there's really no near field in pressure in a spherical wave.

If you did the analysis for any of the other situations, you would find that there was a near field in pressure also. It only has it in velocity. If you look at the expression for velocity, take the j omega-- where the heck is it? I have it upside down, but if I flip this thing upside down and flip this upside down, I would see that-- or just solve this thing for u of omega. Take this over here, bring this over here, and flip it around. You see that u or r and omega is a big function of omega. And it has this stuff in it, which when you get to large omega, it becomes real. But for small omega, it's a mess. It's real and imaginary. That's the near field.

When you're out of the near field in antennas or in acoustics, when the impedance scene looking out there is real, when you get to such big wave fronts that the thing looks real. If you just write the expression solve this, just bring this up, bring this over here, and turn it upside down, you notice that the u has the near field in it. And the near field and the u is caused by this stuff. This all becomes real when omega gets large enough. And, in general, when you look at impedances looking out, radiation impedances, you'll get terms that are higher order in this even.

In other words, you may have an mega over c plus omega over c squared, cubed, whatnot for oddball shapes or real-- not oddball, but real shapes. But not so simple mathematically. And those are all associated with a near field. And they will be in general, in both pressure and velocity. There's no complication in our thought process, but it's something that I don't want you to think just exists always in the velocity and not in the pressure, or in the e field and not in the h field.

So now what you do in general is you have to find from whatever source you have the total ratio of whatever you want the force, the velocity, on this oddball shape, or whatever it is. When you find that, you'll get a function of omega, but you won't be so lucky that you've get a nice rational function, in general. You'll get, god knows, depending on the geometry some function.

And what you do to get even the models that are shown in Baranick for the piston and the piston in a tube, piston in an infinite baffle, what you do then is you take this frequency function, which has a magnitude and a phase as you could get from it. And you make a rational function approximation to that given frequency function. And there are various ways you can do that. You have to sign an error criteria and you might take a mean square or you might take a weighted one. But you make an approximation.

And what you would normally do is, perhaps you'd start with the magnitude of this of what you found. In other words, you found now a very complicated, strange function. And you will say, well, I will allow myself, let's say five poles and maybe five 0's, or three 0's, or whatever. And with those poles and 0's, you would s minus s1, s minus s2, s minus s5 let's say. And upstairs, you might have s minus sa, s minus sd, or whatever. And what you would do then, is you would select the a through d and the 1 through 5 to give you the best, by your criterion, approximation to the actual magnitude, function, whatever you had to, over some frequency band.

And once you make this rational function approximation, then of course, you have a network. A network which all of you can synthesize. Maybe not as one big network, but you can take off this. If you want to get very simple, you can take off this part and you can synthesize that. Then you can cascade it with the current source, an op amp, or whatever with the next set of them, and the next set of them, et cetera. But you can realize that as a lumped parameter network. I mean, they can be shown that you could realize this without any active things as well. Just as r's, l's and c's.

And so, you then take this thing and that gives you some sort of a complicated network. You'll see minimally complicated ones. Here's one that you'll see in Baranick for a piston and an infinite baffle. You'll see something like this.

This is one, and I'll draw one for which velocity is across. The reason you know it's velocity right away is because I put a capacitor in there. And the mass, which all of these impedances are going to see, is an inductor in the analogy that velocity is the through variable. So if I put a capacitor in there, you know right away it's going to be the across variable is velocity and through variable is pressure, or force, whatever you want to make it. And you see a model that looks something like this, for example.

Let's see. I have to be a little bit careful. To the truth, I'm not sure what you see. You have to look it up. This might be one that you'll see.

If this were a model that I saw for some radiation impedance looking out from some device, what resistor in this thing could you tell me the value of right now?

AUDIENCE: If there is a [INAUDIBLE].

DR. AMAR G. BOSE: Parallel combination of this is z0 because this thing becomes a short circuit as you go to high frequencies. The impedance will get real as you go to high frequencies. And so the impedance of this at high frequencies is-- the impedance of this is this in parallel with that.

If I had a different network now, let's see. Let me just look at this. I'll try to make up one here.

This is an impedance now. Since there's an inductor here, this must be-- the through variable must be u, the cross variable is some kind of pressure or force. What can I say about the elements in here at high frequencies?

AUDIENCE: [UNINTELLIGIBLE].

DR. AMAR G. BOSE: Somebody else this time.

AUDIENCE: [INAUDIBLE]

DR. AMAR G. BOSE: Yeah.

If you compute, let's say you have a loudspeaker now and you've found out what the volume-- you've analyzed the thing. You've found out what the velocity of the piston of the cone is. And somebody asks you, well, how many watts acoustical power is being radiated by that thing?

If you have a model, all you have to do is the radiated power from the loudspeaker must be the power that goes into this box. Because average power flowing out of these terminals in terms of the transducer is radiated power. Average power flowing past this border in terms of the model is power dissipated in all the resistors. So at any frequency, if this was a decent model for whatever frequency you're dealing with, to find the radiated power if you knew the velocity, all you have to do is say, well, there's u amperes going through here at that given frequency. Compute the power in here and the power in here, and that's the radiated power. So it's very nice to have a model like this to think in terms of.

You can see, for example, what happens at low frequencies versus what happens at high frequencies. You can see how the power changes just by again, using your insight from electrical networks to tell you now about things like radiated power. So there's a very nice reason to look at the rest of the world like models. Like electrical models if you're electrical engineer.

Again, it's taking a discipline that you got way back in network theory and showing that you can use it for problems that you wouldn't have imagined that it would be applicable. Questions?

I'm just trying to think whether I can launch on to another thing or not. I basically said, I believe, all that I want to say about radiation.

Yeah, let me tell-- I don't know whether I'm going to even bother to do directivity factor and calculations. I don't know, I'll think about it before next time. But I'm calculating something that you're never going to need and the if you ever did need, if you understand what it's all about, you could get it like that. So I'll just I'll just discuss now what directivity factors, or directivity-- there's those two words for it. Directivity factor and there's directivity-- one of them's a logarithm and the other, I don't know what they call them. But it's a measure of directivity.

Now, when you go, as I think I mentioned months before, you go into one of the stores and you buy a TV antenna. They may list it in terms of the directivity index or direcr-- or maybe that's the word. Or directivity factor. And the higher this is, the better the antenna. Providing you're willing to put it on a turning-- on a pole with a motor. Because the more directive the antenna is, what it means is it's very sensitive. A flashlight has a height directivity factor compared to an electric light bulb. Electric light bulb goes equally all directions. So if you put 10 watts of light in there, you will not see that electric light bulb. In the night time, outdoors it won't be a fraction as bright as a flashlight with 10 watts of light directed exactly to you. So the higher the directivity factor, it means that you're going to have a strong on axis signal, and you're going to have an on axis if you wish, that's very directional.

Now, in acoustics what you do is you always use as your-- this is usually expressed as a logarithm directivity factor. So it's relative to something. In acoustics, it's relative to a sphere.

So if you have that as your reference, and you have a loudspeaker, let's say in a box, or a microphone. Thank god we have reciprocity applying both in electromagnetics and in acoustics. And so we can speak of reception or transmission with the same terminology, the directivity factors.

This fellow here, as he gets to be higher and higher frequency, and higher frequency means higher frequency such that the wavelength gets to be less than the circumference of the diaphragm. It gets to be pretty directional.

Now, directivity factor or directivity index, whatever, is a measure of the-- well, it doesn't have an axis exactly. The measure of the intensity of this fellow on axis, the biggest, strongest point to the intensity measured at the same point if you wish. These are all far field. That a sphere would give you, if it were radiating the same power as this. That's the definition. In other words, you put w watts into this thing, you put w watts into this thing. And you go out here in the far field and you say, well, what's the intensity of this fellow right in front of him? And then I measure that. I measure the intensity of the sphere, and I take the ratio of them. So you can see that the higher the directivity factor, the stronger this thing's going to be on axis and the more it's going to be beamed.

In fact, you can think of most acoustical sources, not a sphere but most-- I've never seen a pulsating sphere. But most practical sources, you can think of them as-- for example, all the loudspeakers that you see. You can think of them as a flashlight at high frequencies and as a electric light bulb at low frequencies. When the whole cabinet and everything is small compared to sixth of a wavelength, then the thing is omnidirectional. Which means below a couple hundred Hertz, most loudspeakers. no matter whether they point backwards, forwards, or sideways are omnidirectional sources. Contrary to your intuition that, oh, my god, it's going the wrong way or something. It's omnidirectional.

And at high frequencies, it's quite direction, like a flashlight. And so directivity factors this way.

Now in electromagnetics, your reference for your antennas when you go to buy them is the dipole. Here it's a monopole in acoustics. And the dipole radiation out here on axis in the far field with w watts into the antenna, and then w watts measured at the same point, w watts into one of these things with a lot of reflectors, et cetera, and directors all the way. The kind of huge arrays you see in some people's houses in the suburbs. The ratio on axis of this thing, either sensitivity in reception, or transmission, doesn't matter, to putting out with w watts into it. And think of it as a transmitter. The ratio of the intensity on the axis to the intensity of a dipole is what you call directivity factor.

And the only key to this is how you measure these things. And to measure them, you have to know that this fellow is radiated. Putting w watts into the amplifier that drive this doesn't mean anything. As we just talked about, you're getting maybe 1% of that out.

So if you were going to actually measure this, what'd you have to do is you'd have to measure the power going in each direction of this, and integrate that. You're out here at r or whatever it is, at a constant r. You would go around and you'd measure the power over that whole sphere. And that tells you then what the total power radiated is. That's what you compare to this, which you can calculate. Set those two equal, and then find the intensity on axis versus the intensity of the sphere at the same point. And those are the calculat-- or measurements, actually. In the case of complicated things, you do this by measurement. So I don't know whether I will even bother to calculate it for a simple case.

Can you ask any questions based on what I've said?

AUDIENCE: Can I figure out a different way to answer that question.

DR. AMAR G. BOSE: I can't tell if I'm doing a lousy job or a good job. So I'll think about whether we want to calculate it. If you understand that, I think it may be enough.