You are modeling heat conduction in a solid cylinder. After separating variables in cylindrical coordinates, you obtain a radial ODE whose general solution includes both J₀(r) and Y₀(r). The cylinder extends from r = 0 to r = a. Which solution do you discard, and what is the decisive reason?
AJ₀, because it oscillates and cannot represent a steady-state temperature distribution.
BY₀, because it diverges logarithmically as r → 0, and the physical boundary condition requires a finite temperature at the central axis.
CY₀, because it violates a mathematical theorem about orthogonality of Bessel functions.
DNeither — both solutions are needed for the general solution on the full disk.
Y₀ diverges logarithmically as x → 0, making it physically inadmissible for problems on a domain that includes the origin. A solid cylinder includes the central axis (r = 0), where the temperature must be finite. This physical boundary condition — not a mathematical preference — eliminates Y₀ from the solution. If the domain were an annular region (a ring excluding the origin), Y₀ would contribute because the singularity at r = 0 would lie outside the domain. The geometry and boundary conditions together determine which solutions survive.
Question 2 Multiple Choice
Why does Bessel's equation arise naturally when modeling physical problems with cylindrical symmetry, such as a vibrating circular drumhead?
ABessel's equation is the general form of all second-order linear ODEs and applies universally.
BCylindrical problems are described in cylindrical coordinates; separating variables in those coordinates produces Bessel's equation as the radial component, with the Laplacian's form in r generating the characteristic x²y'' + xy' terms.
CThe boundary conditions at r = a are always homogeneous Dirichlet, which forces solutions to satisfy Bessel's equation.
DBessel functions are a generalization of trigonometric functions, so they arise whenever waves are involved.
The cylindrical Laplacian ∇² contains terms like (1/r)(d/dr)(r dy/dr) = y'' + y'/r, which after multiplying through by r² produces exactly the x²y'' + xy' structure of Bessel's equation. Separating variables in cylindrical coordinates makes this appear naturally for the radial component — it is a consequence of the coordinate geometry, not a choice. In Cartesian coordinates, the Laplacian produces sine and cosine equations; in cylindrical coordinates, it produces Bessel equations.
Question 3 True / False
The Bessel function Y_ν is excluded from solutions on a full disk because it violates a mathematical theorem, independent of any physical interpretation of the problem.
TTrue
FFalse
Answer: False
Y_ν is excluded for physical reasons, not because of a mathematical theorem. Y_ν is a perfectly valid mathematical function and a legitimate solution to Bessel's equation — it is linearly independent from J_ν. The reason it is discarded in full-disk problems is that it diverges at the origin, and the physical boundary condition (finite temperature, pressure, displacement at the central axis) makes that divergence inadmissible. In annular domains that exclude the origin, Y_ν is retained. The choice is always driven by physics and geometry, not mathematics alone.
Question 4 True / False
Bessel functions of the first kind J_ν(x) exhibit decaying oscillatory behavior for large x, analogous to how circular waves decrease in amplitude as they spread outward from a point source.
TTrue
FFalse
Answer: True
For large x, J_ν(x) ≈ √(2/πx) cos(x − νπ/2 − π/4) — a damped cosine whose amplitude decreases as 1/√x. This behavior makes physical sense: in cylindrical geometry, a wave spreading outward from the axis must cover an ever-larger circumference (growing as 2πr), so its amplitude must decay to conserve energy. The analogy to ripples spreading across a pond is exact: they oscillate regularly but decrease in height as the ring expands.
Question 5 Short Answer
Explain why the orthogonality of Bessel functions uses a weighted inner product with a factor of x (or r in physical problems), rather than the standard unweighted inner product used for ordinary Fourier series.
Think about your answer, then reveal below.
Model answer: The weight factor x comes directly from the cylindrical coordinate area element. In Cartesian geometry, equal strips of width dx have equal area, so the standard unweighted integral ∫f(x)g(x)dx reflects equal weighting. In cylindrical geometry, an annular strip at radius r has area proportional to r·dr — a larger ring at greater radius covers more area than a thin ring near the origin. The weighted inner product ∫₀ᵃ x J_ν(λ_m x) J_ν(λ_n x) dx = 0 reflects this geometry: functions are orthogonal in the inner product that respects the cylindrical area element.
This is not an arbitrary mathematical convention — it is the natural inner product for the function space on a disk. The Sturm-Liouville form of Bessel's equation identifies x as the weight function for exactly this reason. When computing Fourier-Bessel series coefficients, you use this weighted orthogonality just as you use the standard Fourier inner product ∫f(x)sin(nπx/L)dx to compute Fourier sine coefficients. Both are orthogonal expansions adapted to the geometry of their domain.