# Calculate the volume integral of divergence over a sphere

1. Apr 18, 2017

### Vitani11

1. The problem statement, all variables and given/known data
For the vector field F(r) = Ar3e-ar2rˆ+Br-3θ^ calculate the volume integral of the divergence over a sphere of radius R, centered at the origin.
2. Relevant equations
Volume of sphere V= ∫∫∫dV = ∫∫∫r2sinθdrdθdφ
Force F(r) = Ar3e-ar2rˆ+Br-3θ^ where ^ denote basis (unit vectors), F is a vector and is a function of r which is also a vector.
D(V) divergence over volume
θ = zenith angle (limits 0 to π)
Φ = azimuthal angle (limits 0 to 2π)
R = radius, r = variable of integration (limits 0 to R)
3. The attempt at a solution
The divergence of this thing is A(3r2e-ar2-2ar4e-ar2)
Here is the integral I set up:

D(V) = A∫∫∫(3r2e-ar2-2ar4e-ar2)r2sinθdrdθdΦ and I got down to the following: ... = 12Aπ∫r4e-ar2dr - 8aAπ∫r6e-ar2dr. My problem is integrating this. I mean is this even possible? If this were complex I am sure I could do it... but the integral of e-ar2 alone is already difficult to integrate.e-r2 I know has the solution of √π when limits are from -∞ to ∞ but these limits are also different.

2. Apr 18, 2017

### Vitani11

Can I do this using a Fourier transform or expansion you think?

3. Apr 18, 2017

### TSny

4. Apr 18, 2017

### Vitani11

Oh... thank you.

5. Apr 18, 2017

### Vitani11

Okay now even fixing that I am having the same issue. I still need to integrate ∫r4e-ar2dr and ∫r6e-ar2dr. Divergence obtained is A(5r2e-ar2-2ar4e-ar2)+Bcosθ/r4sinθ. When doing the integral the limits are from 0 to π and therefore the term with the B coefficient goes away since the integrals limits become 0 to 0, also after integrating the Φ component to 2π (since it doesn't depend on anything in the integral) I have this: 2πA[∫∫5r4e-ar2sinθdrdθ-2πA∫∫2ar6e-ar2sinθdrdθ. Then after integrating sinθ I get 20πA∫r4e-ar2dr-8πAa∫r6e-ar2dr

6. Apr 18, 2017

### TSny

7. Apr 19, 2017

### Staff: Mentor

Are you familiar with the divergence theorem?

8. Apr 19, 2017

### TSny

Following Chestermiller: If you are allowed to use the divergence theorem then that's definitely the way to go!

9. Apr 20, 2017

### Vitani11

Okay awesome, thank you. I'll learn it

10. Apr 25, 2017

### Vitani11

So when I used the divergence theorem I got a value which was not 0, but when I integrated using grit over the volume instead of Gausses theorem I get something which is zero. I had to use the gamma function and what not for the second approach and there is more room for error when solving the problem, but this broke down simply and equaled 0 even after all the algebra which makes it more convincing than the first approach. Would it make sense for this this to be 0? What does this question even mean? I know the maths and what a divergence is... am I just summing up all the bits of outward force over the volume of the sphere? Using gausses theorem and the other method tells me that the divergence in the θ hat direction =0. Is this because of symmetry for θ in terms of the divergence throughout the spheres volume?

11. Apr 25, 2017

### TSny

I let Mathematica do the volume integral of the divergence. It gives a nonzero result which agrees with the result using the divergence theorem.

I'm not sure what "divergence in the θ hat direction" means. Can you clarify that?

12. Apr 25, 2017

### Vitani11

I meant that the divergence throughout the volume of the sphere is symmetric when integrating it throughout the zenith angle from 0 to pi i.e. it just cancels because the sum of the vectors represented by θ hat all sum to 0. Does this make sense? Also what did Mathematica tell you? I don't have software like that to check my work.

13. Apr 25, 2017

### TSny

We are integrating the divergence of the vector field F throughout the volume of the sphere. The divergence at each point is just a number. So, the integral is a sum of numbers rather than a sum of vectors.

However, F does have a $\hat{\theta}$ component. The contribution of this component of F to the divergence of F is Bcosθ/r4. As you noted in post #5, the integral of this expression over the volume of the sphere is zero. The only contribution to the volume integral of ∇ ⋅ F over the sphere come from the $\hat{r}$ component of F. This makes sense when considering the divergence theorem, because only the $\hat{r}$ component of F gives any flux through the surface of the sphere.

I let Mathematica do the integral ∫(5r2e-ar2-2ar4e-ar2)r2dr from r = 0 to r = R and it gave me an answer of R5e-aR2.

14. Apr 26, 2017

### Vitani11

Okay, thank you.