In this study a new numerical method which does not require a regularization parameter is developed for solving inverse and ill-posed problems associated with both steady and unsteady heat conduction applications. Various types of convergence and accuracy criteria, and boundary condition formulations are investigated which confirm that the numerical method developed produces a convergent and stable numerical solution with respect to the number of boundary elements and number of iterations.