Describir: An improved contact modification routine for a computationally efficient CFD simulation of packed beds