An outstanding problem in the field of relativistic heavy-ion collisions is the inability for any simulation model to accurately describe experimental flow data in extremely central collisions -in particular, models always predict either an elliptic flow that is too large or triangular flow that is too small (or both). We reassess the status of this puzzle in light of recent progress in Bayesian parameter estimation, in which a large model parameter space can be efficiently explored to determine what parameters are necessary for a good fit to experimental results, and how well state-of-the-art models are able to describe data. We explore predictions for flow in ultra-central collisions from multiple recent Bayesian models that were tuned to various observables in different collision systems at typical centralities. We find that, while ultra-central data can now be described with better accuracy than in previous calculations, tension with experimental observation remains, and progressively gets worse as one goes to more central collisions. Thus, the physics of ultra-central collisions is still not fully understood. A resolution to this puzzle will be an important step in solidifying our understanding of the physics of the strong interactions in these extreme conditions, as well as increasing our confidence in the results of precision analyses.