Process optimization of selective esterification of glycerol with lauric acid catalyzed by highly uniformed SBA-15 catalysts post impregnated with 12-tungstophosphorus acid (HPW) has been elucidated via response surface methodology (RSM). The significance of factors including reaction time (t), reactants ratio (R) and temperature (T) was investigated and led to the development of response models for lauric acid conversion (R1) and monolaurin yield (R2). The optimized monolaurin yield operating condition (50% yield at t=3.18 h, R=4.6 and T=175 ℃) was identified by using the developed models with high accuracy (98% confidence level, small standard deviation and error) and reproducibility (3 times). Comparison done with reported results showed the good performance of developed catalyst with significant reduction in reaction time needed to achieve similar lauric acid conversion and monolaurin yield.