We derive a simplified model for two-dimensional (2D) channel flows with recirculated regions at moderate Reynolds numbers based on an extension of the boundary layer (BL) theory and averaging across the channel. The model reproduces symmetry-breaking bifurcations and resulting flow structures accurately. Analytical estimates for the decay rates toward the parabolic profile before and after a sudden change in the walls agree well with the full numerical simulations. A seemingly chaotic steady flow is also discovered in a channel with periodic expansions and contractions.