This paper extends grid-free Monte Carlo methods like walk on spheres and walk on stars to solve time-dependent diffusion problems with initial and boundary conditions, eliminating the need for volumetric meshing and sequential time stepping. The method uses random walks with finite time budgets and kernel sampling techniques to directly estimate solutions at any requested time while maintaining parallel and progressive evaluation properties.