An algorithm is presented for the design and simulation of multieffect evap
orators, with forward or backward feed. The mass and enthalpy balances and
the heat transfer rate equations posed around each effect are linearized as
suggested by Lambert et al. [1], but a new iterative scheme is proposed to
ensure convergence and stability when the solution boiling-point rise is c
onsidered Boiling-point rise can reduce the thermal driving force significa
ntly leading to important errors in the calculation.
The algorithm is based on the elimination of an unknown, namely, the produc
t stream flow rate, for simulation in the balances set by assuming a given
product concentration and the reduction of the heat transfer rate equations
set by elimination of the corresponding equation. The simulation problem i
s thus converted into a design that can be sah ed by known numerical method
s rendering the required area. If that is not available, the calculation is
repeated assuming a new product composition, selected by a criterion prese
nted here, until convergence on the actual area is attained.