The temperature of a rod of length L can be described by the 1D non-dimensional heat equation aT at = a(T); 072 + Q (x, t) (1) with time t, temperature T, temperature dependent thermal diffusivity o (7), and a source term describing additional heating/cooling along the rod, Q(x, t). On the left end of the rod, the boundary has a prescribed temperature gradient, aT Or (