Modelling of breakthrough curves in fixed-bed zeolite columns