  The period doubling phenomenon was recently discovered in RR Lyrae stars with
the Kepler space telescope and has been theoretically explained by hydrodynamic
calculations. However, peculiar solutions of the Florida-Budapest turbulent
convective hydrodynamic code suggest that bifurcation cascade may evolve to
chaos in these dynamical systems. We show that chaotic behaviour may be
recovered from the radius variations of the model using the global flow
reconstruction method. The fractal (Lyapunov) dimension of the underlying
dynamical attractor is calculated to be ~2.2. Compared to the radius, the
luminosity variations proved to be less suitable for such investigations due to
their complexity. That suggest that even the continuous Kepler data would
require transformation before conducting a similar analysis.
