Thermal diffusion rate of a Fourier mode ($\beta =\xi k^2$)

With thermal diffusivity $\xi$, a Fourier mode of angular wavenumber $k$ decays at rate $\beta =\xi k^2$. To avoid confusing this $k$ with thermal conductivity, use $\alpha =\xi$ and angular wavenumber $K$ below: the same rate is $\beta =\alpha K^2$, with units ${\mathrm{s}}^{-1}$.

Consider an insulated bar $0\le x\le L$ with constant thermal diffusivity, no sources and no heat flux through either end or the sides. An original initial temperature profile and its exact evolution areThe cosine has zero temperature gradient at both ends, satisfying the insulated Neumann boundary conditions. Its spatial mean is zero, so conservation of energy leaves the mean temperature $\overline{T}$ unchanged while the nonuniform part disappears. A short derivation substitutes $T=\overline{T}+A(t)\cos (K_{n}x)$ into the heat equation: the second spatial derivative supplies $-K_n^2$, hence $A^{\prime }=-\alpha K_n^{2}A$. The exponential decay time is ${\tau }_n=1/(\alpha K_n^2)$; doubling angular wavenumber quarters this time. If cycles-per-length spatial frequency $f_x$ is used instead, $K=2\pi f_x$ and the rate is $4{\pi }^2\alpha f_x^2$.

For the two model materials in thermal diffusivity and $L=0.010\mathrm{m}$, ${\tau }_{1,A}\approx 5.07\mathrm{s}$ and ${\tau }_{1,B}\approx 20.3\mathrm{s}$. After $10\mathrm{s}$, the remaining amplitude fractions are approximately $0.139$ and $0.610$. A narrow hot-and-cold pattern disappears faster than a broad one because its temperature curvature, and hence the imbalance in conductive heat flux, is larger.