Surface Integral: Evaluating with Spherical Coordinates

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
4 replies · 2K views
TheFerruccio
Messages
216
Reaction score
0

Homework Statement



Find
[tex]\iint\limits_S \mathbf{F}\cdot \hat n\, dA[/tex]

Homework Equations



[tex]\mathbf{F} = [1, 1, a][/tex]

[tex]S: s^2+y^2+4z^2 = 4, z \geq 0[/tex]

The Attempt at a Solution



I parameterized in spherical coordinates

[tex]x=4\sin{\phi}\cos{\theta}[/tex]
[tex]y=4\sin{\phi}\sin{\theta}[/tex]
[tex]z=\cos{\phi}[/tex]

Then, I found the surface normal vector, and finding the normal vector is what exploded into something that I couldn't simplify very well. I have a feeling that, because it exploded, that there is a simpler way for me to go about doing this. I thought about using the divergence theorem, but I didn't see how I could use it with an open surface.
 
Physics news on Phys.org
Well, what did you get for that "normal vector"? I, taking the cross product of the the derivatives of [itex]\vec{r}= 4 sin(\phi)cos(\theta)\vec{i}+ 4 sin(\phi)sin(\theta)\vec{j}+ cos(\phi)\vec{k}[/itex], got, giving the "vector differential of surface area",
[tex]4sin^2(\phi)cos(\theta)\vec{i}+ 4sin^2(\phi)sin(\theta)\vec{j}+ 16 sin(\phi)cos(\phi)\vec{k}[/tex]

With [itex]\vec{F}= \vec{i}+ \vec{j}+ a\vec{k}[/itex] then the integrand is
[tex](4 sin^2(\phi)cos(\theta)+ 4sin^2(\phi)sin(\theta)+ 16a sin(\phi)cos(\phi))d\theta d\phi[/tex]

I don't see anything terribly difficult to integrate in that- use the standard identity [itex]sin^2(\phi)= (1/2)(1- cos(2\phi)[/itex].

Or, avoid the surface integral entirely by using the divergence theorem: integrate [itex]\nabla\cdot \vec{F}[/itex] over the interior of the ellipsoid.
 
HallsofIvy said:
Well, what did you get for that "normal vector"? I, taking the cross product of the the derivatives of [itex]\vec{r}= 4 sin(\phi)cos(\theta)\vec{i}+ 4 sin(\phi)sin(\theta)\vec{j}+ cos(\phi)\vec{k}[/itex], got, giving the "vector differential of surface area",
[tex]4sin^2(\phi)cos(\theta)\vec{i}+ 4sin^2(\phi)sin(\theta)\vec{j}+ 16 sin(\phi)cos(\phi)\vec{k}[/tex]

With [itex]\vec{F}= \vec{i}+ \vec{j}+ a\vec{k}[/itex] then the integrand is
[tex](4 sin^2(\phi)cos(\theta)+ 4sin^2(\phi)sin(\theta)+ 16a sin(\phi)cos(\phi))d\theta d\phi[/tex]

I don't see anything terribly difficult to integrate in that- use the standard identity [itex]sin^2(\phi)= (1/2)(1- cos(2\phi)[/itex].

Or, avoid the surface integral entirely by using the divergence theorem: integrate [itex]\nabla\cdot \vec{F}[/itex] over the interior of the ellipsoid.

Yes, that's what I got. However, you putting "normal vector" in quotes hints at me entirely borking up the terminology. I was basing "normal vector" off of: http://en.wikipedia.org/wiki/Surface_normal The whole reason I got stuck was because I thought I had to divide, then, by the magnitude, so I can get a pure direction. Is this not the case?

My confusion over using the divergence theorem was addressed in the first post. The surface isn't closed (since it's only the top half), and I think I'd be getting a net 0 divergence over the entire volume (the vector field uniformly travels in and out of the object). It doesn't seem like it would be helpful.

Maybe my confusion is from mixing up the concept of calculating flux (which is what I'm doing here) and outright finding the surface area. Is finding the surface area itself where I should be focusing on finding the surface normal?

*edit* also, since the integrand is now in spherical coordinates, isn't the _______ (forgot the name) now [tex]\sin{\theta}d\theta d\phi[/tex]?
 
Last edited:
TheFerruccio said:
Yes, that's what I got. However, you putting "normal vector" in quotes hints at me entirely borking up the terminology. I was basing "normal vector" off of: http://en.wikipedia.org/wiki/Surface_normal The whole reason I got stuck was because I thought I had to divide, then, by the magnitude, so I can get a pure direction. Is this not the case?
Yes and no. The vector [itex]\hat{n}[/itex] in the integral indeed refers to the unit vector normal to the surface, but to actually evaluate the integral, you also need dA. It turns out, the cross product method of finding the normal automatically gives you dA as well. In other words, using the notation from the wiki article, you have

[tex]\hat{n}\,dA = \frac{\partial\mathbf{x}}{\partial s}\times\frac{\partial\mathbf{x}}{\partial t}\,ds\,dt[/tex]
*edit* also, since the integrand is now in spherical coordinates, isn't the _______ (forgot the name) now [tex]\sin{\theta}d\theta d\phi[/tex]?
In the normal vector you found, there's a common factor of [itex]\sin\phi[/itex]. That's exactly the factor you're referring to here. (Note you switched the roles of theta and phi from how you were using them in the original post.)