This work is dedicated to analyze the heat removal phenomena during the simulation of the continuous casting process of the steel squared sections; it is a useful method to produce big mounts of steel and some other metals such as aluminum. Radiation, forced convection and thermal conduction are the physical phenomena involved. The method used to solve the heat removal and distribution is a finite difference. Steel is discretized using a regular squared mesh. Here a 3D problem is simplified to a 2D in which every node represents a 3D steel volume, but the assumptions are appropriately justified in order to reduce computing time without sacrifice precision on calculation. The efficiency of the method is also evaluated and alternatives to improve the approaching and results are commented. © 2012 American Institute of Physics.