Digital Library of Mathematical Functions
About the Project
10 Bessel FunctionsApplications

§10.73 Physical Applications


§10.73(i) Bessel and Modified Bessel Functions

Bessel functions first appear in the investigation of a physical problem in Daniel Bernoulli’s analysis of the small oscillations of a uniform heavy flexible chain. For this problem and its further generalizations, see Korenev (2002, Chapter 4, §37) and Gray et al. (1922, Chapter I, §1, Chapter XVI, §4).

Bessel functions of the first kind, Jn(x), arise naturally in applications having cylindrical symmetry in which the physics is described either by Laplace’s equation 2V=0, or by the Helmholtz equation (2+k2)ψ=0.

Laplace’s equation governs problems in heat conduction, in the distribution of potential in an electrostatic field, and in hydrodynamics in the irrotational motion of an incompressible fluid. See Jackson (1999, Chapter 3, §§3.7, 3.8, 3.11, 3.13), Lamb (1932, Chapter V, §§100–102; Chapter VIII, §§186, 191–193; Chapter X, §§303, 304), Happel and Brenner (1973, Chapter 3, §3.3; Chapter 7, §7.3), Korenev (2002, Chapter 4, §43), and Gray et al. (1922, Chapter XI). In cylindrical coordinates r, ϕ, z, (§1.5(ii) we have

10.73.1 2V=1rr(rVr)+1r22Vϕ2+2Vz2=0,

and on separation of variables we obtain solutions of the form ±nϕ±κzJn(κr), from which a solution satisfying prescribed boundary conditions may be constructed.

The Helmholtz equation, (2+k2)ψ=0, follows from the wave equation

10.73.2 2ψ=1c22ψt2,

on assuming a time dependence of the form ±kt. This equation governs problems in acoustic and electromagnetic wave propagation. See Jackson (1999, Chapter 9, §9.6), Jones (1986, Chapters 7, 8), and Lord Rayleigh (1945, Vol. I, Chapter IX, §§200–211, 218, 219, 221a; Vol. II, Chapter XIII, §272a; Chapter XV, §302; Chapter XVIII; Chapter XIX, §350; Chapter XX, §357; Chapter XXI, §369). It is fundamental in the study of electromagnetic wave transmission. Consequently, Bessel functions Jn(x), and modified Bessel functions In(x), are central to the analysis of microwave and optical transmission in waveguides, including coaxial and fiber. See Krivoshlykov (1994, Chapter 2, §2.2.10; Chapter 5, §5.2.2), Kapany and Burke (1972, Chapters 4–6; Chapter 7, §A.1), and Slater (1942, Chapter 4, §§20, 25).

Bessel functions enter in the study of the scattering of light and other electromagnetic radiation, not only from cylindrical surfaces but also in the statistical analysis involved in scattering from rough surfaces. See Smith (1997, Chapter 3, §3.7; Chapter 6, §6.4)Beckmann and Spizzichino (1963, Chapter 4, §§4.2, 4.3; Chapter 5, §§5.2, 5.3; Chapter 6, §6.1; Chapter 7, §7.1.), Kerker (1969, Chapter 5, §5.6.4; Chapter 7, §7.5.6), and Bayvel and Jones (1981, Chapter 1, §§1.6.5, 1.6.6).

More recently, Bessel functions appear in the inverse problem in wave propagation, with applications in medicine, astronomy, and acoustic imaging. See Colton and Kress (1998, Chapter 2, §§2.4, 2.5; Chapter 3, §3.4). In the theory of plates and shells, the oscillations of a circular plate are determined by the differential equation

10.73.3 4W+λ22Wt2=0.

See Korenev (2002). On separation of variables into cylindrical coordinates, the Bessel functions Jn(x), and modified Bessel functions In(x) and Kn(x), all appear.

§10.73(ii) Spherical Bessel Functions

The functions jn(x), yn(x), hn(1)(x), and hn(2)(x) arise in the solution (again by separation of variables) of the Helmholtz equation in spherical coordinates ρ,θ,ϕ1.5(ii)):

10.73.4 (2+k2)f=1ρ2ρ(ρ2fρ)+1ρ2sinθθ(sinθfθ)+1ρ2sin2θ2fϕ2+k2f.

With the spherical harmonic Y,m(θ,ϕ) defined as in §14.30(i), the solutions are of the form f=g(kρ)Y,m(θ,ϕ) with g=j, y, h(1), or h(2), depending on the boundary conditions. Accordingly, the spherical Bessel functions appear in all problems in three dimensions with spherical symmetry involving the scattering of electromagnetic radiation. See Jackson (1999, Chapter 9, §9.6), Bayvel and Jones (1981, Chapter 1, §1.5.1), and Konopinski (1981, Chapter 9, §9.1). In quantum mechanics the spherical Bessel functions arise in the solution of the Schrödinger wave equation for a particle in a central potential. See Messiah (1961, Chapter IX, §§7–10).

§10.73(iii) Kelvin Functions

The analysis of the current distribution in circular conductors leads to the Kelvin functions berx, beix, kerx, and keix. See Relton (1965, Chapter X, §§10.2, 10.3), Bowman (1958, Chapter III, §§51–53), McLachlan (1961, Chapters VIII and IX), and Russell (1909). The McLachlan reference also includes other applications of Kelvin functions.

§10.73(iv) Bickley Functions

See Bickley (1935) and Altaç (1996).

§10.73(v) Rayleigh Function

For applications of the Rayleigh function σn(ν)10.21(xiii)) to problems of heat conduction and diffusion in liquids see Kapitsa (1951a).