A numerical simulation procedure is proposed for analyzing the partial replacement of pulverized coal injection by hydrogen, oxygen, and blast furnace gas (BFG) injections mixed with pulverized coal (PCI) within the tuyeres of large blast furnaces. The massive use of hydrogen-rich gas is extremely interesting for ironmaking blast furnaces in the context of net-zero carbon hot metal production. Likewise, this new approach allows for increasing productivity and for reducing the specific emissions of carbon dioxide toward a net-zero carbon ironmaking technology. Nevertheless, the mixture of pulverized coal injection and gas injection is a complex technology. In addition to the impact on chemical reactions and energy exchange, the internal temperature and gas flow patterns can also change drastically. With a view to assessing the state of the furnace in this complex operation, a comprehensive mathematical model utilizing multiphase theory was developed. The model simultaneously handles bulk solids (sinter, pellets, small coke, granular coke, and also iron ore), gas, liquid metal and slag, and coal powder phases. The associated conservation equations take into account momentum, mass, chemical species, and energy while being discretized and solved using finite volume techniques. The numerical model was validated against the reference operating conditions using 220 kg per ton of pig iron (kg/tHM) of pulverized coal. Therefore, the combined injection of different concentrations of fuel hydrogen, blast furnace gas, and oxygen was simulated for replacing 40, 60, and 80 kg/tHM of coal injection. Theoretical analysis showed that the best scenario with stable operation conditions could be achieved with a productivity increase of 20% corresponding to a CO2 reduction of 15% and 60 kg/tHM of PCI replacement.