Numerical investigation of the Liebau phenomenon