Order is found within the chaotic nonlinear flow between rotating coaxial cylinders. The flow stability analysis is carried out for a pseudoplastic fluid through bifurcation diagram and Lyapunov exponent histogram. The fluid is assumed to follow the Carreau-Bird model, and mixed boundary conditions are imposed. The low-order dynamical system, resulted from Galerkin projection of the conservation of mass and momentum equations, includes additional nonlinear terms in the velocity components originated from the shear-dependent viscosity. It is observed that the base flow loses its radial flow stability to the vortex structure at a lower critical Taylor number, as the shear-thinning effects increase. The emergence of the vortices corresponds to the onset of a supercritical bifurcation, which is also seen in the flow of a linear fluid. However, unlike the Newtonian case, shear-thinning Taylor vortices lose their stability as the Taylor number reaches a second critical number corresponding to the onset of a Hopf bifurcation. Complete flow field together with viscosity maps are given for different scenarios in the bifurcation diagram. 810 N. Ashrafi range of excited spatio-temporal scales. In order to assess the effect of the smaller length scales on the flow, a high resolution of the flow field is needed. It is by now well established that low-order dynamical systems can be a viable alternative to conventional numerical methods for the linear and nonlinear flows [2,3]. Despite the severe degree of truncation in the formulation of these models, some of the basic qualitative elements of the onset of Taylor vortices and destabilization of the cellular structure can be recovered using low-order dynamical systems.A number of studies have been directed on the issue of CCF of Newtonian fluids in the finite gaps between rotating cylinders. A comprehensive glossary of dynamic situations of the flow of the Newtonian fluids between coaxial cylinders is recently collected by Hoffmann et al both theoretically and experimentally [4]. The exchange of stability based on the variation of a control parameter among known flows i.e., CCF, TVF, SPI (spiral TVF), wavy vortex flow (WVF), and wavy spiral vortex flow (WSPI) is explained in their work. It is resulted that the wavy structures appearing via a secondary bifurcation out of Taylor vortex flow and out of spiral vortex flow, respectively, mediate transitions between Taylor and spiral vortices and vice versa. Many of these findings are also predictable by the low-dimensional approach taken in the present study, for the case of non-Newtonian flows. This approach was originally used in the seminal work of Lorenz to investigate the stability of Rayleigh-Bernard convection problem [5]. Later, Kuhlmann [3] extended the method to the TVF of Newtonian fluids in the narrow gap, followed by the work of Berger [6] to cover finite gaps and Ashrafi and Khayat for the narrow gap flow of weakly shear-thinning fluids [2].Flow computations of nonlinear fluid flows, especially the TVF, require extra ca...