Numerous periodic and aperiodic dynamic states obtained in a model for hydrogen peroxide decomposition in the presence of iodate and hydrogen ions (the Bray-Liebhafsky reaction) realized in an open reactor (CSTR), where the flow rate was the control parameter, have been investigated numerically. Between two Hopf bifurcation points, different simple and complex oscillations and different routes to chaos were observed. In the region of the mixed-mode evolution of the system, the transitions between two successive mixed-mode simple states are realized by period-doubling of the initial state leading to a chaotic window in which the next dynamic state emerges mixed with the initial one. It appears in increasing proportions in concatenated patterns until total domination. Thus, with increasing flow rate the period-doubling route to chaos was obtained, whereas with decreasing flow rate the peak-adding route to chaos was obtained. Moreover, in very narrow regions of flow rates, chaotic mixtures of mixed-mode patterns were observed. This evolution of patterns repeats until the end of the mixed-mode region at high flow rates that corresponds to chaotic mixtures of one large and many small amplitude oscillations. Starting from the reverse Hopf bifurcation point and decreasing the flow rate, simple small amplitude sinusoidal oscillations were encountered and then the period-doubling route to chaos. With a further decreasing flow rate, the mixed-mode oscillations emerge inside the chaotic window.