The purpose of the present paper is to predict the grain size of steel during the hot-rolling process. The basis represents a macroscopic simulation system that can cope with temperatures, stresses and strains of steel in a complete continuous rolling mill, including reversible pre-rolling and finishing rolling with several tenths of rolling passes. The grain size models, newly introduced in the present paper, are one-way coupled to the macro-scale calculations performed with the slice model assumption. Macroscale solution is based on a novel radial basis function collocation method. This numerical method is truly meshless by involving the space discretization in arbitrarily distributed nodes without meshing. A new efficient node generation algorithm is implemented in the present paper and demonstrated for irregular domains of the slice as they appear in different rolling passes. Multiple grain size prediction models are considered. Grain size prediction models are based on empirical relations. Austenite grain size at each rolling pass as well as the ferrite grain size at the end of rolling are predicted in this simulation. It is also shown that based on the rolling schedule, it is highly likely that recrystallization takes place at each pass throughout a continuous rolling mill. The simulation system is coded as a user-friendly computer application for industrial use based on programing language C# and an open source developer platform .NET and runs on regular personal computers The computational time for a typical rolling simulation is usually less than one hour and can thus be straightforwardly used to optimize the rolling mill design in a reasonable time.