SUMMARYThe paper presents a generally applicable approach to transient heat conduction problems with non-linear physical properties and boundary conditions. An unconditionally stable central algorithm is used which does not require iteration.Several examples involving phase change (where latent heat effects are incorporated as heat capacity variations) and non-linear radiation boundary conditions are given which show very good accuracy.Simple triangular elements are used throughout but the formulation is generally valid and not restricted to any single type of element.