The transition between non-aerated and aerated regimes in uncovered unbaffled stirred tanks (UUSTs) was investigated using computational fluid dynamics (CFD). The Volume of Fluid (VOF) method was employed to model the free surface dynamics under various operational conditions. The simulations were able to predict the peculiar experimentally observed behavior of UUSTs and revealed that at velocities below a critical threshold (N_<_Ncrit), the system remains non-aerated, while exceeding this threshold (N_≥_Ncrit) induces bubble ingestion, leading to significant changes in power consumption and flow patterns. The CFD simulations accurately predicted the behaviour of the Power Number (Np) as well as the vortex shape inside the tank both in subcritical and supercritical regimes and showed good agreement with original experimental data and correlations from the literature. Additionally, the modeling of the aerated regime successfully predicted the vortex shape, the bubble dispersion within the tank, and the cavities formed behind the blades.