A mathematical model for the adsorption of biomolecules on affinity membranes is presented. The model considers convection, diffusion and adsorption kinetics on the membrane module as well as the influence of dead end volumes and lag times; an analysis of flow distribution on the whole system is also included. The parameters used in the simulations were obtained from equilibrium and dynamic experimental data measured for the adsorption of human IgG on A2P-Sartoepoxy affinity membranes. The identification of a bi-Langmuir kinetic mechanisms for the experimental system investigated was paramount for a correct process description and the simulated breakthrough curves were in good agreement with the experimental data. The proposed model provides a new insight into the phenomena involved in the adsorption on affinity membranes and it is a valuable tool to assess the use of membrane adsorbers in large scale processes.