The mathematical basis functions are spherical harmonics in angle and spherical Bessel functions in radius, selected to satisfy the equipment’s boundary conditions. For axisymmetric spherical modules, the angular dependence is expanded in Legendre polynomials (P_\ell(\cos\theta)), while the radial dependence uses spherical Bessel functions or their equivalent sine-over-radius form. The allowable eigenvalues are determined by whether the sphere has a fixed surface value, zero flux, or convective heat or mass transfer.

The governing principle is to choose an orthogonal eigenfunction set that satisfies both the spherical differential operator and the physical boundary conditions. Dirichlet, Neumann, and Robin conditions produce different characteristic equations, eigenvalues, and expansion coefficients.

How the Spherical Transport Problem Is Formulated

Heat and mass diffusion share the same mathematical structure

For constant properties and no internal generation, temperature or concentration can be represented by a generic field (\Theta(r,\theta,\phi,t)). The transient equation has the diffusion form

[ \frac{\partial \Theta}{\partial t}

D\nabla^2\Theta, ]

where (D) is the relevant diffusivity.

For heat conduction,

[ D=\alpha=\frac{k_T}{\rho \hat C_p}, ]

while for mass diffusion, (D) is the species diffusivity.

Spherical coordinates separate naturally

The Laplacian in spherical coordinates separates into radial and angular parts. A separated solution can be written schematically as

[ \Theta(r,\theta,\phi,t)

R(r)Y(\theta,\phi)G(t). ]

The time-dependent factor generally has the form

[ G(t)=\exp(-D\lambda^2t), ]

so higher spatial modes decay more rapidly.

Axisymmetry reduces the angular basis

If the equipment and loading are independent of the azimuthal angle (\phi), the angular problem depends only on

[ \eta=\cos\theta. ]

The angular equation then becomes Legendre’s equation, whose regular solutions are the Legendre polynomials

[ P_\ell(\eta)=P_\ell(\cos\theta). ]

These polynomials are orthogonal over (-1\leq\eta\leq1), allowing the field to be expanded as

[ \Theta(r,\theta,t)

\sum_{\ell=0}^{\infty} \sum_{n} A_{\ell n} R_{\ell n}(r) P_\ell(\cos\theta) e^{-D\lambda_{\ell n}^2t}. ]

For fully three-dimensional loading, the Legendre polynomials are replaced by spherical harmonics (Y_\ell^m(\theta,\phi)).

Which Radial Basis Functions Are Used?

Spherical Bessel functions are the general radial basis

The radial eigenvalue equation contains the angular separation index (\ell). Its regular solution at the sphere’s center is the spherical Bessel function

[ R_\ell(r)=j_\ell(\lambda r). ]

The second independent solution, the spherical Neumann function (y_\ell(\lambda r)), is singular at (r=0) and is therefore excluded for a solid sphere containing the origin.

For the spherically symmetric mode (\ell=0),

[ j_0(\lambda r)=\frac{\sin(\lambda r)}{\lambda r}. ]

Thus, the commonly used form

[ R_0(r)\propto\frac{\sin(\lambda r)}{r} ]

is specifically the (\ell=0) radial solution, or the result of the transformation

[ \Theta(r,t)=\frac{\Psi(r,t)}{r}. ]

The (\Psi/r) transformation simplifies the radial equation

For a spherically symmetric problem, the radial diffusion equation is

[ \frac{\partial \Theta}{\partial t}

D\left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2\frac{\partial\Theta}{\partial r} \right) \right]. ]

Using

[ \Theta(r,t)=\frac{\Psi(r,t)}{r} ]

converts the spherical radial operator into the one-dimensional diffusion equation

[ \frac{\partial\Psi}{\partial t}

D\frac{\partial^2\Psi}{\partial r^2}. ]

The resulting spatial eigenfunctions are combinations of (\sin(\lambda r)) and (\cos(\lambda r)). Regularity at the center removes the cosine contribution in the usual solid-sphere formulation.

Orthogonality enables series solutions

The radial and angular eigenfunctions form orthogonal sets under the appropriate spherical weight. For radial modes, the volume element introduces an (r^2) weighting:

[ \int_0^l r^2 R_{\ell n}(r)R_{\ell m}(r),dr=0, \qquad n\ne m. ]

This orthogonality permits the initial temperature or concentration distribution to be projected onto the eigenfunctions, producing the coefficients (A_{\ell n}).

How Boundary Conditions Determine the Eigenvalues

The center requires a regularity or symmetry condition

At the center of a solid sphere, the field must remain finite. The physical radial flux must also vanish:

[ \left.\frac{\partial\Theta}{\partial r}\right|_{r=0}=0. ]

This is the spherical symmetry condition. It is not an optional numerical convenience; it enforces continuity through the center and eliminates singular radial solutions.

Dirichlet conditions impose a fixed surface value

A fixed surface temperature or concentration is represented by

[ \Theta(l,t)=\Theta_s. ]

After subtracting the steady boundary value, the homogeneous eigenfunction condition becomes

[ R_{\ell n}(l)=0. ]

For the spherically symmetric mode, this gives

[ \sin(\lambda_n l)=0, ]

and therefore

[ \lambda_n l=n\pi, ]

with the indexing chosen according to whether the zero mode is included.

This condition models a sphere whose surface is maintained by a well-controlled bath, jacket, or external concentration reservoir.

Neumann conditions impose a specified flux

An insulated or impermeable surface has zero normal flux:

[ -k_T\left.\frac{\partial\Theta}{\partial r}\right|_{r=l}=0 ]

for heat transfer, or

[ -D\left.\frac{\partial C}{\partial r}\right|_{r=l}=0 ]

for mass diffusion.

For the transformed spherically symmetric solution, this produces the characteristic equation

[ \lambda_n l=\tan(\lambda_n l). ]

More generally, the derivative of the appropriate spherical Bessel function must vanish:

[ \left.\frac{d}{dr}j_\ell(\lambda_{\ell n}r)\right|_{r=l}=0. ]

An insulated or impermeable boundary conserves the total thermal energy or species mass within the sphere, subject to any internal sources or sinks.

Robin conditions represent convective transfer

A convective surface is described by a mixed condition. For heat transfer,

[ -k_T\left.\frac{\partial\Theta}{\partial r}\right|_{r=l}

h_T\left[\Theta(l,t)-\Theta_\infty(t)\right]. ]

For mass transfer, the analogous condition is

[ -D\left.\frac{\partial C}{\partial r}\right|_{r=l}

k_m\left[C(l,t)-C_\infty(t)\right]. ]

After nondimensionalization, the transfer coefficient appears through a Biot-type parameter. For the spherically symmetric transformed solution, the characteristic equation can be written in a form such as

[ \lambda_n l

\left(1-A l^2\right)\tan(\lambda_n l), ]

depending on the definition of (A), the sign convention, and the transformed dependent variable.

The important physical point is that Robin conditions account for resistance outside the solid, including a fluid boundary layer, jacket-side film, or external mass-transfer film.

Nonhomogeneous conditions require a steady-state adjustment

Eigenfunction expansions are most straightforward when the boundary conditions are homogeneous. If the surface temperature or concentration varies with time, the field is commonly decomposed into a particular or quasi-steady part plus a transient part that satisfies homogeneous conditions.

This step prevents boundary terms from remaining unresolved when the diffusion equation is transformed.

Why Boundary Conditions Must Match the Basis

Boundary conditions define the eigenvalue problem

The differential operator alone does not determine the basis. The operator together with the center condition, surface condition, and any symmetry constraints defines the Sturm-Liouville eigenvalue problem.

Consequently, two spheres with identical material properties can require different eigenvalues if one has a fixed surface temperature and the other has convective heating.

Matching basis functions removes unwanted boundary terms

When finite Fourier or eigenfunction transforms are derived by integration by parts, boundary terms appear. A basis satisfying the same homogeneous boundary conditions makes those terms vanish or convert into known constants.

Using a sine basis for a convective boundary, for example, generally fails to represent the actual surface physics correctly. The resulting transformed problem can have incorrect eigenvalues, coefficients, or both.

Boundary conditions must reflect the actual module

For pilot plant equipment, the mathematical choice should follow the physical interface:

A controlled bath or imposed surface composition suggests Dirichlet behavior.
An insulated wall, symmetry plane, or impermeable surface suggests Neumann behavior.
A jacket, flowing gas, liquid film, or external bulk phase generally suggests Robin behavior.
An inlet, outlet, or reacting interface may require a separate specified-value or flux condition rather than a spherical surface condition.

Understanding the Trade-offs

Closed-form eigenvalues are limited

Dirichlet and some Neumann problems produce simple characteristic equations and explicit eigenvalue sequences. Robin problems usually lead to transcendental equations that require numerical root finding.

The additional numerical work is justified when external transfer resistance materially affects the process.

The sine-over-radius form has limited scope

Writing the radial solution as

[ \frac{\sin(\lambda r)}{r} ]

is appropriate for the spherically symmetric (\ell=0) mode. It is not the complete basis for arbitrary angular variation, where (j_\ell(\lambda r)) must be used.

Using the simpler form for a non-axisymmetric or strongly nonuniform loading can omit physically important angular modes.

Ideal assumptions may not hold in pilot equipment

The standard analytical basis assumes conditions such as constant properties, spherical geometry, and a linear diffusion or conduction equation. Real modules may include variable diffusivity, internal reaction, fluid motion, contact resistance, nonuniform surface transfer, or imperfect spherical shape.

These effects may require numerical methods such as finite-volume, finite-element, or spectral computation.

Experimental boundary conditions are rarely perfect

A supposedly fixed bath temperature may vary, and an estimated film coefficient may change with flow rate, viscosity, or concentration. Model validation should therefore compare measured surface and internal data, rather than treating the assumed boundary condition as exact.

The boundary condition is often one of the largest sources of uncertainty in inferred diffusivity or transfer coefficients.

Making the Right Choice for Your Goal

The basis should be selected from the geometry, symmetry, and interface physics of the pilot module.

If your primary focus is spherically symmetric transient behavior: Use the (\ell=0) spherical Bessel basis, or the equivalent (\sin(\lambda r)/r) representation, with the center symmetry condition and the appropriate surface characteristic equation.
If your primary focus is angularly nonuniform heating or concentration: Expand the solution in Legendre polynomials for axisymmetric systems or spherical harmonics for fully three-dimensional systems, using (j_\ell(\lambda r)) radially.
If your primary focus is a controlled surface temperature or concentration: Apply a Dirichlet condition and use eigenfunctions that vanish at the surface after subtracting the imposed boundary value.
If your primary focus is insulation, impermeability, or symmetry: Apply a Neumann condition with zero normal flux and use derivative-compatible eigenfunctions.
If your primary focus is jacket-side or film-limited transfer: Apply a Robin condition using the measured or estimated heat- or mass-transfer coefficient, then determine the eigenvalues numerically.
If your primary focus is agreement with pilot plant measurements: Verify that the assumed boundary conditions reflect the real inlet, outlet, wall, symmetry, and external-fluid conditions before fitting transport parameters.

Correctly matching spherical Bessel or harmonic basis functions to the physical boundary conditions turns the transient model into a defensible representation of the pilot plant rather than merely a formal mathematical solution.

Summary Table:

Aspect	Key choice	Application
Angular basis	Legendre polynomials or spherical harmonics	Axisymmetric or 3D loading
Radial basis	Spherical Bessel functions (j_\ell)	Regular solutions within solid spheres
Spherical symmetry	(\sin(\lambda r)/r) form	(\ell=0) transient behavior
Dirichlet condition	Fixed surface temperature or concentration	Controlled bath or reservoir
Neumann condition	Zero or specified flux	Insulated or impermeable surfaces
Robin condition	Convective heat or mass transfer	Jackets, films, and flowing phases
Eigenvalues	Roots of characteristic equations	Determined by boundary conditions
Model validation	Compare predicted and measured data	Pilot plant parameter estimation

Need a reliable model for transient heat or mass transfer in your pilot equipment? Contact LABPARK to discuss your application. LABPARK provides Educational and Vocational Unit Operations Pilot Plants for chemical engineering, bioprocess and biotechnology, and environmental and water treatment. Our systems help universities, research institutes, and enterprises study spherical transport behavior, evaluate heat- and mass-transfer performance, and connect mathematical models with practical experiments. We can help you select suitable modules, define realistic boundary conditions, and build training or research setups aligned with your process objectives.