How to calculate divergence of some special fields

  • Level: Graduate 
  • Thread starter Thread starter netheril96
  • Start date Start date
  • Tags Tags
    Divergence Fields
Join the discussion
Registration is free. Start your own thread to ask a follow-up.
4 replies · 2K views
netheril96
Messages
193
Reaction score
0
[tex]\[<br /> \nabla \cdot \frac{{\vec e_r }}{{r^2 }} = 4\pi \delta (\vec r)<br /> \][/tex]
This can be seen from[tex]\[<br /> \nabla \cdot \frac{{\vec e_r }}{{r^2 }} = \frac{1}{{r^2 }}\frac{\partial }{{\partial r}}(r^2 \cdot \frac{1}{{r^2 }}) = \frac{1}{{r^2 }}\frac{\partial }{{\partial r}}(1) = 0(r \ne 0)<br /> \][/tex]
And from Gauss' Theorem[tex]\[<br /> \int_V {(\nabla \cdot \frac{{\vec e_r }}{{r^2 }})dV = \oint_S {\frac{{\vec e_r }}{{r^2 }} \cdot d\vec S} } = 4\pi <br /> \][/tex]
But if I want to directly using the formula of divergence in spherical coordinates,I can only get[tex]\[<br /> \nabla \cdot \frac{{\vec e_r }}{{r^2 }} = \frac{1}{{r^2 }}\frac{\partial }{{\partial r}}(\frac{{r^2 }}{{r^2 }})<br /> \][/tex]
And integrating this over a volume cannot give me the result of 4π[tex]\[<br /> \int_V {(\nabla \cdot \frac{{\vec e_r }}{{r^2 }})dV = } \int_0^\pi {\sin \theta d\theta \int_0^{2\pi } {d\phi \int_0^R {\frac{\partial }{{\partial r}}(\frac{{r^2 }}{{r^2 }})} } } dr = 4\pi \int_0^R {\frac{\partial }{{\partial r}}(\frac{{r^2 }}{{r^2 }})} dr<br /> \][/tex]
(Here V is a sphere with radius of R)
So how can I connect it with Dirac Delta?
By the way,I post this here because this problem arises in the electrostatic field of a point charge and I found nothing about such thing in any book concerning δ(x).
 
Physics news on Phys.org
The problem is that [itex]\frac{r^2}{r^2}=\infty[/itex] at [itex]r=0[/itex].
 
gabbagabbahey said:
The problem is that [itex]\frac{r^2}{r^2}=\infty[/itex] at [itex]r=0[/itex].

So how can I get[tex]\[<br /> \int_0^R {\frac{\partial }{{\partial r}}(\frac{{r^2 }}{{r^2 }})} dr = 1<br /> \][/tex]
Without integration,you cannot conclude some function with a singularity is δ(x)
 
Other than just using Gauss' Law, I suppose an appropriate limiting procedure can be used. I'd start with your expression for [itex]\mathbf{\nabla}\cdot\left(\frac{\textbf{e}_r}{r^2}\right)[/itex] and calculate the limit of it as [itex]r\to 0[/itex]
 
As you have seen [tex]\delta({\vec r})[/tex] is not easily treated in spherical coordinates.
What is wrong with your first two lines? They constitute one of the definitions of the delta function, which is as 'direct' as you can get.