A Dual Reciprocity Boundary Element Method is formulated to solve non-linear heat conduction problems. The approach is based on using the Kirchhoff transform along with lagging of the effective non-linear thermal diffusivity. A posteriori error estimate is used to provide effective estimates of the temporal and spatial error. A numerical example is used to demonstrate the approach