The present study is concerned with breakup models for agglomerates in turbulent flows. First, a brief literature review is provided describing the state of the art concerning the description and modeling of breakup of agglomerates with special emphasis on the role of fluid forces. That comprises turbulent and drag (inertia) stresses. Furthermore, the breakup by rotary stresses is taken into account. Based on the idea to describe these processes using first principles within an Eulerian-Lagrangian approach relying on the large-eddy simulation methodology, modeling ideas from the literature are revisited and extended in order to be applied in a four-way coupled simulation. Building on a deterministic collision and agglomeration model, the breakup process is tackled in a similar manner leading to a highly efficient procedure, which allows to track a huge number of agglomerates and particles but with the restriction to compact dry powder agglomerates. For this kind of agglomerates an enhanced evaluation of the strength is derived relying on an improved determination of the average porosities and coordination number. Furthermore, special emphasis is put on the post-breakup treatment, i.e., models for the sizes and velocities of the generated fragments. Finally, the proposed models are applied to a test case inspired by the experimental study of Weiler [1] for the deagglomeration in a dry powder disperser. The effect of the different physical mechanisms (turbulence, drag, rotation) is analyzed concerning breakup rates and resulting size distributions. Furthermore, the dominant breakup regions are identified separately. Since agglomeration processes are simultaneously taken into account, re-agglomeration of the fragments is found to play a non negligible role. Finally, the effect of the properties of the powder on the breakup and re-agglomeration processes is investigated. (C) 2019 The Author(s). Published by Elsevier B.V.