The mean velocity versus β for different dc-bias, obtained by numerical solving Eqs. (1) and (2) with uniformly distributed initial condition and time interval . The remaining parameters are , a = 1.512, and .