The computational fluid dynamics numerical forecasts for the pulsations of sequences of W Virginis stellar models exhibit two approaches to irregular behavior that are a clear signature of low dimensional chaos. The first indication comes from first return maps in which one plots one maximum radius, or any other suitable variable, versus the next one. The sequence of models shows a period doubling bifurcation, or cascade, leading to chaos. The near quadratic shape of the map is indicative of chaos and implies an underlying horseshoe map. Other sequences of models follow a somewhat different route, but also to chaos, namely the Pommeau–Manneville or tangent bifurcation route.