Most bubble breakage models have been developed for multiphase simulations using Euler-Euler (EE) approaches. Commonly, they are linked with population balance models (PBM) and are validated by making use of Reynolds-averaged Navier-Stokes (RANS) turbulence models. The latter, however, may be replaced by alternate approaches such as Large Eddy simulations (LES) that play a pivotal role in current developments based on lattice Boltzmann (LBM) technologies. Consequently, this study investigates the possibility of transferring promising bubble breakage models from the EE framework into Euler-Lagrange (EL) settings aiming to perform LES. Using our own model, it was possible to reproduce similar bubble size distributions (BSDs) for EL and EE simulations. Therefore, the critical Weber (Wecrit) number served as a threshold value for the occurrence of bubble breakage events. Wecrit depended on the bubble daughter size distribution (DSD) and a set minimum time between two consecutive bubble breakage events. The commercial frameworks Ansys Fluent and M-Star were applied for EE and EL simulations, respectively. The latter enabled the implementation of LES, i.e., the use of a turbulence model with non-time averaged entities. By properly choosing Wecrit, it was possible to successfully transfer two commonly applied bubble breakage models from EE to EL. Based on the mechanism of bubble breakage, Wecrit values of 7 and 11 were determined, respectively. Optimum Wecrit were identified as fitting the shape of DSDs, as this turned out to be a key criterion for reaching optimum prediction quality. Optimum Wecrit values hold true for commonly applied operational conditions in aerated bioreactors, considering water as the matrix.