A mathematical model for the four-phase (gas, powder, liquid, and solids) flow in a two-dimensional ironmaking blast furnace is presented by extending the existing two-fluid flow models. The model describes the motion of gas, solid, and powder phases, based on the continuum approach, and implements the so-called force balance model for the flow of liquids, such as metal and slag in a blast furnace. The model results demonstrate a solid stagnant zone and dense powder hold-up region, as well as a dense liquid flow region that exists in the lower part of a blast furnace, which are consistent with the experimental observations reported in the literature. The simulation is extended to investigate the effects of packing properties and operational conditions on the flow and the volume fraction distribution of each phase in a blast furnace. It is found that solid movement has a significant effect on powder holdup distribution. Small solid particles and low porosity distribution are predicted to affect the fluid flow considerably, and this can cause deterioration in bed permeability. The dynamic powder holdup in a furnace increases significantly with the increase of powder diameter. The findings should be useful to better understand and control blast furnace operations.