In this paper, we develop and analyze the use of the Optimal Homotopy Asymptotic Method (OHAM) to find the approximate analytical solution for partial differential equation involving fuzzy heat equation. OHAM allows the solution of the partial differential equation to be calculated in the form of an infinite series in which the components can be easily computed. The convergence theorem of this method in fuzzy case is presented and proved. This method provides us with a convenient way to control the convergence of approximation series. The method is tested on linear fuzzy heat equations and comparing the exact solution that were made with numerical results showed the effectiveness and accuracy of this method. © 2016 American Scientific Publishers All rights reserved.