Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(BCRs) under unsteady state control is far from enough. In order to study the flow structures in this operation, the volume of fluid (VOF) model and the standard k-ε model to simulate the evolution of gas-liquid flow in BCRs under the start-up state are combined. For both the symmetry and asymmetry flow, the layout of the gas-inlets, the gas-in velocity, the liquid viscosity and the aspect ratio of the BCR all have effects on the liquid velocity distribution. The simulation results could provide some information for the design and scale-up of the BCRs.