In this work, the cohesion-specific inputs for a recent continuum theory for cohesive particles are estimated for moderately cohesive particles that form larger agglomerates via discrete element method (DEM) simulations of an oscillating shear flow. In prior work, these inputs (critical velocities of agglomeration and breakage and collision cylinder diameters) were determined for lightly cohesive particles via the DEM of simple shear flow—i.e., a system dominated by singlets and doublets. Here, the DEM is again used to extract the continuum theory inputs, as experimental measurements are infeasible (i.e., collisions between particles of a diameter of <100 μm). However, simulations of simple shear flow are no longer feasible since the rate of agglomeration grows uncontrollably at higher cohesion levels. Instead, oscillating shear flow DEM simulations are used here to circumvent this issue, allowing for the continuum theory inputs of larger agglomerate sizes to be determined efficiently. The resulting inputs determined from oscillating shear flow are then used as inputs for continuum predictions of an unbounded, gas−solid riser flow. Although the theory has been previously applied to gas−solid flows of lightly cohesive particles, an extension to the theory is needed since moderately cohesive particles give rise to larger agglomerates (that still readily break). Specifically, the wider distribution of agglomerate sizes necessitates the use of polydisperse kinetic-theory-based closures for the terms in the solids momentum and granular energy balances. The corresponding continuum predictions of entrainment rate and agglomerate size distribution were compared against DEM simulations of the same system with good results. The DEM simulations were again used for validation, as it is currently extremely challenging to detect agglomerate sizes and the number of fractions in an experimental riser flow.