Taking the spatial Fourier transform of the diffusion equation gives and . Hence . Fourier inversion and the Gaussian transform from part (a), or direct evaluation of the inverse Gaussian integral, give the heat kernel
Substituting the integral for and interchanging absolutely integrable terms produces the convolution solution
For the stated nonnegative compactly supported data, understood to be nonzero on a set of positive measure, the kernel is strictly positive everywhere. Therefore for every real and every . Property P holds: diffusion has infinite propagation speed. The distant value can be extremely small, but is not zero. For continuous initial data, any nontrivial nonnegative automatically has the required positive-measure set.
Thermal diffusion length 2026-10-14
The thermal diffusion length is the characteristic length over which thermal conduction changes a temperature distribution during time in the constant-thermal diffusivity model. It is a scale, not a sharp front or distance travelled at a constant speed: the heat equation has infinite propagation speed mathematically. Different definitions of penetration depth introduce numerical factors. In particular, the normalized one-dimensional heat kernel has spatial variance , so its root-mean-square width is .
Figure 1. . Original predictions for the two stipulated materials in thermal diffusivity. Left: , with four times the time needed to double length. Right: on logarithmic axes. Equal thermal conductivity with four times the volumetric heat capacity halves the reached length at a given time and quadruples the time at a given length. The scale differs from the insulated-bar decay time by the factor . Details.
At , the model materials have thermal diffusion lengths and . For their thermal diffusion times are and . These values describe the same heat equation as the thermal diffusion rate of a Fourier mode; the different numerical factors come from the definition of the scale.
Thermal diffusion wave 2026-10-18
A thermal diffusion wave is a periodically forced solution of the heat equation, with spatially decreasing amplitude and varying wave phase. In a semi-infinite rod , assume constant thermal diffusivity, negligible side losses, a nearly uniform temperature across each cross-section, and the boundary conditions and as . A steady periodic thermal state then satisfies
Choose the square root in the complex numbers with so that the amplitude decays into the rod. Substitution of the phasor into the heat equation proves the formula: . Taking the real part gives the displayed cosine. In this convention is the spatial attenuation slope and is the spatial lag slope, both in ; combines them into one complex number. If a complex wavenumber is defined by instead, .
With , let and define a positive downstream lag . Then
and either measurement determines thermal diffusivity:
The logarithm uses an amplitude ratio, not a ratio of squared amplitudes. The lag is in radians and must be tracked continuously across cycles: a measured phase difference modulo leaves an integer-cycle ambiguity. Use closer sensor spacings or vary frequency continuously to resolve that ambiguity. Equal attenuation and lag slopes provide a useful diagnostic of the no-side-loss model.
The apparent phase velocity is , which depends on frequency. A thermal diffusion wave therefore does not give a material-independent travel speed for heat; it is a forced pattern in a dissipative diffusion equation. The ideal heat equation still has infinite propagation speed, and the observed lag is not the arrival time of a sharp front.