Prediction of axial liquid velocity profile in bubble columns