Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Physics / Transmission Lines / Notes By Chapter and Appendix / Appendix H 3D props

saved old Appendix H REVIEWED

DOCX · 24.7 KB
Open DOCX file

Archived superseded appendix from Phil's transmission lines notes, dated 11.14.13 and replaced by new Appendices H and I. It proves that the Laplacian of 1/r is -4πδ(r) using the divergence theorem, extends this to h(r)/r, and gets the Helmholtz and Poisson Green's functions. A 2D version shows the Laplacian of ln(1/r) is -2πδ(r).

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
Saved old Appendix H PhL 11.14.13 I replaced this old App H with a pair of new appendices, App H and App I. They subsume all the results developed in this old App H. So this is just for archive, can be tossed eventually. Appendix H : Calculation of 2[1/r] , 2[h(r)/r] and 2[e-jkr/r], and 22D[ln(1/r)] __________________________________________________________________________________ Theorem 1: -2[1/r] = 4πδ(r) (H.1) Proof: Let volume V be all of 3D space. Carve out from V a small spherical cavity of radius a centered at r = 0. If we call this spherical volume Va and then V' = V - Va is the original volume with the spherical cavity. In order to show that some function g(r) = δ(r), one has to show that lima→0 ∫V' dV g(r) = 0 (H.2a) lima→0 ∫Va dV g(r) = 1 (H.2b) This is basically the definition of δ(r). Since δ(r) has units L-3, g(r) has units L-3. Our candidate function of interest is g(r) = - (1/4π) 2[1/r] (H.3) Using 2 in spherical coordinates acting on a function of r, one finds that, since ∂r(1) = 0, 2[1/r] = (1/r2)∂r(r2∂r) [1/r] = 0 r > 0 (H.4) so that g(r) = - (1/4π) 2[1/r] = 0 r > 0 . (H.5) Thus, condition (H.2a) is trivially satisfied since r > 0 everywhere in volume V'. It remains to verify condition (H.2b). Consider the integral appearing in the left side of (H.2b) ∫Va dV g(r) = - (1/4π) ∫Va dV 2[1/r] = - (1/4π) ∫Va dV [1/r] . (H.6) The divergence theorem says ∫V dV div F = ∫S dS F (H.7) where V is any closed volume whose surface is S, and dS points out. Using V = Va and F = [1/r] = ∂r(1/r) = -r-2 we find that LHS (H.7) = ∫Va dV div [1/r] = ∫Va dV 2[1/r] = ∫Va dV [-4πg(r)] = -4π ∫Va dV g(r) RHS (H.7) = ∫S dS [1/r] = ∫dΩ [a2 ] [1/r]|r=a = ∫dΩ[a2 ] [-a-2] = -4π which tells us that ∫Va dV g(r) = 1 for any a. Thus, lima→0 ∫Va dV g(r) = 1 and we have then verified (H.2b). Therefore we conclude that - (1/4π)2[1/r] = δ(r) (H.8) or 2[1/r] = - 4πδ(r) which is (H.1). QED __________________________________________________________________________________ Theorem 2: 2[h(r)/r] = - 4π h(0) δ(r) + h"(r)/ r (H.9) Proof: Start with this vector identity 2(φψ) = φ2ψ + ψ2φ + 2 φ ψ (H.10) so that 2(h r-1) = h2(r-1) + r-12h + 2 h (r-1) = - h 4π δ(r) + r-12h + 2 [ h' (-r-2) = - 4π h(0) δ(r) + r-12h - 2 r-2 h'(r) (H.11) Algebra shows that, 2h = (1/r2)∂r(r2∂r)h(r) = h"(r) + (2/r)h'(r) (H.12) so then 2(h r-1) = - 4π h(0) δ(r) + r-1 [h"(r) + (2/r)h'(r) ] - 2 r-2 h'(r) = - 4π h(0) δ(r) + h"(r)/ r which is the claim of (H.9). QED __________________________________________________________________________________ Example of Theorem 2: h(r) = e-jkr h = e-jkr h(0) = 1 h' = -jk e-jkr h" = -k2 e-jkr 2(e-jkr/r) = - 4π 1 δ(r) + [-k2 e-jkr ] / r = -4πδ(r) - k2(e-jkr/r) Thus, ( 2+k2) (e-jkr/r) = - 4πδ(r) (H.13) or - (2+k2) (e-jkr/4πr) = δ(r) (H.14) If we translate the origin from 0 to x' and define R = |r - x'| and rename r to be x, (H.13) becomes - (2+k2) (e-jkR/4πR) = δ(x-x') (H.15) Thus we may identify e-jkR/4πR as the Helmholtz equation Green's Function. In the case k = 0 we get - 2(1/4πR) = δ(x-x') (H.16) and we may identify (1/4πR) as the Poisson equation Green's Function as shown in (A.0.3) and (A.0.4). __________________________________________________________________________________ Theorem 3: 2[ln(1/r)] = - 2πδ(r) in two dimensions (H.1)' Our 2D proof here will be an "copy, paste and edit" version of the 3D proof of Theorem 1. We shall use the word "volume" with symbol V to denote area. A "spherical cavity" now means a circular hole. It should be understood that 2 = ∂x2+∂y2 and δ(r)= δ(x)δ(y). We indicate all edits in red and equation numbers are the same as for Theorem 1 with a prime added. Proof: Let volume V be all of 2D space. Carve out from V a small spherical cavity of radius a centered at r = 0. If we call this spherical volume Va and then V' = V - Va is the original volume with the spherical cavity. In order to show that some function g(r) = δ(r), one has to show that lima→0 ∫V' dV g(r) = 0 (H.2a)' lima→0 ∫Va dV g(r) = 1 (H.2b)' This is basically the definition of δ(r). Since δ(r) has units L-2, g(r) has units L-2. Comment: When any differential operator like or 2 is applied to ln(r0/r), the result is independent of r0 so we can always take r0 = 1. For example, ∂x [ln(r0/r)] = ∂x [ lnr0 + ln(1/r)] = ∂x ln(1/r). In what follows, ln(r) and ln(1/r) are always acted upon by differential operators, so we can interpret these objects as dimensionless quantities ln(r/r0) and ln(r0/r). Then it is clear below that dim [g(r)] = L-2. Our candidate function of interest is g(r) = - (1/2π) 2[ln(1/r)] = +(1/2π) 2 [ ln(r)] (H.3)' Using 2 in polar coordinates acting on a function of r, one finds that, since ∂r(1) = 0, 2[ln(r)] = (1/r)∂r(r∂r) [ln(r)] = 0 r > 0 (H.4)' so that g(r) = - (1/2π) 2[ln(1/r)] = 0 r > 0 . (H.5)' Thus, condition (H.2a)' is trivially satisfied since r > 0 everywhere in volume V'. It remains to verify condition (H.2b)'. Consider the integral appearing in the left side of (H.2b)' ∫Va dV g(r) = - (1/2π) ∫Va dV 2[ln(1/r)] = - (1/2π) ∫Va dV [ln(1/r)] . (H.6)' The divergence theorem says ∫V dV div F = ∫S dS F (H.7)' where V is any closed volume whose surface is S, and dS points out. Using V = Va and F = [ln(1/r)] = ∂r(ln(1/r)) = - ∂r(lnr) = [ -r-1] we find that [ here dS = dS where dS = adθ = a piece of circumference of the circle bounding Va ] LHS (H.7) = ∫Va dV div [ln(1/r)] = ∫Va dV 2[ln(1/r)] = ∫Va dV [-2πg(r)] = -2π ∫Va dV g(r) RHS (H.7) = ∫S dS [ ln(1/r)] = ∫dθ [a ] [ ln(1/r)]|r=a = ∫dθ [a ] [-a-1] = -2π which tells us that ∫Va dV g(r) = 1 for any a. Thus, lima→0 ∫Va dV g(r) = 1 and we have then verified (H.2b)'. Therefore we conclude that - (1/2π)2[ln(1/r)] = δ(r) (H.8)' or 2[ln(1/r)] = - 2πδ(r) which is (H.1)'. QED