Numerical simulation of extraction columns using a drop population model