The steady-state behaviour of binary and multicomponent packed distillation columns is simulated using two film plug flow and back-mixing models. Orthogonal collocation on finite elements is demonstrated to be a powerful and efficient method of solution of these models. Close agreement is obtained with experimental results reported by other workers suggesting that the procedures described form a good basis for column design.