A generalised cascade algorithm for steady state simulation of multiple effect evaporator systems is discussed. The algorithm, capable of handling any feed arrangement, includes heat recovery features such as liquor and condensate flash units and feed preheating. This is made possible by the use of composite flow fractions which fully describe the internal flow connections. The user-provided models for these units and all effects may be of varying degrees of sophistication.