AbstractIn this paper, we develop a CFD model for describing a bubble column reactor for carrying out a consecutive first‐order reaction sequence A → B → C. Three reactor configurations, all operating in the homogeneous bubbly regime, were investigated: (I) column diameter DT = 0.1 m, column height HT = 1.1 m, (II) DT = 0.1 m, HT = 2 m, and (III) DT = 1 m, HT = 5 m. Eulerian simulations were carried out for superficial gas velocities UG in the range of 0.005–0.04 m/s, assuming cylindrical axisymmetry. Additionally, for configurations I and III fully three‐dimensional transient simulations were carried out for checking the assumption of cylindrical axisymmetry. For the 0.1 m diameter column (configuration I), 2‐D axisymmetric and 3‐D transient simulations yield nearly the same results for gas holdup ϵG, centerline liquid velocity VL(0), conversion of A, χA, and selectivity to B, SB. In sharp contrast, for the 1 m diameter column (configuration III), there are significant differences in the CFD predictions of ϵG, VL(0), χA, and SB using 2‐D and 3‐D simulations; the 2‐D strategies tend to exaggerate VL(0), and underpredict ϵG, χA, and SB. The transient 3‐D simulation results appear to be more realistic. The CFD simulation results for χA and SB are also compared with a simple analytic model, often employed in practice, in which the gas phase is assumed to be in plug flow and the liquid phase is well mixed. For the smaller diameter columns (configurations I and II) the CFD simulation results for χA are in excellent agreement with the analytic model, but for the larger diameter column the analytic model is somewhat optimistic. There are two reasons for this deviation. Firstly, the gas phase is not in perfect plug flow and secondly, the liquid phase is not perfectly mixed. The computational results obtained in this paper demonstrate the power of CFD for predicting the performance of bubble column reactors. Of particular use is the ability of CFD to describe scale effects.