A study has been conducted of the modeling and control of an industrial multiple-effect evaporator system used for black liquor recovery. The developed model is a combination of phenomenological and empirical (neutral networks) approaches, based on industrial data. An advance model-predictive control strategy, allowing for manipulated variable constraints, compares favorably with a classical PID-based scheme.