Turbulent airflow in channels of rectangular cross section with symmetric centerbodies is studied numerically. Shock wave configurations formed in the channel and in front of the entrance are examined. Solutions of the unsteady Reynolds-averaged Navier-Stokes equations are obtained with finite-volume solvers of second-order accuracy. The solutions demonstrate an expulsion/swallowing of the shocks with variations of the free-stream Mach number or angle of attack. Effects of the centerbody length and thickness on the shock wave stability and flow bifurcation are examined. Bands of the Mach number and angle of attack, in which there exist non-unique flow fields, are identified.