Leukocyte filters are applied clinically to remove leukocytes from blood. In order to optimize leukocyte filters, a mathematical model to describe the leukocyte filtration process was developed by modification of a general theoretical model for depth filtration. The model presented here can be used to predict the time-dependent leukocyte filtration as a function of cell-cell interaction in the filter, filter efficiency, filter capacity, filter dimensions, and leukocyte concentration in the suspension applied to the filter. The results of different leukocyte filtration experiments previously reported in the literature could be well described by the present model.