Register to reply

Find the power radiated using the Poyting vector

by blueyellow
Tags: power, poyting, radiated, vector
Share this thread:
blueyellow
#1
Dec21-11, 08:53 AM
P: 241
1. The problem statement, all variables and given/known data

As a model for electric quadrupole radiation, consider two oppositely oriented oscillating electric dipoles, +p[itex]_{0}[/itex] cos [itex]\omega[/itex] t at z=+d/2, and -p[itex]_{0}[/itex] cos [itex]\omega[/itex] t at z=-d/2, hence separated by a distance d along the z-axis. In the radiation zone, the scalar and vector potentials are given by:

[itex]\phi[/itex]=-([itex]\mu[/itex](subscript0) p(subscript 0) omega^2 d/4pi)((cos^2 theta)/r) cos [[itex]\omega[/itex] (t-r/c)]

A=-(mu(subscript 0) p(subscript 0) omega^2 d)/(4pi cr) cos [itex]\theta[/itex] cos [[itex]\omega[/itex] (t-r/c)] z-hat

a) Find the electric and magnetic fields, E and B.
b) Find the Poynting vector, N, and the power radiated, P= loop integral (subscript S) N.n-hat da

3. The attempt at a solution

a) I have done this but I am not writing all of it out.
b) N= E cross H
= [([itex]\mu[/itex](subscript 0) p(subscript0)^2 [itex]\omega[/itex]^4 d^2)/([itex]\mu[/itex] *16 pi^2 c r^2)] cos^3 [itex]\theta[/itex] [[([itex]\omega[/itex]^2)/(c^2)] sin^2 ([itex]\omega[/itex](t-(r/c))) - ([itex]\omega[/itex]/c) (sin [itex]\theta[/itex])/([itex]\mu[/itex][itex]\mu_{0}[/itex] r) cos ([itex]\omega[/itex](t-(r/c)))
-2sin[itex]\theta[/itex] cos([itex]\omega[/itex](t-(r/c)))((p[itex]_{0}[/itex]/[itex]\mu[/itex]) [([itex]\omega[/itex]^2 d)/(4pi c r)] cos[itex]\theta[/itex](([itex]\omega[/itex]/c) sin([itex]\omega[/itex](t-(r/c)) - (sin[itex]\theta[/itex]/[itex]\mu[/itex][itex]\mu_{0}[/itex] r) cos([itex]\omega[/itex](t-(r/c)))

But how does one work out the power radiated from this? My books and notes don't provide a simple answer. I am having trouble visualising the electric quadrupole and don't know what the unit vector is. And it has to be integrated over the area. The area of what? How can there be an area when it is two electric dipoles? Please help.
Phys.Org News Partner Science news on Phys.org
New type of solar concentrator desn't block the view
Researchers demonstrate ultra low-field nuclear magnetic resonance using Earth's magnetic field
Asian inventions dominate energy storage systems
vela
#2
Dec21-11, 02:15 PM
Emeritus
Sci Advisor
HW Helper
Thanks
PF Gold
P: 11,722
You want to integrate over a closed surface, e.g. a sphere, to find the total power being radiated away. This is just like when you used the integral form of Gauss's law except instead of integrating the electric field, you're integrating the Poynting vector.
blueyellow
#3
Dec24-11, 06:38 AM
P: 241
I found an easier way to do this. From my notes it says:

P=[itex]\oint_{S}[/itex] N.n-hat da
=(1/μ[itex]_{0}[/itex]) [itex]\int^{2pi}_{0}[/itex][itex]\int^{pi}_{0}[/itex]B[itex]_{\phi}[/itex]E[itex]_{\theta}[/itex] r[itex]^{2}[/itex] sin[itex]\theta[/itex] d[itex]\theta[/itex] d[itex]\phi[/itex]
=(1/μ[itex]_{0}[/itex]) [itex]\int^{2pi}_{0}[/itex][itex]\int^{pi}_{0}[/itex] [itex]\frac{\mu_{0}I_{0}\delta l}{4\pi}[/itex] sin[itex]\theta[/itex] ([itex]\frac{-i\omega}{rc}[/itex]+[itex]\frac{1}{r^{2}}[/itex]exp(i(kr - omega*t)) [itex]\frac{\mu_{0}I_{0}\delta l}{4\pi}[/itex] sin[itex]\theta[/itex] sin [itex]\theta[/itex] (-i*omega/r) exp(i(kr - omega*t)) sin [itex]\theta[/itex] d[itex]\theta[/itex] d[itex]\phi[/itex]

One problem: what do I do with the [itex]\delta[/itex]l's? How are they to be gotten rid of? Do I just ignore them?

Thaakisfox
#4
Dec24-11, 08:23 AM
P: 263
Find the power radiated using the Poyting vector

There is no dl at all. You didnt plug the proper fields in the calculation of the poynting vector. You should be calculating electric quadrupole radiation not magnetic dipole radiation. The fields you wrote for the poynting vector are the fields of a magnetic dipole.

Write the fields which u calculated for the electric quadrupole..
blueyellow
#5
Dec24-11, 09:09 AM
P: 241
"Write the fields which u calculated for the electric quadrupole.."

I don't know how to do that. How do I go from a dipole field to a quadrupole field then? Multiply everything by 2? That doesn't seem right. Yes, I will look in my notes and books and see if they say anything, but I am asking just in case this turns out to be one of those complicated questions that require 'thinking' rather than bookwork.
blueyellow
#6
Dec24-11, 09:11 AM
P: 241
Never mind, I just realised I copied down the equations for B and E from the question above the question I was supposed to be doing.
Thaakisfox
#7
Dec24-11, 09:16 AM
P: 263
no, the quadrupole is the next term in the multipole expansion of the fields. After the dipole comes the quadrupole in the expansion.
blueyellow
#8
Dec24-11, 10:53 AM
P: 241
So

P=[itex]\oint[/itex]N.n-hat da
=[itex]\int[/itex] div N d(tau)
=[itex]\int[/itex] [itex]\frac{\partial N_{z}}{\partial z}[/itex] d(tau)

but N has no z component so the integral goes to zero.

That the power seems zero doesn't seem right.
blueyellow
#9
Dec24-11, 10:57 AM
P: 241
I sort of ignored the n-hat because I didn't see any other way to go about doing the question. What would be the direction of n-hat anyway? All we know about the quadrupole by considering two dipoles is:

+p[itex]_{0}[/itex] cos ωt at z=d/2
-p[itex]_{0}[/itex] cos ωt at z=-d/2

so n-hat could be pointing in either the r-hat or theta-hat direction. Does the n-hat need to be taken into account to do the question properly?
blueyellow
#10
Dec24-11, 11:13 AM
P: 241
Previously I made a mistake while posting.

N= E cross H
= [([itex]\mu[/itex](subscript 0) p(subscript0)^2 [itex]\omega[/itex]^4 d^2)/([itex]\mu[/itex] *16 pi^2 c r^2)] cos^3 [itex]\theta[/itex] [[([itex]\omega[/itex]^2)/(c^2)] sin^2 ([itex]\omega[/itex](t-(r/c))) - ([itex]\omega[/itex]/c) (sin [itex]\theta[/itex])/([itex]\mu[/itex][itex]\mu_{0}[/itex] r) cos ([itex]\omega[/itex](t-(r/c)))
-2sin[itex]\theta[/itex] cos([itex]\omega[/itex](t-(r/c)))((p[itex]_{0}[/itex]/[itex]\mu[/itex]) [([itex]\omega[/itex]^2 d)/(4pi c r)] cos[itex]\theta[/itex](([itex]\omega[/itex]/c) sin([itex]\omega[/itex](t-(r/c)) - (sin[itex]\theta[/itex]/[itex]\mu[/itex][itex]\mu_{0}[/itex] r) cos([itex]\omega[/itex](t-(r/c)))] z-hat

I had missed the z-hat out at the end.
blueyellow
#11
Dec28-11, 03:30 PM
P: 241
I changed it from cylindrical to spherical coordinates. I found that in spherical coordinates it only has a rho-component. So only [itex]\frac{1}{\rho^{2}}[/itex] [itex]\frac{\partial}{\partial\rho}[/itex] ([itex]\rho^{2}[/itex] N[itex]_{\rho}[/itex]) matters when trying to work out div N.

But N[itex]_{\rho}[/itex] has no [itex]\rho[/itex] terms in it. So div N=0. This is the same as the result I got when I did it in cylindrical coordinates.

I spoke to someone else who is doing the same homework and he also doesn't think it makes sense for the power radiated to be zero. Is the power radiated really zero? If so, why? If not, is there some sort of trick that will help one to answer the question? Thanks if anyone replies.
vela
#12
Dec29-11, 03:13 AM
Emeritus
Sci Advisor
HW Helper
Thanks
PF Gold
P: 11,722
What did you get for E and B?
blueyellow
#13
Dec29-11, 06:19 AM
P: 241
E=-grad[itex]\phi[/itex]
=[itex]\frac{\partial\phi}{\partial r}[/itex]r-hat+[itex]\frac{1}{r}[/itex][itex]\frac{\partial\phi}{\partial\theta}[/itex]θ-hat +[itex]\frac{\partial\phi}{\partial z}[/itex]z-hat
=[[itex]\frac{-\mu_{0}p_{0}\omega^{2}d}{4\pi}[/itex][itex]\frac{cos^{2}θ}{r}[/itex]([itex]\frac{\omega}{c}[/itex]sin([itex]\omega[/itex](t-[itex]\frac{r}{c}[/itex]))+[itex]\frac{-\mu_{0}p_{0}\omega^{2}d}{4\pi}[/itex](-2sinθ [itex]\frac{cosθ}{r}[/itex])cos([itex]\omega[/itex](t-[itex]\frac{r}{c}[/itex]))]r-hat+[[itex]\frac{\mu_{0}p_{0}\omega^{2}d}{4r^{2}\pi}[/itex]sinθ cosθ cos(ω(t-[itex]\frac{r}{c}[/itex]))]θ-hat
=[itex]\frac{\mu_{0}p_{0}\omega^{2}d cosθ}{4\pi r}(-cosθ \frac{\omega}{c}sin(\omega[/itex](t-[itex]\frac{r}{c}[/itex])))r-hat+([itex]\frac{\mu_{0}p_{0}\omega^{2}d}{4\pi r^{2}}[/itex]sinθcosθcos([itex]\omega(t-\frac{r}{c}[/itex])))θ-hat

B=[itex]\nabla[/itex][itex]\times[/itex]A

=[itex]\frac{1}{r}[/itex][itex]\frac{\partial A_{z}}{\partialθ}[/itex]r-hat -[itex]\frac{\partial A_{z}}{\partial r}[/itex]theta-hat
=[itex]\frac{\mu_{0}p_{0}\omega^{2}d sinθ}{4\pi cr^{2}}[/itex]sinθcos(ω(t-r/c)))r-hat+([itex]\frac{\mu_{0}p_{0}\omega^{2}d sinθ}{4\pi cr}[/itex]cosθ([itex]\frac{ω}{c}[/itex]sin(ω(t-r/c))-[itex]\frac{sinθ}{r}[/itex]cos(ω(t-r/c)))θ-hat

Have I gone wrong somewhere?
blueyellow
#14
Dec29-11, 12:22 PM
P: 241
Sorry, I think I know where I have gone wrong now. E is not -grad phi in this case.
vela
#15
Dec29-11, 12:30 PM
Emeritus
Sci Advisor
HW Helper
Thanks
PF Gold
P: 11,722
What were your results in spherical coordinates? You probably don't want to use cylindrical coordinates here.

Even if you do, you have to remember to rewrite the potentials in terms of the cylindrical coordinate ##(r, \theta, z)## first because ##r## and ##\theta## aren't defined the same way in cylindrical and spherical coordinates.

Since you're looking for the fields in the radiation zone, where r is large, you can drop contributions proportional to 1/r2. The 1/r terms will dominate.
blueyellow
#16
Dec29-11, 12:54 PM
P: 241
So I have to find E and B in spherical coordinates? I can't just find them in cylindrical coordinates, and then convert them to spherical when I find N? Because as I said, when I convert the result I got for N into spherical coordinates, and then I do the integral, it just goes to zero.

I have realised more now that the question is about a quadrupole, but it is asking you to consider two dipoles. Does that mean I have to multiply everything I did by two?

What is a 'radiation zone'?
Thanks if anyone explains.
vela
#17
Dec29-11, 01:04 PM
Emeritus
Sci Advisor
HW Helper
Thanks
PF Gold
P: 11,722
Despite the presence of ##\hat{z}##, the expressions for the potentials are, I'm pretty sure, in spherical coordinates. (You should check that.) So using spherical coordinates is the path of least resistance. Also, the end result won't and, more important, can't matter on which coordinate system you choose.

What textbook are you using? For terms like radiation zone, you should look those up in your text rather than asking here. You want to learn how to find that kind of basic information on your own.
blueyellow
#18
Dec29-11, 01:42 PM
P: 241
I am using Griffiths. It does not explain what radiation zone is. It just highlights it in bold. It says 'we are interested in fields that survive at large distances from the source, in the so called radiation zone'.

'large distances'?
That means the radiation zone is everywhere? If the radiation zone is everywhere, and there is no difference between one place and another with regards to whether it is in the radiation zone or not, then why is it mentioned in the question?


Register to reply

Related Discussions
Total Power Radiated Introductory Physics Homework 2
Power radiated by the Sun Introductory Physics Homework 11
Power radiated by the sun Introductory Physics Homework 0
Radiated power, what does this even mean? Introductory Physics Homework 9
Net Power vs. Net Radiated Introductory Physics Homework 1