The design and control of blast furnace (BF) ironmaking must be optimized in order to be competitive and sustainable, particularly under the more and more demanding and tough economic and environmental conditions. To achieve this, it is necessary to understand the complex multiphase flow, heat and mass transfer, and global performance of a BF under different conditions. Mathematical modeling, often coupled with physical modeling, plays an important role in this area. This paper reviews the recent developments in this direction. The emphasis is given to mathematical models for different BF regions from the top charging system, body, and finally down to raceway and hearth. The needs for the further research and developments are also discussed.