Abstract A mathematical model for the catalytic autothermal reforming (ATR) reaction of synthetic crude glycerol to hydrogen in a fixed bed tubular reactor (FBTR) and over an in-house developed metal oxide catalyst is presented in this work. The heterogeneous model equations account for a two-phase system of solid catalyst and bulk feed gas. Also, the ATR of crude glycerol reaction scheme and intrinsic kinetic rate model over an active, selective, and stable nickel-based catalyst were integrated in the developed model. Also, the model was validated using experimental data generated in our labs for the ATR of synthetic crude glycerol. The modelling results adequately described the detailed gas product composition and distribution, temperature profiles, and conversion propagation in the axial direction of the fixed bed reactor over a wide range of reaction temperature (773−923 K) and weight-time (12.71−158.23 g cat·min·molC−1). The crude glycerol conversion predicted with the model showing a close resemblance to those obtained experimentally with an average absolute deviation (AAD) of less than 8%. The maximum crude glycerol conversion and hydrogen yield were found to be 92% and 3 mol hydrogen/mol crude glycerol, respectively. Also, the gas product concentration profile in the reactor was adequately described (90%) accuracy with a hydrogen concentration of 39 vol%.