Why Are Bessel Functions Of The Second Kind Excluded? Ensure Regular Radial Transport In Solid Cylinders

Bessel functions of the second kind are excluded because they become singular at the cylinder’s centerline. In a solid cylindrical apparatus whose domain includes (r=0), the radial solution generally has the form (A J_0(\lambda r) + B Y_0(\lambda r)). Since (Y_0(\lambda r)) approaches negative infinity as (r) approaches zero, any nonzero (B) would predict an infinite temperature or concentration at the center, which is physically unacceptable.

For radial transport in a solid cylinder, the finite-value condition at (r=0) forces the coefficient of (Y_0) to zero, leaving (J_0(\lambda r)) as the physically admissible radial function.

Why the Radial Equation Produces Two Function Families

The cylindrical radial equation

After separating variables in a cylindrical heat-conduction or mass-diffusion problem, the radial component commonly satisfies a Bessel equation of order zero:

[ r^2 R'' + rR' + \lambda^2 r^2 R = 0. ]

Its general solution is

[ R(r) = A J_0(\lambda r) + B Y_0(\lambda r), ]

where (J_0) is the Bessel function of the first kind and (Y_0) is the Bessel function of the second kind.

The two solutions have different behavior

The function (J_0(\lambda r)) remains finite at the origin:

[ J_0(0)=1. ]

The second-kind function is singular there:

[ Y_0(\lambda r)\rightarrow -\infty \qquad \text{as } r\rightarrow 0. ]

Mathematically, both functions solve the differential equation away from the origin, but only one is compatible with a regular centerline.

Why the Centerline Condition Matters

Temperature and concentration must remain finite

In a conventional solid cylindrical rod, catalytic bed, or pipe wall, the axis (r=0) is part of the physical material. Temperature, concentration, and other transported quantities cannot become infinite at that location under an ordinary finite-source model.

If (B\neq 0), the term (B Y_0(\lambda r)) produces an unbounded field at the centerline. That violates the required physical condition:

[ |R(0)|<\infty. ]

Therefore,

[ B=0, ]

and the radial solution reduces to

[ R(r)=A J_0(\lambda r). ]

Symmetry also supports the same conclusion

A regular centerline requires zero radial flux by symmetry:

[ \left.\frac{dR}{dr}\right|_{r=0}=0. ]

The first-kind solution satisfies this regularity condition. The second-kind solution does not: its derivative is also singular as (r) approaches zero.

Thus, both finite field values and finite, symmetric flux behavior eliminate (Y_0).

What This Means for Radial Transport Models

Solid cylindrical domains use (J_0) only

For a solid cylinder extending from the centerline to an outer radius (R),

[ 0\leq r\leq R, ]

the admissible radial eigenfunctions have the form

[ J_0(\lambda_n r). ]

The outer boundary condition then determines the allowed eigenvalues (\lambda_n). For example, a specified wall temperature, convective boundary, or zero-flux condition produces a corresponding eigenvalue equation involving (J_0) and often (J_1).

The exclusion is a boundary-condition requirement

The coefficient of (Y_0) is not removed because second-kind functions are mathematically invalid. They are valid solutions for (r>0).

They are excluded because the physical domain includes a singular point, and the solution must be regular throughout that domain.

When (Y_0) Is Appropriate

Annular domains do not include the origin

The conclusion changes when the domain is annular:

[ r_{\mathrm{in}}\leq r\leq r_{\mathrm{out}}, \qquad r_{\mathrm{in}}>0. ]

Examples include double-pipe heat exchangers, concentric reactor channels, and shell-side or tube-side flow regions surrounding a central solid.

Because (r=0) lies outside the modeled region, (Y_0(\lambda r)) is finite everywhere within the domain.

Both functions may be needed

For an annulus, the radial solution generally remains

[ R(r)=A J_0(\lambda r)+B Y_0(\lambda r). ]

The two independent functions provide the flexibility needed to satisfy boundary conditions at both the inner and outer walls.

In this setting, excluding (Y_0) would usually remove a necessary degree of freedom and could prevent the model from satisfying both wall conditions.

Understanding the Trade-offs

The rule depends on geometry

The statement that (Y_0) must be excluded applies to domains containing the axis. It does not apply universally to every cylindrical transport problem.

The decisive question is whether the radial interval reaches (r=0), not whether the equipment is generally described as cylindrical.

A singularity may represent a different physical model

An infinite mathematical solution is usually unacceptable for ordinary heat or mass transport in a solid cylinder. However, singular behavior can arise in specialized models involving idealized line sources, point sources, or concentrated forcing at the axis.

Those cases require separate source terms and an appropriate interpretation of the singularity. They are not the standard regular-center problem represented by a solid cylindrical pilot apparatus.

Numerical implementations must preserve regularity

A numerical solver should impose the centerline condition explicitly or use a basis composed only of regular eigenfunctions. Retaining (Y_0) in a discretized solution can create severe mesh-dependent values near the axis and distort predicted temperatures, concentrations, or fluxes.

The analytical restriction therefore also protects numerical stability and physical interpretability.

Making the Right Choice for Your Goal

Choose the radial basis according to the actual transport domain and its boundary conditions.

If your primary focus is transport in a solid cylinder: Use only (J_0(\lambda_n r)) for the radial eigenfunctions and enforce finite value and symmetric zero-flux behavior at (r=0).
If your primary focus is transport in an annular apparatus: Retain both (J_0(\lambda r)) and (Y_0(\lambda r)), because the origin is outside the domain and both functions remain finite between the inner and outer radii.
If your primary focus is modeling a singular central source: Do not discard (Y_0) automatically, but define the source and its mathematical boundary condition explicitly before interpreting the solution physically.

The correct criterion is simple: exclude (Y_0) when the physical domain includes (r=0) and the transported quantity must remain finite there.

Summary Table:

Transport domain	Radial basis	Key condition
Solid cylinder, 0 ≤ r ≤ R	J₀(λₙr) only	Field remains finite and flux is zero at r = 0
Annular domain, rᵢ > 0	A J₀(λr) + B Y₀(λr)	Both functions are finite within the modeled region
Singular central source	Model-dependent	Define the source term and singular boundary condition explicitly
Numerical implementation	Regular basis or explicit centerline condition	Prevent mesh-dependent values and preserve physical stability

Need reliable pilot-scale radial transport equipment? Contact LABPARK to discuss Educational and Vocational Unit Operations Pilot Plants tailored to your application. LABPARK supports universities, research institutes, and enterprises with pilot plants for chemical engineering, bioprocess and biotechnology, and environmental and water treatment. Our team can help you select the right solid-cylinder or annular configuration, align the design with centerline and wall boundary conditions, and build a practical system for teaching, research, process development, and scale-up.