A scalable new mathematical model based on the principles of Process Integration has been developed for the analysis of multiple effect evaporator (MEE) systems. It uses the concepts of stream analysis, temperature path and internal heat exchange for the formulation of the model equations. In addition to the above concepts, the model also takes into account the variable physico-thermal properties of steam/vapor, condensate and liquor, while simulating the MEE system. The present model consists of a set of linear equations and does not present any stability or oscillation problems during solution as is generally seen in the case of models that are based on sets of nonlinear equations. The model equations are automatically generated through a computer program and the model was run for three different liquor and flow sequences to prove its utility. The results obtained are compared with published models.