A multicomponent mathematical model is presented for a large-scale slurry bubble column reactor operating in the heterogeneous or churn-turbulent flow regime. The model accounts for both the Fischer-Tropsch reaction as well as the water gas shift reaction and the individual paraffin and olefin formation rates. It provides all the data necessary for reliable scale up, process optimization and prediction of the performance of industrial scale Fischer-Tropsch bubble column reactors.