Abstract. This paper presents a mathematical model based on EulerlLagrange approach for time-dependent calculations of two-(gas-liquid) and three-(gasliquid-solid) phase flows in a bubble column. The fluid phase flow was calculated based on the Euler approach by solving the three-dimensional Reynolds-averaged conservation equations in a time dependent way. The set of equations was closed using the standard k-e turbulence model. Two-way coupling was accounted for by adding dispersed phase source terms in all conservation equations of the continuous phase and additionally considering wake-induced turbulence. Bubble motion was calculated by solving the equation of motion taking into account drag force, pressure, added mass force, transverse lift force, buoyancy and gravity. The tracking of solid particles was based on considering drag force, pressure, added mass force, Saffman force, Magnus force and gravity. The introduction of the effective density in the continuous phase conservation equations permitted to perform the calculations up to relatively high volume fractions of the dispersed phase. The interaction between bubbles and particles was considered by adding simple correction terms to the drag coefficient of both bubbles and particles. Alternatively, the interaction was described by collisions between bubbles and particles using a stochastic collision model. The comparison of the predicted results with experimental data from the literature showed reasonable good agreement for the cases of two-phase and three-phase flows.