Statement

For smooth rapidly decaying initial data u0u_0 and a smooth source f(t,x)f(t,x) with a common compact spatial support on a finite time interval, a solution of

tuνΔu=f,u(0)=u0,\partial_tu-\nu\Delta u=f,\qquad u(0)=u_0,

is given by the heat Duhamel formula

u(t)=Ttu0+0tTtsf(s)ds,u(t)=T_tu_0+\int_0^t T_{t-s}f(s)\,ds,

where TtT_t is the .

Verification and scope

Differentiating the integral gives the endpoint value f(t)f(t); its remaining terms give νΔ0tTtsf(s)ds\nu\Delta\int_0^t T_{t-s}f(s)\,ds. For these smooth data, spatial derivatives may instead fall on ff, justifying differentiation near s=ts=t. The integral vanishes at time zero. Rougher sources require the corresponding convergence and regularity estimates before the formula is interpreted as a classical solution.