The derivation for Heat Equation I am reading starts by stating
Net change of heat inside $[x,x+\Delta x]$ = Net flux of heat across boundaries + Total heat generated inside $[x,x+\Delta x]$ and writes the conservation equation
$$\textit{Total Heat Inside} [x,x+\Delta x]= cpA \int _{ x}^{x+\Delta x}u(s,t) ds $$
and later equates net flux across the boundary with
$u_x(x+\Delta x,t) - u_x(x,t)$,
so the heat flow is proportional to the spatial gradient of temperature at both ends. What I dont understand is why is net flux proportional to the spatial gradient, why cant we use the temparature variable $u$ itself to measure the difference? If the temparature difference at both ends is high, then we would expect higher flow from high temparature to low, right? What would be the physical intuition behind this? Was this model arrived at by experimentation?

