In the present paper, the possibilities to cover all fluidization states by a single drag law in steady-state CFD multiphase flow modeling are evaluated. The time-averaged drag force is expressed as the product of the drag force calculated from the traditional drag laws for homogeneous conditions, and a correction function. Closure correlations for the correction function are developed by nonlinear regression modeling based on data collected from 69 transient 2D simulations of bubbling, turbulent and circulating fluidized beds. The correlations are given as functions of eight variables: the solid volume fraction, the distance from the nearest wall, the height above the air distributor, the slip velocity between the phases, the gas velocity, the particle size, the solid density and the gas viscosity. The results indicate that covering all fluidized bed conditions in a single drag correlation is feasible, although fully satisfactory results were not obtained for the surface and freeboard regions of a bubbling fluidized bed (BFB) with correlations that were acceptable in circulating fluidized bed (CFB) conditions. A correlation that covers the whole range of fluidized states is complicated and thus the modeling task could be divided into development of separate correlations for different regions that could be combined into a single correlation by means of blending functions. The validity of this approach was demonstrated by developing a separate correlation for the dilute conditions above a height of 1.5 m in CFB risers. Results show that the accuracy of the predictions significantly improved in dilute CFB conditions where a much simpler correlation with six input variables could be used. The modeling approach is a good starting point for the development of a general drag law for CFD simulations of fluidized beds. (C) 2014 Elsevier B.V. All rights reserved.