Spatially averaged film thickness models have been the preferred approach for simulation of falling film reactors. This paper will discuss both the general and specific limitations of available models. An improved model is proposed, which has been found to perform well in industrially relevant conditions. For its specific application to sulphonation/sulphation reactions yielding anionic surfactants, the improved model provides enhanced predictions about reactor performance, including relative colour intensity at different operational conditions.