Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Math / Sneddon / Sned p203 Ku=mu u and the Lk=delta problem Jim et al / Lg=d done with Jim and sent to Jim

jim delta 1_18_11

DOCX · 19.6 KB
Open DOCX file

A short letter in Word format from Phil to his correspondent Jim, dated 1/18/11 by file name. It uses a regulated 1/r, the 3D Laplacian identity and the relation between radial and n-dimensional delta functions to derive the first solution. It then shows the second solution agrees by subtracting the two, with a Maple check, and ends by crediting Jim's contributions.

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
Jim, I'm sure you are pretty bored by now with the problem I proposed, but here is a simple way to verify the two solutions you came up with, where there is no need for Stakgold's Theorem and all that stuff. If you don't want to look at equations, at least jump down to item 3. 1. Use this shorthand notation: [1/r] ≡ limε→0 1/ . (0) We have this 3D known fact that - 23D,radial ([1/r]) = - r-2∂r(r2∂r[1/r]) = 4π δ(3)(r) (1) and we know that δ(3)(r) = δ(r)/[4πr2] so the above says - r-2∂r(r2∂r[1/r]) = δ(r)/r2 (2) Multiply both sides by r - r-1∂r(r2∂r[1/r]) = δ(r)/r (3) Use the fact that δ(2)(r) = δ(r)/[2πr] so the above says - r-1∂r(r2∂r[1/r]) = 2π δ(2)(r) (4) Now use "Jim's Theorem" which says this r-2∂r(r2∂rf) = r-1∂r2(rf) // = 2r-1∂rf + ∂r2f (5) and multiply it by r to get r-1∂r(r2∂rf) = ∂r2(rf) (6) and install this into (4) to get - ∂r2(r[1/r]) = 2π δ(2)(r) (7) and this then gives your first solution to my posed problem L([1/r]) = 2π δ(2)(r) 2. Your second solution is this: (1/r)∂r(r[1/r]) = 2π δ(2)(r) (8) If this is correct, then if we subtract (7) from (8) we should get zero, (1/r)∂r(r[1/r]) + ∂r2(r[1/r]) = 0 (9) We shall now show that this is in fact zero and in this way we verify your second solution. For the first step we have Maple show this fact, (1/r)∂r(rf(r)) + ∂r2(rf(r)) = (1/r) [ r2∂r2+ 3r∂r + 1]f(r) (10) If we take f(r) = (1/r), this thing is trivially true: [ r2∂r2+ 3r∂r + 1](1/r) = 2/r - 3/r + 1/r = 0 (11) Not surprisingly, if we take f(r) = [1/r], it is still true. Maple gives [ r2∂r2+ 3r∂r + 1][1/r] = -ε2 (2r2-ε2)/(r2+ε2)5/2 → ε2 2/r3 → 0 (12) Setting f(r) = [1/r] in (10) and using (12) we conclude that (1/r)∂r(r[1/r]) + ∂r2(r[1/r]) = 0 which is the fact (9) we sought to show, and your second solution is then verified. 3. I wish to give proper acknowledgement to Jim for solving this problem and, in doing so, pointing out these useful facts: (a) The fact that δ(n)(r) = δ(r)/[Sn(1)rn-1] where Sn(1) is the surface area of a sphere in n dimensions. (b) That the proper way to regulate 1/r is to say [1/r] ≡ limε→0 1/ . Other methods I tried did not work, such as [1/r] ≡ limε→0 1/(r+ε) or [1/r] ≡ limε→0 1/r1+ε . (c) The fact that ∂r(r2∂rf) = r-1∂r2(rf) which I used above. (d) for coming up with both solutions noted above. When I asked Jim about this problem, I had no solutions at all.