A modified two dimensional (2-D) Monte Carlo (MC) technique was used to simulate primary recrystallization in automotive steels containing fine particles. In order to consider anisotropic properties of grain boundary energy and grain boundary mobility, functions of boundary misorientation were introduced. Orientation-dependent stored energy developed in 80% cold-rolled interstitial free (IF) sheet steel was evaluated by reconstructing of data measured using electron back-scattered diffraction (EBSD) analysis. A subgrain method based on subgrain structure is used for quantitative analysis of the stored energy. The simulation reveals that particles affect evolution of microstructure during recrystallization. The simulation provided a theoretical foundation for understanding effect of particles on the final microstructures and crystallographic textures.