We performed numerical and experimental studies on the viscous folding in diverging microchannel flows which were recently reported by Cubaud and Mason (Phys Rev Lett 96:114501, 2006a). We categorized the flow patterns as ''stable'', ''folding,'' and ''chaotic'' depending on channel shape, flow ratio, and viscosity ratio between two fluids. We focused on the effect of kinematic history on viscous folding, in particular, by changing the shape of diverging channels: 90°, 45°, and hyperbolic channel. In experiments, the proposed power-law relation (f $ _ c 1 ; where f is the folding frequency, and _ c is the characteristic shear rate) by Cubaud and Mason (Phys Rev Lett 96:114501, 2006a) was found to be valid even for hyperbolic channel. The hyperbolic channel generated moderate flows with smaller folding frequency, amplitude, and a delay of onset of the folding compared with other two cases, which is considered to be affected by compressive stress when compared to the simulation results. In each channel, the folding frequency increases and the amplitude decreases as the thread width decreases since higher compressive stress is applied along the thin thread. The secondary folding was also reproduced in the simulation, which was attributed to locally heterogeneous development of compressive stresses along the thread. This study proves that the viscous folding can be controlled by the design of flow kinematics and of the compressive stresses at the diverging region.