In the present study a mathematical model has been applied to interpret the permeate flux decay that occurs during the process of concentrating skim milk by ultrafiltration using a commercial membrane module. The effects on membrane fouling of two operational variables, temperature and transmembrane pressure, have been studied using technical parameters. An energy analysis has demonstrated that the major energy consumption takes place in the thermal process and not in the mechanical pumping of the fluid. In addition, higher increments in permeate volume can be achieved by increasing transmembrane pressure, not temperature. The mathematical analysis presented here permits the evaluation of optimum values of the engineering parameters necessary to design and operate skim milk ultrafiltration units. (c) 2008 Elsevier Ltd. All rights reserved.