A method has been developed for the solution of inverse heat diffusion problems to find the initial condition, boundary condition, and the source and sink function in the heat diffusion equation. The method has been used in the development of a source‐and‐sink method to find the boundary conditions in inverse Stefan problems. Green's functions have been used in the solution, and the problems are solved by using two approaches: a series solution approach, and a time incremental approach. Both can be used to find the boundary conditions without reliance on the flux information to be supplied at both sides of the interface. The methods are efficient in that they require less equations to be solved for the conditions. The numerical results have shown to be accurate, convergent, and stable. Most of all, the results do not degrade with time as in other time marching schemes reported in the literature. Algorithms can also be easily developed for the solution of the conditions.