Bridgman directional solidification of investment castings is a key technology for the production of reliable and highly efficient gas turbine blades. In this paper, a mathematical model for three-dimensional (3D) simulation of solidification process of single crystal investment castings was developed based on basic heat transfer equations. Complex heat radiation among the multiple blade castings and the furnace wall was considered in the model. Temperature distribution and temperature gradient in superalloy investment castings of single blade and multiple ones were investigated, respectively. The calculated cooling curves were compared with the experimental results and agreed well with the latter. It is indicated that the unsymmetrical temperature distribution and curved liquid-solid interface caused by the circle distribution of multiple turbine blades are probably main reasons why the stray grain and other casting defects occur in the turbine blade.