1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Calculate the volume integral of divergence over a sphere

  1. Apr 18, 2017 #1
    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. jcsd
  3. Apr 18, 2017 #2
    Can I do this using a Fourier transform or expansion you think?
     
  4. Apr 18, 2017 #3

    TSny

    User Avatar
    Homework Helper
    Gold Member

  5. Apr 18, 2017 #4
    Oh... thank you.
     
  6. Apr 18, 2017 #5
    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
     
  7. Apr 18, 2017 #6

    TSny

    User Avatar
    Homework Helper
    Gold Member

  8. Apr 19, 2017 #7
    Are you familiar with the divergence theorem?
     
  9. Apr 19, 2017 #8

    TSny

    User Avatar
    Homework Helper
    Gold Member

    Following Chestermiller: If you are allowed to use the divergence theorem then that's definitely the way to go!
     
  10. Apr 20, 2017 #9
    Okay awesome, thank you. I'll learn it
     
  11. Apr 25, 2017 #10
    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?
     
  12. Apr 25, 2017 #11

    TSny

    User Avatar
    Homework Helper
    Gold Member

    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?
     
  13. Apr 25, 2017 #12
    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.
     
  14. Apr 25, 2017 #13

    TSny

    User Avatar
    Homework Helper
    Gold Member

    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.
     
  15. Apr 26, 2017 #14
    Okay, thank you.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Calculate the volume integral of divergence over a sphere
  1. Volume integral (Replies: 10)

Loading...