This article presents a novel model to simulate melting of a phase change material, with natural convection and radiation. For the phase change problem, the enthalpy formulation is used. Energy equation is solved by finite differences, whereas fluid flow equations are solved by the lattice Boltzmann method. For radiation intensity, the radiative transfer equation is solved by the discrete ordinates method, and then radiation flux is added into the energy equation. The model is first validated with literature results. Then, a glass brick wall filled with a fatty acid is simulated to evaluate the heat transfer processes. The results show that (1) natural convection plays an important role in the transitional behaviour of the global heat transfer process and (2) the long wave radiation has few impacts on the melting process. (c) 2012 Elsevier Ltd. All rights reserved.