<p>Oxygen-enriched bottom-blown smelting is a core technology in non-ferrous metal smelting due to its advantages of wide feedstock adaptability, high smelting intensity, and low energy consumption. However, the extremely high temperature, airtight operating conditions, and multiphase environment inside the furnace make it difficult for conventional experiments to fully reveal the physicochemical mechanisms of the molten bath. As a low-cost, repeatable tool, numerical simulation effectively overcomes this limitation. This paper systematically reviews the research progress of multiphase flow numerical simulation and its influencing factors in bottom-blown furnaces (BBFs). First, it reviews the development history of smelting processes and numerical simulation technologies for BBFs, before elaborating on the furnace structure, working principles, and commonly used multiphase flow and turbulence models. On this basis, the influence mechanisms of key parameters such as bottom-blowing injector structure, inclination angle, arrangement mode, spacing, molten bath depth, and gas flow rate on molten bath flow, mixing efficiency, splashing behavior, gas holdup, and reaction kinetics are analyzed in detail. Finally, this paper discusses the limitations of current isothermal models and highlights future research needs in non-isothermal multi-physics coupling, high-temperature validation, and refractory lining erosion prediction.</p>