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!

Solving the Helmholtz Equation for a Point Source

  1. Feb 7, 2013 #1
    1. The problem statement, all variables and given/known data

    By integrating (2-55), over a small volume containing the origin, substituting ψ = Ce-jβr/r, and letting r approach zero, show that C = 1/4π, thus proving (2-58).


    2. Relevant equations

    (2-55): ∇2ψ + β2ψ = -δ(x)δ(y)δ(z)
    (2-58): ψ = e-jβr/(4πr)

    3. The attempt at a solution

    I am SO lost on this.

    I integrated the RHS first. I would assume I have to do this in spherical coordinates. I know from the internet that to convert the delta function between coordinate systems you need to incorporate the jacobian, but because it's at the origin (i.e., x,y,z = 0) I just assumed that:

    ∫∫∫(-δ(x)δ(y)δ(z))r2sin(θ)dr dθ d∅
    = -∫∫∫r2sin(θ)dr dθ d∅
    = -(4π/3)r03

    Where I've integrated from r = 0 to r0, θ from 0 to pi, and ∅ from 0 to 2pi. Is this even correct?

    For the next part, the LHS, I'm totally lost. I'm trying to follow the instructions chronologically, thus

    ∫∫∫(∇2 + β2)ψ r2sin(θ)dr dθ d∅

    Where I'm using the same integration limits as before. I'm really, really unsure how to even do this, but I assumed that psi is regarded as scalar (or otherwise does not have dependency on r, theta or phi). Thus, starting with the r component:

    ∫∫∫(1/r2)(∂/∂r)[r2(∂/∂r)] (ψ r2sin(θ)dr dθ d∅)
    = ∫∫∫(1/r2)(∂/∂r)[r2 (ψ 2r sin(θ)dr dθ d∅)]
    = ∫∫∫(1/r2)(∂/∂r)[2r3 (ψ sin(θ)dr dθ d∅)]
    = ∫∫∫(1/r2)[6r2 (ψ sin(θ)dr dθ d∅)]
    = ∫∫∫[6 (ψ sin(θ)dr dθ d∅)]
    = 4ψ∏r0

    For the next part (θ component):
    ∫∫∫(1/[r2 sin(θ)])(∂/∂θ)[sin(θ)(∂/∂θ)] (ψ r2sin(θ)dr dθ d∅)
    =∫∫∫(1/[r2 sin(θ)])(∂/∂θ)[sin(θ)cos(θ)] (ψ r2dr dθ d∅)
    =∫∫∫(1/[r2 sin(θ)])(∂/∂θ)[(1/2)sin(2θ)] (ψ r2dr dθ d∅)
    =∫∫∫(1/[r2 sin(θ)])[cos(2θ)] (ψ r2dr dθ d∅)
    =∫∫∫(cos(2θ)/ sin(θ)]) (ψ dr dθ d∅)

    But I can't integrate this because (cos(2θ)/ sin(θ)]) doesn't converge. I'm really lost on where to go from here and don't fully understand what I'm being asked to do.

    ANY help is appreciated. Thank you. Please note that I'm an electrical eng so our notation may be different, and I may have some difficulty with more of the hardcore physics notation. I have posted this on the intro physics board but nothing doing, so I'm not sure if maybe this is more advanced. If it isn't, my apologies.
     
  2. jcsd
  3. Feb 7, 2013 #2

    fzero

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    No, it's not correct. You've dropped the delta function in the second line ([itex]\delta(\vec{x})\neq 1[/itex] everywhere). There's no particular reason to integrate over spherical variables here on the RHS, but if you do, then you will actually need to use the Jacobian.


    [itex]\psi[/itex] is a scalar function, but it definitely depends on [itex]r[/itex]. You are told in the problem to substitute the expression for it in 2-58 into 2-55. Why don't you actually try doing that?
     
  4. Feb 7, 2013 #3
    Do the right hand side in Cartesian coordinates (you don't need to worry about limits as long as origin is inside the volume).


    That's not correct...
    The radial part of the Laplacian is [tex] \frac{1}{r^2} \frac{\partial}{\partial r} \left(r^2 \frac{\partial \psi}{\partial r}\right) [/tex] The derivatives hit [itex] \psi [/itex].

    The key here is to notice that ψ diverges in the origin, so you can't just plug it into the formulas and differentiate. The most intuitively clear way to proceed is to use the divergence theorem on the first term. You should have something like [tex] \int dV \nabla^2 \psi = \int dS \nabla \psi \cdot \hat{n}. [/tex] Take the surface to be a sphere, assume that ψ does not depend on the angles and you can evaluate the integral easily.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Solving the Helmholtz Equation for a Point Source
Loading...