Modeling of breakthrough curves in fixed-bed zeolite columns