Gas-solid multiphase flows are prone to develop an instability known as clustering.Two-fluid models, which treat the particulate phase as a continuum, are known to reproduce the qualitative features of this instability, producing highly-dynamic, spatiotemporal patterns. However, it is unknown whether such simulations are truly aperiodic or a type of complex periodic behavior. By showing that the system possesses a sensitive dependence on initial conditions and a positive largest Lyapunov exponent, λ 1 ≈ 1/τ , we provide a tentative answer: continuum predictions of clustering are chaotic. We further demonstrate that the chaotic behavior is dimensionally dependent, a conclusion which unifies previous results and strongly suggests that the chaotic behavior is not a result of the fundamental kinematic instability, but of the secondary (inherently multidimensional) instability.a) E-mail: william.fullmer@colorado.edu b) Corresponding Author, E-mail: hrenya@colorado.edu 1 Granular matter is a collection of discrete, interacting solid particles which, like classic (molecular) matter, can be generally classified into one of three states 1 : i) under static conditions granular heaps or piles can sustain gravity-induced stress and behave like a solid 2 ; ii) dense granular flows characterized by enduring and multi-particle contacts behave similarly to a fluid 3 ; and iii) rapid granular flows characterized by instantaneous contacts described as a granular gas 4 . However, it is worth noting that all three granular states are only superficially similar to their molecular counterparts due to the dissipative nature of particle-particle contacts (inelasticity, friction, etc.) 1,4 . In industrial operations, rapid (collision-dominated) granular flows are often encountered in devices in which the particles are fluidized by a gas 5 ; such gas-solid flows are the focus of this work. Rapid gas-solid flows are prone to an instability termed clustering in which particles tend to form spatially inhomogeneous patterns of high and low concentrations 6 . Continuum or two-fluid models have long been known to be able to predict the qualitative nature of the clustering instability 7 and more recent quantitative assessments have also shown promising results 8-10 However,it is yet unknown whether or not such predictions are chaotic. We take a first step in answering this question by simulating fluidization in an unbounded domain and calculating a positive largest Lyapunov exponent, thereby indicating that continuum predictions of clustering are in fact chaotic. Further, we show that chaotic behavior may be reduced to periodic behavior by constraining the dimensionality of the system.