In the present investigation, a numerical three-dimensional (3-D) heat flow model for friction stir welding (FSW) has been developed, based on the method of finite differences. The algorithm, which is implemented in MATLAB 5.2, is provided with a separate module for calculation of the microstructure evolution and the resulting hardness distribution. The process model is validated by comparison with in-situ thermocouple measurements and experimental hardness profiles measured at specific time intervals after welding to unravel the strength recovery during natural aging. Furthermore, the grain structure within the plastically deformed region of the as-welded materials has been characterized by means of the electron backscattered diffraction (EBSD) technique in the scanning electron microscope (SEM). Some practical applications of the process model are described toward the end of the article.