Abstract
A new robust Monte Carlo approach is developed and proposed for solving heat conduction problems. This method is very beneficial when implemented within integrated multiscale physical systems where the use of linear equations solvers is limited due to low scalability and the requirement of large memory for data storage. The developed method is based on introducing virtual “thermal particles” and has the potential of wide capabilities in various complex multidimensional simulations in plasma physics and other applications. The numerical results show good agreement with analytical methods as well as other numerical solutions.