A mathematical model of an isothermal semibatch bubble reactor has been developed to describe the esterification reaction of free fatty acids with superheated alcohol vapor. The proposed model accounts for the effects of mass transfer followed by chemical reaction in the liquid phase. The fluid physical properties are calculated by published correlations, and the partition coefficient of alcohol vapor was estimated based on thermodynamic models of vapor-liquid equilibria. Experimental data of acid-catalyzed esterification of oleic acid with superheated ethanol vapor at different conditions of temperature and gas superficial velocities were used to estimate the liquid side mass transfer coefficients and kinetic parameters. The results obtained showed that the model satisfactorily fits the experimental data for all operating conditions and was also able to simulate and predict experimental results for intermediate conditions with a coefficient of determination R-2 of 0.987. In addition, the estimated values of the mass transfer coefficient agreed with data reported in the literature for gas-liquid reactors.