Prediction of heterogeneous reactive azeotropes in esterification systems