In the lower zone of blast furnace (BF), slag trickles down through a packed bed of carbonaceous particles in the form of films, rivulets or droplets. During its downward flow, there are significant interactions occurring between slag and other phases, ie gas and packing particles. In particular, the interaction between gas and slag is closely associated with the key flow phenomena such as loading, flooding or channelling in the lower zone of the BF. With gas introduced laterally through tuyeres and slag flowing downward from the cohesive zone to the hearth, this interaction varies throughout the lower zone and in the cohesive zone:  strong cross-flow of gas relative to the slag occurs in front of the raceway  counter-current flow of gas and slag occurs in the upper parts of the lower zone  both strong cross- and counter-current flow of gas and slag around and within the cohesive zone. These interactions can be linked to furnace irregularities, affecting smooth operation and limiting production. Therefore, understanding the gas-slag interaction and its influencing factors remains critical for BF process control. In the current review, the physical behaviour and key influencing factors of gas and slag flow are initially examined. Numerical models describing gas and slag flows in a packed bed at micro-, mesoand macro-scopic levels, are reviewed. Each simulation approach has its advantages and limitations, as well as unique assumptions when used in BF process modelling. Applications of these approaches to gas-slag flow in key identified functional zones will highlight the differences between approaches and modelling outcomes. As no single approach can fully address all aspects of gasslag interactions, an integration of different approaches is considered an effective and beneficial method in providing a more complete understanding of gas-slag flow in the BF process.