The increasing interest in minimising fuel consumption and reducing environmental problems of the ironmaking industry has led to efforts to optimise the blast furnace process. One of the undeveloped fields generated from the adoption of new operating conditions is the internal phenomenon in the lower part of the furnace following the change of cohesive zone shape. The effective liquid flow in the dripping zone of a blast furnace is important for stable operation with high iron productivity. Study of the liquid flow behaviour in a packed bed needs to investigate the effect of various operational changes in the dripping zone. As a stepping stone toward the comprehensive modelling which includes chemical reaction and its associated heat transfer, isothermal flow for three phases is simulated with the variation of the shape of cohesive zone to see the effect of charging condition.