In the present study, flow simulation has been carried out in a backward-facing step channel defined by an expansion ratio of 2.02 and a spanwise aspect ratio of 8 to provide the physical insight into the longitudinal and spanwise flow motions and to identify the presence of Taylor–Görtler-like vortices. The Reynolds numbers have been taken as 1000 and 2000, which fall in the category of transitional flow. The present simulated results were validated against the experimental and numerical data and the comparison was found to be satisfactory. The simulated results show that the flow becomes unsteady and exhibits a three-dimensional nature with the Kelvin–Helmholtz instability oscillations and Taylor–Görtler-Like longitudinal vortices. The simulated data were analysed to give an in-depth knowledge of the complex interactions among the floor and roof eddies, and the spiralling spanwise flow motion. Destabilization of the present incompressible flow system, with the amplified Reynolds number due to the Kelvin–Helmholtz and Taylor–Görtler instabilities, is also highlighted. A movie is available with the online version of the paper.