Codex @codex  0

Side losses need not destroy the method. Let $\theta$ be the small temperature disturbance about a steady background in surroundings at fixed temperature. For a rod of cross-sectional area $S$ and perimeter $\mathcal{P}$, a uniform real surface heat transfer coefficient $h_s$ produces a loss per unit volume of $H\theta$, with $H=h_s\mathcal{P}/S$. Assume constant $k,\rho ,c_p,H$, a cross-section nearly uniform in temperature, and no distributed heat source. The balance becomesHere $H$ has units $\mathrm{W}{\mathrm{m}}^{-3}{\mathrm{K}}^{-1}$, while $\lambda$ has units ${\mathrm{s}}^{-1}$. This is a distributed version of Newton's law of cooling; a constant linearized contribution from thermal radiation may be included over a small temperature range. It does not represent arbitrary time-dependent coupling to a surrounding material with its own thermal storage.

For the decaying exponential branch, the phasor equation givesso thatReal positive side loss increases $a$ and decreases $b$ at fixed angular frequency and thermal diffusivity. Their product eliminates $\lambda$, yielding the combined inferenceThis result requires that attenuation and lag belong to the same exponential complex wavenumber at the same frequency, and that the loss term is real. It does not generally survive a superposition caused by a nearby finite end, spatially varying coefficients, transverse temperature gradients, or unequal sensor delays. The inferred $\lambda =\alpha (a^2-b^2)$ supplies an additional model check; a negative value signals inconsistency with this passive-loss model rather than negative physical cooling.

Consider an original stipulated model, not measured material data: thermal diffusivity $\alpha =10^{-4}{\mathrm{m}}^2{\mathrm{s}}^{-1}$, period $P=100\mathrm{s}$, driven temperature amplitude ${\Theta }_0=4\mathrm{K}$, and thermometers at $x_1=0.020\mathrm{m}$ and $x_2=0.050\mathrm{m}$. Without side losses, $a=b=17.7245{\mathrm{m}}^{-1}$ and thermal penetration depth $\delta =0.05642\mathrm{m}$. ThusThe two no-loss estimates and the combined estimate all return $10^{-4}{\mathrm{m}}^2{\mathrm{s}}^{-1}$.

Now retain the same thermal diffusivity and drive but introduce $\lambda =0.010{\mathrm{s}}^{-1}$. The exact complex wavenumber gives $a=19.1863{\mathrm{m}}^{-1}$ and $b=16.3742{\mathrm{m}}^{-1}$, henceUsing attenuation alone with the no-loss formula would give $8.5343\times 10^{-5}{\mathrm{m}}^2{\mathrm{s}}^{-1}$; using lag alone would give $1.17174\times 10^{-4}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. Their disagreement is informative. If an independent volumetric heat capacity measurement gives $\rho c_p=2.0\times 10^6\mathrm{J}{\mathrm{m}}^{-3}{\mathrm{K}}^{-1}$, the combined estimate implies thermal conductivity $k=200\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$ and $H=2.0\times 10^4\mathrm{W}{\mathrm{m}}^{-3}{\mathrm{K}}^{-1}$.

The exact complex amplitudes can be checked independently by discretizing $\alpha {\hat{\theta }}^{\prime \prime }-(\lambda +i\omega )\hat{\theta }=0$ and solving the resulting linear system with a tridiagonal matrix. Halving the spatial step quarters the discrepancy, as expected for centered finite differences. Separately fitting the real time traces to a constant plus sine and cosine terms recovers both amplitudes and their lag; this also checks the sign convention used in the combined estimate.

The Ångström periodic heating method measures thermal diffusivity by periodically heating one end of a rod and comparing the resulting temperature oscillations at two separated thermometers. The farther signal has a smaller amplitude and a later wave phase. Their combination measures how quickly the material redistributes heat, without requiring the absolute heater power. The method bears the name of Anders Jonas Ångström, the physicist who became a professor at Uppsala in 1858.^[ref]

An actual teaching apparatus is PASCO's TD-8513: four metal bars share a heater/cooler using the thermoelectric effect, and each has two embedded thermistors about $3\mathrm{cm}$ apart. Foam covers reduce environmental heat transfer. Its manual instructs students to measure sensor spacing and connect the thermistors through a temperature interface to a computer. For the aluminium bar, a $0$–$8\mathrm{V}$ triangular drive is applied with a $30\mathrm{s}$ period. Recording begins after the amplitudes stabilize; two cycles provide amplitudes and the delay between corresponding peaks. Repeated runs at $40$ and $50\mathrm{s}$ test consistency. The manual converts the inferred thermal diffusivity to thermal conductivity using mass density and specific heat capacity.^[ref]

For a general rod experiment, attach small thermocouples, thermistors or other thermometers at measured positions $x_1<x_2$, record both against the same clock, and establish a steady periodic thermal state. A sinusoid is especially convenient. For a switched or triangular drive, isolate the same Fourier series component at both sensors: unrelated peaks in distorted waveforms need not measure that component's phase difference. Subtract each sensor's mean temperature; a steady spatial background does not belong in the oscillating amplitude ratio. Small excursions help keep thermal conductivity, mass density and specific heat capacity constant.

The appeal is that the rod supplies two complementary clocks: decreasing amplitude measures spatial damping, while increasing wave phase lag measures delay. Unlike laser flash analysis, which follows a transient after one pulse, this method repeatedly samples a steady periodic thermal state. The timing determines thermal diffusivity $\alpha$ in ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. To infer thermal conductivity $k$ in $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$, also measure or otherwise justify the mass density $\rho$ and specific heat capacity at constant pressure $c_p$:Here the stationary-solid approximation neglects thermal expansion and mechanical work, as in thermal diffusivity. A timing measurement alone cannot distinguish two materials with the same $\alpha$ and different volumetric heat capacities.

Let a homogeneous plate $0<x<L$ begin at uniform temperature $T_0$, cool with zero-temperature Dirichlet boundary conditions until time $t_b$, then switch to insulated Neumann boundary conditions. For thermal diffusivity $\kappa >0$, the Fourier sine series during cooling is $T=(4T_0/\pi ){\sum }_{n\ge 1,n\text{ odd}}\sin (n\pi x/L)e^{-\kappa (n\pi /L{)}^{2}t}/n$. Integrating the insulated heat equation shows that its integral average is conserved. All nonconstant Fourier cosine series modes then decay, proving that the final uniform temperature is $T_{\infty }=(8T_0/{\pi }^2){\sum }_{n\ge 1,n\text{ odd}}{e}^{-\kappa (n\pi /L{)}^2{t}_b}/n^2$. Its limits are $T_0$ at zero bath time and zero at infinite bath time. Instantaneous changes of boundary condition can create initial-boundary corner incompatibilities without invalidating the positive-time solution.

A two-phase diffusion model with constant properties has similarity solutions built from the complementary error function in each half-line. Equilibrium gives $T_a=-m{C}_a$ and $C_s(a)=k_D{C}_a$. The Stefan condition and interfacial solute balance determine $\lambda$ and $C_a$; finite diffusion on both sides must be retained when constructing the narrow transition near $\lambda =0$.

In a radially stratified shearing sheet, one can normalize the scalar buoyancy perturbation by writing $b_x=-N^2\theta$ and $D_t\theta =u_x+\xi {\nabla }^2\theta$. With this convention $\theta$ has dimensions of length, rather than temperature; $\xi$ is thermal diffusivity.

For an imposed field $B_0$ in a layer of depth $d$, the Chandrasekhar number compares magnetic tension with viscous and magnetic-diffusive effects. Here $\eta$ is magnetic diffusivity and $\nu$ is kinematic viscosity. In thermal-diffusion units the linear magnetic-force coefficient is $Q\zeta$, where $\zeta =\eta /\kappa$ and $\kappa$ is thermal diffusivity.

With permeability $K$, depth $d$, kinematic viscosity $\nu$ and the relevant effective thermal diffusivity ${\kappa }_T$, the Darcy Rayleigh number compares buoyancy-driven Darcy advection to thermal diffusion. It differs from the clear-fluid Rayleigh number by replacing $d^3$ with $Kd$. Effective heat capacity and porosity factors may be absorbed into the definition of ${\kappa }_T$ and the thermal-time scale, and must be kept consistent with the chosen Darcy law convention.

Diffusion is net transport associated with random microscopic motion. A localized patch of a dilute solute can spread through a stationary solvent: individual molecules move in both directions, but more cross from the crowded side of an imaginary surface than from the less crowded side. This imbalance gives a net diffusive flux down the concentration gradient; the molecules do not need to sense which direction has lower concentration. At uniform concentration, microscopic exchanges continue even though the net diffusive flux vanishes. ^[ref]

The distinction from advection is easiest to see in a dye patch. Uniform bulk motion carries its centre downstream; diffusion broadens it relative to that centre. Both can happen together. A moving or spreading patch by itself does not identify the mechanism: a nonuniform velocity field can also stretch a patch. ^[ref]

For a dilute nonreacting solute in a uniform medium at fixed temperature, Fick's first law models the relative diffusive flux as $\mathbf{j}=-D\nabla C$, with scalar diffusion coefficient $D>0$ of units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Here $C$ is the amount of substance per unit volume, so $\mathbf{j}$ has units $\mathrm{mol}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$. With carrier velocity field $\mathbf{u}$, the total flux is $C\mathbf{u}+\mathbf{j}$. The continuity equation, with no sources or chemical reactions, givesFor constant $D$ and incompressible flow, this reduces to $C_t+\mathbf{u}\cdot \nabla C=D{\nabla }^{2}C$, where ${\nabla }^2$ is the Laplacian. In a stationary carrier it is the diffusion equation. These assumptions matter: variable $D$ must remain inside the divergence, and changing sources require extra terms. ^[ref]

The useful surprise is the common mathematics: thermal conduction smooths temperature through the analogous heat equation, while molecular diffusion smooths composition. This analogy already motivated the original formulation of Fick's laws. ^[ref] However, the thermal diffusivity and the solutal diffusivity need not be equal. A diffusive time scale $t_D\sim L^2/D$ means that crossing ten times the distance requires roughly one hundred times the time. This is why molecular diffusion is especially effective across thin layers, and why chaotic advection can assist mixing by making those layers thinner. Physical diffusion is a transport mechanism, distinct from the broader probabilistic class of Markov diffusion.

A finite rod generally contains both exponential branches. With prescribed oscillating temperature at $x=0$ and an insulated far end, ${\hat{\theta }}^{\prime }(L)=0$, the exact solution isThe correction near a sensor is controlled by $e^{-2a(L-x)}$. The important separation is between the farther sensor and the end, not merely between the heater and the end. For a far end clamped at the mean temperature, $\hat{\theta }(L)=0$, replace both hyperbolic cosines by hyperbolic sines; other end boundary conditions require their corresponding solution.

In the lossy Ångström method worked example, an insulated rod of length $L=0.080\mathrm{m}$ gives $D=0.354866$ and $\Delta \varphi =0.615651\mathrm{rad}$. Incorrectly applying the semi-infinite formula returns $1.29418\times 10^{-4}{\mathrm{m}}^2{\mathrm{s}}^{-1}$, about $29.4\%$ high. At $L=0.300\mathrm{m}$, it returns $9.99877\times 10^{-5}{\mathrm{m}}^2{\mathrm{s}}^{-1}$, about $0.0123\%$ low. Side loss remains real in both cases; the error arises because the measured spatial response no longer has one exponential complex wavenumber.

Choose the period and sensor spacing so that the lag and attenuation are resolvable while the downstream signal remains larger than sensor noise. Keep transverse equilibration fast compared with the period; compare the transverse thermal diffusion time with $P$ and check that surface exchange does not create large transverse temperature gradients. A short period reduces far-end influence but can expose sensor response times; a long period increases penetration and finite-end bias. Calibrate both sensors together, exchange their positions to check channel bias, compare successive cycles, and repeat at several spacings and periods. In a valid single-exponential region, linear regression of $\log A$ against $x$ and unwrapped lag against $x$ should give consistent slopes. These are proposed diagnostic controls. If they fail, fit the finite-rod or coupled-sensor model rather than interpreting every delay as intrinsic thermal diffusivity.

Fourier analysis describes a function by the oscillations that compose it. A graph records values at successive positions or times; a frequency spectrum records how much of each frequency is present, together with its wave phase. The same complicated shape can become a short list of simple components. For periodic functions, the components form a Fourier series; on the whole real line, the Fourier transform uses a continuous range of frequencies, and the Fourier inversion theorem puts them back together.

For a concrete example, consider the trigonometric polynomialWith $t$ measured in seconds, its three Fourier harmonics have frequencies $1$, $3$, and $5$ cycles per second and amplitudes $1$, $1/2$, and $1/4$. Each cosine has zero wave phase. Adding them creates the sharper peaks in the figure. The frequency spectrum makes the recipe immediately visible; changing an amplitude changes the shape without changing the component's frequency. Fourier orthogonality lets us recover each Fourier coefficient by averaging against the corresponding sine or cosine, because the other Fourier harmonics cancel over a full period of a function.

The deeper reason these components are so useful is that they simplify operations. By Euler's formula, the exponential function $e^{i\omega x}$ packages a cosine and a sine into one Fourier mode. A translation or a derivative changes only its multiplier:Thus Fourier modes are simultaneous eigenfunctions of these operations. In Fourier series, this gives a literal diagonalization in the Fourier basis; on the real line, the Fourier transform gives the continuous counterpart, with each frequency handled separately. The Translation property of the Fourier transform and the Fourier transform of a derivative turn shifting and differentiating a whole function into multiplication of its Fourier transform. A convolution likewise becomes multiplication by the convolution theorem. This is beautiful because one description simplifies several apparently different operations at once.

For example, the heat equation $u_t=\kappa u_{xx}$, with thermal diffusivity $\kappa >0$, evolves a spatial Fourier mode asFine ripples have larger angular frequencies and decay faster: the formula explains why heating smooths a jagged temperature profile. For a general initial profile with sufficient regularity, the Fourier transform gives $\hat{u}(\omega ,t)=e^{-\kappa {\omega }^{2}t}\hat{u}(\omega ,0)$, and Fourier inversion recovers the temperature. Similarly, multiplying selected parts of a frequency spectrum suppresses unwanted components. Claude E. Shannon's Communication in the Presence of Noise (1949), Section III, describes ideal filtering as a projection that retains some frequency coordinates and removes others.^[ref]

The physical motivation is already visible in Joseph Fourier's Théorie analytique de la chaleur (1822). Chapter III constructs heat profiles using sums of sines and cosines; Chapter IX, articles 355–359, passes from Fourier series to integrals as the spatial interval becomes unbounded. These heat problems explain the motivation for the discrete and continuous descriptions. Their modern formulations specify which functions admit reconstruction and in what sense the expansions converge.^[ref][ref]

For a diffusion equation with fixed homogeneous boundary conditions, a fundamental diffusion mode is an eigenfunction of its nonnegative decay operator belonging to the smallest decay eigenvalue ${\nu }_0$, when such an eigenvalue exists. It is the slowest decaying separated normal mode. With a discrete spectrum and an isolated lowest eigenvalue, initial data having a nonzero projection onto its eigenspace approach this decay law at late times. Strongly vanishing thermal diffusivity can destroy existence of a lowest mode, as illustrated by the loss of a spectral gap under strongly degenerate radial diffusion.

The half-space cooling model solves the one-dimensional heat equation for initially hot material cooled at a fixed-temperature surface. For constant thermal diffusivity $\kappa$, a fixed interior-temperature isotherm deepens proportional to $\sqrt{\tau }$. Applied to an oceanic lithosphere, it explains the leading age dependence of cooling; finite plate thickness and basal heat supply limit its extrapolation to old ages.

Laser flash analysis measures thermal diffusivity by heating one face of a specimen briefly and recording the opposite face's temperature rise. The original 1961 experiment used a flash tube rather than a laser: a blackened thin metal specimen sat in an opaque ceramic holder, with thermocouple wires pressed against its rear face. An amplifier, oscilloscope and camera recorded the signal. Uniform front illumination and the black coating promoted uniform absorption; the opaque low-thermal conductivity holder blocked direct rear illumination and reduced edge losses. The delayed rear response reveals internal thermal conduction rather than direct heating of the sensor. ^[ref]

A modern version uses a laser or flash lamp, a specimen of measured thickness $L$, controlled initial temperature, and a rear-facing infrared detector with a synchronized recorder. The readout is a temperature rise versus time, starting from a pre-pulse baseline. Thermal conductivity additionally requires mass density and specific heat capacity; it is not obtained from timing alone. ^[ref]

Under ideal one-dimensional, insulated, constant-property assumptions, instantaneous uniform absorption in a very thin front layer gives the normalized rear riseHere $T_0$ is the initial uniform temperature and $T_{\infty }$ is the final uniform temperature in the insulated model. Fitting this curve estimates thermal diffusivity; each term has the same $\alpha K_n^2$ decay as the insulated-bar Fourier modes. Equivalently, the half-rise time $t_{1/2}$ solves $R(t_{1/2})=1/2$. Writing $s=\alpha t/L^2$, the series gives $s_{1/2}\approx 0.13879$, henceThis coefficient comes from the stated ideal heat equation solution, not from an empirical sensor calibration. It has been independently checked by solving the series equation. Absolute absorbed pulse energy cancels from the normalized curve; determining heat capacity from its absolute height needs additional information. ^[ref]

Parker and colleagues identify finite flash duration, heat loss, nonuniform illumination, recording bandwidth and electrical interference as practical limits. A specimen too thin responds during the pulse; one too thick allows losses before the rear rise finishes. They shielded the electrical chain and recorded a baseline; their reported tests were preliminary metal measurements near ambient and moderately elevated temperatures. ^[ref]

Useful additional controls are to measure thickness at several points, verify the recorder's time scale and detector response independently, record a no-pulse baseline, and repeat at different pulse energies while avoiding phase transitions. Compare different thicknesses of the same homogeneous material: ideal half-rise times scale as $L^2$. Failure of these checks can expose finite-pulse effects, sensor delay, optical transmission, coating storage or edge losses; a thickness error is especially consequential because the inferred thermal diffusivity depends on $L^2$. These are proposed diagnostic controls, not invented measurements. When the pulse is not short or the sample exchanges appreciable heat with its surroundings, fit a model that includes those effects instead of applying the ideal half-rise coefficient unchanged.

The ratio of magnetic diffusivity $\eta$ to thermal diffusivity $\kappa$ controls their relative relaxation times. In nonrotating vertical-field magnetoconvection, oscillatory linear onset in the simple stress-free layer requires $0<\zeta <1$, so thermal perturbations diffuse faster than magnetic perturbations. This is distinct from the magnetic Prandtl number, $\nu /\eta$.

For a wall shear $u\simeq sz$, heating over streamwise length $L$ gives an outer temperature rise proportional to $L/(sz)$. Its transverse thermal conduction grows as $z^{-3}$ relative to the local source. If $\chi$ is thermal diffusivity, comparing conduction with heating gives $\delta \sim (\chi L/s{)}^{1/3}$, up to a numerical constant depending on the precise comparison. The same scale follows by comparing transverse diffusion time $z^2/\chi$ with residence time $L/(sz)$.

Cooling oceanic lithosphere contracts and increases its column mass, causing isostatic deepening. Integrating a half-space cooling profile gives an excess mass proportional to $\sqrt{t}$. Converting this to seafloor depth includes the density contrast between mantle and seawater. Finite-plate cooling and basal heat supply limit indefinite square-root subsidence at old ages.

Use positive jumps $\Delta T=T-T_{\infty }$ and $\Delta S=S-S_{\infty }$ for the hot salty lower reservoir. The Boussinesq density excess below the interface is proportional to $\beta \Delta S-\alpha \Delta T$. Static stability therefore requireswith one the neutral limiting value. For $R_{\rho }<1$ the lower layer is lighter overall and ordinary overturning convection replaces the statically stable interface.

In the stable hot-below, salty-below configuration, thermal diffusivity exceeds salt diffusivity. Heat transfer can destabilize thin adjacent fluid regions before the stabilizing salt difference is erased, causing the diffusive regime of double-diffusive convection. As $R_{\rho }\to \infty$, stabilization by salt dominates, interfacial convective exchange is strongly suppressed, and the specified heat-flux model tends to zero at fixed $\Delta T$. This limit is not a claim that molecular heat conduction itself ceases.