A generalized mathematical model has been developed to simulate the non-isothermal performance of bubble column reactors. Steady-state reactor behavior is rigorously described without limitations on the extent of macro-mixing, flow configuration, absorption-reaction regime, and nature of reactive species. A parametric sensitivity analysis shows that the idealized lumped-distributed models do not offer a realistic portrait of the non-isothermal behavior of bubble column reactors. Implications on the possible occurrence of static bifurcations are discussed.