In this study, a solid acid catalyst SO42−/TiO2/ZSM-5 (STZ) was prepared. Biodiesel was synthesized from tall oil fatty acid (TOFA) by the ethanol esterification. The fixed bed and pervaporation (PV) coupling increased the esterification conversion of TOFA and ethanol by about 16 %. The spherical STZ catalyst’s appearance, composition, and catalytic activity were assessed through TGA, XRD, SEM, BET, FTIR, NH3-TPD and XPS analyses. The reaction parameters of ethanol/TOFA molar ratio (10:1–20:1), reaction temperature (75–85 ℃), catalyst dosage (6–12 g) and flow rate (0.3–0.9 mL/min) were optimized by the response surface method (RSM). Under RSM optimization, the TOFA esterification reaction achieved a conversion rate of 99.56 ± 0.23 %. After five cycles of esterification, TOFA conversion remained above 87 %. A computational fluid dynamics (CFD) model was established to predict the esterification reaction behavior in the reactor, and the TOFA conversion rate was found to be close to the experimental results through systematic simulation, by reaction kinetics equation, continuity equation, Navier-Stokes equation, and Brinkman equation. In addition, the TOFA biodiesel prepared by this study complies with EN 14214 and ASTM D6751 standards.