Many key forging components of heavy equipment are manufactured by large steel ingots. Macrosegregation in steel ingots is a key defect formed during the solidification process. Over the past few decades, numerical modeling has played a more and more important role in the study of macrosegregation. Various models have been developed and applied to different ingot casting processes. This paper focused on the application of macrosegregation models to the steel ingot. Firstly, the formation mechanism and influencing factors of macrosegregation were introduced. Then, the existing macrosegregation models and their recent development were summarized. Macrosegregation models accounting for such mechanisms as solidification shrinkage-induced flow and mushy zone deformation were analyzed, respectfully. To model macrosegregation due to solidification shrinkage, the key was to solve the free surface. A simple derivation showed that the multi-phase (including gas phase) models were equivalent to the VOF-based segregation models in dealing with the shrinkage-induced flow. Finally, our recent research work on numerical modeling of macrosegregation in steel ingots was illustrated, including application of the developed multi-component and multi-phase macrosegregation model to a 36 t steel ingot, and numerical simulation of multiple pouring process. The carbon and sulphur concentrations at about 1800 sampling points, covering the full section of a 36 t ingot, were measured. By detailed temperature recording, accurate heat transfer conditions between the ingot and mould were obtained. Typical macrosegregation patterns, including the bottom-located negative segregation and the pushpin-like positive segregation zone in the top riser, have been reproduced both in the measurements and the predictions. The carbon and sulphur concentrations predicted by the three dimensional multi-component and multi-phase macrosegregation models agreed well with the measurements, thus proving that the model can well predict the macrosegregation formation in steel ingots. As for the multi-pouring process simulation, the results show a high concentration of carbon at the bottom and a low concentration of carbon at the top of the mould after the multi-pouring process with carbon content high in the first ladle and low in the last ladle. Therefore, the multiple pouring process could get the initial solute distribution with the opposite form of segregation. Such carbon concentration distribution would reduce the negative segregation at the bottom and the positive segregation at the top of the solidified ingot, thus proving the ability of the multiple pouring process for the control of macrosegregation.