A valveless pump consisting of a pumping chamber with an elastic tube was simulated using an immersed boundary (IB) method. The interaction between the motion of the elastic tube and the pumping chamber generated a net flow toward the outlet throughout a full cycle of the pump. The net flow rate of the valveless pump was examined by varying the stretching coefficient (phi), bending coefficient (gamma), the aspect ratio (l/d) of the elastic tube, and the frequency (f) of the pumping chamber. As the stretching and bending coefficients of the elastic tube increased, the net flow through the valveless pump decreased. Elastic tubes with aspect ratios in the range of 2 <= l/d <= 3 generated a higher flow rate than that generated for tubes with aspect rations of l/d = 1 or 4. As the frequency of the pumping chamber increased, the net flow rate of the pump for l/d = 2 increased. However, the net flow rate for l/d = 3 was nonlinearly related to the pumping frequency due to the complexity of the wave motions. Snapshots of the fluid velocity vectors and the wave motions of the elastic tube were examined over one cycle of the pump to gain a better understanding of the mechanism underlying the valveless pump. The relationship between the average gap in the elastic tube and the average flow rate of the pump was analyzed. A smaller gap in the elastic tube during the expansion mode and a wider gap in the elastic tube during the contraction mode played a dominant role in generating a high average flow rate in the pump, regardless of the stretching coefficient (phi), the aspect ratio (l/d) of the elastic tube, or the pumping frequency of the pumping chamber (f). (C) 2012 Elsevier Inc. All rights reserved.