1540470959-Boundary_Value_Problems_and_Partial_Differential_Equations__Powers

(jair2018) #1

5.8 Some Applications of Bessel Functions 331


Thea’s andb’s are determined from Eqs. (4) and (5) by using the orthogonality
relation ∫
a
0


J 0 (λnr)J 0 (λmr)rdr= 0 , n=m.

B. Spherical Waves


In spherical(ρ,θ,φ)coordinates (see Section 5.9), the Laplacian operator∇^2
becomes


∇^2 u=ρ^12 ∂ρ∂

(

ρ^2 ∂∂ρu

)

+ρ (^2) sin^1 (φ)∂φ∂


(

sin(φ)∂φ∂u

)

+^1

ρ^2 sin^2 (φ)

∂^2 u
∂θ^2.

Consider a wave problem in a sphere when the initial conditions depend only
on the radial coordinateρ:


1
ρ^2


∂ρ

(

ρ^2 ∂u
∂ρ

)

=^1

c^2

∂^2 u
∂t^2

, 0 <ρ<a, 0 <t,( 13 )

u(a,t)= 0 , 0 <t,( 14 )
u(ρ, 0 )=f(ρ), 0 <ρ<a,( 15 )
∂u
∂t

(ρ, 0 )=g(ρ), 0 <ρ<a.( 16 )

Assumingu(ρ,t)=R(ρ)T(t), we separate variables and find


T′′+λ^2 c^2 T= 0 , (17)
(
ρ^2 R′

)′

+λ^2 ρ^2 R= 0 , 0 <ρ<a, (18)
R(a)= 0 , (19)
∣∣
R( 0 )

∣∣

bounded. (20)

Again, the condition (20) has been added becauseρ=0 is a singular point.
Equation (18) may be put into the form


R′′+ρ^2 R′+λ^2 R= 0 ,

and comparison with Eq. (1) shows thatα=− 1 /2,γ=1, andρ= 1 /2; thus
the general solution of Eq. (18) is


R(ρ)=ρ−^1 /^2

[

AJ 1 / 2 (λρ)+BY 1 / 2 (λρ)

]

.

We know that nearρ=0,


J 1 / 2 (λρ)∼const×ρ^1 /^2 ,
Y 1 / 2 (λρ)∼const×ρ−^1 /^2.
Free download pdf