Due to demands of lower costs and higher productivity in the steel industry, the volume of operating blast furnaces has grown during the last decades. As the height is limited by the allowable pressure drop, the hearth diameter has grown considerably and, along with this, also draining-related problems. In this paper a mathematical model is developed for simulating the drainage in the case where an impermeable region exists in the blast furnace hearth. The model describes the quasi-stationary drainage process of a hearth with two operating tapholes, where the communication between the two pools of molten slag and iron can be controlled by parameterized expressions. The model also considers the case where the buoyancy of the liquids is sufficient for lifting the coke bed. The implications of different size of the liquid pools, communication between the pools, bed porosity, etc. are studied by simulation, and conclusions concerning their effect on the drainage behavior and evolution of the liquid levels in the hearth are drawn. The simulated liquid levels are finally demonstrated to give rise to a pressure profile acting on the hearth which agrees qualitatively with signals from strain gauges mounted in the hearth wall of an industrial ironmaking process.