A high-resolution, quasi-geostrophic numerical model is utilized to examine two-layer baroclinic flow in a cylinder. Particular attention is given to the role of horizontal shear of the basic state induced by viscosity near the cylinder wall, and to the desymmetrization brought about by the cylindrical geometry, in the transition to baroclinic chaos. Solutions are computed for both f-plane and β-plane situations, and the results are compared to previous laboratory experiments. Agreement in the former case is found to be good, although the onset of chaos occurs at slightly lower forcing in the laboratory when its basic flow is prograde, and at higher forcing amplitude when the experimental basic azimuthal currents are retrograde. This suggests that the modest discrepancies may be attributable to computationally neglected ageostrophic effects in the interior fluid and Ekman boundary layers. When β ≠ 0, the numerical and laboratory results are in excellent agreement. The computational simulations indicate that the viscous sidewall boundary layer plays a pivotal role in the dynamics. Moreover, in contrast to previous studies performed in a periodic, rectilinear channel, the route to chaos is largely temporal and involves relatively few spatial modes. The reduction in symmetries upon going from f-plane channel to either f-plane or β-plane cylinder models leads to fewer secondary instabilities and fewer spatial modes that are active in the dynamics.
Absolute/convective secondary instabilities and the role of confinement in free shear layers