Simulation of the grain growth process, as a function of steel heat transfer conditions, is helpful for predicting grain structures of continuous cast steel products. Many authors have developed models based on numerical methods to simulate grain growth during metal solidification. Nevertheless, the anisotropic nature of grain structures makes necessary the employment of new mathematical methods such as chaos theory, fractals, and probabilistic and stochastic theories of simulation. The problem is significant for steelmakers to avoid defects in products and to control the steel microstructure during the continuous casting process. This work discusses the influence of nodal solidification times and computer algorithms on the dynamic formation of the chill, columnar, and equiaxed zones including physical phenomena such as nucleation and grain growth. Moreover, the model incorporates pre-nucleation and pre-growth routines in the original algorithm. There is a description of the influence of the mathematical parameter criteria and probabilities over the grain morphology obtained after solidification. Finally, an analysis of these algorithms elucidates the differences between these structures and those obtained from models considering only the solidification.