Researchers have compared two key computational approaches, the "fluid Ansatz" and the "WallGo" code, to determine the velocity of bubble walls in cosmological phase transitions. These transitions are crucial for understanding the early universe, and the bubble velocity ($v_w$) is a fundamental parameter for predicting gravitational wave signals. The study reveals that both methods agree closely in the regime of reasonably mild phase transitions, with a strength parameter $\alpha \lesssim 0.01$.

The agreement is particularly good when only top-quark annihilation is considered. However, a noticeable discrepancy appears once scattering processes are included. The work also investigates the limitations of linearizing the Boltzmann equation when applying the fluid Ansatz to stronger phase transitions. It is observed that non-linear contributions induce significant shifts in the predicted terminal velocity as $\alpha \to 1$, even though the non-linear contribution to the wall pressure remains quantitatively small compared to the equilibrium and linearized non-equilibrium parts.

These findings have important implications for the interpretation of future gravitational wave observations. Strong phase transitions are precisely the primary targets for future gravitational wave observatories. Therefore, the study emphasizes the need for higher precision computations of $v_w$ in the semi-classical approach and suggests that a treatment beyond the WKB approximation may be needed to adequately understand these extreme early universe events.